[R-sig-ME] REML estimates not given
Renwick, A. R.
a.renwick at abdn.ac.uk
Wed Sep 17 14:11:48 CEST 2008
I am running an lmer model using a poisson distribution and log link. However, when I ask REML=TRUE I do not get a REML estimate only the ML estimate.
mix3<-lmer(RoundedOverlap~Sex+Margin+sess+Sex:Margin+(1|Site), family=poisson, data=both,REML=TRUE)
mix3 at deviance
mix4<-lmer(RoundedOverlap~Sex+Margin+sess+Sex:Margin+(1|Site), family=poisson, data=both,REML=FALSE)
mix4 at deviance
I am using both the latest R (R version 2.7.2 ) and lme4 package
I wasw ondering if anybody knows why REML does not seem to be being used and if this is an issue with poisson distributed models in general.
School of Biological Sciences
University of Aberdeen
The University of Aberdeen is a charity registered in Scotland, No SC013683.
More information about the R-sig-mixed-models