R中结合for循环、paste与aov做批量方差分析出现contrasts报错求助

1. 调用aov()时直接传入了存储变量名的字符串，未从main数据框中提取对应列的实际值，aov()将字符串判定为单个常量，导致自变量被识别为仅1个水平，触发contrast报错
2. 赋值语句左侧直接使用paste()的返回值属于无效语法，无法生成对应名称的变量存储结果

基础修复（匹配你原本生成全局变量的需求） #

for(i in 1:length(DependentVariablesAnova)){
  for(j in 1:length(IndependentVariablesAnova)){
    # 生成合法的公式对象
    cur_formula <- as.formula(paste0(DependentVariablesAnova[i], " ~ ", IndependentVariablesAnova[j]))
    # 执行ANOVA检验
    cur_aov <- aov(cur_formula, data = main)
    # 动态生成变量名并赋值到全局环境
    res_name <- paste("Anova", DependentVariablesAnova[i], IndependentVariablesAnova[j], sep = "_")
    assign(res_name, cur_aov)
  }
}

优化方案（推荐） #

用列表统一存储所有检验结果，避免生成大量零散全局变量，更方便后续批量处理：

anova_res_list <- list()
for(i in 1:length(DependentVariablesAnova)){
  for(j in 1:length(IndependentVariablesAnova)){
    cur_formula <- as.formula(paste0(DependentVariablesAnova[i], " ~ ", IndependentVariablesAnova[j]))
    res_name <- paste0(DependentVariablesAnova[i], "_", IndependentVariablesAnova[j])
    anova_res_list[[res_name]] <- aov(cur_formula, data = main)
  }
}
# 调用示例：提取因变量为y、自变量为x的检验结果
# anova_res_list[["y_x"]]

如果修复代码后仍出现相同报错，说明对应自变量本身的有效值只有1个，可通过以下代码提前筛选符合要求的自变量：

# 计算每个自变量的去重有效值数量，筛选出≥2的变量
valid_ind_vars <- names(which(sapply(main[,IndependentVariablesAnova], function(col) length(unique(na.omit(col)))) >= 2))
# 后续循环直接用valid_ind_vars替代原IndependentVariablesAnova即可

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