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@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.
