#加载包
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](https://i-blog.csdnimg.cn/blog_migrate/2923d3a6d4fbc403e0dab9e4447e9c76.png)
p2 <- p + stat_compare_means(method = 'wilcox.test', ref.group = 'setosa', label = 'p.signif') #第二三组与第一组比较
p2
![p2.png](https://i-blog.csdnimg.cn/blog_migrate/d64ca7c449f174cf2b8ad2f5af0c2290.png)
cowplot::plot_grid(p1, p2, labels = c('p1', 'p2'), label_fontface = 4)#组合图形
![p3.png](https://i-blog.csdnimg.cn/blog_migrate/07a2943ea60a726094f3281263673505.png)