如何为Tukey HSD图添加p值？代码报错排查求助

错误原因 #

• contrasts can be applied only to factors with 2 or more levels 错误直接源于Type列仅含单一水平（"CONTROLE"）：Tukey HSD用于多组间差异检验，必须至少2个不同分组才能计算对比p值。
• 额外隐患：value列用逗号作为小数点分隔符，直接转数值类型会失败，导致后续统计分析无法正常运行。

分步解决方案 #

1. 数据预处理修正 #

先解决数据格式问题，确保数值列和因子列符合分析要求：

# 替换逗号小数点为点，转换为数值类型
df$value <- as.numeric(gsub(",", ".", df$value))
# 将Type和Roof转换为因子类型
df$Type <- factor(df$Type)
df$Roof <- factor(df$Roof)

2. 过滤无效分组 #

提前过滤掉Type仅含单一水平的Roof分组，避免后续分析报错：

library(dplyr)
# 筛选出Type至少有2个水平的Roof
valid_roofs <- df %>%
  group_by(Roof) %>%
  summarize(n_types = n_distinct(Type)) %>%
  filter(n_types >= 2) %>%
  pull(Roof)
# 保留有效分组数据
df_filtered <- df %>% filter(Roof %in% valid_roofs)

3. 按Roof分组计算Tukey HSD并绘图 #

使用agricolae包计算Tukey HSD结果，再用ggsignif添加到分面箱线图中：

library(agricolae)
library(ggplot2)
library(ggsignif)

# 按Roof分组计算Tukey HSD对比结果
tukey_results <- df_filtered %>%
  group_by(Roof) %>%
  do({
    # 拟合ANOVA模型
    aov_mod <- aov(value ~ Type, data = .)
    # 计算Tukey HSD
    hsd_res <- HSD.test(aov_mod, "Type", group = FALSE)
    # 整理成ggsignif所需格式：指定对比组、p值、标签位置
    hsd_res$comparison %>%
      mutate(y_pos = max(.$value) * 1.08) # 标签位置可根据数据调整
  }) %>%
  ungroup()

# 绘制带p值的分组箱线图
ggplot(df_filtered, aes(x = Type, y = value)) +
  geom_boxplot(fill = "#619CFF", alpha = 0.7) +
  facet_wrap(~Roof, scales = "free_y") + # 分面展示不同Roof
  geom_signif(
    data = tukey_results,
    aes(xmin = group1, xmax = group2, annotations = paste0("p = ", round(p, 3)), y_position = y_pos),
    manual = TRUE, tip_length = 0.01
  ) +
  theme_minimal()

关键注意事项 #

• 确保每个Roof分组下的Type至少包含2个不同水平（如CONTROLE、soprema），否则该分组无法进行组间差异检验。
• 如果部分Roof分组的Type水平数不足，代码会自动过滤这些分组，避免报错。

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