R语言如何访问动态命名创建的线性模型lm的系数？

你原来的代码报错是因为paste0("lm_", q[[i]])返回的是字符串，而coef()需要接收线性模型对象本身，可通过以下两种方式解决：

方法1：直接使用get()函数获取动态变量对象 #

get()可以通过变量名字符串获取对应环境中的变量实体，修改后的可用代码如下：

q = c(0.01, 0.02)
for(i in seq_along(q)){
  assign(paste0("lm_", q[[i]]), lm(mtcars$mpg ~ mtcars$disp)) 
}

# 正确访问系数
for(i in seq_along(q)){
  # 先通过get获取对应的lm模型对象
  current_lm = get(paste0("lm_", q[[i]]))
  intercept = coef(current_lm)[[1]]
  slope = coef(current_lm)[[2]]
  cat("Write the coefficients to some file...")
}

方法2：更推荐的最佳实践：用命名列表存储所有模型 #

动态生成独立变量的方式维护成本高，更建议把所有模型存在同一个命名列表中，访问更方便，不需要操作全局环境变量：

q = c(0.01, 0.02)
# 初始化空列表存储模型
lm_list = list()
for(i in seq_along(q)){
  model_name = paste0("lm_", q[[i]])
  lm_list[[model_name]] = lm(mtcars$mpg ~ mtcars$disp)
}

# 访问系数
for(model_name in names(lm_list)){
  intercept = coef(lm_list[[model_name]])[[1]]
  slope = coef(lm_list[[model_name]])[[2]]
  cat("Write the coefficients to some file...")
}

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