1

次の例では、各モデルの残差プロットをファイルに書き込みます。ディスプレイでそれらを見る必要はありません。

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()
}

このコードをループなしで実行すると、必要なものが得られます。ただし、ループ内に空のファイルが作成されます。

何か案は？

ありがとうございました。

4

1 に答える 1

3

答えは FAQ にあります。実際にはFAQ 7.22です。ただし、これは、nlme パッケージの関数が実際のプロットを行うために格子/トレリス グラフィックスを使用していることに気付くまで明らかではありませんplot.lme(ヘルプ ページには の参照がありますがplot.lme、明らかではありません)。

解決策の短い形式 (ただし、問題を完全に理解するために FAQ と他のドキュメントを読むことをお勧めします) はplot、printコマンドで をラップすることです。

于 2014-06-19T19:59:10.640 に答える