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等格子函数允许对格子对象应用修改。您首先将绘图对象分配给命名变量，为其创建mod，然后再次绘制（）它。

help(package='lattice')
?Lattice

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

将abline添加到augPred图

1 个答案: