r - 从列表创建和调用线性模型

Question

所以我试图比较不同的线性模型，以确定一个模型是否比另一个更好。但是我有几个模型，所以我想创建一个模型列表，然后调用它们。那可能吗？

 Models <- list(lm(y~a),lm(y~b),lm(y~c)
 Models2 <- list(lm(y~a+b),lm(y~a+c),lm(y~b+c))

 anova(Models2[1],Models[1])

感谢您的帮助!

Answer 1

如果您有两个模型列表，并且想要比较每对模型，那么您需要 Map :

models1 <- list(lm(y ~ a), lm(y ~ b), lm(y ~ c)
models2 <- list(lm(y ~ a + b), lm(y ~ a + c), lm(y ~ b + c))

Map(anova, models1, models2)

这基本上等同于以下 for 循环:

out <- vector("list", length(models1))
for (i in seq_along(out) {
  out[[i]] <- anova(models1[[i]], models2[[i]])
}

Map 是函数的一个例子，你可以在 https://github.com/hadley/devtools/wiki/Functionals 找到更多关于它们的信息。