R语言——ggplot2 + ggpubr包添加显著性

#加载包
library(ggplot2)
library(ggpubr)

p <- iris %>% 
  ggplot(aes(Species, Sepal.Length)) + 
  geom_violin(aes(fill = Species),trim = F) + 
  stat_boxplot(geom = 'errorbar', width = 0.1) + 
  geom_boxplot(aes(fill = Species), width = 0.2, show.legend = F) +
  geom_jitter(aes(shape = Species)) + 
  stat_summary(fun.y = median, geom = 'point', color = 'white', size = rel(2)) + 
  guides(shape = F, color = F) + 
  theme_bw() + 
  scale_fill_brewer(palette = 'Set2')

mycomparision <- list(c('setosa', 'versicolor'), c('setosa', 'virginica'), c('versicolor', 'virginica'))  #设置比较组
bartlett.test(Sepal.Length ~ Species)
                #Bartlett test of homogeneity of variances

 #data:  Sepal.Length by Species
 #Bartlett's K-squared = 16.006, df &