[R-sig-ME] Using Observations as Random Effect in GLMM?

Daniel Hocking dhocking at unh.edu
Sat Jan 21 21:44:13 CET 2012


Hi everyone,

I am having trouble with overdispersion when trying to model count  
data using a GLMM. Beyond going to a negative binomial or Poisson- 
lognormal distribution, I have seen the suggestion (from Ben Bolker I  
believe) to include observation as a random effect. For example using  
the lme4 package my code would look something like this:

glmer(count ~ SoilT + SoilT2 + RH + rain24 + drought +
rain24*SoilT + drought*rain24 + (1 | plot) + (1 | obs), data = Data,
family = poisson)

When I try this I get a fitted vs. residual plot with large residuals  
at low fitted values funneling down to small residuals as the fitted  
values get larger. This indicates heterogeneity. I was wondering if  
that is expected for some reason with observation-level random effects  
or if this model just doesn't meet the assumptions of GLMM for my data?

Thanks,
Dan
------------------------------------------------------------------------------------
Daniel J. Hocking
122 James  Hall
Department of Natural Resources & the Environment
University of New Hampshire
Durham, NH 03824

dhocking at unh.edu
http://sites.google.com/site/danieljhocking/
http://quantitativeecology.blogspot.com/
http://richnessoflife.blogspot.com/

"Without data, you are just another person with an opinion."




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