Propensity score matching and subclassification in observational studies with multi‐level treatments

In this article, we develop new methods for estimating average treatment effects in observational studies, in settings with more than two treatment levels, assuming unconfoundedness given pretreatment variables. We emphasize propensity score subclassification and matching methods which have been among the most popular methods in the binary treatment literature. Whereas the literature has suggested that these particular propensity-based methods do not naturally extend to the multi-level treatment case, we show, using the concept of weak unconfoundedness and the notion of the generalized propensity score, that adjusting for a scalar function of the pretreatment variables removes all biases associated with observed pretreatment variables. We apply the proposed methods to an analysis of the effect of treatments for fibromyalgia. We also carry out a simulation study to assess the finite sample performance of the methods relative to previously proposed methods.

[1]  Marie Davidian,et al.  Comment: Demystifying Double Robustness: A Comparison of Alternative Strategies for Estimating a Population Mean from Incomplete Data. , 2008, Statistical science : a review journal of the Institute of Mathematical Statistics.

[2]  F. Thoemmes,et al.  Propensity Score Analysis , 2015 .

[3]  Joseph Kang,et al.  Demystifying Double Robustness: A Comparison of Alternative Strategies for Estimating a Population Mean from Incomplete Data , 2007, 0804.2958.

[4]  G. Imbens,et al.  Large Sample Properties of Matching Estimators for Average Treatment Effects , 2004 .

[5]  M. Frölich.  Programme Evaluation with Multiple Treatments , 2002, SSRN Electronic Journal.

[6]  G. Imbens,et al.  Efficient Estimation of Average Treatment Effects Using the Estimated Propensity Score , 2000 .

[7]  Stephen R Cole,et al.  The consistency statement in causal inference: a definition or an assumption? , 2009, Epidemiology.

[8]  G. Imbens.  The Role of the Propensity Score in Estimating Dose-Response Functions , 1999 .

[9]  David A. Williams,et al.  Burden of illness and treatment patterns for patients with fibromyalgia. , 2012, Pain medicine.

[10]  S. Schneeweiss,et al.  Evaluating uses of data mining techniques in propensity score estimation: a simulation study , 2008, Pharmacoepidemiology and drug safety.

[11]  W. G. Cochran.  The effectiveness of adjustment by subclassification in removing bias in observational studies. , 1968, Biometrics.

[12]  J. Hahn.  On the Role of the Propensity Score in Efficient Semiparametric Estimation of Average Treatment Effects , 1998 .

[13]  Kosuke Imai,et al.  Causal Inference With General Treatment Regimes , 2004 .

[14]  Lane F Burgette,et al.  A tutorial on propensity score estimation for multiple treatments using generalized boosted models , 2013, Statistics in medicine.

[15]  Guido W. Imbens,et al.  A Martingale Representation for Matching Estimators , 2009, SSRN Electronic Journal.

[16]  Matias D. Cattaneo,et al.  Efficient semiparametric estimation of multi-valued treatment effects under ignorability , 2010 .

[17]  R. Bennett,et al.  Minimal Clinically Important Difference in the Fibromyalgia Impact Questionnaire , 2009, The Journal of Rheumatology.

[18]  A. Dawid.  Some Misleading Arguments Involving Conditional Independence , 1979 .

[19]  Markus Frlich,et al.  Finite-Sample Properties of Propensity-Score Matching and Weighting Estimators , 2004, Review of Economics and Statistics.

[20]  J. Robins,et al.  Marginal Structural Models and Causal Inference in Epidemiology , 2000, Epidemiology.

[21]  David A. Williams,et al.  Long-term Evaluation of Opioid Treatment in Fibromyalgia , 2015, The Clinical journal of pain.

[22]  Kenneth B Wells,et al.  Burden of Illness , 1968 .

[23]  K. Imai,et al.  Covariate balancing propensity score , 2014 .

[24]  P. Rosenbaum,et al.  Invited commentary: propensity scores. , 1999, American journal of epidemiology.

[25]  Donald B. Rubin,et al.  Bayesian Inference for Causal Effects: The Role of Randomization , 1978 .

[26]  B. Shepherd,et al.  GUIDO IMBENS, DONALD RUBIN, Causal Inference for Statistics, Social, and Biomedical Sciences: An Introduction. New York: Cambridge University Press. , 2016, Biometrics.

[27]  Jeremy A Rassen,et al.  Matching by Propensity Score in Cohort Studies with Three Treatment Groups , 2013, Epidemiology.

[28]  Elaine L. Zanutto,et al.  Matching With Doses in an Observational Study of a Media Campaign Against Drug Abuse , 2001, Journal of the American Statistical Association.

[29]  G. Imbens,et al.  The Propensity Score with Continuous Treatments , 2005 .

[30]  Elizabeth A Stuart,et al.  Improving propensity score weighting using machine learning , 2010, Statistics in medicine.

[31]  D. Rubin,et al.  The central role of the propensity score in observational studies for causal effects , 1983 .

[32]  Richard K. Crump,et al.  Dealing with limited overlap in estimation of average treatment effects , 2009 .

[33]  E. Foster,et al.  Propensity Score Matching: An Illustrative Analysis of Dose Response , 2003, Medical care.

[34]  M. Lechner.  Identification and Estimation of Causal Effects of Multiple Treatments Under the Conditional Independence Assumption , 1999, SSRN Electronic Journal.

[35]  Til Stürmer,et al.  Confounder summary scores when comparing the effects of multiple drug exposures , 2010, Pharmacoepidemiology and drug safety.

[36]  Elaine L. Zanutto,et al.  Using Propensity Score Subclassification for Multiple Treatment Doses to Evaluate a National Antidrug Media Campaign , 2005 .

[37]  Michael Lechner,et al.  The performance of estimators based on the propensity score , 2013 .

[38]  G. Gurtner,et al.  Statistics in medicine. , 2011, Plastic and reconstructive surgery.

[39]  G. Imbens,et al.  Efficient Estimation of Average Treatment Effects Using the Estimated Propensity Score , 2002 .

[40]  D. Rubin.  Estimating causal effects of treatments in randomized and nonrandomized studies. , 1974 .