[英]Adding abline to augPred plot
对于可能是一个非常基本的问题,我深表歉意,我对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(我的组)中的特定时间段之间放置线条并对该时间段之间的值取平均值来读取图表的方法...
当然,任何其他方式/帮助或指导将不胜感激!
提前致谢
我认为我们不能说出它是否“完全不正确”,因为您没有描述问题,也没有提供任何数据。 (作图问题是紧挨是完全不正确,虽然)。我可以告诉你,答案是不使用abline
,因为augPred对象与绘制plot.augPred
返回(和图形)格对象。 abline
是基本图形功能,不与点阵设备共享坐标系。 格对象是可以修改的列表。 您的绘图在WSZ_Code的不同级别上可能具有不同的面板,但是所需线条的位置完全不清楚,尤其是因为您使用省略号进行了尾随。 您引用“ times”,但是没有“ times”变量。
存在诸如trellis.focus
和update.trellis
类的晶格函数,这些函数允许将更改应用于晶格对象。 您将首先将绘图对象分配给一个命名变量,对其进行修改,然后再次对其进行plot()。
help(package='lattice')
?Lattice
(如果这是急事,那么最好手动进行任何计算并使用ImageMagick编辑pdf或png输出。)
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.