我正在绘制来自nlme的一些分组数据,我有120个面板。默认绘图plot(dataG)
将它们分为2行,60列,填满了屏幕,但是难以阅读。当我指定布局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))
答案 0 :(得分:0)
我能够使用lattice
包生成12 x 10的面板图。希望它可以让你朝着正确的方向前进。
> roots <-read.table("COBS Roots 2008-2013noCNfake08.txt", header = TRUE)
> roots$EU <- with(roots, factor(plot):factor(depth))
> library(lattice)
> xyplot(mass ~ year | EU, data = roots, layout = c(12,10))