Sensitivity Analysis for Multiple Comparisons in Matched Observational Studies Through Quadratically Constrained Linear Programming

ABSTRACT A sensitivity analysis in an observational study assesses the robustness of significant findings to unmeasured confounding. While sensitivity analyses in matched observational studies have been well addressed when there is a single outcome variable, accounting for multiple comparisons through the existing methods yields overly conservative results when there are multiple outcome variables of interest. This stems from the fact that unmeasured confounding cannot affect the probability of assignment to treatment differently depending on the outcome being analyzed. Existing methods allow this to occur by combining the results of individual sensitivity analyses to assess whether at least one hypothesis is significant, which in turn results in an overly pessimistic assessment of a study's sensitivity to unobserved biases. By solving a quadratically constrained linear program, we are able to perform a sensitivity analysis while enforcing that unmeasured confounding must have the same impact on the treatment assignment probabilities across outcomes for each individual in the study. We show that this allows for uniform improvements in the power of a sensitivity analysis not only for testing the overall null of no effect, but also for null hypotheses on specific outcome variables while strongly controlling the familywise error rate. We illustrate our method through an observational study on the effect of smoking on naphthalene exposure. Supplementary materials for this article are available online.

[1]  K. Gabriel,et al.  On closed testing procedures with special reference to ordered analysis of variance , 1976 .

[2]  Leslie T. Stayner,et al.  IARC Monographs on the evaluation of carcinogenic risks to humans: Some traditional herbal medicines, some mycotoxins, naphthalene and styrene , 2002 .

[3]  P. Rosenbaum Impact of Multiple Matched Controls on Design Sensitivity in Observational Studies , 2013, Biometrics.

[4]  K. Santucci,et al.  Association of naphthalene with acute hemolytic anemia. , 2000, Academic emergency medicine : official journal of the Society for Academic Emergency Medicine.

[5]  D. Rubin,et al.  Reducing Bias in Observational Studies Using Subclassification on the Propensity Score , 1984 .

[6]  J. Goeman,et al.  The Sequential Rejection Principle of Familywise Error Control , 2010, 1211.3313.

[7]  J. Angerer,et al.  Pilot study on the naphthalene exposure of German adults and children by means of urinary 1- and 2-naphthol levels. , 2004, International journal of hygiene and environmental health.

[8]  Onyebuchi A Arah,et al.  Bias Formulas for Sensitivity Analysis of Unmeasured Confounding for General Outcomes, Treatments, and Confounders , 2011, Epidemiology.

[9]  P. Rosenbaum A Characterization of Optimal Designs for Observational Studies , 1991 .

[10]  Paul R. Rosenbaum,et al.  Design sensitivity in observational studies , 2004 .

[11]  J. L. Hodges,et al.  Estimates of Location Based on Rank Tests , 1963 .

[12]  W. G. Cochran The Planning of Observational Studies of Human Populations , 1965 .

[13]  B. Hansen Full Matching in an Observational Study of Coaching for the SAT , 2004 .

[14]  E. Stuart,et al.  Using full matching to estimate causal effects in nonexperimental studies: examining the relationship between adolescent marijuana use and adult outcomes. , 2008, Developmental psychology.

[15]  Y. D. Kim,et al.  Effects of occupation, lifestyle and genetic polymorphisms of CYP1A1, CYP2E1, GSTM1 and GSTT1 on urinary 1-hydroxypyrene and 2-naphthol concentrations. , 2001, Carcinogenesis.

[16]  P. Strickland,et al.  Involuntary Tobacco Smoke Exposure and Urinary Levels of Polycyclic Aromatic Hydrocarbons in the United States, 1999 to 2002 , 2009, Cancer Epidemiology Biomarkers & Prevention.

[17]  E. Lehmann Elements of large-sample theory , 1998 .

[18]  Paul R Rosenbaum,et al.  Sensitivity Analysis for m‐Estimates, Tests, and Confidence Intervals in Matched Observational Studies , 2007, Biometrics.

[19]  Tyler J. VanderWeele,et al.  On the definition of a confounder , 2013, Annals of statistics.

[20]  Joseph J. Sabia Does sex education affect adolescent sexual behaviors and health? , 2006, Journal of policy analysis and management : [the journal of the Association for Public Policy Analysis and Management].

[21]  E. L. Lehmann,et al.  Nonparametric Confidence Intervals for a Shift Parameter , 1963 .

[22]  Iarc Monographs Some traditional herbal medicines, some mycotoxins, naphthalene and styrene. , 2002, IARC monographs on the evaluation of carcinogenic risks to humans.

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

[24]  T. Kawamoto,et al.  A Study for the Proper Application of Urinary Naphthols, New Biomarkers for Airborne Polycyclic Aromatic Hydrocarbons , 1999, Archives of environmental contamination and toxicology.

[25]  T. Speed,et al.  On the Application of Probability Theory to Agricultural Experiments. Essay on Principles. Section 9 , 1990 .

[26]  P. Rosenbaum,et al.  Substantial Gains in Bias Reduction from Matching with a Variable Number of Controls , 2000, Biometrics.

[27]  S. Zionts,et al.  Programming with linear fractional functionals , 1968 .

[28]  J. Booth,et al.  Resampling-Based Multiple Testing. , 1994 .

[29]  Outcome Measures in Clinical Trials , 2005 .

[30]  Chunrong Jia,et al.  A Critical Review of Naphthalene Sources and Exposures Relevant to Indoor and Outdoor Air , 2010, International journal of environmental research and public health.

[31]  J. Goeman,et al.  The Inheritance Procedure: Multiple Testing of Tree-structured Hypotheses , 2012, Statistical applications in genetics and molecular biology.

