如何使用 anova() 循环比较两个模型列表答案

【问题标题】：How to loop over comparison of two model lists with anova()如何使用 anova() 循环比较两个模型列表
【发布时间】：2022-07-19 22:41:57
【问题描述】：

我是 R Studio 的新手，不知道如何制作一个可重现的示例，但我希望你能掌握它：

我有两个列表，model1_list 和 model2_list，它们都包含 47 个 ERGM 模型输出。我想像这样将每个模型与另一个列表中的对应模型进行比较

anova(model1_list[[1]], model2_list[[1]])

但是，我无法想出一个解决方案来一次循环所有 anova() 测试。我试过了

anova_m1_m2 <- list
for(i in c(1:47)){
  m1 <- model1_list[[i]]
  m2 <- model2_list[[i]]
  a_m1_m2 <- anova(m1, m2)
  anova_m1_m2[[(i)]] <- a_m1_m2
}

但是得到错误信息 Error in anova_m1_m2[[(i)]] <- a_m1_m2 : object of type 'builtin' is not subsettable

有什么想法吗？

【问题讨论】：

标签： r loops anova

【解决方案1】：

修复它

anova_m1_m2 <- list()

【讨论】：

更好的做法是anova_m1_m2 <- vector(mode = "list", length = length(model1_list))。请参阅here 了解原因。