Re: PROC COUNTREG - Scoring with Zero Inflate Conway Maxwell Poisson (ZICMP)

Psueri — Mon, 29 Apr 2024 10:43:48 GMT

I found that the value of Lambda is Lambda = mu^nu
But the values of P(y=0), P(y!=0) and Predicted are still wrong

PROC COUNTREG - Scoring with Zero Inflate Conway Maxwell Poisson (ZICMP)

Psueri — Fri, 26 Apr 2024 16:13:50 GMT

Hello Sas Community,

I am currently working on a project where I need to utilize PROC COUNTREG in SAS, particularly focusing on utilizing ZICMP distribution with spatialeffects.

I'm trying to manually score based on the parameters, and from the documentation, I've reconstructed the following steps.

Let:

n be the number of observations
m be the number of model regressors
s be the number of spatial effects regressors
z be the number of zero model regressors
d be the number of dispersion model regressors
X1 (n x m) be the matrix containing the values of m model regressors
X2 (n x s) be the matrix containing the values of s spatial effects regressors
Z1 (n x z) be the matrix containing the values of z zero model regressors
D1 (n x d) be the matrix containing the values of d dispersion model regressors
Wmat (n x n) be the spatial weights matrix
Beta1 (m x 1) be the vector containing the coefficients of m model regressors
Beta2 (s x 1) be the vector containing the coefficients of s spatial effects regressors
Gamma (z x 1) be the vector containing the coefficients of z zero model regressors
Delta (d x 1) be the vector containing the coefficients of d dispersion model regressors
Intercept be the value of the intercept of the model
Inf_Intercept be the value of the intercept of the zero model
Disp_Intercept be the value of the intercept of the dispersion model

, or

But the values of lambda vector obtained in the output of the countreg differs from the ones computed as exp(xbeta), that is instead the value mu in the countreg output.

(in documentation there is a typo in the ZICMP, but it's correct in CMP and the results are the same of countreg output)

Now there isn't a finite form for this Normalization Factor, following the paper:

Approximating the Conway–Maxwell–Poisson distribution normalization constant

Steven B. Gillispie ∗ and Christopher G. Green

University of Washington, Seattle, USA — Department of Statistics

Technical Report no. 615

June 6, 2013

I found this approximation formula

And then

Now, the predicted results obtained do not match those obtained from the output of the countreg. The differing values are: Lambda (as mentioned above), P(y=0), P(y!=0) and Predicted (those three latter are probably due to the Lambda error)

Any examples, tips, or resources you could share would be greatly appreciated. Thank you in advance for your assistance!

Best regards,

Piero

Psueri Tracker

Re: PROC COUNTREG - Scoring with Zero Inflate Conway Maxwell Poisson (ZICMP)

PROC COUNTREG - Scoring with Zero Inflate Conway Maxwell Poisson (ZICMP)