I am a new user and I am happy to be now part of this community ^^ I am having issues with how to add the adjusted pvalues (after the Bonferroni correction) on top of the boxplot (using the packages ggplot2) for statistical significance.

I want to add " for p< 0.001, for p < 0.0001, for p < 0.00001" on top of each bar plot, but don't know which function should I use for the same purpose.

Thanking you.

boxplot(MEs_tdf1[,-1], boxfill = NA, border = NA) #invisible boxes - only axes and plot area
boxplot(MEs_tdf1[MEs_tdf1$cluster=="1", -1], xaxt = "n", add = TRUE, boxfill="red",
        boxwex=0.25, at = 1:ncol(MEs_tdf1[,-1]) - 0.15) #shift these left by -0.15
boxplot(MEs_tdf1[MEs_tdf1$cluster=="2", -1], xaxt = "n", add = TRUE, boxfill="blue",
        boxwex=0.25, at = 1:ncol(MEs_tdf1[,-1]) + 0.15)