[R] lme: random effects of a quadratic term

From: Montse Rue <montse.rue_at_gmail.com>
Date: Sat 26 Mar 2005 - 20:20:43 EST


Hello,

I am estimating the following model:

so2.lme<-lme(so2~1+I(alcadakm^2)+dia,data=subjectes2,na.action=na.omit)

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?

Thanks!

Montse Rue



R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Sat Mar 26 20:28:08 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:30:56 EST