R语言嵌套循环生成唯一组合及执行单因素ANOVA问题求助

原有代码问题排查 #

• 循环内未按当前迭代的Feature、ID筛选数据集，每次都对全量数据执行ANOVA，结果重复无效
• 列表存储仅用ID作为索引，不同Feature的计算结果会被覆盖
• 直接存入ANOVA的summary对象，未提取结构化的统计量，无法直接合并为数据框

修正后嵌套循环实现 #

# 提前提取唯一值，避免循环内重复计算
unique_feat <- unique(data1$Feature)
unique_id <- unique(data1$ID)
# 初始化结果列表
res_list <- list()
idx <- 1

for (feat in unique_feat) {
  for (id in unique_id) {
    # 筛选当前分组的子集
    cur_df <- data1[data1$Feature == feat & data1$ID == id, ]
    # 跳过样本量不足的分组，避免运行报错
    if (nrow(cur_df) < 2 || length(unique(cur_df$Sub)) < 2) next
    # 执行单因素ANOVA
    aov_fit <- summary(aov(Value ~ as.factor(Sub), data = cur_df))
    # 提取组间因子的统计量（自由度、平方和、均方、F值、P值）
    group_res <- aov_fit[[1]][1, ]
    # 补充分组标识
    group_res$Feature <- feat
    group_res$ID <- id
    # 存入结果列表
    res_list[[idx]] <- group_res
    idx <- idx + 1
    
    # 如需单独输出每个分组的结果，取消注释下行即可
    # write.csv(group_res, paste0("ANOVA_", feat, "_", id, ".csv"), row.names = F)
  }
}

# 合并为统一结果表
anova_res <- do.call(rbind, res_list)
# 查看输出
print(anova_res)

结果使用说明 #

输出的anova_res已经包含所有你需要的组间因子统计参数，可直接用于图表生成，以下是F值柱状图的示例代码：

library(ggplot2)
ggplot(anova_res, aes(x = paste0(Feature, "-", ID), y = `F value`, fill = Feature)) +
  geom_col() +
  labs(x = "Feature-ID组合", y = "ANOVA组间F值") +
  theme(axis.text.x = element_text(angle = 45, hjust = 1))

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