Sage Journals: Discover world-class research

Abstract

Clustered and multivariate survival times, such as times to recurrent events, commonly arise in biomedical and health research, and marginal survival models are often used to model such data. When a large number of predictors are available, variable selection is always an important issue when modeling such data with a survival model. We consider a Cox’s proportional hazards model for a marginal survival model. Under the sparsity assumption, we propose a penalized generalized estimating equation approach to select important variables and to estimate regression coefficients simultaneously in the marginal model. The proposed method explicitly models the correlation structure within clusters or correlated variables by using a prespecified working correlation matrix. The asymptotic properties of the estimators from the penalized generalized estimating equations are established and the number of candidate covariates is allowed to increase in the same order as the number of clusters does. We evaluate the performance of the proposed method through a simulation study and analyze two real datasets for the application.

Keywords

Clustered failure time correlation structure diverging number of predictors generalized estimating equations marginal Cox’s proportional hazards model multivariate survival time

Get full access to this article

View all access options for this article.

References

Cox

DR.

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

Tibshirani

The lasso method for variable selection in the Cox model. Stat Med 1997; 16: 385–395.

Tibshirani

Regression shrinkage and selection via the LASSO. J Royal Stat Soc Ser B 1996; 58:267–288.

Faraggi

Simon

Bayesian variable selection method for censored survival data. Biometrics 1998; 54: 1475–1485.

Fan

Variable selection for Cox’s proportional hazards model and frailty model. Ann Stat 2002; 30: 74–99.

Fan

Variable selection via nonconcave penalized likelihood and its oracle properties. J Am Stat Assoc 2001; 96: 1348–1360.

Huang

Harrington

Penalized partial likelihood regression for right-censored data with bootstrap selection of the penalty parameter.

Biometrics 2002; 58: 781–791.

Gui

Penalized Cox regression analysis in the high-dimensional and low-sample size settings, with applications to microarray gene expression data.

Bioinformatics 2005; 21: 3001–3008.

Efron

Hastie

Johnstone

, et al. Least angle regression. Ann Stat 2004; 32: 407–499.

10.

Zhang

Adaptive Lasso for Cox’s proportional hazards model. Biometrika 2007; 94: 691–703.

11.

Zou

A note on path-based variable selection in the penalized proportional hazards model. Biometrika 2008; 95: 241–247.

12.

Bickel

Tsybakov

, et al. Regularization in statistics. Test 2006; 15: 271–344.

13.

Tibshirani

Regression shrinkage and selection via the lasso: a retrospective. J Royal Stat Soc: Ser B 2011; 73: 273–282.

14.

McGilchrist

Aisbett

CW.

Regression with frailty in survival analysis.

Biometrics 1991; 47: 461–466.

15.

Hougaard

Analysis of multivariate survival data. New York: Springer, 2000.

16.

Niu

Peng

A new estimating equation approach for marginal hazard ratio estimation. Computat Stat Data Analys 2015; 87:46–56.

17.

Liang

Zeger

SL.

Longitudinal data analysis using generalized linear models. Biometrika 1986; 73: 13–22.

18.

Cai

Hypothesis testing of hazard ratio parameters in marginal models for multivariate failure time data. Lifetime Data Analysis 1999; 5: 39–53.

19.

Cai

Fan

, et al.

Variable selection for multivariate failure time data.

Biometrika 2005; 92: 303–316.

20.

Wang

Consistent model selection and data-driven smooth tests for longitudinal data in the estimating equations approach. J Royal Stat Soc: Series B 2009; 71: 177–190.

21.

Wang

Zhou

Penalized generalized estimating equations for high-dimensional longitudinal data analysis.

Biometrics 2012; 68: 353–360.

22.

Wang

GEE analysis of clustered binary data with diverging number of covariates. Ann Stat 2011; 39: 389–417.

23.

Kalbfleisch

Prentice

RL.

The statistical analysis of failure time data. 2nd ed. New York: John Wiley & Sons, 2002.

24.

Segal

Neuhaus

JM.

Robust inference for multivariate survival data.

Stat Med 1993; 12: 1019–1031.

25.

Shih

Louis

TA.

Assessing gamma frailty models for clustered failure time data.

Lifetime Data Analysis 1995; 1: 205–220.

26.

Hunter

Variable selection using MM algorithms.

Ann Stat 2005; 33: 1617–1642.

27.

Craven

Smoothing noisy data with spline functions: estimating the correct degree of smoothing by the method of generalized cross-validation. Proc Nat Acad Sci USA 1979; 31: 377–403.

28.

Wang

Tsai

CL.

Tuning parameter selectors for the smoothly clipped absolute deviation method.

Biometrika 2007; 94: 553–568.

29.

WJ.

Penalized estimating equations.

Biometrics 2003; 59: 126–132.

30.

Cho

Model selection for correlated data with diverging number of parameters. Statistica Sinica 2013; 23: 901–927.

31.

Monaco

Cai

Grizzle

Bootstrap analysis of multivariate failure time data.

Stat Med 2005; 24: 3387–3400.

32.

Andersen

Gill

RD.

Cox’s regression model for counting processes: a large sample study. Ann Stat 1982; 10: 1100–1120.

33.

Xie

Yang

, et al. Asymptotics for generalized estimating equations with large cluster sizes. Ann Stat 2003; 31: 310–347.

34.

Balan

Schiopu-Kratina

, et al. Asymptotic results with generalized estimating equations for longitudinal data. Ann Stat 2005; 33: 522–541.

35.

Clayton

DG.

A model for association in bivariate life tables and its application in epidemiological studies of familial tendency in chronic disease incidence. Biometrika 1978; 65: 141–152.

36.

Grambsch

Therneau

TM.

Proportional hazards tests and diagnostics based on weighted residuals. Biometrika 1994; 81: 515–526.

37.

Pan

Akaike’s information criterion in generalized estimating equations.

Biometrics 2001; 57: 120–125.

38.

Levy

RI.

The Framingham study: The epidemiology of atherosclerotic disease. J Am Med Assoc 1981; 245: 512–512.

39.

Freedman

DA.

On the so-called “Huber sandwich estimator” and “robust standard errors”. Am Statistician 2006; 60: 299–302.

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.27 MB

Variable selection via penalized generalized estimating equations for a marginal survival model

Abstract

Keywords

Get full access to this article

References

Supplementary Material