How to randomly generate Y from zero-inflated negative binomial dbn fit with proc genmod?

[RobertWF1]

I've recently fit a zero-inflated negative binomial distribution to hospital inpatient visit data using proc genmod:

ods trace on;
proc genmod data = IP_count_data;
  	model Y = / dist=zinb;
  	zeromodel;
	output out=ZINBFit p=p pzero=pzero;
	ods output ZeroParameterEstimates = ZeroParameterEstimates; 
	ods output ParameterEstimates = ParameterEstimates;
run;

I'd like to randomly generate observations Y from this same zero-inflated negative binomial distribution. This can be done using the following expression, pulled from an earlier question in the forum:

Y_simulated = rand('BERNOULLI', 1-pzero)*rand('NEGBINOMIAL', k, p);

I can use the pzero = zero inflation probability parameter generated from proc genmod. 

However, proc genmod does not produce the parameters k (= # of successes) and p (probability of success) which are necessary in the rand('NEGBINOMIAL', k, p) function. Can k and p be manually calculated from the proc genmod output?

[WeiChen]

@Rick_SAS might know how to do this.

[RobertWF1]

Maybe I should repost my question in the Statistical Procedures forum?

https://communities.sas.com/t5/Statistical-Procedures/bd-p/statistical_procedures

[StatDave · Accepted solution]

See this note that shows how to estimate distribution parameters for use in DATA step functions, including the negative binomial. 

[RobertWF1]

Fantastic, thank you!