R循环回归中如何为每个i匹配对应专属变量并加入回归模型

核心结论：完全可以实现，你之前尝试失败的核心原因是混淆了Stata的宏拼接逻辑和R的公式解析规则，整个实现完全可以复用你现有写的as.formula(paste())框架，不需要额外引入复杂包。

实现逻辑 #

你需要的专属变量命名规则是「前缀var4 + 当前循环i对应的取值」，本质是动态生成对应名称的字符串嵌入公式，和你之前拼接被解释变量、核心解释变量的逻辑完全一致：

• Stata中写var4[i]本质是宏自动拼接，把i的取值嵌到变量名里；
• R中要主动把这个拼接动作写在字符串生成环节，不能直接把var4[i]写进公式字符串——因为R公式里的[]会被解析为子集提取操作，不是字符串拼接，自然识别不到你要的专属变量。

可直接复用的代码框架 #

在你原有代码基础上修改，只需要加一步动态生成当前i对应的专属变量名即可：

# 提前定义基础向量
depvars <- c("dep1", "dep2")  # 存储所有被解释变量名的向量
indepvars <- c("ind1", "ind2")# 存储所有核心处理变量名的向量
regression <- list()          # 存储所有回归结果的空列表

# 双层循环遍历组合
for (i in depvars) {
  for (j in indepvars) {
    # 动态拼接当前i对应的专属变量名：var4前缀 + i的当前取值
    special_var_for_i <- paste0("var4", i)
    # 把专属变量名嵌入完整公式字符串
    full_formula <- paste(
      i, "~", j,
      "+ var1 + var2 + var3",
      "+", special_var_for_i,
      "+ factor(city) + factor(year)"
    )
    # 执行回归并存入列表
    regression[[paste0(i, "_", j)]] <- glm(
      as.formula(full_formula),
      data = 你对应的数据集名
    )
  }
}

常见踩坑提醒 #

• 不要直接在公式字符串里写var4[i]，这是最容易触发报错的点，R会把这段解析为「提取var4向量中i位置的元素」，和你要的变量名拼接逻辑完全无关；
• 如果你是按数值索引做循环（比如写法是for(i in 1:length(depvars))），拼接专属变量的时候不能直接用数值i，要先取对应位置的被解释变量字符串，比如current_depvar <- depvars[i]，再拼paste0("var4", current_depvar)，不然会拼出var41、var42这类和实际变量名不匹配的结果；
• 确保paste0("var4", 对应变量后缀)生成的字符串，和你数据集里实际存在的专属变量名100%匹配，否则会触发object 'xxx' not found的报错；
• 如果i的取值和专属变量的后缀不是完全对应关系，只要提前存一个命名映射向量，取的时候按i值匹配对应变量名即可，整体拼接逻辑不变。

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