Abstract
BACKGROUND: Availability of detailed data from electronic health records (EHRs) has increased the potential to examine the comparative effectiveness of dynamic treatment strategies using observational data. Inverse probability (IP) weighting of dynamic marginal structural models can control for time-varying confounders. However, IP weights for continuous treatments may be sensitive to model choice. METHODS: We describe a target trial comparing strategies for treating anemia with darbepoetin in hemodialysis patients using EHR data from the UK Renal Registry 2004 to 2016. Patients received a specified dose (microgram/week) or did not receive darbepoetin. We compared 4 methods for modeling time-varying treatment: (A) logistic regression for zero dose, standard linear regression for log dose; (B) logistic regression for zero dose, heteroscedastic linear regression for log dose; (C) logistic regression for zero dose, heteroscedastic linear regression for log dose, multinomial regression for patients who recently received very low or high doses; and (D) ordinal logistic regression. RESULTS: For this dataset, method (C) was the only approach that provided a robust estimate of the mortality hazard ratio (HR), with less-extreme weights in a fully weighted analysis and no substantial change of the HR point estimate after weight truncation. After truncating IP weights at the 95th percentile, estimates were similar across the methods. CONCLUSIONS: EHR data can be used to emulate target trials estimating the comparative effectiveness of dynamic strategies adjusting treatment to evolving patient characteristics. However, model checking, monitoring of large weights, and adaptation of model strategies to account for these is essential if an aspect of treatment is continuous.