[32]  Joseph P. Romano,et al.  Exact and Approximate Stepdown Methods for Multiple Hypothesis Testing , 2003 .

[33]  F. Dominici,et al.  Association between the Medicare hospice benefit and health care utilization and costs for patients with poor-prognosis cancer. , 2014, JAMA.

[34]  S. Hecht,et al.  Human urinary carcinogen metabolites: biomarkers for investigating tobacco and cancer. , 2002, Carcinogenesis.

[35]  B. Ripley,et al.  Robust Statistics , 2018, Wiley Series in Probability and Statistics.

[36]  P. Rosenbaum Design of Observational Studies , 2009, Springer Series in Statistics.

[37]  A. Conn,et al.  An Efficient Method to Solve the Minimax Problem Directly , 1978 .

[38]  N. Meinshausen Hierarchical testing of variable importance , 2008 .

[39]  T. Shakespeare,et al.  Observational Studies , 2003 .

[40]  J. S. Maritz,et al.  A note on exact robust confidence intervals for location , 1979 .

[41]  Elizabeth A. Stuart,et al.  An Introduction to Sensitivity Analysis for Unobserved Confounding in Nonexperimental Prevention Research , 2013, Prevention Science.

[42]  P. Rosenbaum Weighted M-statistics With Superior Design Sensitivity in Matched Observational Studies With Multiple Controls , 2014 .

[43]  C. Dolea,et al.  World Health Organization , 1949, International Organization.

[44]  J. Lamour,et al.  Hemolysis from exposure to naphthalene mothballs. , 1991, The New England journal of medicine.

[45]  Paul R. Rosenbaum,et al.  A Note on Optimal Matching With Variable Controls Using the Assignment Algorithm , 2001 .

[46]  J. I The Design of Experiments , 1936, Nature.

[47]  K E Warner,et al.  Smoking and lung cancer: an overview. , 1984, Cancer research.

[48]  E. Mackenzie,et al.  On Estimation of the Survivor Average Causal Effect in Observational Studies When Important Confounders Are Missing Due to Death , 2009, Biometrics.

[49]  S. Marcus,et al.  Using Omitted Variable Bias to Assess Uncertainty in the Estimation of an AIDS Education Treatment Effect , 1997 .

[50]  S. Holm A Simple Sequentially Rejective Multiple Test Procedure , 1979 .

[51]  Samuel Karlin,et al.  Mathematical Methods and Theory in Games, Programming, and Economics , 1961 .

[52]  Binbing Yu,et al.  Sensitivity analysis for trend tests: application to the risk of radiation exposure. , 2005, Biostatistics.

[53]  I. Shpitser,et al.  A New Criterion for Confounder Selection , 2011, Biometrics.

[54]  Magdalena Cerdá,et al.  Effect of the 2010 Chilean Earthquake on Posttraumatic Stress: Reducing Sensitivity to Unmeasured Bias Through Study Design , 2013, Epidemiology.

[55]  Joseph L. Gastwirth,et al.  Asymptotic separability in sensitivity analysis , 2000 .

[56]  Elizabeth A Stuart,et al.  Matching methods for causal inference: A review and a look forward. , 2010, Statistical science : a review journal of the Institute of Mathematical Statistics.

[57]  Paul W. Holland,et al.  The sensitivity of linear regression coefficients' confidence limits to the omission of a confounder , 2009, 0905.3463.

[58]  J. L. Hodges,et al.  Rank Methods for Combination of Independent Experiments in Analysis of Variance , 1962 .

[59]  Amirreza Ziai TRADITIONAL HERBAL MEDICINE , 2002 .

[60]  Samuel Karlin,et al.  Mathematical Methods and Theory in Games, Programming, and Economics , 1961 .

[61]  E. C. Hammond,et al.  Smoking and lung cancer: recent evidence and a discussion of some questions. , 1959, Journal of the National Cancer Institute.

[62]  K. Strassburger,et al.  The partitioning principle: a powerful tool in multiple decision theory , 2002 .

[63]  T. VanderWeele,et al.  Generalized Cornfield conditions for the risk difference , 2014, 1404.7175.

[64]  Yogendra P. Chaubey Resampling-Based Multiple Testing: Examples and Methods for p-Value Adjustment , 1993 .

[65]  X Zhang,et al.  Multivariate Matching and Bias Reduction in the Surgical Outcomes Study , 2001, Medical care.

[66]  Nico Voigtländer,et al.  Persecution Perpetuated: The Medieval Origins of Anti-Semitic Violence in Nazi Germany , 2012 .

[67]  Peter J. Huber,et al.  Robust Statistics , 2005, Wiley Series in Probability and Statistics.

[68]  S. Daniels,et al.  Tobacco Smoke Exposure Is Associated With the Metabolic Syndrome in Adolescents , 2005, Circulation.

[69]  Paul R. Rosenbaum,et al.  Sensitivity of Two-Sample Permutation Inferences in Observational Studies , 1990 .

[70]  Abba M Krieger,et al.  Causal conclusions are most sensitive to unobserved binary covariates , 2006, Statistics in medicine.

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

[72]  T D WILSON,et al.  Smoking and Lung Cancer , 1960, Journal of the Irish Medical Association.

[73]  G. W. Imbens Sensitivity to Exogeneity Assumptions in Program Evaluation , 2003 .

[74]  Paul R. Rosenbaum,et al.  Sensitivity Analysis for Equivalence and Difference in an Observational Study of Neonatal Intensive Care Units , 2009 .

[75]  D. Rubin [On the Application of Probability Theory to Agricultural Experiments. Essay on Principles. Section 9.] Comment: Neyman (1923) and Causal Inference in Experiments and Observational Studies , 1990 .