Displaying 2 results from an estimated 2 matches for "lmeres".
Did you mean:
limeres
2005 Mar 26
1
lme: random effects of a quadratic term
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
2011 Jan 17
1
Using anova() with glmmPQL()
...quot;))
mnth=factor(c(rep("Nov",times=40),rep("Jan",times=40)),levels=c("Nov","Jan"))
subjectnum=c(1:10,1:10,1:10,1:10,1:10,1:10,1:10,1:10)
subject=factor(paste(roostsitu,subjectnum,sep=""))
dataf=data.frame(mnth,roostsitu,subjectnum,subject,stmass)
lmeres=lme(fixed=stmass~mnth*roostsitu,random=~1|subject/mnth,na.action=na.exclude)
anova(object=lmeres,test="Chisq")
numDF denDF F-value p-value
(Intercept) 1 36 31143.552 <.0001
mnth 1 36 95.458 <.0001
roostsitu 3 36 10.614...