Is there such a package that provides for zero-inflated negative binomial mixed-effects model estimation in R?

By that I mean:

  • Zero-inflation where you can specify the binomial model for zero inflation, like in function zeroinfl in package pscl:

    zeroinfl(y~X|Z, dist = "negbin")
    where Z is the formula for the zero inflation model;

  • Negative binomial distribution for the count part of the model;

  • Random effects specified similar to function lmer of package lme4.

I understand glmmADMB can do all that, except the formula for zero inflation cannot be specified (it is just an intercept, i.e. Z is just 1). But are there any other packages that can do it all?

I will be very thankful for your help!

  • $\begingroup$ I'm looking for the same thing. I think glmmADMB is what I want, but I can't get it running. $\endgroup$ – gregmacfarlane Dec 7 '12 at 20:41
  • $\begingroup$ I'm wondering whether or not the ZIM or aod packages can do what you want to do? $\endgroup$ – Graeme Walsh Jun 26 '13 at 10:28
  • 2
    $\begingroup$ As an update, the glmmTMB package by Ben Bolker supports a zero-inflated generalized linear mixed model (ZIGLMM). $\endgroup$ – JWilliman Feb 22 '18 at 21:55

I think this is the package you need: glmmADMB. I downloaded it here: http://otter-rsch.com/admbre/examples/glmmadmb/glmmADMB.html

But I still had some problems to get it to run, so I followed the instructions provided in this link and now it works fine http://glmmadmb.r-forge.r-project.org/

Hope this helps!

  • $\begingroup$ Note that this package only allows the fitting of a constant term for the zero-inflation part of the model. $\endgroup$ – user2390246 Apr 19 '18 at 16:19

The pscl package provides for a zero inflated Poisson model. I don't think that it can do a negative binomial model, but it might be a place to start. The linked JSS article also discusses related packages, which may lead you to what you're looking for.

  • 1
    $\begingroup$ The pscl package does (now) allow negative binomial models by using zeroinfl(..., dist = "negbin", ...) $\endgroup$ – user2390246 Apr 19 '18 at 16:20

Depending on what you're trying to do, you might want to look at the aster package. Aster models allow joint analysis of multiple variables that have different probability distributions, and recently have been updated to allow for random effects. They were designed for life history analysis and will work in situations where you can split your response into distinct parts with different distributions, (e.g. survival = Bernoulli, reproduction = Poisson). They can handle "zero-inflation" by modeling the majority of zeroes as bernoulli, and the remainder of the response as negative binomial.

You'll find plenty of documentation here:



Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.