nlme からグループ化されたデータをプロットしており、120 のパネルがあります。デフォルトのプロットplot(dataG)
では、2 行 60 列で画面いっぱいに表示されますが、読みにくいです。layout を指定するとplot(dataG), layout= c(12,10))
、適切な数の行と列が得られますが、列はすべて滑らかになります。
すべてがnlme内で起こっていることが問題なのかどうかはわかりませんが、nlmeの本で解決策を見つけていません.
ここでデータを見つけることができます
https://www.dropbox.com/s/79tssi252ai0ez8/COBS%20Roots%202008-2013noCNfake08.txt .
そして、グループ化を定義するコード:
roots<-read.table("COBS Roots 2008-2013noCNfake08.txt", header = TRUE)
library(nlme)
roots$EU<- with(roots, factor(plot):factor(depth))
rootsG<-groupedData(mass ~ year | EU, data=roots)
plot(rootsG, layout = c(12, 10))