r - 使用 lapply 拟合多个模型——如何在 lm 对象中保持模型公式自包含

Question

以下代码将 4 个不同的模型公式拟合到 mtcars数据集，使用 for 循环或 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 的调用中然后计算表达式。