Displaying 1 result from an estimated 1 matches for "lightintensity".
2009 Jul 10
1
problems with contrast matrix
...either with lmer (for PE) or with lme (for LA).
All models tested converge but the by-default contrast matrix is really hard
to interpret because all the contrasts always involve all the
wavelengths.
Here is an example of the syntax for lmer and lme, for the complete model:
model1<-lmer(PE~Sex+LightIntensity+ordered(WA)+(1|Animal),data=datamic,na.action=na.omit,family="poisson",REML=FALSE)
model2<-lme(LA~Sex+LightIntensity+ordered(WA),random=~1|Animal,data=datamic,na.action=na.omit,method="ML")
So I built a specific contrast matrix:
cmat.r<-matrix(c(-3/6,-2/6,-1/6,+1/6,+2/6...