Dylan Small Department of Statistics, Wharton School, University of Pennsylvania. Based on joint work with Paul Rosenbaum

Transcription

1 Instrumental variables and their sensitivity to unobserved biases Dylan Small Department of Statistics, Wharton School, University of Pennsylvania Based on joint work with Paul Rosenbaum

2 Overview Instrumental Variable (IV) Method: Method of controlling for unmeasured confounding Example: Effect of World War II military service on future earnings. Sensitivity to unobserved biases for IV method. How strength of IVs affects sensitivity to unobserved biases: Implications for designing studies with IVs.

3 Example: WWII Veteran Status and Earnings Does military service raise or lower earnings? Angrist and Krueger (1994) studied this in context of WWII military service and 1980 earnings (using 5% public use sample of US Census). Lower earnings? Military service in WWII interrupts education or career. Higher earnings? Labor market might favor veterans, GI Bill increases education.

4 This is association not causation: WWII Vets might not be comparable to Non-Vets in terms of health, criminal behavior WWII Vets (76% of men) earned on average $4500 more in 1980 than Non-Vets.

5 Addressing Confounding Confounding Variable: Variable that is (i) not comparable between treatment and control groups. (ii) affects outcome. e.g., health, criminal behavior. If all confounders measured, they can be adjusted for by regression, propensity scores, matching methods But health, criminal behavior not measured in the Census

