I am trying to simulate from observed data that I have fit to a zero-inflated Poisson regression model. I fit the data in R using
zeroinfl() from the package
pscl, but I am having trouble figuring out how to derive the ZIP distribution from the coefficient estimates.
I know how to derive the predicted counts from these coefficient estimates (more information here: http://www.ats.ucla.edu/stat/stata/faq/predict_zip.htm), but can anyone help me understand how to find/derive estimates for my distribution parameters (i.e. lambda for the Poisson distribution, p for the Bernoulli distribution) that I can then sample from?