如何在R中为双向ANOVA添加基于Tukey HSD的事后检验标签

针对双向ANOVA的Tukey HSD事后检验标签添加需求，推荐使用emmeans包自动生成分组字母，彻底避免手动标注错误，具体实现步骤如下：

1. 安装并加载所需包 #

install.packages("emmeans")
library(emmeans)
library(ggplot2)

2. 处理模型并提取分组标签 #

使用emmeans获取交互项的均值，再通过cld()函数基于Tukey检验规则自动生成分组字母：

# 你的原始数据
df = data.frame(variable = c('2021','2021','2021','2021',
                             '2022','2022','2022','2022'), 
                value = c(6,7.3,6.5,7.2,6.3,7.2,6.5,7.4),
                type = c('lawn','meadow','lawn','meadow',
                         'lawn','meadow','lawn','meadow'))

# 拟合双向ANOVA模型
M1 <- aov(value ~ type * variable, data = df)

# 提取交互项的均值并生成Tukey分组字母
emm <- emmeans(M1, ~ type * variable)
tukey_groups <- cld(emm, adjust = "tukey", Letters = letters)

# 转换为数据框并整理格式（去除分组字母中的空格）
tukey_df <- as.data.frame(tukey_groups)
tukey_df$.group <- gsub(" ", "", tukey_df$.group)

3. 绘制带标签的ggplot图形 #

将自动生成的分组标签添加到均值误差棒上方：

ggplot(df, aes(x = variable, y = value))+
  stat_summary(fun = mean, geom = 'point', na.rm = T, size = 5)+
  stat_summary(fun.data = mean_cl_normal, geom = 'errorbar', width = .5, na.rm = T)+
  # 添加事后检验标签（位置偏移至误差棒上方0.1处）
  geom_text(data = tukey_df, aes(x = variable, y = upper.CL + 0.1, label = .group), 
            size = 5, fontface = "bold")+
  facet_wrap(~type)+
  theme_bw()

方案核心优势 #

• 零手动错误风险：cld()函数直接根据Tukey检验的p值生成分组字母，相同字母代表组间无显著差异，不同字母代表差异显著；
• 适配双向交互项：完美处理type * variable这类交互效应的分组逻辑，无需手动整理两两比较结果；
• ggplot无缝集成：生成的分组数据框可直接用于geom_text/geom_label，无需额外格式转换。

内容的提问来源于stack exchange，提问作者Trisha Miller