Question

我有一些在循环中运行模型的代码。循环的每次迭代运行稍微不同的模型，结果存储在变量中。存储这些对象的好方法是什么，以便在循环终止后可以访问它们？我想到了这样的事情：

fit.list <- list(n)
for (i in 1:n) {
    fit <- glm(......)
    fit.list[i] <- fit
}

但后来我想访问每个模型结果，例如summary(fit.list[4])或plot(fit.list[15])，但这似乎不起作用。

Answer 1

尝试

plot(fit.list[[15]])

单个[函数提取包含所请求组件的列表，即使该列表长度为1。

双[[函数提取单个声明的组件并返回它但不在列表中;即你得到的组件本身不是包含该组件的列表。

以下是插图：

> mylist <- list(a = 1, b = "A", c = data.frame(X = 1:5, Y = 6:10))
> str(mylist)
List of 3
 $ a: num 1
 $ b: chr "A"
 $ c:'data.frame':  5 obs. of  2 variables:
  ..$ X: int [1:5] 1 2 3 4 5
  ..$ Y: int [1:5] 6 7 8 9 10
> str(mylist["c"])
List of 1
 $ c:'data.frame':  5 obs. of  2 variables:
  ..$ X: int [1:5] 1 2 3 4 5
  ..$ Y: int [1:5] 6 7 8 9 10
> str(mylist[["c"]])
'data.frame':   5 obs. of  2 variables:
 $ X: int  1 2 3 4 5
 $ Y: int  6 7 8 9 10

注意最后两个命令输出的区别。 str(mylist["c"])说“List of 1”而str(mylist[["c"]])说“'data.frame':”。

使用plot(fit.list[15])，您要求R绘制列表对象，而不是列表中该元素中包含的模型。

Answer 2

也许可以试试

fit.list <- list()
for (i in 1:5) {
counts <- c(18,17,15,20,10,20,25,13,12)
outcome <- gl(3,1,9)
treatment <- gl(3,3)
print(d.AD <- data.frame(treatment, outcome, counts))
glm.D93 <- glm(counts ~ outcome + treatment, family=poisson())
fit.list[[i]] <-glm.D93
}

请注意fit.list[[i]]而不是fit.list[i]

在R中存储和访问模型结果

2 个答案: