Adjusted restricted mean survival times in observational studies

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.