使用lapply来拟合多个模型--如何保持lm对象中模型公式的自包含性

Question

以下代码适用于mtcars数据集的4种不同的模型公式，可用于循环或lapply。在这两种情况下，结果中存储的公式称为formulas[[1]]、formulas[[2]]等，而不是人类可读的公式。

formulas <- list(
  mpg ~ disp,
  mpg ~ I(1 / disp),
  mpg ~ disp + wt,
  mpg ~ I(1 / disp) + wt
)
res <- vector("list", length=length(formulas))
for (i in seq_along(formulas)) {
  res[[i]] <- lm(formulas[[i]], data=mtcars)
}
res
lapply(formulas, lm, data=mtcars)

是否有办法在结果中显示完整、可读的公式？

Answer 1

这应该能行

lapply(formulas, function(x, data) eval(bquote(lm(.(x),data))), data=mtcars)

它反覆了

[[1]]

Call:
lm(formula = mpg ~ disp, data = data)

Coefficients:
(Intercept)         disp  
   29.59985     -0.04122  


[[2]]

Call:
lm(formula = mpg ~ I(1/disp), data = data)

Coefficients:
(Intercept)    I(1/disp)  
      10.75      1557.67  

....etc

我们使用bquote将公式插入到对lm的调用中，然后计算表达式。

Answer 2

为什么不只是：

lapply( formulas, function(frm) lm( frm, data=mtcars))
#------------------
[[1]]

Call:
lm(formula = frm, data = mtcars)

Coefficients:
(Intercept)         disp  
   29.59985     -0.04122  

[[2]]

Call:
lm(formula = frm, data = mtcars)

Coefficients:
(Intercept)    I(1/disp)  
      10.75      1557.67  
snpped....

如果您希望结果的名称具有公式的“字符”化版本，那么它就是“

 names(res) <- as.character(formulas)
 res[1]
#-----
$`mpg ~ disp`

Call:
lm(formula = frm, data = mtcars)

Coefficients:
(Intercept)         disp  
   29.59985     -0.04122 

Answer 3

你也可以尝试这样的方法

library(purrr)
library(tibble)
models <- map(formulas, lm, data = mtcars)
models