In observational studies with censored data, exposure-outcome associations are commonly measured with adjusted hazard ratios from multivariable Cox proportional hazards models. The difference in restricted mean survival times (RMSTs) up to a pre-specified time point is an alternative measure that offers a clinically meaningful interpretation. Several regression-based methods exist to estimate an adjusted difference in RMSTs, but they digress from the model-free method of taking the area under the survival function. We derive the adjusted RMST by integrating an adjusted Kaplan-Meier estimator with inverse probability weighting (IPW). The adjusted difference in RMSTs is the area between the two IPW-adjusted survival functions. In a Monte Carlo-type simulation study, we demonstrate that the proposed estimator performs as well as two regression-based approaches: the ANCOVA-type method of Tian et al and the pseudo-observation method of Andersen et al. We illustrate the methods by reexamining the association between total cholesterol and the 10-year risk of coronary heart disease in the Framingham Heart Study. KEYWORDS inverse probability weighting, observational studies, propensity score, restricted mean survival time, survival analysis, time-to-event data
BACKGROUNDIn observational studies with time-to-event outcomes, the adjusted hazard ratio (HR) has become the effect measure of choice to quantify exposure-outcome associations when adjustment for confounders is required. The adjusted HR is commonly estimated by a multivariable Cox proportional hazards model. There are, however, two limitations of the adjusted HR. First, the interpretation of survival benefit using HRs is challenging. 1,2 The HR is a relative measure and does not communicate any information about the absolute effect. It is well established that absolute measures are valuable for public health decision-making. [3][4][5] Second, the HR depends on follow-up duration if the proportional hazards assumption is violated. In the case of non-proportional hazards, reporting only the HR may result in incorrect conclusions. 2,6 3832