Incorporating prognostic factors into causal estimators: a comparison of methods for randomised controlled trials with a time-to-event outcome.
In randomised controlled trials, the effect of treatment on those who comply with allocation to active treatment can be estimated by comparing their outcome to those in the comparison group who would have complied with active treatment had they been allocated to it. We compare three estimators of the causal effect of treatment on compliers when this is a parameter in a proportional hazards model and quantify the bias due to omitting baseline prognostic factors. Causal estimates are found directly by maximising a novel partial likelihood; based on a structural proportional hazards model; and based on a 'corrected dataset' derived after fitting a rank-preserving structural failure time model. Where necessary, we extend these methods to incorporate baseline covariates. Comparisons use simulated data and a real data example. Analysing the simulated data, we found that all three methods are accurate when an important covariate was included in the proportional hazards model (maximum bias 5.4%). However, failure to adjust for this prognostic factor meant that causal treatment effects were underestimated (maximum bias 11.4%), because estimators were based on a misspecified marginal proportional hazards model. Analysing the real data example, we found that adjusting causal estimators is important to correct for residual imbalances in prognostic factors present between trial arms after randomisation. Our results show that methods of estimating causal treatment effects for time-to-event outcomes should be extended to incorporate covariates, thus providing an informative compliment to the corresponding intention-to-treat analysis.
Medical and Pharmaceutical Statistics Research Unit, Department of Mathematics and Statistics, Lancaster University, Lancaster, U.K. l.v.hampson@lancaster.ac.uk
Lisa V
LV
Hampson
Chris
C
Metcalfe
D015982:Bias; D056808:Biostatistics; D015984:Causality; D006801:Humans; D016013:Likelihood Functions; D008297:Male; D009203:Myocardial Infarction / Q000178:diet therapy / Q000401:mortality; D010349:Patient Compliance / Q000706:statistics & numerical data; D011379:Prognosis; D016016:Proportional Hazards Models; D011471:Prostatic Neoplasms / Q000401:mortality / Q000601:surgery; D016032:Randomized Controlled Trials as Topic / Q000706:statistics & numerical data*; D055502:Secondary Prevention; D013997:Time Factors; D016896:Treatment Outcome; D057832:Watchful Waiting