Comparing treatment effects after adjustment with multivariable Cox proportional hazards regression and propensity score methods

PURPOSE: To compare adjusted effects of drug treatment for hypertension on the risk of stroke from propensity score (PS) methods with a multivariable Cox proportional hazards (Cox PH) regression in an observational study with censored data. METHODS: From two prospective population-based cohort studies in The Netherlands a selection of subjects was used who either received drug treatment for hypertension (n = 1293) or were untreated 'candidates' for treatment (n = 954). A multivariable Cox PH was performed on the risk of stroke using eight covariates along with three PS methods. RESULTS: In multivariable Cox PH regression the adjusted hazard ratio (HR) for treatment was 0.64 (CI(95%): 0.42, 0.98). After stratification on the PS the HR was 0.58 (CI(95%): 0.38, 0.89). Matching on the PS yielded a HR of 0.49 (CI(95%): 0.27, 0.88), whereas adjustment with a continuous PS gave similar results as Cox regression. When more covariates were added (not possible in multivariable Cox model) a similar reduction in HR was reached by all PS methods. The inclusion of a simulated balanced covariate gave largest changes in HR using the multivariable Cox model and matching on the PS. CONCLUSIONS: In PS methods in general a larger number of confounders can be used. In this data set matching on the PS is sensitive to small changes in the model, probably because of the small number of events. Stratification, and covariate adjustment, were less sensitive to the inclusion of a non-confounder than multivariable Cox PH regression. Attention should be paid to PS model building and balance checking