Cross Validated is a question and answer site for people interested in statistics, machine learning, data analysis, data mining, and data visualization. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

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!

share|improve this question
I'm looking for the same thing. I think glmmADMB is what I want, but I can't get it running. – gregmacfarlane Dec 7 '12 at 20:41
I'm wondering whether or not the ZIM or aod packages can do what you want to do? – Graeme Walsh Jun 26 '13 at 10:28

I think this is the package you need: glmmADMB. I downloaded it here:

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

Hope this helps!

share|improve this answer

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.

share|improve this answer

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:

share|improve this answer

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

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