2014-02-24 57 views
0

我正在繪製一些來自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

您對在'lattice'包使用'xyplot'?我能夠生產出接近你想要的東西。 –

+0

我不反對它,特別是如果它可以使用分組數據。我在nlme之外沒有任何'晶格'的體驗。 – Nazer

+0

'lattice'對分組數據非常有用。我在下面添加了一個答案。 –

回答

0

我是能夠產生使用lattice包12由面板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)) 

enter image description here