Cause‐Specific Cumulative Incidence Estimation and the Fine and Gray Model Under Both Left Truncation and Right Censoring

Summary The standard estimator for the cause-specific cumulative incidence function in a competing risks setting with left truncated and/or right censored data can be written in two alternative forms. One is a weighted empirical cumulative distribution function and the other a product-limit estimator. This equivalence suggests an alternative view of the analysis of time-to-event data with left truncation and right censoring: individuals who are still at risk or experienced an earlier competing event receive weights from the censoring and truncation mechanisms. As a consequence, inference on the cumulative scale can be performed using weighted versions of standard procedures. This holds for estimation of the cause-specific cumulative incidence function as well as for estimation of the regression parameters in the Fine and Gray proportional subdistribution hazards model. We show that, with the appropriate filtration, a martingale property holds that allows deriving asymptotic results for the proportional subdistribution hazards model in the same way as for the standard Cox proportional hazards model. Estimation of the cause-specific cumulative incidence function and regression on the subdistribution hazard can be performed using standard software for survival analysis if the software allows for inclusion of time-dependent weights. We show the implementation in the R statistical package. The proportional subdistribution hazards model is used to investigate the effect of calendar period as a deterministic external time varying covariate, which can be seen as a special case of left truncation, on AIDS related and non-AIDS related cumulative mortality.

[1]  M. Mandel Censoring and Truncation—Highlighting the Differences , 2007 .

[2]  Mei-Jie Zhang,et al.  A Mass Redistribution Algorithm for Right Censored and Left Truncated Time to Event Data. , 2009, Journal of statistical planning and inference.

[3]  Thomas M Braun,et al.  Comparing the small sample performance of several variance estimators under competing risks , 2007, Statistics in medicine.

[4]  Somnath Datta,et al.  The Kaplan–Meier Estimator as an Inverse-Probability-of-Censoring Weighted Average , 2001, The American statistician.

[5]  S Chevret,et al.  A Note on Including Time‐dependent Covariate in Regression Model for Competing Risks Data , 2005, Biometrical journal. Biometrische Zeitschrift.

[6]  Martin Schumacher,et al.  Misspecified regression model for the subdistribution hazard of a competing risk by Latouche A, Boisson V, Chevret S and Porcher R. Statistics in Medicine 2006; DOI: 10.2002/sim.2600 , 2007, Statistics in medicine.

[7]  Kholoud Porter,et al.  Effective therapy has altered the spectrum of cause-specific mortality following HIV seroconversion , 2006, AIDS.

[8]  Mei-Cheng Wang,et al.  Estimating the Occurrence Rate for Prevalent Survival Data in Competing Risks Models , 1995 .

[9]  Robert Gray,et al.  A Proportional Hazards Model for the Subdistribution of a Competing Risk , 1999 .

[10]  R. Gill,et al.  Cox's regression model for counting processes: a large sample study : (preprint) , 1982 .

[11]  David A. Binder,et al.  Fitting Cox's proportional hazards models from survey data , 1992 .

[12]  H Putter,et al.  Tutorial in biostatistics: competing risks and multi‐state models , 2007, Statistics in medicine.

[13]  R. Gray A Class of $K$-Sample Tests for Comparing the Cumulative Incidence of a Competing Risk , 1988 .

[14]  Cascade Changes in the uptake of antiretroviral therapy and survival in people with known duration of HIV infection in Europe: results from CASCADE. , 2000 .

[15]  Grace L. Yang,et al.  Estimation of the truncation probability in the random truncation model , 1998 .

[16]  Danyu Lin,et al.  On fitting Cox's proportional hazards models to survey data , 2000 .

[17]  Pao-Sheng Shen,et al.  The Product-Limit Estimate as an Inverse-Probability-Weighted Average , 2003 .