这真的很奇怪。假设我们有一个lmer
模型。生成诊断图的代码可以正常工作,但是当我将它放入for循环中时,它将停止工作!它只是产生一个空图:
require(lme4)
m <- lmer(y ~ (1|x), data = data.frame(y = c(1,1,1,1), x = c("a", "a", "b", "b")))
for (i in 1:2) {
# this code works only outside the for-loop!
png(filename = paste0(i, ".png"))
plot(m)
dev.off()
}
如果代码在for循环外运行,则代码通常会生成正确的PNG文件!这真是奇怪。问题出在哪里?