保存多個線性回歸模型

Question

我模擬了一些數據，我想循環模擬 1000 次。 我每次都需要將數據擬合為線性模型並保存所有信息。

我的一位朋友告訴我，我可以每次提取標准誤差的系數並將它們保存起來，使它們成為矩陣或數據框。 但是，我想知道是否可以保存模型。 我想模擬 1000 次並保存所有模型。

但是我如何從 1000 次分析整體模型？ 就像模型是否合適一樣，我想使用 summary(model) 並查看每個變量是否顯着以及 f 統計量是否顯着？

Answer 1

您可以將所有模型存儲為列表，然后使用save函數將它們寫入文件：

lst <- lapply(1:1000, function(repetition) {
  mod <- lm(...)  # Replace this with the code you need to train your model
  return(mod)
})
save(lst, file="myfile.RData")

稍后，您可以使用load("myfile.Rdata")再次加載數據，然后使用lst[[1]] 、 lst[[2]] 、...