R语言线性回归报错variable lengths differ，批量循环执行失败求助

错误原因 #

报错Error in model.frame.default : variable lengths differ的核心原因：

• 你构造id_n的方式错误：c(id_1, id_2,...)会把15列的所有元素拼接成一个超长向量（长度=15×13=195），循环时每次取出的是这个向量里的单个元素（长度1），而线性回归中自变量（如Vegetariana等）的长度是13，两者长度不匹配，导致报错。
• 单独执行linear_model(id_1)正常，是因为id_1是数据框中完整的一列（长度13），和自变量列长度一致。

修正步骤 #

1. 正确构造目标列名向量 #

不要直接拼接列，而是存储列名：

id_n <- paste0("id_", 1:15)  # 生成id_1到id_15的列名

2. 修改线性回归函数 #

让函数接受列名，并指定data参数，确保所有变量都来自同一个数据集：

linear_model <- function(y_col, data) {
  # 动态生成公式
  formula_str <- paste(y_col, "~ Vegetariana + Vegana + una_comida + dos_comidas + Low + Medium + Delivery")
  # 执行回归并返回摘要
  summary(lm(formula = formula_str, data = data))
}

3. 批量执行循环（或用更简洁的lapply） #

方法1：for循环

# 创建列表存储结果
results <- list()

for(col_name in id_n) {
  results[[col_name]] <- linear_model(col_name, data = data)
}

# 查看某个结果，比如id_1的回归摘要
results[["id_1"]]

方法2：lapply（更简洁）

results <- lapply(id_n, linear_model, data = data)
names(results) <- id_n  # 给列表元素命名，方便查看

额外说明 #

• 如果你的id_1等列不在data数据框中（而是全局环境的向量），需要调整函数逻辑，确保传入的因变量和自变量长度一致。但根据你提供的数据集结构，更合理的做法是将所有列放在同一个数据框中，避免全局环境变量混乱。
• 循环中直接打印a不会保存结果，建议用列表存储所有回归摘要，方便后续查看和分析。

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