【原】aov()和ggplot()函数实现单因素方差分析多重比较字母标记

ggpubr包适合针对两个组的差异比较可视化，这个下一期给大家介绍。本期主要介绍使用aov()和ggplot()函数实现单因素方差分析多重比较字母标记。

方差分析aov()函数的代码，见以下两个文章：

用glht()和plot()函数实现方差分析主效应多重比较，附字母标记

用R语言做单因素方差分析及多重比较

选用R语言自带案例数据集PlantGrowth，研究两个处理和一个对照组对植物产量的影响，每组10例共3记录，主要考察处理对提高植物产量有无影响。

数据构成：因变量weight，因子变量group，三个水平依次为ctrl、trt1、trt2。

sig1 <- aov(weight ~ group,data=PlantGrowth)summary(sig1)

方差分析表：

三组间的weight平均水平时有统计学差异的（P＜0.05）。接下来lsd多重比较。

mylsd <- LSD.test(sig1,"group",p.adj="bonferroni")mylsd

多重比较的字母标记结果如下表：

接下来用ggplot函数绘制箱图并添加多重比较的字母标记结果。

绘图代码见我的陪伴统计社群。

===