将 abline 添加到 augPred 图

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

help(package='lattice')
?Lattice

（如果这是一项紧急工作，您最好手动进行任何计算并使用 ImageMagick 编辑 pdf 或 png 输出。）