如何在R中用grep筛选变量后批量运行ANOVA并存储结果？

原代码的核心问题有两个：

• aov(i ~ treatment)里的i是字符串，无法被公式解析，必须动态构造符合要求的公式对象
• 循环过程中没有存储分析结果，所有输出仅临时打印，执行完后无留存

方案1：修复原生for循环 #

先构造存储结果的列表，再在循环中动态生成公式，把每个变量的模型、方差分析摘要、Tukey事后检验结果都存进去：

# 筛选目标变量
varlist <- grep("num_weeks_", names(crao2), value=TRUE)
# 创建空列表存储结果
anova_results <- vector("list", length(varlist))
names(anova_results) <- varlist

for (i in varlist) {
  # 动态构造公式
  formula <- as.formula(paste(i, "~ treatment"))
  # 拟合模型（注意原代码里的data参数：如果你的数据对象是df就替换crao2）
  model <- aov(formula, data = crao2)
  # 存储结果：每个元素是包含三个部分的子列表
  anova_results[[i]] <- list(
    model = model,
    summary = summary(model),
    tukey = TukeyHSD(model)
  )
}

# 查看某个变量的结果示例
anova_results[[varlist[1]]]$summary
anova_results[[varlist[1]]]$tukey

方案2：用lapply实现（基础R风格） #

lapply可以直接遍历varlist，自动返回结果列表，代码更简洁：

varlist <- grep("num_weeks_", names(crao2), value=TRUE)

anova_results_lapply <- lapply(varlist, function(var) {
  formula <- as.formula(paste(var, "~ treatment"))
  model <- aov(formula, data = crao2)
  list(
    model = model,
    summary = summary(model),
    tukey = TukeyHSD(model)
  )
})
# 给列表元素命名，方便后续索引
names(anova_results_lapply) <- varlist

方案3：用purrr::map实现（tidyverse风格） #

如果习惯tidyverse生态，map的写法更直观，支持简洁的公式风格语法：

library(purrr)

varlist <- grep("num_weeks_", names(crao2), value=TRUE)

anova_results_map <- map(varlist, ~{
  formula <- as.formula(paste(.x, "~ treatment"))
  model <- aov(formula, data = crao2)
  list(
    model = model,
    summary = summary(model),
    tukey = TukeyHSD(model)
  )
})
names(anova_results_map) <- varlist

补充说明 #

• 注意原代码中循环用的data = df，但变量是从crao2筛选的，确保数据对象一致，否则会出现变量找不到的错误
• 存储的列表每个元素对应一个目标变量的完整分析结果，后续可以直接提取需要的部分进行查看或进一步处理

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