Parametric Spatial Cure Rate Models for Interval‐Censored Time‐to‐Relapse Data

Several recent papers (e.g., Chen, Ibrahim, and Sinha, 1999, Journal of the American Statistical Association 94, 909-919; Ibrahim, Chen, and Sinha, 2001a, Biometrics 57, 383-388) have described statistical methods for use with time-to-event data featuring a surviving fraction (i.e., a proportion of the population that never experiences the event). Such cure rate models and their multivariate generalizations are quite useful in studies of multiple diseases to which an individual may never succumb, or from which an individual may reasonably be expected to recover following treatment (e.g., various types of cancer). In this article we extend these models to allow for spatial correlation (estimable via zip code identifiers for the subjects) as well as interval censoring. Our approach is Bayesian, where posterior summaries are obtained via a hybrid Markov chain Monte Carlo algorithm. We compare across a broad collection of rather high-dimensional hierarchical models using the deviance information criterion, a tool recently developed for just this purpose. We apply our approach to the analysis of a smoking cessation study where the subjects reside in 53 southeastern Minnesota zip codes. In addition to the usual posterior estimates, our approach yields smoothed zip code level maps of model parameters related to the relapse rates over time and the ultimate proportion of quitters (the cure rates).

[1]  Chap T. Le,et al.  Applied Survival Analysis , 1998 .

[2]  D. Cox,et al.  Analysis of Survival Data. , 1985 .

[3]  C. R. Rao,et al.  Linear Statistical Inference and its Applications , 1968 .

[4]  J G Ibrahim,et al.  Bayesian Semiparametric Models for Survival Data with a Cure Fraction , 2001, Biometrics.

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

[6]  B. Carlin,et al.  Robust Bayesian Analysis in Medical and Epidemiological Settings , 2000 .

[7]  A I Goldman,et al.  Survivorship analysis when cure is a possibility: a Monte Carlo study. , 1984, Statistics in medicine.

[8]  V. Farewell,et al.  The use of mixture models for the analysis of survival data with long-term survivors. , 1982, Biometrics.

[9]  N. Anthonisen,et al.  Effects of multiple attempts to quit smoking and relapses to smoking on pulmonary function. Lung Health Study Research Group. , 1998, Journal of clinical epidemiology.

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

[11]  B. Carlin,et al.  Identifiability and convergence issues for Markov chain Monte Carlo fitting of spatial models. , 2000, Statistics in medicine.

[12]  L. Wasserman,et al.  Computing Bayes Factors by Combining Simulation and Asymptotic Approximations , 1997 .

[13]  D. Clayton,et al.  Empirical Bayes estimates of age-standardized relative risks for use in disease mapping. , 1987, Biometrics.

[14]  J. Ibrahim,et al.  Bayesian Inference for Multivariate Survival Data with a Cure Fraction , 2002 .

[15]  Vernon T. Farewell,et al.  Mixture models in survival analysis: Are they worth the risk? , 1986 .

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

[17]  A Yu Yakovlev,et al.  Stochastic Models of Tumor Latency and Their Biostatistical Applications , 1996 .

[18]  D. Rubin,et al.  Inference from Iterative Simulation Using Multiple Sequences , 1992 .

[19]  Joseph Berkson,et al.  Survival Curve for Cancer Patients Following Treatment , 1952 .

[20]  Bradley P. Carlin,et al.  BAYES AND EMPIRICAL BAYES METHODS FOR DATA ANALYSIS , 1996, Stat. Comput..

[21]  Bradley P. Carlin,et al.  Hierarchical Spatio-Temporal Mapping of Disease Rates , 1997 .

[22]  B. Carlin,et al.  Semiparametric spatio‐temporal frailty modeling , 2003 .

[23]  B. Carlin,et al.  Hierarchical Proportional Hazards Regression Models for Highly Stratified Data , 1999, Biometrics.

[24]  B. Carlin,et al.  Frailty modeling for spatially correlated survival data, with application to infant mortality in Minnesota. , 2003, Biostatistics.

[25]  J G Ibrahim,et al.  The Large Sample Distribution of the Weighted Log Rank Statistic Under General Local Alternatives , 1997, Lifetime data analysis.

[26]  Joseph G. Ibrahim,et al.  A New Bayesian Model For Survival Data With a Surviving Fraction , 1999 .

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

[28]  A. Gelfand,et al.  Proper multivariate conditional autoregressive models for spatial data analysis. , 2003, Biostatistics.

[29]  J. Besag Spatial Interaction and the Statistical Analysis of Lattice Systems , 1974 .

[30]  Louise Ryan,et al.  Modeling Spatial Survival Data Using Semiparametric Frailty Models , 2002, Biometrics.

[31]  Purushottam W. Laud,et al.  Predictive Model Selection , 1995 .

[32]  Alan E. Gelfand,et al.  Model choice: A minimum posterior predictive loss approach , 1998, AISTATS.