在augPred图中添加abline

Question

很抱歉，这可能是一个非常基本的问题，我对R.

我希望读取我的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 (我的组)中放置特定时间段之间的线，并在这段时间内对值进行平均.

当然，任何其他方式/帮助或指导都将不胜感激！

提前感谢

Answer 1

我认为我们无法判断这是否“完全不正确”，因为你没有描述这个问题，也没有包括任何数据。(不过，绘图问题几乎是完全不正确的。)我可以告诉您，答案是不使用abline，因为augPred对象是用返回(并绘制)格子对象的plot.augPred绘制的。abline是一个基本的图形函数，不与点阵设备共享坐标系统。格型对象是可以修改的列表。您的绘图可能在不同级别的WSZ_Code上有不同的面板，但所需行的位置完全不清楚，特别是因为您使用省略号进行跟踪。你指的是“时间”，但没有“时间”变量。

有一些格函数，如trellis.focus和update.trellis，允许将修改应用于格对象。首先，将plot对象赋值给一个命名的变量，然后将mods分配给它，然后将它重新分配()。

help(package='lattice')
?Lattice

(如果这是一项匆忙的工作，您最好手工计算并使用ImageMagick编辑pdf或png输出。)