[R] plot.augPred with grouping factors

From: Dan Bebber <danbebber_at_forestecology.co.uk>
Date: Sun, 30 Dec 2007 10:14:27 -0000


I fitted the following model:
dgr.lme2 <- lme(Average ~ poly(Time, 2) * Treatment, random = ~1|Replicate, data = dgr, correlation = corAR1(), method = "ML")

I plotted the data in separate panels for each Treatment: xyplot(Average ~ Time|Treatment, groups = Replicate, data = dgr, type = c("p","l"), ylab = expression(paste(Delta, G)))

I would now like to overlay the population-level predictions from dgr.lme2 onto each Treatment-level panel.

How can I tell plot.augPred to do this? Or do I need a bespoke panel function?
(plot.augPred gives panels for each Replicate by default)

Many thanks,
Dan Bebber

Checked by AVG Free Edition.

11:51



R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Mon 31 Dec 2007 - 08:49:21 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Mon 31 Dec 2007 - 09:30:21 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.