[R-sig-ME] Overdispersion lme4 binomial

Jarrod Hadfield j.hadfield at ed.ac.uk
Sun Aug 1 17:11:06 CEST 2010


Dear Chris,

Over-dispersion does not occur with a binary response variable so you  
don't need to test for it.

This does not mean that between-datum heterogeneity in the probability  
of success is absent, only that it cannot be observed. For example,  
take 1000 random draws from a binomial distribution with constant  
probability (0.5):

table(rbinom(1000, 1, 0.5))

and compare the frequency of outcomes with a 1000 draws from 1000  
binomial distributions with different probabilities of success (but  
with mean = 0.5)

table(rbinom(1000, 1, runif(1000)))

The data look the same, and so the between-datum heterogeneity  
(residual variance if you like) although it may exist cannot be  
estimated from the data.

Cheers,

Jarrod


Quoting Chris Mcowen <cm744 at st-andrews.ac.uk>:

> Dear List,
>
> I am wanting to test for overdispersion in my model and am unsure   
> how for my specific case.
>
> I have 2 random factors, 7 fixed factors that have multiple levels   
> and are categorical and then i have a binary response (True or False).
>
> model1 <- lmer(threattf~1+(1|order/family) + geophyte + seasonality   
> + pollendispersal + breedingsystem*fruit + habit + lifeform +   
> woodyness, family=binomial)
>
> I would be very grateful if somebody could point me in the right   
> direction for testing for overdispersion under such scenarios?
>
> Please see part of the output below -
>
> Thanks for any help, and if more data is required feel free to ask.
>
> Chris
>
> Generalized linear mixed model fit by the Laplace approximation
> Formula: threattf ~ 1 + (1 | order/family) + geophyte + seasonality   
> +      pollendispersal + breedingsystem * fruit + habit + lifeform +  
>       woodyness
>   AIC  BIC logLik deviance
>  1562 1649 -764.2     1528
> Random effects:
>  Groups       Name        Variance Std.Dev.
>  family:order (Intercept) 0.26932  0.51896
>  order        (Intercept) 0.00000  0.00000
> Number of obs: 1242, groups: family:order, 43; order, 9
>
> Fixed effects:
>                        Estimate Std. Error z value Pr(>|z|)
> (Intercept)            -0.10413    0.98004  -0.106  0.91538
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>
>



-- 
The University of Edinburgh is a charitable body, registered in
Scotland, with registration number SC005336.




More information about the R-sig-mixed-models mailing list