You want an 'augPred' plot with multiple lines per panel, from
multiple factors in the model? Because you provided such a simple,
self-contained example, I can offer much more substantive comments than
I might have otherwise. First the bad news: I don't see an easy way to
get what you
All,
I have a question RE plotting the prediction lines of a random effects
model via augPred. I'll illustrate via the Pixel dataset within the
nlme package:
library(nlme)
attach(Pixel)
fm1Pixel = lme(pixel ~ day + I(day^2), data = Pixel, random = list(Dog =
~ 1))
plot(augPred(fm1Pixel)) ###