4

假设我有一些模型存储在列表中：

mods <- list()
mods[[1]] <- lm(mpg ~ disp, data = mtcars)
mods[[2]] <- lm(mpg ~ disp + factor(cyl), data = mtcars)
mods[[3]] <- lm(mpg ~ disp * factor(cyl), data = mtcars)

我想用stats::AIC. 我正在寻找我会从中得到的输出AIC(mods[[1]], mods[[2]], mods[[3]])，但我希望它能够概括为任意长的列表。我以为

do.call(AIC, mods)

会起作用，但它会返回一些非常冗长且无益的东西。（如果列表已命名，则会导致错误，除非其中一个名称是object，对应于 的第一个参数AIC，但随后您将再次获得详细输出。）

失败后do.call，我开始思考eval(parse())解决方案，但我想我应该先在这里问。

4

2 回答 2

4 
summary(do.call(AIC, mods))
       df         AIC       
 Min.   :3   Min.   :153.4  
 1st Qu.:4   1st Qu.:159.6  
 Median :5   Median :165.8  
 Mean   :5   Mean   :163.1  
 3rd Qu.:6   3rd Qu.:168.0  
 Max.   :7   Max.   :170.2

但这可能不是你想要的。巴蒂斯特给出了答案：

my.aic <- function(x) {
  x <- do.call(AIC, x)
  rownames(x) <- NULL
  return(x)
}
my.aic(mods)
##   df      AIC
## 1  3 170.2094
## 2  5 165.7680
## 3  7 153.4352

这相当接近：

AIC(mods[[1]], mods[[2]], mods[[3]])
##           df      AIC
## mods[[1]]  3 170.2094
## mods[[2]]  5 165.7680
## mods[[3]]  7 153.4352
于 2013-05-19T02:10:30.027 回答
2

这是一个使用eval(parse()). 幸运的是，马修在我输入之前给出了更好的答案。

AIC_l <- function(L, FUN = "AIC") {
    args <- paste0("mods[[", seq_along(L), "]]", collapse = ", ")
    my_call <- paste0(FUN, "(", args, ")")
    eval(parse(text = my_call))
}
AIC_l(L = 1:3, FUN = "AIC")
于 2013-05-19T02:14:00.853 回答