Добавление abline к графику augPred

извините за то, что, вероятно, будет очень простым вопросом, я очень новичок в 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

This is my model: 

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 (моя группа), либо усреднив значения между этим периодом...

Конечно, любой другой способ/помощь или руководство будут высоко оценены!

заранее спасибо


r multi-level mixed-models
person Rbeginner    schedule 04.05.2013    source источник
comment
также - если это совершенно неверно, пожалуйста, дайте мне знать!   -  person Rbeginner    schedule 04.05.2013

Ответы (1)


arrow_upward
0
arrow_downward

Я не думаю, что мы можем сказать, является ли это «совершенно неверным», поскольку вы не описали вопрос и не включили никаких данных. (Однако вопрос построения графика близок к тому, чтобы быть совершенно неверным.) Я могу сказать вам, что ответ НЕ использовать abline, поскольку объекты augPred отображаются с помощью plot.augPred, который возвращает (и отображает) объект решетки. abline является базовой графической функцией и не использует общую систему координат с устройством решетки. Решетчатые объекты представляют собой списки, которые можно изменять. На вашем графике, вероятно, были разные панели на разных уровнях WSZ_Code, но расположение нужных строк совершенно неясно, тем более что вы заканчиваете многоточием. Вы ссылаетесь на «раз», но нет переменной «раз».

Существуют функции решетки, такие как trellis.focus и update.trellis, которые позволяют применять модификации к объектам решетки. Сначала вы должны назначить объект графика именованной переменной, сделать для него модификации, а затем снова использовать функцию plot().

help(package='lattice')
?Lattice

(Если это срочная работа, возможно, вам лучше выполнить какие-либо расчеты вручную и использовать ImageMagick для редактирования вывода в формате pdf или png.)

person IRTFM    schedule 04.05.2013
comment
Ок, спасибо - посмотрю. Извините за неясность, я сказал, что я очень новичок в R! - person Rbeginner; 05.05.2013