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 = 2, p-value = 0.0003345,方差不齐用非参数检验。

p1 <- p + stat_compare_means(method = 'wilcox.test', comparisons = mycomparision, label = 'p.signif') #以显著性（星星）显示结果
p1

p1.png

p2 <- p + stat_compare_means(method = 'wilcox.test', ref.group = 'setosa', label = 'p.signif') #第二三组与第一组比较
p2

p2.png

cowplot::plot_grid(p1, p2, labels = c('p1', 'p2'), label_fontface = 4)#组合图形

p3.png