Sage Journals: Discover world-class research

Abstract

We introduce a supervised learning mixture model for censored durations (C-mix) to simultaneously detect subgroups of patients with different prognosis and order them based on their risk. Our method is applicable in a high-dimensional setting, i.e. with a large number of biomedical covariates. Indeed, we penalize the negative log-likelihood by the Elastic-Net, which leads to a sparse parameterization of the model and automatically pinpoints the relevant covariates for the survival prediction. Inference is achieved using an efficient Quasi-Newton Expectation Maximization algorithm, for which we provide convergence properties. The statistical performance of the method is examined on an extensive Monte Carlo simulation study and finally illustrated on three publicly available genetic cancer datasets with high-dimensional covariates. We show that our approach outperforms the state-of-the-art survival models in this context, namely both the CURE and Cox proportional hazards models penalized by the Elastic-Net, in terms of C-index, AUC(t) and survival prediction. Thus, we propose a powerful tool for personalized medicine in cancerology.

Keywords

Coxs proportional hazards model CURE model Elastic-Net regularization high-dimensional estimation mixture duration model survival analysis

Get full access to this article

View all access options for this article.

References

Alizadeh

Eisen

Davis

et al.

Distinct types of diffuse large b-cell lymphoma identified by gene expression profiling. Nature 2000; 403: 503–511.

Rosenwald

Wright

Chan

et al.

The use of molecular profiling to predict survival after chemotherapy for diffuse large-b-cell lymphoma. New Engl J Med 2002; 346: 1937–1947.

Golub

Slonim

Tamayo

et al.

Molecular classification of cancer: class discovery and class prediction by gene expression monitoring. Science 1999; 286: 531–537.

Hastie

Tibshirani

Botstein

et al.

Supervised harvesting of expression trees. Genome Biol 2001; 2: research0003.1–research0003.12

Tibshirani

Hastie

Narasimhan

et al.

Diagnosis of multiple cancer types by shrunken centroids of gene expression. Proc Natl Acad Sci 2002; 99: 6567–6572.

Beer

Kardia

Huang

et al.

Gene-expression profiles predict survival of patients with lung adenocarcinoma. Nat Med 2002; 8: 816–824.

Bhattacharjee

Richards

Staunton

et al.

Classification of human lung carcinomas by mRNA expression profiling reveals distinct adenocarcinoma subclasses. Proc Natl Acad Sci 2001; 98: 13790–13795.

Sørlie

Perou

Tibshirani

et al.

Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications. Proc Natl Acad Sci 2001; 98: 10869–10874.

Shipp

Ross

Tamayo

et al.

Diffuse large b-cell lymphoma outcome prediction by gene-expression profiling and supervised machine learning. Nat Med 2002; 8: 68–74.

10.

Van’t Veer

Dai

Van De Vijver

et al.

Gene expression profiling predicts clinical outcome of breast cancer. Nature 2002; 415: 530–536.

11.

Banfield

Raftery

. Model-based gaussian and non-gaussian clustering. Biometrics 1993; 49: 803–821.

12.

Cox

. Regression models and life-tables. J Roy Stat Soc B 1972; 34: 187–220.

13.

Simon

Friedman

Hastie

et al.

Regularization paths for Coxs proportional hazards model via coordinate descent. J Stat Softw 2011; 39: 1–13.

14.

Farewell

. The use of mixture models for the analysis of survival data with long-term survivors. Biometrics 1982; 38: 1041–1046.

15.

Kuk

Chen

. A mixture model combining logistic regression with proportional hazards regression. Biometrika 1992; 79: 531–541.

16.

Kuo

Peng

. A mixture-model approach to the analysis of survival data. Biostat – Basel 2000; 5: 255–272.

17.

De Angelis

Capocaccia

Hakulinen

et al.

Mixture models for cancer survival analysis: application to population-based data with covariates. Stat Med 1999; 18: 441–454.

18.

Klein

Moeschberger

. Survival analysis: techniques for censored and truncated data, New York: Springer Science & Business Media, 2005.

19.

Zou

Hastie

. Regularization and variable selection via the elastic net. J R Stat Soc Series B Stat Methodol 2005; 67: 301–320.

20.

Dempster

Laird

Rubin

. Maximum likelihood from incomplete data via the EM algorithm. J Roy Stat Soc B 1977; 39: 1–38.

21.

Zhu

Byrd

et al.

Algorithm 778: L-BFGS-B: Fortran subroutines for large-scale bound-constrained optimization. ACM Trans Math Soft 1997; 23: 550–560.

22.

Andrew G and Gao J. Scalable training of L1-regularized log-linear models. In: International conference on machine learning. New York: ACM, 2007, pp. 33–40

23.

. On the convergence properties of the EM algorithm. Ann Stat 1983; 11: 95–103.

24.

Breslow

. Contribution to the discussion of the paper by DR Cox. J Roy Stat Soc B 1972; 34: 216–217.

25.

Fan J, Feng Y, Wu Y, et al. High-dimensional variable selection for Coxs proportional hazards model. In: Borrowing strength: theory powering applications–a festschrift for Lawrence D. Brown. Institute of Mathematical Statistics, Beachwood, OH, 2010, pp. 70–86

26.

Wei

. The accelerated failure time model: a useful alternative to the cox regression model in survival analysis. Stat Med 1992; 11: 1871–1879.

27.

Mukherjee

Maiti

. On some properties of positive definite Toeplitz matrices and their possible applications. Linear Algebra Appl 1988; 102: 211–240.

28.

Heagerty

Lumley

Pepe

. Time-dependent roc curves for censored survival data and a diagnostic marker. Biometrics 2000; 56: 337–344.

29.

Blanche

Latouche

Viallon

. Time-dependent AUC with right-censored data: a survey. Risk Assess Eval Pred 2013; 215: 239–251.

30.

Harrell

Lee

Mark

. Tutorial in biostatistics multivariable prognostic models: issues in developing models, evaluating assumptions and adequacy, and measuring and reducing errors. Stat Med 1996; 15: 361–387.

31.

Heagerty

Zheng

. Survival model predictive accuracy and roc curves. Biometrics 2005; 61: 92–105.

32.

Uno

Cai

Pencina

et al.

On the C-statistics for evaluating overall adequacy of risk prediction procedures with censored survival data. Stat Med 2011; 30: 1105–1117.

33.

Kaplan

Meier

. Nonparametric estimation from incomplete observations. J Am Stat Assoc 1958; 53: 457–481.

34.

Zhou

Yau

Gray

et al.

Enhanced NFκB and AP-1 transcriptional activity associated with antiestrogen resistant breast cancer. BMC Cancer 2007; 7: 1–1.

35.

Oskarsson

Acharyya

Zhang

et al.

Breast cancer cells produce tenascin c as a metastatic niche component to colonize the lungs. Nat Med 2011; 17: 867–874.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

0.38 MB

C-mix: A high-dimensional mixture model for censored durations,with applications to genetic data

Abstract

Keywords

Get full access to this article

References

Supplementary Material