添加abline到augPred情節

道歉什麼可能是一個很基本的問題，我是很新的R.添加abline到augPred情節

我期待以平均出的值提供給讀出我的augPred情節一段時間之間的預測。

> head(tthm.groupeddata) 
Grouped Data: TTHM ~ Yearmon | WSZ_Code 
    WSZ_Code Treatment_Code Year Month TTHM CL2_FREE BrO3 Colour PH TURB Yearmon 
1  2    3 1996  1 30.7  0.35 0.00030 0.75 7.4 0.055 Jan 1996 
2  6    1 1996  2 24.8  0.25 0.00055 0.75 6.9 0.200 Feb 1996 
3  7    4 1996  2 60.4  0.05 0.00055 0.75 7.1 0.055 Feb 1996 
4  7    4 1996  2 58.1  0.15  NA 0.75 7.5 0.055 Feb 1996 
5  7    4 1996  3 62.2  0.20  NA 2.00 7.6 0.055 Mar 1996 
6  5    2 1996  3 40.3  0.15 0.00140 2.00 7.7 0.055 Mar 1996

這是我的模型：

modellme<- lme(TTHM ~ Yearmon, random = ~ 1|WSZ_Code, data=tthm.groupeddata)

和我目前的劇情：

plot(augPred(modellme, order.groups=T),xlab="Date", ylab="TTHM concentration", main="TTHM Concentration with Time for all Water Supply Zones")

我想辦法通過向讀出圖要麼之間放置線在一個特定的WSZ_Code（我的組）中的特定時間段，並在此期間取平均值...

當然，任何其他方式/幫助或指導將不勝感激！

在此先感謝

來源

2013-05-04 Rbeginner

也 - 如果這是完全不正確，請不要讓我知道！ – Rbeginner 2013-05-04 14:06:09

我不認爲我們可以判斷它是否是「完全不正確」，因爲你還沒有描述的問題，並沒有包括任何數據。（雖然繪圖問題接近完全不正確）。我可以告訴你，答案是不使用abline，因爲augPred對象用plot.augPred繪製，它返回（並繪製）格子對象。 abline是一個基本的圖形功能，並不與網格設備共享一個座標系統。格對象是可以修改的列表。您的繪圖可能在WSZ_Code的不同層次上有不同的面板，但所需線條的位置完全不清楚，特別是因爲您使用省略號進行了跟蹤。你提到「時代」，但沒有「時代」變量。

有格子函數，如trellis.focus和update.trellis允許將修改應用於格子對象。您首先將繪圖對象分配給一個命名變量，對其進行修改，然後再繪製（）它。

help(package='lattice') 
?Lattice

（如果這是一個繁忙的工作，你可能會更好用手工做任何計算和使用ImageMagick編輯PDF或PNG輸出。）

來源

2013-05-04 15:48:29

好的，謝謝 - 看看。抱歉不清楚，我確實說過我對R很新！ – Rbeginner 2013-05-05 13:02:04

添加abline到augPred情節

回答

相關問題