[R] negative binomial lmer

[Gregor Gorjanc]

Ben Bolker <bolker <at> ufl.edu> writes:

...
>  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:
> 
> www.cefe.cnrs.fr/esp/TBElston_Parasitology2001.pdf

Actually it is very simple

lmer(y ~ effA + (1 | effB), family=quasipoisson)

i.e. this fits the following model

y_ijk ~ Poisson(\lambda_ijk)
log(lambda_ijk) = \mu + effaA_i + effB_ij + e_ijk
effB_i ~ Normal(0, \sigma^2_b)
e_ijk ~ Normal(0, \sigma^2_e)

Gregor