[R-sig-ME] single argument anova for GLMMs (really, glmer, or dispersion?)

Andrew Robinson A.Robinson at ms.unimelb.edu.au
Sat Dec 13 22:14:16 CET 2008


On Sat, Dec 13, 2008 at 12:00:07PM -0500, Ben Bolker wrote:
> > 
> > I think that this is fair enough and well put, John, but I'm going to
> > push back in the other direction with a hypothetical example.  Let's
> > say that you have your over-dispersed count data.  What do you lose if
> > you simply take some convenient and credible transformation of the
> > response variable and then use lme, paying close attention to your
> > conditional distribution plots?  
> 
>   Besides the aesthetic preference for fully specified models etc.
> (although there's also the danger of forgetting that "all models
> are wrong etc." and believing the model too much), the most common
> reason in ecological contexts for not being able to get away with
> transformation is that the data are zero-rich (someone mentioned
> zero-inflated/hurdle models earlier in this discussion, which
> basically amounts to modeling presence/absence [either of
> "structural" zeros or of all zero values] and conditional
> density separately).  There's nothing you can do to transform
> a spike in the data (at zero or elsewhere) into anything
> other than a spike ...

Yes indeed :), you are correct.  But, that's a problem that bedevils
glms too.  I was more thinking along the lines: when would it be
inadvisable to use a transformation and explicit variance model
instead of a glm?

Cheers,

Andrew

-- 
Andrew Robinson  
Department of Mathematics and Statistics            Tel: +61-3-8344-6410
University of Melbourne, VIC 3010 Australia         Fax: +61-3-8344-4599
http://www.ms.unimelb.edu.au/~andrewpr
http://blogs.mbs.edu/fishing-in-the-bay/




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