Re: [R] negative binomial lmer

From: Ben Bolker <>
Date: Sat 29 Jul 2006 - 00:52:27 EST

Tracy Feldman <tracysfeldman <at>> writes:

> To whom it may concern:
> I have a question about how to appropriately conduct an lmer analysis for
negative binomially distributed
> data. I am using R 2.2.1 on a windows machine.
> I am trying to conduct an analysis using lmer (for non-normally distributed
data and both random and fixed
> effects) for negative binomially distributed data. To do this, I have been
using maximum likelihood,
> comparing the full model to reduced models (containing all but one effect, for
all effects). However, for
> negative binomially distributed data, I need to estimate the parameter theta.
 I have been doing this by
> using a negative binomial glm of the same model (except that all the effects
are fixed), and estimating mu
> as the fitted model like so:

 I haven't tried it, but you could also consider using a Poisson-lognormal (rather than neg binomial, which is Poisson-gamma) distribution, which might make this all work rather well in lmer:

    Ben Bolker mailing list PLEASE do read the posting guide and provide commented, minimal, self-contained, reproducible code. Received on Sat Jul 29 01:05:27 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Sat 29 Jul 2006 - 08:17:20 EST.

Mailing list information is available at Please read the posting guide before posting to the list.