JOINT LONGITUDINAL-SURVIVAL-CURE MODELS AND THEIR APPLICATION TO PROSTATE CANCER

Many scientific investigations generate both longitudinal data and sur- vival data. Methods for the combined analysis of both kinds of data have been developed in recent years, with the main emphasis being on modeling and estima- tion. In cancer research it is common for there to be long term survivors or cured patients and methods have been developed to analyze such data. In this article, we review both joint models for the analysis of longitudinal and survival data and cure models. We then present a joint longitudinal-survival-cure model to analyze data from a study of prostate cancer patients treated with radiation therapy. In this model each patient is assumed to be either cured or susceptible to clinical recur- rence. The cured fraction is modeled as a logistic function of baseline covariates. The longitudinal PSA data is modeled as a non-linear hierarchical mixed model, with different models for the cured and susceptible groups. The clinical recurrences are modeled as a time-dependent proportional hazards model for those in the sus- ceptible group. The baseline variables are covariates in both the failure time and longitudinal models. We use both a Monte Carlo EM algorithm and Markov chain Monte Carlo techniques to fit the model. The results from the two estimation methods are compared. We focus on both selected parameters of the model and derived interpretable quantities.

[1]  J. Ibrahim,et al.  A Bayesian semiparametric joint hierarchical model for longitudinal and survival data. , 2003, Biometrics.

[2]  K. Yamaguchi.  Accelerated Failure-Time Regression Models with a Regression Model of Surviving Fraction: An Application to the Analysis of “Permanent Employment” in Japan , 1992 .

[3]  D. Schoenfeld,et al.  A model for markers and latent health status , 2000 .

[4]  Chin-Shang Li,et al.  A semi‐parametric accelerated failure time cure model , 2002, Statistics in medicine.

[5]  K. Dear,et al.  A Nonparametric Mixture Model for Cure Rate Estimation , 2000, Biometrics.

[6]  Marie Davidian,et al.  An estimator for the proportional hazards model with multiple longitudinal covariates measured with error. , 2002, Biostatistics.

[7]  R. Little.  Pattern-Mixture Models for Multivariate Incomplete Data , 1993 .

[8]  S. Zeger,et al.  Joint analysis of longitudinal data comprising repeated measures and times to events , 2001 .

[9]  Chin-Shang Li,et al.  Identifiability of cure models , 2001 .

[10]  M. Wulfsohn,et al.  A joint model for survival and longitudinal data measured with error. , 1997, Biometrics.

[11]  D. Thomas,et al.  Simultaneously modelling censored survival data and repeatedly measured covariates: a Gibbs sampling approach. , 1996, Statistics in medicine.

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

[13]  J. P. Sy,et al.  Estimation in a Cox Proportional Hazards Cure Model , 2000, Biometrics.

[14]  V. Farewell.  A model for a binary variable with time-censored observations , 1977 .

[15]  R Henderson,et al.  Joint modelling of longitudinal measurements and event time data. , 2000, Biostatistics.

[16]  Jeremy M G Taylor,et al.  The joint modeling of a longitudinal disease progression marker and the failure time process in the presence of cure. , 2002, Biostatistics.

[17]  I. Kaplan,et al.  A model of prostatic carcinoma tumor kinetics based on prostate specific antigen levels after radiation therapy , 1991, Cancer.

[18]  I. Kaplan,et al.  Prostate-specific antigen kinetics after external beam irradiation for carcinoma of the prostate. , 1994, International journal of radiation oncology, biology, physics.

[19]  G. A. Whitmore,et al.  Failure Inference From a Marker Process Based on a Bivariate Wiener Model , 1998, Lifetime data analysis.

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

[21]  J. M. Taylor,et al.  A comparison of smoothing techniques for CD4 data measured with error in a time-dependent Cox proportional hazards model. , 1998, Statistics in medicine.

[22]  A Pollack,et al.  The fall and rise of prostate‐specific antigen: Kinetics of serum prostate‐specific antigen levels after radiation therapy for prostate cancer , 1993, Cancer.

[23]  A A Tsiatis,et al.  Evaluating surrogate markers of clinical outcome when measured with error. , 1998, Biometrics.

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

[25]  M H Gail,et al.  Evaluating serial cancer marker studies in patients at risk of recurrent disease. , 1981, Biometrics.

[26]  R. Prentice.  Covariate measurement errors and parameter estimation in a failure time regression model , 1982 .

[27]  Yan Wang,et al.  Jointly Modeling Longitudinal and Event Time Data With Application to Acquired Immunodeficiency Syndrome , 2001 .

[28]  M. Wulfsohn,et al.  Modeling the Relationship of Survival to Longitudinal Data Measured with Error. Applications to Survival and CD4 Counts in Patients with AIDS , 1995 .

[29]  Anthony Y. C. Kuk,et al.  A mixture model combining logistic regression with proportional hazards regression , 1992 .

[30]  K. Doksum,et al.  Gaussian models for degradation processes-part I: Methods for the analysis of biomarker data , 1995, Lifetime data analysis.

[31]  J. M. Taylor,et al.  Survival Analysis Using Auxiliary Variables Via Multiple Imputation, with Application to AIDS Clinical Trial Data , 2002, Biometrics.

[32]  S. Berman.  A stochastic model for the distribution of HIV latency time based on T4 counts , 1990 .

[33]  J M Taylor,et al.  Semi-parametric estimation in failure time mixture models. , 1995, Biometrics.

[34]  N M Laird,et al.  Mixture models for the joint distribution of repeated measures and event times. , 1997, Statistics in medicine.

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

[36]  Anastasios A. Tsiatis,et al.  A semiparametric estimator for the proportional hazards model with longitudinal covariates measured with error , 2001 .

[37]  D. Harrington,et al.  Counting Processes and Survival Analysis , 1991 .

[38]  Hong Chang,et al.  Model Determination Using Predictive Distributions with Implementation via Sampling-Based Methods , 1992 .

[39]  Nathaniel Schenker,et al.  Analysis of Censored Survival Data with Intermittently Observed Time-Dependent Binary Covariates , 1998 .

[40]  J. Ware,et al.  Random-effects models for longitudinal data. , 1982, Biometrics.

[41]  D. Rubin.  INFERENCE AND MISSING DATA , 1975 .

[42]  N M Laird,et al.  Model-based approaches to analysing incomplete longitudinal and failure time data. , 1997, Statistics in medicine.

[43]  Seymour Geisser,et al.  8. Predictive Inference: An Introduction , 1995 .

[44]  A Tsodikov,et al.  A proportional hazards model taking account of long-term survivors. , 1998, Biometrics.

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

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