Abstract
Count data with extra zeros are common in many biomedical applications and the zero-inflated generalized power series (ZIGPS) distribution may be appropriate, in which the baseline discrete distribution is a generalized power series distribution, which is a natural extension of power series distribution. In this article, a Monte Carlo EM (MCEM) algorithm is proposed to obtain maximum likelihood estimates and standard errors for ZIGPS distribution with missing responses. A Monte Carlo approximation method combining the Gibbs sampler and M-H algorithm is used to implement E-step, whereas the M-step is completed via the method of conditional maximization. As classical model selection procedure such as Akaike's information criterion (AIC) becomes problematic for our considered ZIGPS distribution with incomplete data, some variations on AIC are presented under two different missingness mechanisms, namely, missing at random (MAR) and missing not at random (MNAR), respectively. The most attractive point is that our methods cannot only be used to select distributions and variables, but also can be used to find out whether there exists zero-inflation or not, which is conducted via score test or Bayesian test in previous literature. Finally, a simulation study and a real example are used to illustrate the proposed methodology.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.