[R] One inflated Poisson or Negative Binomal regression

Paul Johnson pauljohn at ku.edu
Mon Apr 19 19:15:53 CEST 2004


Dear Peter:

I notice there is a R code for a Zero-inflated Poisson/NB process on the 
Stanford Political Science Computational Lab (Prof. Simon Jackman) web 
page.  If I were wanting to do  a one-inflated model, I would start with 
that because, at least to my eye, it is very easy to  follow.  Mind you, 
I did not try this myself, but I bet you could make it go.  In the file 
zeroinfl.r, look at the function:

 zeroinflNegBin <- function(parms){

it is pretty clear you'd have to supply a probability model for the 
outcomes valued 1 and then fit them into the overall likelihood.

pj


http://pscl.stanford.edu/content.html
Peter Flom wrote:

>Hello
>
>I am interested in Poisson or (ideally) Negative Binomial regression
>with an inflated number of 1  responses
>
>I have seen JK Lindsey's fmr function in the gnlm library, which fits
>zero inflated Poisson (ZIP) or zero inflated negative binomial
>regression, but the help file states that for ' Poisson or related
>distributions  the mixture involves the zero category'.
>
>I had thought of perhaps subtracting 1 from all the counts and then
>fitting the ZIP or ZINB models, and then adding 1, but am not sure if
>this is legitimate, or if there is some better method.
>
>Contextual details:
>The dependent variable is number of primary sexual partners in the last
>year.  The independent variables include a) Being married or in a
>committed relationship  b) using hard drugs  c) sex  d) age
>
>N is c. 500
>
>Not surprisingly, there are a large number of 1 responses, especially
>for those who are married or in a relationship.  More surprisingly, the
>mean number of partners is the same (1.05 vs. 1.02) for people in and
>not in relationships, but the variances are very different, mostly
>because those in a relationhsip are much more likely to say exactly 1.
>
>Thanks in advance
>
>Peter
>
>Peter L. Flom, PhD
>Assistant Director, Statistics and Data Analysis Core
>Center for Drug Use and HIV Research
>  
>

-- 
Paul E. Johnson                       email: pauljohn at ku.edu
Dept. of Political Science            http://lark.cc.ku.edu/~pauljohn
1541 Lilac Lane, Rm 504                              
University of Kansas                  Office: (785) 864-9086
Lawrence, Kansas 66044-3177           FAX: (785) 864-5700




More information about the R-help mailing list