我正在绘制来自 nlme 的一些分组数据,并且我有 120 个面板。默认绘图plot(dataG)
将它们放在 2 行 60 列中,这会填满屏幕,但很难阅读。当我指定 layoutplot(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))