6 Unmeasured Confounding Graph is conditional on measured confounders (race, education up to 8 years, location of birth) Earnings Veteran Status Unobserved Confounders (Health, criminal behavior, etc.

7 Instrumental Variables Strategy Y=Outcome W=Treatment Z=IV Graph is conditional on measured confounders (race, education up to 8 years, location of birth) W: Veteran Status Y:Earnings Z: Year of Birth Unobserved Confounders (Health etc.) Extract variation in W from Z that is free of unobserved confounders and use this variation to estimate the causal effect of W on Y. Key IV Assumptions: (1) Z independent of unobserved variables; (2) Z does not have direct effect on outcome.

8

9 IV Applications in Health Research Outcome (Y ) Treatment (W ) IV ( Z ) Reference Mortality More intensive vs. less intensive treatment for heart attack patients Distance lived from cardiac care center Mortality Mortality Conventional vs. atypical antipsychotics Premature baby delivered at high level NICU vs. local hospital Prescribing physician s preference Mother s differential distance between high level NICU and local hospital McLellan, McNeil and Newhouse (1994) Wang, Schneeweiss et al. (2005) Lorch, Baiocchi, Ahlberg and Small (2012) Birth weight Maternal smoking State cigarette taxes Evans and Ringel (1999) Birth weight Maternal smoking Random assignment of free smoker s counseling Permutt and Hebel (1989) Heart attack HDL cholesterol Genes that affect HDL Voight et al. (2012)

10 Prototype IV Design: Matched Pair Encouragement Design Consider a matched pair design in which there are I matched pairs (say matched for measured confounders) and one unit j in each pair i is encouraged to receive treatment ( Z ij = 1) and the other unit j is not encouraged to receive treatment ( Z = 0). ij ' In this context, the encouragement variable Z is said to be a valid instrumental variable (IV) if Z is effectively randomly assigned: 1 1 PZ ( i1 = 1, Zi2 = 0) =, PZ ( i1 = 0, Zi2 = 1) = 2 2 (i.e., Z is not related to any unmeasured confounders). Inference can be based on two stage least squares or permutation inference. 95% CI for effect of military service on earnings using 1926 vs as IV: (-$1445, -$500)

11 A picture of the IV argument -- We created matched triples: men matched on quarter of birth, race, age, education up to 8 years and location of birth. -- This figure provides reason to doubt military service increases earnings by $ From 1924 to 1926, the proportion of veterans stayed about constant and the earnings stayed about the same. From 1926 to 1928, the proportion of veterans decreased by 50% but earnings increased, suggesting military service decreases earnings.

12 Sensitivity Analysis IV method assumes that the IV (encouragement) is effectively randomly assigned: 1 1 PZ ( i1 = 1, Zi2 = 0) =, PZ ( i1 = 0, Zi2 = 1) = 2 2 There is often concern about whether this is true. In WWII Study, there are gradual long term trends in apprenticeship, education, employment and nutrition that might bias comparisons of workers born two years apart. A sensitivity analysis asks how departures from random assignment of the IV of various magnitudes might alter a study s conclusion.

13 Model for Sensitivity Analysis For subject ij, let π ij denote the probability that ij is encouraged, π = PZ ( = 1). ij ij Suppose that two subjects ij and ik may differ in their odds of being encouraged by at most a factor of Γ 1 because they differ in terms of an unobserved covariate, uij uik, 1 πij (1 πik ) Γ i, jk,. Γ π (1 π ) ik ij If Γ= 1, IV is randomly assigned. If Γ> 1, then distribution of treatment assignments is unknown but magnitude of departure from random assignment controlled by Γ.

14 Sensitivity Analysis for WWII Study Upper Bound on One-Sided Significance Level for 1926 vs IV Γ H : β 0 H : β 4, β = causal effect of military service on earnings

15 Strength of IV We d like our study to be as insensitive to bias as possible, i.e., finding is significant for as large a Γ as possible. How does strength of IV affect sensitivity to bias. An IV is strong if encouragement has a strong effect on treatment received; An IV is weak if encouragement has only a weak effect on treatment received. Study Strong IV Weak IV World War II Study 1926 vs vs Maternal Smoking Study Random assignment of State cigarette taxes free counseling Effects of Weak IVs 1. Increased Variance 2. Increased Sensitivity to Bias

16

17 Effect of Weak IVs I: Increased Variance Y W X Z X Unobserved Variables If Z is a weak IV, then the variance of the IV estimate will be higher because less variation in W from Z can be extracted. 95% CI for effect of military service using 1926 vs IV: (-$1,445, -$500). 95% CI for effect of military service using 1924 vs IV: (-$10,130, $10,750)

18 Effect of Weak IVs II: Increased Sensitivity to Bias Power of a Sensitivity Analysis (Rosenbaum, 2004) 1 Suppose Z were in fact a valid IV so that PZ ( i1 = 1, Zi2 = 0) = PZ ( i1 = 0, Zi2 = 1) =, 2 but we didn t know this and wanted to allow for some sensitivity to bias measured by Γ Suppose also that β β0 (true causal effect minus null hypothesis causal effect) was large, so that H0 : β = β0 was substantially in error. We would like to be able to reject H0 : β = β0 when the bias could be up to some Γ (e.g., Γ =1.5). Power of a sensitivity analysis at Γ : Probability that we will reject H0 : β = β0 for Γ assuming that Z is a valid IV and a given value of β β0.

19 Effect size: ( β β0)/ σ = 1 Number of pairs I Strength of IV: Γ , ,000 lim P(Treat IV=1)- I P(Treat IV=0) When the IV is valid ( Γ= 1), the power is of course greater for stronger IVs but there is good power for all cases with sample size of 10,000 pairs. Valid but weak IVs eventually get it right. But when 1 Γ>, the power can tend to 1 or 0 depending on the strength of the IV. Weak IVs are quite sensitive to small biases.

20 Practical Consequences 1. Weak IVs that might have small bias are dangerous to use. Weak IVs are sensitive to quite small biases ( Γ> 1 yet Γ close to 1), even when the effect size ( β β0)/ σ is quite large. Unless one is confident that a weak IV is perfectly valid ( Γ= 1), its extreme sensitivity to small biases is likely to limit its usefulness to the study of enormous effects, ( β β0)/ σ >> Strong IVs that might be moderately biased are useful. A strong IV may provide useful information even if moderate biases are plausible. Consider two studies, a small study with a strong IV and a large study with a weak IV, which would have the same power if both IVs are unbiased. When there is concern that the IVs might be biased, the small study with a strong IV has considerable advantages.

21 Potential IVs in Health Outcomes Research Potential IV Differential Distance to Nearest Provider of Treatment A vs. Treatment B Geographic or Hospital Preference for Treatment A vs. B Physician Preference for Treatment A vs. B Calendar Time (one treatment may become more common over time) Genetic Variants Timing of Admission to Hospital Insurance Plan Coverage for Treatment A vs. B Randomized Encouragement at Point of Care for Treatment A vs. B When No Clear Cut Choice Strength Weak or Strong Weak or Strong Weak or Strong Weak or Strong Usually Weak Weak or Strong Weak or Strong Potentially Strong Reference for this talk: Small, D. and Rosenbaum, P. (2008). War and wages: the strength of instrumental variables and their sensitivity to unobserved biases. Journal of the American Statisical Association,