我模拟了一些数据,我想循环模拟1000次。我需要每次都将数据拟合到线性模型并保存所有信息。
我的一位朋友告诉我,我每次都可以提取标准误差系数并保存它们,使它们成为矩阵或数据帧。但是,我想知道我是否可以保存模型。我想模拟1000次并保存所有模型。
但我如何从1000次分析整体模型? 如果模型适合,我想使用摘要(模型) 并查看每个变量是否显着以及f统计量是否显着?
答案 0 :(得分: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]]
,...