【发布时间】:2014-03-26 09:03:33
【问题描述】:
我正在绘制来自 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))
【问题讨论】:
-
您反对在
lattice包中使用xyplot吗?我能够生产出接近您想要的东西。 -
我不反对,尤其是如果它可以处理分组数据。我只是在
nlme之外没有任何lattice经验。 -
lattice非常适合分组数据。我在下面添加了一个答案。