R使用循环批量生成多变量对的单因素ANOVA、统计摘要及绘图

核心修正说明 #

• 原伪代码的错误根源是直接传入变量对象而非变量名字符串，且没有适配dplyr、ggplot的非标准求值规则，本方案用.data代词即可无额外依赖解决该问题
• 变量对统一用字符串格式存储，后续新增分析配对只需在列表中新增一行即可，无需修改核心逻辑

完整可运行代码 #

# 提前加载所需包，不要放在循环内重复加载
library(dplyr)
library(ggplot2)

# 定义变量对列表：每个元素格式为 c("自变量列名", "因变量列名")
var_pairs <- list(
  c("Gender", "Attr_scores"),
  c("Gender", "Fear_scores"), 
  c("Type", "Appr_scores"), 
  c("Emotion", "Avoid_scores")
)

# 可选：初始化列表存储所有结果，方便后续查看导出
all_desc <- list()
all_aov <- list()

# 循环处理每一组变量对
for (pair in var_pairs) {
  iv <- pair[1] # 提取自变量名
  dv <- pair[2] # 提取因变量名
  
  # 1. 分组描述统计
  desc_res <- group_by(df1, .data[[iv]]) %>%
    summarise(
      N = n(),
      Mean = mean(.data[[dv]], na.rm = TRUE),
      Sd = sd(.data[[dv]], na.rm = TRUE)
    )
  # 打印结果，也可写入csv保存
  cat("===== 分组描述统计：", iv, "→", dv, "=====\n")
  print(desc_res)
  all_desc[[paste0(iv, "_", dv)]] <- desc_res
  
  # 2. 单因素ANOVA
  aov_formula <- as.formula(paste0(dv, " ~ ", iv))
  aov_res <- aov(aov_formula, data = df1)
  cat("\n===== ANOVA结果：", iv, "→", dv, "=====\n")
  print(summary(aov_res))
  all_aov[[paste0(iv, "_", dv)]] <- aov_res
  
  # 3. 绘制带误差棒的抖动散点图并保存
  p <- ggplot(df1, aes(x = .data[[iv]], y = .data[[dv]])) + 
    geom_jitter(position = position_jitter(0.2), alpha = 0.7) + # 加透明度优化重叠点显示
    stat_summary(fun.data = mean_sdl, fun.args = list(mult = 1), 
                 geom = "pointrange", color = "red", linewidth = 1, size = 1.2) +
    labs(x = iv, y = dv) + # 自动设置坐标轴标签
    theme_bw(base_size = 14) # 可选优化主题
  # 自动生成文件名
  file_name <- paste0(iv, "_", dv, "_plot.png")
  ggsave(file_name, p, width = 1600, height = 900, units = "px", dpi = 150)
  cat("\n图片已保存为：", file_name, "\n\n")
}

额外使用说明 #

• 如需新增分析配对，仅需在var_pairs列表中按照格式添加新的c("自变量名", "因变量名")条目即可
• 存储在all_desc和all_aov中的结果可以按需导出，例如导出某个描述统计结果为csv：write.csv(all_desc[["Gender_Attr_scores"]], "gender_attr_desc.csv", row.names = FALSE)

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