title: Risk Prediction for Prostate Cancer Recurrence Through Regularized Estimation with Simultaneous Adjustment for Nonlinear Clinical Effects.
abstract: In biomedical studies, it is of substantial interest to develop risk prediction scores using high-dimensional data such as gene expression data for clinical endpoints that are subject to censoring. In the presence of well-established clinical risk factors, investigators often prefer a procedure that also adjusts for these clinical variables. While accelerated failure time (AFT) models are a useful tool for the analysis of censored outcome data, it assumes that covariate effects on the logarithm of time-to-event are linear, which is often unrealistic in practice. We propose to build risk prediction scores through regularized rank estimation in partly linear AFT models, where high-dimensional data such as gene expression data are modeled linearly and important clinical variables are modeled nonlinearly using penalized regression splines. We show through simulation studies that our model has better operating characteristics compared to several existing models. In particular, we show that there is a non-negligible effect on prediction as well as feature selection when nonlinear clinical effects are misspecified as linear. This work is motivated by a recent prostate cancer study, where investigators collected gene expression data along with established prognostic clinical variables and the primary endpoint is time to prostate cancer recurrence. We analyzed the prostate cancer data and evaluated prediction performance of several models based on the extended c statistic for censored data, showing that 1) the relationship between the clinical variable, prostate specific antigen, and the prostate cancer recurrence is likely nonlinear, i.e., the time to recurrence decreases as PSA increases and it starts to level off when PSA becomes greater than 11; 2) correct specification of this nonlinear effect improves performance in prediction and feature selection; and 3) addition of gene expression data does not seem to further improve the performance of the resultant risk prediction scores.
authors.affiliation: Department of Biostatistics and Bioinformatics Emory University Atlanta, GA 30322, USA.
authors.forename: Qi
authors.initials: Q
authors.lastname: Long
authors.forename: Matthias
authors.initials: M
authors.lastname: Chung
authors.forename: Carlos S
authors.initials: CS
authors.lastname: Moreno
authors.forename: Brent A
authors.initials: BA
authors.lastname: Johnson
references.citation: Stat Methods Med Res. 1995 Sep;4(3):187-96
references.pmid: 8548102
references.citation: Am J Pathol. 2004 Nov;165(5):1799-807
references.pmid: 15509548
references.citation: Biostatistics. 2009 Oct;10(4):659-66
references.pmid: 19553356
references.citation: Epidemiology. 2010 Jan;21(1):128-38
references.pmid: 20010215
references.citation: Stat Med. 1997 Feb 28;16(4):385-95
references.pmid: 9044528
references.citation: Biometrika. 1965 Jun;52:203-23
references.pmid: 14341275
references.citation: Biometrics. 2009 Jun;65(2):394-404
references.pmid: 18573133
references.citation: Biotechniques. 2008 Mar;44(3):417-23
references.pmid: 18361796
references.citation: J Natl Cancer Inst. 2003 May 7;95(9):634-5
references.pmid: 12734304
references.citation: J Urol. 2003 Dec;170(6 Pt 2):S6-9; discussion S10
references.pmid: 14610404
references.citation: Stat Med. 2000 Aug 15;19(15):1997-2014
references.pmid: 10900448
references.citation: Biom J. 2010 Feb;52(1):70-84
references.pmid: 19937997