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

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 。 但是，直到您意识到nlme软件包中的plot.lme函数使用晶格/网格图形进行实际绘图（在帮助页面上有plot.lme ，但并不明显）之后， plot.lme显而易见。

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