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
