使用 lapply 擬合多個模型——如何在 lm 對象中保持模型公式自包含

Question

以下代碼使用 for 循環或 lapply 將 4 個不同的模型公式擬合到mtcars數據集。 在這兩種情況下，存儲在結果中的formulas[[1]]被稱為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