Bayesian regression analysis of panel count data under frailty nonhomogeneous Poisson process model with an unknown frailty distribution

Abstract

Panel count data arise when recurrent events are observed periodically in a study. The response variable of interest is the number of recurrent events within different time windows instead of the exact onset times of the events. The gamma frailty Poisson process model has been proposed to accommodate the within-subject correlation and overdispersion in panel count data. Although the existing methods based on the gamma frailty Poisson process model have shown some robustness against frailty distribution misspecifications, they are also found to produce biased estimates in some other cases when the gamma frailty assumption is violated. In this paper, we generalize the gamma frailty Poisson process model to allow an unknown frailty distribution for analyzing panel count data. Specifically the frailty distribution is modeled nonparametrically by assigning a Dirichlet Process Gamma Mixture prior. An efficient Gibbs sampler is developed to facilitate the Bayesian computation. Extensive simulation results suggest that the proposed Bayesian approach has an excellent performance in estimating the regression parameters and the baseline mean function and outperforms the corresponding Bayesian method based on the gamma frailty Poisson model when the gamma frailty distribution is misspecified. The proposed method is applied to a skin cancer dataset for an illustration.