R语言循环中调用lm函数使用动态变量名构建线性模型问题求助

lm函数接收公式对象作为第一个参数，直接写test_par[i]时，R会将其识别为从数据集all_data_plant中提取名为test_par的变量的第i个元素，而非将test_par[i]存储的字符串作为因变量名调用，因此运行报错。

方法1：使用reformulate构造公式（最推荐，后续模型调用无兼容问题） #

reformulate专门用于构造公式对象，正确用法如下：

# test_par为存储所有因变量名的字符向量，例如c("weight_dry", "height", "leaf_area")
current_form <- reformulate(termlabels = "weeks", response = test_par[i])
model <- lm(current_form, data = all_data_plant)

如果此前使用该函数未生效，基本为参数传递错误，按上述写法调整即可。

方法2：拼接字符串后转为公式 #

current_form <- as.formula(paste(test_par[i], "~ weeks"))
model <- lm(current_form, data = all_data_plant)

方法3：直接按列名索引数据集 #

不需要构造公式，直接从数据集提取对应因变量列传入：

model <- lm(all_data_plant[[test_par[i]]] ~ weeks, data = all_data_plant)

# 定义所有需要作为因变量的变量名
test_par <- c("weight_dry", "height", "leaf_num", "root_length")
# 定义空列表存储所有模型结果
model_list <- list()
for (i in seq_along(test_par)) {
  current_form <- reformulate("weeks", response = test_par[i])
  model_list[[i]] <- lm(current_form, data = all_data_plant)
  # 给模型命名，方便后续提取结果
  names(model_list)[i] <- test_par[i]
}
# 调用示例：提取weight_dry对应的模型结果
summary(model_list[["weight_dry"]])

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