Providing software solutions since 1976

• Log in
• Create profile
• Welcome
  [Log out]
• Edit profile
• My SAS

• Print  |
• Feedback  |
• Share

  • E-mail
  

  E-mail This Page

  Your Name:*
  Your E-mail:*
  Recipient E-mail:*
  Send a copy to me
  Personal Message:
  
  

  * indicates a required field

The example titled "Modeling Zero-Inflation: Is it Better to Fish Poorly or Not to Have Fished At All?" in the FMM procedure documentation discusses zero-inflated and hurdle models for modeling count data containing excessive zeros. As noted there, the hurdle model supposes two processes at work — one that generates zeros with some probability, and the other that generates events. For instance, the Poisson hurdle model is a mixture of a degenerate distribution at zero and a truncated Poisson distribution. The zero-inflated Poisson (ZIP) model also uses a degenerate zero distribution, but the second process is a regular Poisson distribution which can generate both zeros and events. For the ZIP model, the first process therefore generates only extra zeros beyond those of the regular Poisson distribution. For the hurdle model, the first process generates all of the zeros.

Truncated Poisson and negative binomial distributions are discussed and illustrated in more detail in this note.

The example in the FMM documentation illustrates the ZIP model. The Poisson hurdle model is just as easily fit, but uses the DIST=TRUNCPOISSON option instead of the DIST=POISSON option.

      proc fmm data=catch;
         class gender;
         model count = gender*age / dist=TruncPoisson;
         model       +            / dist=Constant;
         run;

Notice that the results are similar to the ZIP model shown in the FMM documentation.

The above Poisson hurdle model can also be fit using PROC NLMIXED. The following statements fit the model and may help to clarify how the model is fit. A logistic model containing only an intercept is used for the zeros process as can be seen by the statements defining the linear predictor for the zeros model (LINPZERO) and the probability of zero (PI). The events process is defined by its linear predictor (LINPNOZERO) and mean count (MUNOZERO) statements. The next two statements define the log likelihood for the Poisson hurdle model.

      proc nlmixed data=catch;
         parameters a0=0 a1=0 a2=0 b0=0;
         linpzero   = b0;
         pi         = 1/(1+exp(-linpzero));
         linpnozero = a0 + a1*(gender='F')*age + a2*(gender='M')*age;
         munozero   = exp(linpnozero);
         logpnozero = log(pi) - log(1-exp(-munozero)) - munozero -
                      lgamma(count+1) + count*log(munozero);
         if count=0 then ll=log(1-pi); else ll=logpnozero;
         model count ~ general(ll);
         run;

Beginning in SAS 9.3 TS1M2, hurdle models using the negative binomial distribution can also be fit using the DIST=TRUNCNEGBIN option to specify the truncated version of the distribution.

Operating System and Release Information