R-有没有一种方法可以输出和导入模型，而不是每次会议重新运行它们

Question

我必须运行许多用于工作的计算密集型模型（例如，在其他随机效果中嵌套了很多随机效果），从它们中进行预测，然后将其绘制为图。有时，我的经理希望我对预测或图表进行更改，但是这些模型有时需要几个小时才能运行。有什么办法可以通过导出和导入模型对象来节省时间，而不是每次都必须在脚本中重新运行它们？

Answer 1

counts <- c(18,17,15,20,10,20,25,13,12)
outcome <- gl(3,1,9)
treatment <- gl(3,3)
d.AD <- data.frame(treatment, outcome, counts)
glm.D93 <- glm(counts ~ outcome + treatment, family = poisson())

保存模型对象

saveRDS(glm.D93, file="glm.rds")
rm(glm.D93)

检索模型对象

glm.D93 <- readRDS("glm.rds")
anova(glm.D93)
# Analysis of Deviance Table
# Model: poisson, link: log
# Response: counts
# Terms added sequentially (first to last)
          # Df Deviance Resid. Df Resid. Dev
# NULL                          8    10.5814
# outcome    2   5.4523         6     5.1291
# treatment  2   0.0000         4     5.1291

Answer 2

您可以使用saveRDS(object, filename)保存任何对象，然后使用readRDS(filename)读回R。参见?saveRDS。