使用facet_grid绘制箱线图添加Tukey HSD结果仅显示旱季问题求助

核心原因 #

你遇到的问题大概率是Tukey检验未按「季节+站点」的分组维度拆分计算，或者在添加显著性标记时，没有将检验结果与每个子图的分组（季节）做对应绑定，导致只有旱季的结果被渲染。

解决步骤与代码示例 #

以下是基于ggplot2+emmeans+ggsignif的完整解决方案，以3站点×2季节的鱼类丰度数据为例：

1. 模拟/准备数据

确保数据包含站点、季节、丰度三个核心字段：

set.seed(123)
# 构造模拟数据
dat <- expand.grid(
  站点 = c("站点A", "站点B", "站点C"),
  季节 = c("旱季", "雨季"),
  重复 = 1:10
)
dat$丰度 <- rnorm(
  nrow(dat),
  mean = rep(c(5, 8, 6, 7, 9, 5), each=10), # 给不同站点-季节组合设置差异均值
  sd = 1.2
)

2. 按季节分组计算Tukey HSD检验

使用emmeans按季节拆分，计算每个季节内站点间的差异显著性：

library(emmeans)
library(multcomp)

# 拟合包含季节与站点交互项的方差分析模型
model <- aov(丰度 ~ 季节*站点, data = dat)

# 按季节分组，对站点间的丰度做Tukey多重比较
tukey_res <- emmeans(model, pairwise ~ 站点 | 季节)

# 将检验结果整理为数据框，方便后续绘图调用
tukey_df <- as.data.frame(tukey_res$contrasts)
# 提取每个比较对应的季节信息
tukey_df$季节 <- gsub(".*\\| ", "", tukey_df$contrast)
# 整理比较组名称（如"站点A vs 站点B"）
tukey_df$站点对比组 <- sapply(
  strsplit(as.character(tukey_df$contrast), " - "),
  function(x) paste(x, collapse = " vs ")
)

3. 绘制带显著性标记的分面子图

用ggplot2+ggsignif将检验结果匹配到对应子图：

library(ggplot2)
library(ggsignif)

ggplot(dat, aes(x = 站点, y = 丰度)) +
  geom_boxplot(fill = "#87CEEB", alpha = 0.7) +
  # 按季节分面，若要按站点行、季节列排列，可改为facet_grid(站点~季节)
  facet_grid(. ~ 季节) +
  # 为每个子图添加Tukey检验的显著性标记
  geom_signif(
    data = tukey_df,
    aes(
      xmin = as.numeric(factor(sapply(strsplit(as.character(contrast), " - "), "[[", 1))),
      xmax = as.numeric(factor(sapply(strsplit(as.character(contrast), " - "), "[[", 2))),
      annotations = ifelse(p.value < 0.05, paste("p=", round(p.value, 3), sep=""), "ns"),
      y_position = max(dat$丰度) + 0.5
    ),
    manual = TRUE,
    inherit.aes = FALSE
  ) +
  labs(x = "站点", y = "鱼类丰度") +
  theme_bw()

关键注意事项 #

• 必须在检验时通过| 季节指定分组维度，确保每个季节的站点间比较都被计算。
• 添加标记时，要保证tukey_df中的季节字段与分面的季节因子完全匹配（包括大小写、字符格式）。
• 若使用ggpubr的stat_compare_means，需通过group.by = "季节"指定分组，避免仅计算整体检验结果。

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