A Bayesian proportional hazards model for general interval-censored data

The proportional hazards (PH) model is the most widely used semiparametric regression model for analyzing right-censored survival data based on the partial likelihood method. However, the partial likelihood does not exist for interval-censored data due to the complexity of the data structure. In this paper, we focus on general interval-censored data, which is a mixture of left-, right-, and interval-censored observations. We propose an efficient and easy-to-implement Bayesian estimation approach for analyzing such data under the PH model. The proposed approach adopts monotone splines to model the baseline cumulative hazard function and allows to estimate the regression parameters and the baseline survival function simultaneously. A novel two-stage data augmentation with Poisson latent variables is developed for the efficient computation. The developed Gibbs sampler is easy to execute as it does not require imputing any unobserved failure times or contain any complicated Metropolis-Hastings steps. Our approach is evaluated through extensive simulation studies and illustrated with two real-life data sets.

[1]  Somnath Datta,et al.  Inference Based on Imputed Failure Times for the Proportional Hazards Model with Interval-Censored Data , 1998 .

[2]  Xiaoyan Lin,et al.  Bayesian proportional hazards model for current status data with monotone splines , 2011, Comput. Stat. Data Anal..

[3]  G. Casella,et al.  The Bayesian Lasso , 2008 .

[4]  Donglin Zeng,et al.  Semiparametric additive risks model for interval-censored data , 2006 .

[5]  Donald Geman,et al.  Stochastic relaxation, Gibbs distributions, and the Bayesian restoration of images , 1984 .

[6]  Xiaoyan Lin,et al.  Bayesian Proportional Odds Models for Analyzing Current Status Data: Univariate, Clustered, and Multivariate , 2011, Commun. Stat. Simul. Comput..

[7]  J. Goedert,et al.  A prospective study of human immunodeficiency virus type 1 infection and the development of AIDS in subjects with hemophilia. , 1989, The New England journal of medicine.

[8]  Xiaoyan Lin,et al.  A Bayesian approach for analyzing case 2 interval-censored data under the semiparametric proportional odds model , 2011 .

[9]  Grace Wahba.  [Monotone Regression Splines in Action]: Comment , 1988 .

[10]  survBayes : A introduction into the package , 2009 .

[11]  Ming-Hui Chen,et al.  Bayesian analysis of generalized odds-rate hazards models for survival data , 2007, Lifetime data analysis.

[12]  A M Zaslavsky,et al.  A Markov chain Monte Carlo EM algorithm for analyzing interval-censored data under the Cox proportional hazards model. , 1998, Biometrics.

[13]  Jianguo Sun,et al.  Interval Censoring , 2003 .

[14]  Xiaoyan Lin,et al.  A semiparametric probit model for case 2 interval‐censored failure time data , 2010, Statistics in medicine.

[15]  R A Betensky,et al.  Using Conditional Logistic Regression to Fit Proportional Odds Models to Interval Censored Data , 2000, Biometrics.

[16]  Zhigang Zhang,et al.  Regression analysis of interval‐censored failure time data with linear transformation models , 2005 .

[17]  W. Gilks,et al.  Adaptive Rejection Sampling for Gibbs Sampling , 1992 .

[18]  Kurt Hornik,et al.  The Comprehensive R Archive Network , 2012 .

[19]  G. Satten.  Rank-based inference in the proportional hazards model for interval censored data , 1996 .

[20]  Philippe Lambert,et al.  Smooth estimation of survival functions and hazard ratios from interval‐censored data using Bayesian penalized B‐splines , 2011, Statistics in medicine.

[21]  Glenn Heller,et al.  Proportional hazards regression with interval censored data using an inverse probability weight , 2011, Lifetime data analysis.

[22]  Jianguo Sun,et al.  The Statistical Analysis of Interval-censored Failure Time Data , 2006 .

[23]  M. Plummer,et al.  CODA: convergence diagnosis and output analysis for MCMC , 2006 .

[24]  D Commenges,et al.  A penalized likelihood approach for arbitrarily censored and truncated data: application to age-specific incidence of dementia. , 1998, Biometrics.

[25]  W Pan,et al.  A Multiple Imputation Approach to Cox Regression with Interval‐Censored Data , 2000, Biometrics.

[26]  David R. Cox,et al.  Regression models and life tables (with discussion , 1972 .

[27]  Christian P. Robert,et al.  Monte Carlo Statistical Methods , 2005, Springer Texts in Statistics.

[28]  Wei Pan.  Extending the Iterative Convex Minorant Algorithm to the Cox Model for Interval-Censored Data , 1999 .

[29]  Joseph G. Ibrahim,et al.  Bayesian Survival Analysis , 2004 .

[30]  D Sinha,et al.  Bayesian Analysis and Model Selection for Interval‐Censored Survival Data , 1999, Biometrics.

[31]  David B Dunson,et al.  Semiparametric bayes' proportional odds models for current status data with underreporting. , 2011, Biometrics.

[32]  Adrian F. M. Smith,et al.  Sampling-Based Approaches to Calculating Marginal Densities , 1990 .

[33]  D. Finkelstein,et al.  A proportional hazards model for interval-censored failure time data. , 1986, Biometrics.

[34]  Klaus Langohr,et al.  Tutorial on methods for interval-censored data and their implementation in R , 2009 .

[35]  Ming-Hui Chen,et al.  Bayesian dynamic regression models for interval censored survival data with application to children dental health , 2013, Lifetime data analysis.

[36]  T. Cai,et al.  Hazard Regression for Interval‐Censored Data with Penalized Spline , 2003, Biometrics.

[37]  Y. Laurian,et al.  AIDS in subjects with hemophilia. , 1990, The New England journal of medicine.

[38]  Jian Huang,et al.  A Spline‐Based Semiparametric Maximum Likelihood Estimation Method for the Cox Model with Interval‐Censored Data , 2010 .

[39]  Bradley P. Carlin,et al.  Bayesian measures of model complexity and fit , 2002 .

[40]  Donald Geman,et al.  Stochastic Relaxation, Gibbs Distributions, and the Bayesian Restoration of Images , 1984, IEEE Transactions on Pattern Analysis and Machine Intelligence.