[R] plot.augPred with grouping factors

Dan Bebber
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