Remarks on the method of propensity score

Dear Editor, I read with great interest Donald Rubin’s paper ‘The design versus the analysis of observational studies for causal effects: parallels with the design of randomized trials’ (2007; 26(1):20–36) [1], as well as the Letter To The Editor by Ian Shrier (2008; 27(14):2740–2741) [2], and Author Reply by Don Rubin (2008; 27(14):2741–2742) [3]. Shrier’s Letter posed an important question that remains unanswered in Rubin’s reply. I here venture to answer this question and to clarify related issues concerning the interpretation of propensity scores (PS) and their role in causal inference. Shrier’s question was whether, asymptotically, the use of PS methods as described by Rubin may actually increase, not decrease, bias over and above a crude, unadjusted comparison between treated and untreated subjects. The answer is: Yes, and the M-graph cited by Shrier (see also [4, 5]) provides a simple example; the crude estimate is bias-free, while PS methods introduce new bias. This occurs when treatment is strongly ignorable to begin with and becomes non-ignorable at some levels of ei . In other words, although treated and untreated units are balanced in each stratum of ei , the balance only holds relative to the covariates measured; unobserved confounders may be highly unbalanced in each stratum of ei , capable of producing significant bias. Moreover, such imbalance may be dormant in the crude estimate and awakened through the use of PS methods. There are other features of PS methods that are worth emphasizing [6]. First, the PS ei is a probabilistic, not a causal concept. Therefore, in the limit of very large sample, PS methods are bound to produce the same bias as straight stratification on the same set of measured covariates. They merely offer an effective way of approaching the asymptotic estimate which, due to the high dimensionality of X , is practically unattainable with straight stratification. Still, the asymptotic estimate is the same in both cases, and may or may not be biased, depending on the set of covariates chosen. Second, the task of choosing a sufficient (i.e. bias-eliminating) set of covariates for PS analysis requires qualitative knowledge of the causal relationships among both observed and unobserved covariates. Given such knowledge, finding a sufficient set of covariates or deciding whether a sufficient set exists are two problems that can readily be solved by graphical methods [4, 6, 7]. Finally, experimental assessments of the bias-reducing potential of PS methods (such as those described in Rubin [1]) can only be generalized to cases where the causal relationships among covariates, treatment, outcome and unobserved confounders are the same as in the experimental study. Thus, a study that proves bias reduction through the use of covariate set X does not justify the use of X in problems where the influence of unobserved confounders may be different.