I am using ggpubr and ggplot2. First I did a anova test on my 3 clusters and now I would like to include my p-adj values of the TukeyHSD test in my box plot by using
stat_compare_means or stat_pvalue_manually
since I found out that Tukey test cannot be shown with stat_compare_means unfortunately..
Thanks