如何在R Studio中循环对两个ERGM模型列表的对应模型执行anova()比较？报错问题求助

嗨，我来帮你搞定这个问题！你遇到的报错其实是个小疏忽导致的，咱们一步步来解决：

错误原因分析 #

你写的第一行代码anova_m1_m2 <- list有问题——这里你把变量赋值成了R的内置**list()函数**本身，而不是一个空的列表对象。当你后续尝试用[[i]]给它赋值时，R会报错，因为函数是不能被子集赋值的。

修正方案1：初始化空列表后循环 #

先创建一个和你的模型列表长度一致的空列表，再执行循环：

# 初始化空列表，长度匹配模型列表（不用硬写47，更灵活）
anova_m1_m2 <- vector("list", length = length(model1_list))

# 循环遍历每个位置的模型对
for(i in seq_along(model1_list)){
  m1 <- model1_list[[i]]
  m2 <- model2_list[[i]]
  anova_m1_m2[[i]] <- anova(m1, m2)
}

这里用seq_along(model1_list)代替1:47，能自动适配列表长度，万一后续模型数量变化也不用改代码。

更简洁的方案2：用purrr包的map2函数 #

如果你习惯用tidyverse风格的代码，purrr包的map2()函数可以一行搞定对应列表元素的操作，省去手动写循环的麻烦：

# 先加载purrr包（如果没装的话先运行install.packages("purrr")）
library(purrr)

# 批量执行ANOVA比较，结果直接存为列表
anova_m1_m2 <- map2(model1_list, model2_list, ~anova(.x, .y))

后续使用 #

两种方法都会生成一个包含47个ANOVA结果的列表，你可以通过anova_m1_m2[[1]]、anova_m1_m2[[2]]来查看对应位置模型的比较结果。

内容的提问来源于stack exchange，提问作者caritajohanna