添加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情节

回答

相关问题