Essentials of survival analysis

1. Essentials of survival analysis How to practice evidence based oncology European School of Oncology July 2004 Antwerp, Belgium Dr. Iztok Hozo Professor of Mathematics Indiana University NW www.iun.edu/~mathiho

2. Time-to-Event • Time-to-event data are generated when the measure of interest is the amount of time to occurrence of an event of interest. • For Example: • – Time from randomization to death in clinical trial • – Time from randomization to recurrence in a cancer clinical trial • – Time from diagnosis of cancer to death due to the cancer • – Time from diagnosis of cancer to death due to any causes • – Time from remission to relapse of leukemia • – Time from HIV infection to AIDS • – Time from exposure to cancer incidence in an epidemiological cohort study

3. Censoring • Censoring occurs when we have some information, but we don’t know the exact time-to-event measure. • For example, patients typically enter a clinical study at the time randomization (or the time of diagnosis, or treatment) and are followed up until the event of interest is observed. However, censoring may occur for the following reasons: • a person does not experience the event before the study ends; • death due to a cause not considered to be the event of interest (traffic accident, adverse drug reaction,…); and • loss to follow-up, for example, if the person moves. We say that the survival time is censored. These are examples of right censoring, which is the most common form of censoring in medical studies. For these patients, the complete time-to-event measure is unknown; we only know that the true time-to-event measure is greater than the observed measurement.

4. Example: X means an event occurred; O means that the subject was censored.

5. Example 2 (from Kleinbaum: “Survival Analysis”) Consider data from a retrospective study of 13 women who had surgery for breast cancer. The survival times are: 23, 47, 69, 70+, 71+, 100+, 101+, 148, 181, 198+, 208+, 212+, 224+ (the “+” means that that particular patient was censored)

6. Survival Curve - Calculus • S(t) = cumulative survival function = proportion that survive until time tf(t) = frequency distribution of age at deathh(t) = hazard function (i.e. death rate at age t) = event rate • Relationships:

7. Probability Distribution function Probability Density function Survival function Distribution Function, Survival Function and Density Function

8. Creating a Kaplan-Meier curve • For each non-censored failure timetj (time-to-event time) evaluate: • nj = number at risk before time tj • dj = number of deaths from tj-1 to tj • Fraction = estimated probability of surviving past tj-1given that you are at risk at time The Product Limit Formula:

9. Kaplan-Meier Product Limit Estimate Consider data from a retrospective study of 45 women who had surgery for breast cancer. The survival times are: 23, 47, 69, 70+, 71+, 100+, 101+, 148, 181, 198+, 208+, 212+, 224+

11. Survival Curves – more examples

12. Log-Rank test for two groups • Suppose we have two groups, each with a different treatment. Usually, we represent this kind of situation in a 2x2 table.

13. Log-Rank test for two groups • If the data are given through time, we have a series of 2x2 tables. • Expected number of events If the two groups were the same – what would the expected number of events be? • Observed minus expectedThis is a measure of deviation of one treatment from their average (the expected) • Log-rank statistic measures whether the data in the two groups are statistically “different”.

14. Comparing Survival Functions • Question: Did the treatment make a difference in the survival experience of the two groups? • Hypothesis: H0: S1(t)=S2(t) for all t ≥ 0. • Three often used tests: • Log-rank test (aka Mantel-Haenszel Test); • Wilcoxon Test; • Likelihood ratio test.

15. Log-rank example (from Kleinbaum: “Survival Analysis”)

17. Time n1 m1 q1 S1 n2 m2 q2 S2 0 21 0 100% 21 0 100% 1 21 0 100% 21 2 90% 2 21 0 100% 19 2 81% 3 21 0 100% 17 1 76% 4 21 0 100% 16 2 67% 5 21 0 100% 14 2 57% 6 21 3 1 86% 12 0 57% 7 17 1 1 81% 12 0 57% 8 16 0 81% 12 4 38% 10 15 1 2 75% 8 0 38% 11 13 0 75% 8 2 29% 12 12 0 75% 6 2 19% 13 12 1 69% 4 0 19% 15 11 0 69% 4 1 14% 16 11 1 3 63% 3 0 14% 17 10 0 63% 3 1 10% 22 7 1 54% 2 1 5% 23 6 1 5 45% 1 1 0% Total: 9 21 Survival data vs. two-by-two table = different

18. Log-Rank test for several groups • The null hypothesis is that all the survival curves are the same. • Log-rank statistic is given by the sum: • This statistic has Chi-square distribution with (# of groups – 1) degrees of freedom.

19. Cox Proportional Hazards Regression • Most interesting survival-analysis research examines the relationship between survival — typically in the form of the hazard function — and one or more explanatory variables (or covariates). • Most common are linear-like models for the log hazard. • For example, a parametric regression model based on the exponential distribution, • Needed to assess effect of multiple covariates on survival • Cox-proportional hazards is the most commonly used multivariate survival method • Easy to implement in SPSS, Stata, or SAS • Parametric approaches are an alternative, but they require stronger assumptions about h(t).

20. Multivariate methods: Cox proportional hazards • Assumes multiplicative risk—this is the proportional hazard assumption • Conveniently separates baseline hazard function from covariates • Baseline hazard function over time • Covariates are time independent • Nonparametric • Can handle both continuous and categorical predictor variables (think: logistic, linear regression) • Without knowing baseline hazard ho(t), can still calculate coefficients for each covariate, and therefore hazard ratio

21. Limitations of Cox PH model • Covariates normally do not vary over time • True with respect to gender, ethnicity, or congenital condition • One can program time-dependent variables • Baseline hazard function, ho(t), is never specified, but Cox PH models known hazard functions • You can estimate ho(t) accurately if you need to estimate S(t).

22. Hazard Ratio • Interesting to interpretFor example, if HR = 0.70, we can deduce the following: • Relative effect on survival isor 30% reduction of the risk of death • Absolute Difference in survival is given as so, if S = 60%, which represents a 10% difference. • Difference in median survival is given as the difference between the median/HR and the median. For example, if the median is months, then the difference is given as or 10.71 months increase in median survival.