[R] How to tell if there is overdispersion in a GLMM?
    maggy yan 
    kiotoqq at gmail.com
       
    Tue Feb  3 16:27:09 CET 2015
    
    
  
I read something on http://glmm.wikidot.com/faq, under "How can I deal with
overdispersion in GLMMs?":
library(lme4)  ## 1.0-4set.seed(101)
d <- data.frame(y=rpois(1000,lambda=3),x=runif(1000),
            f=factor(sample(1:10,size=1000,replace=TRUE)))
m1 <- glmer(y~x+(1|f),data=d,family=poisson)
overdisp_fun(m1)##        chisq        ratio          rdf            p
## 1026.7780815    1.0298677  997.0000000    0.2497659
library(glmmADMB)  ## 0.7.7
m2 <- glmmadmb(y~x+(1|f),data=d,family="poisson")
overdisp_fun(m2)##        chisq        ratio          rdf            p
## 1026.7585031    1.0298480  997.0000000    0.2499024
In both case, the chisq is > rdf, does it mean there is over dispersion?
thanks for any help
	[[alternative HTML version deleted]]
    
    
More information about the R-help
mailing list