An approach to estimation in relative survival regression

• Published in: Biostatistics (Oxford, England) Vol. 10; no. 1; pp. 136 - 146
• Main Authors: Perme, Maja Pohar, Henderson, Robin, Stare, Janez
• Format: Journal Article
• Language: English
• Published: England Oxford University Press 01.01.2009; Oxford Publishing Limited (England)
• Subjects: Algorithms; Bias; Biometry - methods; Estimating techniques; Follow-Up Studies; Humans; Likelihood Functions; Myocardial Infarction - mortality; Proportional Hazards Models; Risk Factors; Statistical analysis; Survival Analysis; Relative survival; Additive model; EM algorithm
• ISSN: 1465-4644; 1468-4357; 1468-4357
• DOI: 10.1093/biostatistics/kxn021

The goal of relative survival methodology is to compare the survival experience of a cohort with that of the background population. Most often an additive excess hazard model is employed, which assumes that each person's hazard is a sum of 2 components-the population hazard obtained from life tables and an excess hazard attributable to the specific condition. Usually covariate effects on the excess hazard are assumed to have a proportional hazards structure with parametrically modelled baseline. In this paper, we introduce a new fitting procedure using the expectation-maximization algorithm, treating the cause of death as missing data. The method requires no assumptions about the baseline excess hazard thus reducing the risk of bias through misspecification. It accommodates the possibility of knowledge of cause of death for some patients, and as a side effect, the method yields an estimate of the ratio between the excess and the population hazard for each subject. More importantly, it estimates the baseline excess hazard flexibly with no additional degrees of freedom spent. Finally, it is a generalization of the Cox model, meaning that all the wealth of options in existing software for the Cox model can be used in relative survival. The method is applied to a data set on survival after myocardial infarction, where it shows how a particular form of the hazard function could be missed using the existing methods.