Understanding Dependent Competing Risks: A Simulation Study to Illustrate the Relationship Between Cause-Specific Hazard and Marginal Hazard

In the analysis of competing risk data, the observed effect of a covariate can be obtained via a Fine and Gray sub-distribution hazard ratio. Sometimes, it is also desirable to obtain the virtual effect of a covariate as if the competing risks were non-existent. Under the latent failure time scenario, when the event of interest and the competing risk event are independent, the cause-specific hazard ratio obtained from the Cox model where the competing events are censored represents the ratio of the marginal hazards and can be interpreted as the virtual effect of the covariate. However, when the two events are not independent, the cause-specific hazard ratio is not the ratio of the marginal hazards as the ratio depends not only on the marginal hazards but also on the correlation between the competing risk and the event of interest. Using simulation, we investigated the degree to which the cause-specific hazard ratio changes relative to the marginal hazard with this correlation. It was found that the discrepancy between the cause-specific hazard ratio and the theoretical marginal hazard ratio increased as the proportion of competing risk events and the correlation between the events increased (>0.2). Depending on the direction of the correlation, the cause-specific hazard ratio can overor under-estimate the marginal hazard ratio. Using real-life datasets, we show how these results can be used to make inferences on the virtual effects.

[1]  V T Farewell,et al.  The analysis of failure times in the presence of competing risks. , 1978, Biometrics.

[2]  S. Cole,et al.  Competing risk regression models for epidemiologic data. , 2009, American journal of epidemiology.

[3]  R. Tsang,et al.  Long-term outcome in stage I and II follicular lymphoma following treatment with involved field radiation therapy alone. , 2004, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[4]  V. Gebski,et al.  First‐event or marginal estimation of cause‐specific hazards for analysing correlated multivariate failure‐time data? , 2008, Statistics in medicine.

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

[6]  D. Manjunath,et al.  BIVARIATE EXPONENTIAL DISTRIBUTIONS USING LINEAR STRUCTURES , 2002 .

[7]  R. Gilbert,et al.  A dose escalation study of hyperfractionated accelerated radiation delivered with integrated neck surgery (HARDWINS) for the management of advanced head and neck cancer. , 2008, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[8]  M. Crowder Classical Competing Risks , 2001 .

[9]  J Crowley,et al.  Estimation of failure probabilities in the presence of competing risks: new representations of old estimators. , 1999, Statistics in medicine.

[10]  G. Heinze,et al.  Competing risks analyses: objectives and approaches , 2014, European heart journal.

[11]  F. Downton Bivariate Exponential Distributions in Reliability Theory , 1970 .

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

[13]  Merlin C. Thomas,et al.  Competing-risk analysis of ESRD and death among patients with type 1 diabetes and macroalbuminuria. , 2011, Journal of the American Society of Nephrology : JASN.

[14]  Ewout W Steyerberg,et al.  Competing risks and the clinical community: irrelevance or ignorance? , 2011, Statistics in medicine.

[15]  Martin Schumacher,et al.  A competing risks analysis of bloodstream infection after stem‐cell transplantation using subdistribution hazards and cause‐specific hazards , 2007, Statistics in medicine.

[16]  Susanne Rosthøj,et al.  Competing risks as a multi-state model , 2002, Statistical methods in medical research.

[17]  Ramesh C. Gupta,et al.  Competing Risk Analysis , 2014 .