estimated growth curves from lme4

All,

The archive email below mentions the lack of general predict methods in
lme4.  I would like to plot the estimated growth curves for
treatment/placebo
groups (or the estimated patient growth curves that augment the former
w/ the patient estimated
random effects).

I obtained the plot manually by constructing the appropriate functions
from the fixed effects output and plotting using curve() (sample code
below).

For those that have transitioned from lme, do you suggest sticking w/
such an approach
(and possibly automating it w/ additional code, but this would be
dependent on the type
of model that is estimated each time), or going back to lme where I
recall that this was relatively
simple given the nesting mentioned by Douglas Bates below.

Thanks,
David

sample code:
## assume growth curve is quadratic, and coefficeints for placego group
are 10, 50, -20:
quad.fun.Placebo <- function(t) {
y = 100 + 50*t - 20*t^2
y}
## plot from t=0 to t=5
curve(quad.fun.Placebo, 0,5, xlab="Time", ylab="Dependent Variable",
col="red")

ps - why don't the methods show up below?

estimated growth curves from lme4

Thread (2 messages)