Bayesian Functional Mixed Models for Survival Responses with Application to Prostate Cancer

In this chapter, we propose a flexible approach to model functional measurements for survival outcomes. Often the class of models for functional observations are assumed to be linear, which may be too restrictive in some cases. We propose an alternative model, in which the simple linear mixed model has been modified by a more flexible semiparametric spline-based functional mixed model, wherein the usage of splines simplifies parameterizations and the joint modeling framework allows synergistic benefit between the regression of functional predictors and the modeling of survival data. We explicitly model the number and location of change points such that our formulation allows for an unknown set of basis functions characterizing the population-averaged and patient-specific trajectories. In addition, we propose a novel auxiliary variable scheme for a fully Bayesian estimation of our model, which not only allows dimension reduction of the parameter space but also allows efficient sampling from the conditional distributions. We illustrate our approach with a recent prostate cancer clinical trial study.

[1]  Adrian F. M. Smith,et al.  Automatic Bayesian curve fitting , 1998 .

[2]  R. Millikan,et al.  Phase II trial of alternating weekly chemohormonal therapy for patients with androgen-independent prostate cancer. , 1997, Clinical cancer research : an official journal of the American Association for Cancer Research.

[3]  J. Kalbfleisch,et al.  The Statistical Analysis of Failure Time Data: Kalbfleisch/The Statistical , 2002 .

[4]  B. Mallick,et al.  Bayesian Hierarchical Spatially Correlated Functional Data Analysis with Application to Colon Carcinogenesis , 2008, Biometrics.

[5]  J. Kalbfleisch,et al.  The Statistical Analysis of Failure Time Data , 1980 .

[6]  Xihong Lin,et al.  Semiparametric Modeling of Longitudinal Measurements and Time‐to‐Event Data–A Two‐Stage Regression Calibration Approach , 2008, Biometrics.

[7]  Peter Müller,et al.  A Bayesian Semiparametric Survival Model with Longitudinal Markers , 2010, Biometrics.

[8]  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.

[9]  F. Hsieh,et al.  Joint Modeling of Survival and Longitudinal Data: Likelihood Approach Revisited , 2006, Biometrics.

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

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

[12]  D. Cox Regression Models and Life-Tables , 1972 .

[13]  J. Ibrahim,et al.  A Flexible B‐Spline Model for Multiple Longitudinal Biomarkers and Survival , 2005, Biometrics.

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

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

[16]  M. Wand,et al.  Smoothing with Mixed Model Software , 2004 .

[17]  D. Osoba,et al.  Chemotherapy with mitoxantrone plus prednisone or prednisone alone for symptomatic hormone-resistant prostate cancer: a Canadian randomized trial with palliative end points. , 1996, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[18]  V. De Gruttola,et al.  Modelling progression of CD4-lymphocyte count and its relationship to survival time. , 1994, Biometrics.

[19]  R. Kohn,et al.  Nonparametric regression using Bayesian variable selection , 1996 .

[20]  Jane-Ling Wang,et al.  Modeling Longitudinal Data with Nonparametric Multiplicative Random Effects Jointly with Survival Data , 2008, Biometrics.

[21]  P. O'dwyer,et al.  Gonadotropin hormone-releasing hormone analogues: a new therapeutic approach for prostatic carcinoma. , 1986, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[22]  M. Pourahmadi,et al.  Bayesian analysis of covariance matrices and dynamic models for longitudinal data , 2002 .

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

[24]  Nan Laird,et al.  Joint models for efficient estimation in proportional hazards regression models. , 2003, Statistics in medicine.

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

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

[27]  Kim-Anh Do,et al.  Phase III trial of androgen ablation with or without three cycles of systemic chemotherapy for advanced prostate cancer. , 2008, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

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

[29]  W. Catalona,et al.  Measurement of prostate-specific antigen in serum as a screening test for prostate cancer. , 1991, The New England journal of medicine.

[30]  Fang Yao,et al.  FUNCTIONAL PRINCIPAL COMPONENT ANALYSIS FOR LONGITUDINAL AND SURVIVAL DATA , 2007 .

[31]  B. Mallick,et al.  Generalized Nonlinear Modeling With Multivariate Free-Knot Regression Splines , 2003 .

[32]  S R Shi,et al.  DEVELOPMENT OF AN OPTIMAL PROTOCOL FOR ANTIGEN RETRIEVAL: A ‘TEST BATTERY’ APPROACH EXEMPLIFIED WITH REFERENCE TO THE STAINING OF RETINOBLASTOMA PROTEIN (pRB) IN FORMALIN‐FIXED PARAFFIN SECTIONS , 1996, The Journal of pathology.

[33]  N. Reid,et al.  Estimating Risks of Progressing to Aids when Covariates are Measured , 1993 .

[34]  Robert Kohn,et al.  Nonparametric regression using linear combinations of basis functions , 2001, Stat. Comput..

[35]  Bradley P Carlin,et al.  Separate and Joint Modeling of Longitudinal and Event Time Data Using Standard Computer Packages , 2004 .

[36]  Ori Rosen,et al.  A Bayesian Model for Sparse Functional Data , 2008, Biometrics.

[37]  Joseph G. Ibrahim,et al.  Monte Carlo Methods in Bayesian Computation , 2000 .

[38]  W. K. Hastings,et al.  Monte Carlo Sampling Methods Using Markov Chains and Their Applications , 1970 .

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

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

[41]  Taylor Murray,et al.  Cancer statistics, 2000 , 2000, CA: a cancer journal for clinicians.

[42]  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 .

[43]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[44]  D. Pauler,et al.  Predicting time to prostate cancer recurrence based on joint models for non‐linear longitudinal biomarkers and event time outcomes , 2002, Statistics in medicine.