0

我正在绘制来自 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))
4

1 回答 1

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

在此处输入图像描述

于 2014-02-24T16:04:33.880 回答