The GENMOD Procedure

Zero-Inflated Models

Count data that have an incidence of zeros greater than expected for the underlying probability distribution of counts can be modeled with a zero-inflated distribution. In GENMOD, the underlying distribution can be either Poisson or negative binomial. See Lambert (1992), Long (1997) and Cameron and Trivedi (1998) for more information about zero-inflated models. The population is considered to consist of two types of individuals. The first type gives Poisson or negative binomial distributed counts, which might contain zeros. The second type always gives a zero count. Let $\lambda $ be the underlying distribution mean and $\omega $ be the probability of an individual being of the second type. The parameter $\omega $ is called here the zero-inflation probability, and is the probability of zero counts in excess of the frequency predicted by the underlying distribution. You can request that the zero inflation probability be displayed in an output data set with the PZERO keyword. The probability distribution of a zero-inflated Poisson random variable Y is given by

\[  \mr {Pr}(Y=y) = \left\{  \begin{array}{ll} \omega + (1-\omega )\mr {e}^{-\lambda } &  \mbox{for } y=0 \\ (1-\omega )\frac{\lambda ^ y \mr {e}^{-\lambda }}{y!} &  \mbox{for } y = 1,2,\ldots \\ \end{array} \right. \\  \]

and the probability distribution of a zero-inflated negative binomial random variable Y is given by

\[  \mr {Pr}(Y=y) = \left\{  \begin{array}{ll} \omega + (1-\omega )(1+k\lambda )^{-\frac{1}{k}} &  \mbox{for } y=0 \\ (1-\omega ) \frac{\Gamma (y+1/k)}{\Gamma (y+1)\Gamma (1/k)} \frac{(k\lambda )^ y}{(1+k\lambda )^{y+1/k}} &  \mbox{for } y = 1,2,\ldots \\ \end{array} \right. \\  \]

where k is the negative binomial dispersion parameter.

You can model the parameters $\omega $ and $\lambda $ in GENMOD with the regression models:

\begin{eqnarray*}  h(\omega _ i) &  = &  \mb {z}_ i^\prime \bgamma \\ g(\lambda _ i) &  = &  \mb {x}_ i^\prime \bbeta \\ \end{eqnarray*}

where h is one of the binary link functions: logit, probit, or complementary log-log. The link function h is the logit link by default, or the link function option specified in the ZEROMODEL statement. The link function g is the log link function by default, or the link function specified in the MODEL statement, for both the Poisson and the negative binomial. The covariates $\mb {z}_ i$ for observation i are determined by the model specified in the ZEROMODEL statement, and the covariates $\mb {x}_ i$ are determined by the model specified in the MODEL statement. The regression parameters $\bgamma $ and $\bbeta $ are estimated by maximum likelihood.

The mean and variance of Y for the zero-inflated Poisson are given by

\begin{eqnarray*}  \mr {E}(Y) &  = &  \mu = (1-\omega )\lambda \\ \mr {Var}(Y)&  = &  \mu + \frac{\omega }{1-\omega }\mu ^2 \\ \end{eqnarray*}

and for the zero-inflated negative binomial by

\begin{eqnarray*}  \mr {E}(Y) &  = &  \mu = (1-\omega )\lambda \\ \mr {Var}(Y)&  = &  \mu + \left(\frac{\omega }{1-\omega }+\frac{k}{1-\omega }\right)\mu ^2 \\ \end{eqnarray*}

You can request that the mean of Y be displayed for each observation in an output data set with the PRED keyword.