如何在循环中为各abs分组关联Kruskal-Dunn检验结果至对应ggplot图

核心问题出在循环中未正确隔离分组数据与检验结果，导致跨分组的图层或变量污染。以下是针对性的修复代码和关键说明：

修复代码示例 #

# 加载必要包
library(tidyverse)
library(rstatix)
library(ggsignif)

# 获取abs列的唯一分组值
abs_groups <- unique(df.tidy.1$abs)

# 循环处理每个分组
for (abs_val in abs_groups) {
  # 1. 严格子集化当前分组数据，完全隔离其他组
  current_subset <- df.tidy.1 %>% filter(abs == abs_val)
  
  # 2. 针对当前子集执行Kruskal-Wallis检验
  kw_test <- current_subset %>% kruskal_test(your_value_col ~ your_group_col)
  
  # 3. 执行Dunn多重比较（替换为你需要的p值校正方法）
  dunn_test <- current_subset %>% 
    dunn_test(your_value_col ~ your_group_col, p.adjust.method = "bonferroni")
  
  # 4. 整理当前分组的显著性标记数据
  sig_marks <- dunn_test %>%
    filter(p.adjusted < 0.05) %>% # 筛选显著结果
    mutate(
      y_pos = max(current_subset$your_value_col) * 1.1, # 标记位置高于箱线图
      sig_label = paste("p =", round(p.adjusted, 4))
    )
  
  # 5. 生成独立的箱线图，仅使用当前分组的数据和标记
  plot <- ggplot(current_subset, aes(x = your_group_col, y = your_value_col)) +
    geom_boxplot(fill = "#4292c6") +
    ggtitle(paste("Abs:", abs_val)) +
    theme_minimal() +
    # 添加当前分组的显著性标记
    geom_signif(
      data = sig_marks,
      aes(xmin = group1, xmax = group2, annotations = sig_label),
      y_position = sig_marks$y_pos,
      tip_length = 0.02
    )
  
  # 6. 打印/保存当前图，避免图层累积
  print(plot)
  # ggsave(paste0("boxplot_abs_", abs_val, ".png"), plot, width = 7, height = 5)
}

关键修复点 #

• 数据隔离：每次循环内生成独立的current_subset，确保检验和绘图仅基于当前abs分组的 data
• 检验独立：Kruskal-Wallis和Dunn检验均绑定当前子集，不共享全局变量
• 图层独立：每次循环重新定义plot对象，避免在同一个ggplot实例上叠加不同分组的标记
• 标记专属：显著性标记的位置和内容完全基于当前分组的统计结果，不会携带上一次循环的残留数据

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