[R] lme: random effects of a quadratic term

Montse Rue montse.rue at gmail.com
Sat Mar 26 10:20:43 CET 2005


I am estimating the following model:


And when I try to plot the random effects of the quadratic term with
respect to a covariate (mam) I get an error:

> so2.lmeRE<-ranef(so2.lme,augFrame=T)
> plot(so2.lmeRE,form=I(alcadakm^2)~mam)
Error in plot.ranef.lme(so2.lmeRE, form = I(alcadakm^2) ~ mam ) : 
	Only single effects allowed in left side of form.

Any suggestion?


Montse Rue

More information about the R-help mailing list