Parameter Estimation and Hypothesis Testing of Multivariate Poisson Inverse Gaussian Regression

Multivariate Poisson regression is used in order to model two or more count response variables. The Poisson regression has a strict assumption, that is the mean and the variance of response variables are equal (equidispersion). Practically, the variance can be larger than the mean (overdispersion). Thus, a suitable method for modelling these kind of data needs to be developed. One alternative model to overcome the overdispersion issue in the multi-count response variables is the Multivariate Poisson Inverse Gaussian Regression (MPIGR) model, which is extended with an exposure variable. Additionally, a modification of Bessel function that contain factorial functions is proposed in this work to make it computable. The objective of this study is to develop the parameter estimation and hypothesis testing of the MPIGR model. The parameter estimation uses the Maximum Likelihood Estimation (MLE) method, followed by the Newton–Raphson iteration. The hypothesis testing is constructed using the Maximum Likelihood Ratio Test (MLRT) method. The MPIGR model that has been developed is then applied to regress three response variables, i.e., the number of infant mortality, the number of under-five children mortality, and the number of maternal mortality on eight predictors. The unit observation is the cities and municipalities in Java Island, Indonesia. The empirical results show that three response variables that are previously mentioned are significantly affected by all predictors.

[1]  Dimitris Karlis,et al.  A Simulation Comparison of Several Procedures for Testing the Poisson Assumption , 2000 .

[2]  Feng-Chang Xie,et al.  Influence analysis for Poisson inverse Gaussian regression models based on the EM algorithm , 2007 .

[3]  M. E. Ghitany,et al.  An EM Algorithm for Multivariate Mixed Poisson Regression Models and its Application , 2012 .

[4]  F. Famoye,et al.  Generalized poisson regression model , 1992 .

[5]  Gordon E. Willmot,et al.  A mixed poisson–inverse‐gaussian regression model , 1989 .

[6]  Samuel Musili Mwalili,et al.  Poisson Inverse Gaussian (PIG) Model for Infectious Disease Count Data , 2016 .

[7]  Jennifer Leaning,et al.  Consequences of maternal mortality on infant and child survival: a 25-year longitudinal analysis in Butajira Ethiopia (1987-2011) , 2015, BMC Genomics.

[8]  Gordon E. Willmot,et al.  The Poisson-Inverse Gaussian distribution as an alternative to the negative binomial , 1987 .

[9]  Shaw-Pin Miaou,et al.  Modeling Traffic Crash-Flow Relationships for Intersections: Dispersion Parameter, Functional Form, and Bayes Versus Empirical Bayes Methods , 2003 .

[10]  J. Juritz,et al.  BIVARIATE COMPOUND POISSON DISTRIBUTIONS , 1987 .

[11]  Yajie Zou,et al.  The Poisson inverse Gaussian (PIG) generalized linear regression model for analyzing motor vehicle crash data , 2016 .

[12]  Walter Zucchini,et al.  Parameter Estimation for the Sichel Distribution and its Multivariate Extension , 1987 .

[13]  M. S. Holla,et al.  On a poisson-inverse gaussian distribution , 1967 .