R 如何在循环中绘制多个残差图?
在下面的示例中,我想将每个模型的残差图写入一个文件中。我不需要在显示器上看到它们R 如何在循环中绘制多个残差图?,r,plot,anova,mixed-models,R,Plot,Anova,Mixed Models,在下面的示例中,我想将每个模型的残差图写入一个文件中。我不需要在显示器上看到它们 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") pl
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()
}
for(1:500中的i){
temp.model事实上,答案在常见问题解答中。但是,直到您意识到nlme软件包中的plot.lme
函数使用晶格/网格图形来进行实际绘图(在帮助页面上有关于plot.lme
的参考,但不明显),这一点才明显
解决方案的简短形式(但我仍然建议阅读常见问题解答和其他文档以完全理解该问题)是将plot
包装在print
命令中。尝试执行print(plot(temp.model))
在循环内部。不确定绘图函数是如何工作的,但如果它使用网格图形,那应该会有帮助。它确实使用了晶格(网格)函数。是的,它有帮助。谢谢。