r - 如何在循环中绘制多个残差图？

Question

在以下示例中，我想将每个模型的残差图写入文件中。我不需要在我的显示器中看到它们。

for (i in 1:500){
     temp.model<-lme(as.formula(paste("Var",i) ~ X1*X2, sep=""), data = example, random=~1| Exp/Person)
     jpeg(paste("C:/Myfolder", i, ".jpg", sep = ""), quality=50, bg="white")
     plot(temp.model)
     dev.off ()
     graphics.off()
}

当我在没有循环的情况下运行这段代码时，我得到了我想要的。但是，它会在循环中创建空白文件。

有任何想法吗？

谢谢你。

Answer 1

答案在FAQ中，实际上是FAQ 7.22 。然而，在您意识到plot.lmenlme 包中的函数使用 lattice/trellis 图形来进行实际绘图之前，这并不明显（在帮助页面上有参考plot.lme，但并不明显）。

解决方案的简短形式（但我仍然建议阅读常见问题解答和其他文档以完全理解该问题）是将其包装plot在print命令中。