Posterior propriety and computation for the Cox regression model with applications to missing covariates

In this paper, we carry out an in-depth theoretical investigation of Bayesian inference for the Cox regression model. We establish necessary and sufficient conditions for posterior propriety of the regression coefficient, β, in Cox's partial likelihood, which can be obtained as the limiting marginal posterior distribution of β through the specification of a gamma process prior for the cumulative baseline hazard and a uniform improper prior for β. We also examine necessary and sufficient conditions for posterior propriety of the regression coefficients, β, using full likelihood Bayesian approaches in which a gamma process prior is specified for the cumulative baseline hazard. We examine characterisation of posterior propriety under completely observed data settings as well as for settings involving missing covariates. Latent variables are introduced to facilitate a straightforward Gibbs sampling scheme in the Bayesian computation. A real dataset is presented to illustrate the proposed methodology. Copyright 2006, Oxford University Press.

[1]  Margaret S. Pepe,et al.  A mean score method for missing and auxiliary covariate data in regression models , 1995 .

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

[3]  Joseph G. Ibrahim,et al.  A Bayesian justification of Cox's partial likelihood , 2003 .

[4]  J. Ibrahim,et al.  Power prior distributions for regression models , 2000 .

[5]  N. Hjort.  Nonparametric Bayes Estimators Based on Beta Processes in Models for Life History Data , 1990 .

[6]  Adrian E. Raftery,et al.  Bayesian model averaging: development of an improved multi-class, gene selection and classification tool for microarray data , 2005, Bioinform..

[7]  Z. Ying,et al.  Cox Regression with Incomplete Covariate Measurements , 1993 .

[8]  D. Madigan,et al.  Bayesian Model Averaging in Proportional Hazard Models: Assessing the Risk of a Stroke , 1997 .

[9]  W. Tsai,et al.  On using the Cox proportional hazards model with missing covariates , 1997 .

[10]  Bradley P. Carlin,et al.  Bayesian approaches for monitoring clinical trials with an application to toxoplasmic encephalitis prophylaxis , 1993 .

[11]  Hemant Ishwaran,et al.  Computational Methods for Multiplicative Intensity Models Using Weighted Gamma Processes , 2004 .

[12]  Adrian E. Raftery,et al.  Accounting for Model Uncertainty in Survival Analysis Improves Predictive Performance , 1995 .

[13]  丸山 徹.  Convex Analysisの二,三の進展について , 1977 .

[14]  J. Kalbfleisch.  Non‐Parametric Bayesian Analysis of Survival Time Data , 1978 .

[15]  M. May.  Bayesian Survival Analysis. , 2002 .

[16]  Jerald F. Lawless,et al.  Semiparametric methods for response‐selective and missing data problems in regression , 1999 .

[17]  J G Ibrahim,et al.  Maximum Likelihood Methods for Cure Rate Models with Missing Covariates , 2001, Biometrics.

[18]  V. Sondak,et al.  High- and low-dose interferon alfa-2b in high-risk melanoma: first analysis of intergroup trial E1690/S9111/C9190. , 2000, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[19]  N. Breslow.  Covariance analysis of censored survival data. , 1974, Biometrics.

[20]  Myunghee Cho Paik.  Multiple Imputation for the Cox Proportional Hazards Model with Missing Covariates , 1997, Lifetime data analysis.

[21]  Kirby L. Jackson,et al.  Log-linear analysis of censored survival data with partially observed covariates , 1989 .

[22]  Rupert G. Miller,et al.  Survival Analysis , 2022, The SAGE Encyclopedia of Research Design.

[23]  J G Ibrahim,et al.  Estimating equations with incomplete categorical covariates in the Cox model. , 1998, Biometrics.

[24]  M. Pepe,et al.  Auxiliary covariate data in failure time regression , 1995 .

[25]  J G Ibrahim,et al.  Using the EM-algorithm for survival data with incomplete categorical covariates , 1996, Lifetime data analysis.

[26]  J. Ibrahim,et al.  Likelihood-Based Methods for Missing Covariates in the Cox Proportional Hazards Model , 2001 .

[27]  Jun S. Liu,et al.  The Collapsed Gibbs Sampler in Bayesian Computations with Applications to a Gene Regulation Problem , 1994 .

[28]  Ming-Hui Chen,et al.  Propriety of posterior distribution for dichotomous quantal response models , 2000 .

[29]  Joseph G. Ibrahim,et al.  On propriety of the posterior distribution and existence of the maximum likelihood estimator for regression models with covariates missing at random , 2004 .

[30]  M. Pagano,et al.  Survival analysis. , 1996, Nutrition.

[31]  W. Michael Conklin,et al.  Monte Carlo Methods in Bayesian Computation , 2001, Technometrics.

[32]  Ming-Hui Chen,et al.  Monte Carlo Estimation of Bayesian Credible and HPD Intervals , 1999 .

[33]  D J Sargent,et al.  A general framework for random effects survival analysis in the Cox proportional hazards setting. , 1998, Biometrics.

[34]  P Gustafson,et al.  Large hierarchical Bayesian analysis of multivariate survival data. , 1997, Biometrics.

[35]  Joseph G. Ibrahim,et al.  A conditional model for incomplete covariates in parametric regression models , 1996 .