Time to event analysis in the presence of competing risks

Time to event (e.g., death or hospitalization) is a typical outcome of interest in medical research. It is not uncommon to see researchers use all-cause mortality as event of interest, consequently, treating various causes of death the same way. Standard methods are applied to analyze such outcomes: Kaplan-Meier curve to display the estimated survival curve, log rank test to compare survival curves of independent groups, and Cox proportional hazards to model the relationship between outcome and a set of predictor variables and covariates. Time to event for some subjects may not be observed due to end of study or because the subject dropped out. These observations are known to be right-censored and contain the information that the unknown event time is past the observed censored time. Even if death is the outcome of interest, there may be cases where all-cause mortality may not be an appropriate outcome. Instead researchers are only interested in certain causes of death. The paper by Romero-Farina et al published in this issue is a good illustration. The authors aim to determine if there is a difference between men and women with normal myocardial perfusion gated SPECT with regard to the prognoses of these patients using stress test results. The events of interest are cardiac death, non-fatal myocardial infarction (MI), and cardiac revascularization (CR). Non-cardiac death is considered as a competing risk and not a censoring event because it precludes the occurrence of cardiac death, MI and CR while censoring event only precludes observing the time of the event. Figure 1 depicts the competing risk for this particular example. Figure 1 can easily be expanded to add more competing events. Survival function and hazard function are important quantities in the analysis of time-to-event data. Survival function quantifies the probability of a person being event free at a given point in time while hazard function quantifies the risk that a person who is event-free at a given point in time will experience the event in the next instant. Kaplan-Meier curve is commonly used for estimating the survival function, and Cox proportional hazards model is commonly used to model the association of the hazard function with a set of predictor variables or covariates. In the presence of competing risks, each competing event has an associated hazard function known as the cause-specific hazards. Akin to the regular hazard function, a cause-specific hazard quantifies the risk of experiencing an event from a particular cause. Associated with this cause-specific hazard is the cumulative incidence function (CIF), also known as subdistribution function, which quantifies the probability of experiencing an event on or before a given point in time due to a particular cause. How then do we analyze the time to event data when competing events are present? One approach is to treat the competing event as a censoring event and use standard time-to-event methods for analyses. However, standard methods assume that the time to the primary event and the time to censoring are independent. In other words, the hazard at a given time point of subjects who have not experienced the primary event is the same as the hazard of those who have experienced the competing event. This is clearly violated if the competing event prevents the subject from experiencing the primary event. In the case where the competing event is not death (for example, stroke) and the primary event may still occur, the risk of the primary event may be altered after experiencing the competing event (e.g., experiencing a stroke may increase the risk of experiencing MI). As a consequence, the Kaplan-Meier curve which Reprint requests: Inmaculada Aban, PhD, Department of Biostatistics, University of Alabama at Birmingham, 1720 Second Avenue South, Birmingham, AL 35294-0022; caban@uab.edu J Nucl Cardiol 2015;22:466–7. 1071-3581/$34.00 Copyright 2014 American Society of Nuclear Cardiology.