Displaying 1 result from an estimated 1 matches for "glmm3".
Did you mean:
glmm
2004 Nov 09
1
Some questions to GLMM
...ilar
results, but different values for AIC BIC and logLik.
I read in this thread:
http://maths.newcastle.edu.au/~rking/R/help/03b/6849.html
That it should be the same value due to the same algorithm
Maybe as additional comment, I specified a "NULL-model"
> dummy<-rep(1,469)
> glmm3<-GLMM(count~1,random=~1|dummy,poisson)
resp
> Dataset$dummy<-1
> glmm3<-glmmPQL(count~1,random=~1|dummy,poisson)
and GLMM gave the correct value (logLik = (Null deviance/2) from glm),whereas
glmmPQL calculated a logLik almost twice as high.
=========(MASS)=====
> glmm1<-glm...