Instrumental variables methods are often used in experimental or observational settings where we cannot credibly assume the treatment of interest is randomly assigned, but can assume the treatment is assigned my some third variable that’s unrelated to the outcome of interest.
Take two examples of causal relationships researchers might be interested in observing: the effect of distributive benefits on voter turnout and the effect of fertilizer use on farmers’ crop yield. In both cases, the distribution of resources or the use of fertilizer may depend on the outcome we wish to study — on the one hand, politicians may disburse funds according to local electoral characteristics, and on the other, famers with higher crop yield might be both more likely to use fertilizers. Instrumental variables can be a useful approach in these contexts where a treatment \(X\) is endogenous to an outcome \(Y\).
Strong instrumental variables satisfy four key assumptions:
Using the examples above, researchers could choose to instrument distributive benefits with a naturally occurring disaster — the occurence of a hurricane is correlated with the disbursement of relief aid by the government but may be independent of voter turnout. Similarly, researchers might randomly assign some farmers fertilizer vouchers as an instrument for fertilizer use.
The binary_iv_designer()
function allows researchers to declare a design of a single instrument \(Z\) to estimate the effect of an endogenous treatment \(X\) on an outcome \(Y\), all of which are binary variables. Below we declare a design using the fertilizer example.
Model:
We assume our sample has three out of four possible compliance types in a sample of \(N\) farmers: always-takers, never-takers, and compliers. In order for the monotonicy assumption to hold, we assume there are no defiers in our sample (no farmers who will adopt fertilizer use if not given the voucher, and who will not adopt it if given the voucher). Below is how we define \(X(Z)\) with \(Z \in {0,1}\) for each type:
Type | \(X_i(0)\) | \(X_i(1)\) |
---|---|---|
Always-taker | 1 | 1 |
Never-taker | 0 | 0 |
Complier | 0 | 1 |
Defier | 1 | 0 |
We consider a population of size \(N\) and we define potential outcome \(Y\) (in this case crop yield) as a linear function of the form \(Y_{i,type} = \alpha_{type} + \beta_{type} * X_i + \delta_{type} * Z_i + u_y\), where \(u_y\) is a normally distributed individual-level shock. We consider the treatment to have heterogeneous effects among different complier types.
Inquiry:
Firstly, we define the first-stage effect of \(Z\) on \(X\) as \(E[X(Z=1) - X(Z=0)]\).
Our second estimand of interest is the average effect of fertilizer use (our endogenous treatment, \(X\)), on crop productivity: \(ATE = E[Y(X = 1) - Y(X = 0)]\).
Thirdly, we are interested in the average treatment effect among compliers (farmers who use fertilizers when awarded the vouchers but who will not use it if not given voucher): \(LATE = E[Y(X = 1) - Y(X = 0) | type_i = complier]\).
Lastly, when monotonicity does not hold, but independence and exclusion restriction do, we define the estimand as a weighted average of treatment effects for compliers and defiers: \(LATE_{het} = \frac{Pr(type = c)}{Pr(type = c) - Pr(type = d)}*E[Y_i(1)-Y_i(0)| type = c] - \frac{Pr(type = d)}{Pr(type = c) - Pr(type = d)}*E[Y_i(1)-Y_i(0)| type = d]\), where \(c\) and \(d\) are compliers and defiers, respectively (see Imbens 2014: 34).
Data strategy:
We collect survey data on a representative sample of 400 farmers and randomly assign fertilizer vouchers with a probability of 0.5.
Answer strategy:
We estimate the first-stage effect by calculating the difference-in-means estimator of the function \(X ~ Z\). For an estimation of the \(ATE\), \(LATE\), and \(LATE_{het}\) we estimate the effect of endogenous treatment \(X\) on \(Y\) and use a two-stage least squares (2SLS) estimation.
N <- 100
type_probs <- c(0, 0, 1, 0)
assignment_probs <- c(0.5, 0.5, 0.5, 0.5)
outcome_sd <- 1
a <- c(1, 0, 0, 0)
b <- c(0.5, 0.5, 0.5, 0.5)
d <- c(0, 0, 0, 0)
population <- declare_population(N = N, type = sample(1:4,
N, replace = TRUE, prob = type_probs), type_label = c("Always",
"Never", "Complier", "Defier")[type], u_Z = runif(N),
u_Y = rnorm(N) * outcome_sd, Z = (u_Z < assignment_probs[type]),
X = (type == 1) + (type == 3) * Z + (type == 4) * (1 -
Z))
potential_outcomes <- declare_potential_outcomes(Y ~ a[type] +
b[type] * X + d[type] * Z + u_Y, assignment_variables = "X")
reveal <- declare_reveal(outcome_variables = Y, assignment_variables = "X")
estimand <- declare_estimand(first_stage = mean((type ==
3) - (type == 4)), ate = mean(Y_X_1 - Y_X_0), late = mean(Y_X_1[type ==
3] - Y_X_0[type == 3]), late_het = (mean(type == 3) *
mean(Y_X_1[type == 3] - Y_X_0[type == 3]) - mean(type ==
4) * mean(Y_X_1[type == 4] - Y_X_0[type == 4]))/(mean(type ==
3) - mean(type == 4)))
estimator_1 <- declare_estimator(X ~ Z, estimand = "first_stage",
label = "d-i-m")
estimator_2 <- declare_estimator(Y ~ X, estimand = c("ate",
"late", "late_het"), model = lm_robust, label = "lm_robust")
estimator_3 <- declare_estimator(Y ~ X | Z, estimand = c("ate",
"late", "late_het"), model = iv_robust, label = "iv_robust")
binary_iv_design <- population + potential_outcomes + reveal +
estimand + estimator_1 + estimator_2 + estimator_3
We now diagnose the design above and compare our estimates with those of designs that violate the first three assumptions one by one.
# Violate relevance of instrument for treatment (no first-stage effects)
iv_nonrelevant <- binary_iv_designer(N = 500, assignment_probs = c(.2, .3, .7, .5), b_Y = .5)
# Violate relevance of instrument for treatment (no first-stage effects)
iv_nonrandom <- binary_iv_designer(type_probs = c(.5,.5, 0, 0), b_Y = .5)
# Violate exclusion restriction
iv_nonexcl <- binary_iv_designer(d_Y = .1)
We also compare estimates across designs in which we assume violation of the monotonicity assumption with homogeneous and heterogeneous treatment effects across complier types.
# Violates monotonicity with homogeneous treatment effects
iv_defiers <- binary_iv_designer(type_probs = c(0,0,.8,.2))
# Violates monotonicity with heterogeneous treatment effects
iv_het <- binary_iv_designer(type_probs = c(0,0,.8,.2),
b = c(0.5, 0.5, 0.5, 0))
diagnoses <- diagnose_designs(binary_iv_design, iv_nonrelevant, iv_nonrandom,
iv_nonexcl, iv_defiers, iv_het, sims = 25)
We highlight a few takeaways. Firstly, our estimates are unbiased when all assumption hold. We show that violating any of the first three assumptions yields biased estimates.
Design Label | Estimand Label | Estimator Label | Bias | Power | Mean Estimand | |
---|---|---|---|---|---|---|
1 | binary_iv_design | ate | iv_robust | 0.02 | 0.76 | 0.50 |
3 | binary_iv_design | ate | lm_robust | 0.02 | 0.76 | 0.50 |
7 | binary_iv_design | late | iv_robust | 0.02 | 0.76 | 0.50 |
9 | binary_iv_design | late | lm_robust | 0.02 | 0.76 | 0.50 |
15 | iv_nonrelevant | ate | iv_robust | -0.90 | 0.16 | 0.50 |
17 | iv_nonrelevant | ate | lm_robust | 0.57 | 1.00 | 0.50 |
21 | iv_nonrelevant | late | iv_robust | -0.90 | 0.16 | 0.50 |
23 | iv_nonrelevant | late | lm_robust | 0.57 | 1.00 | 0.50 |
29 | iv_nonrandom | ate | iv_robust | -11.50 | 0.04 | 0.50 |
31 | iv_nonrandom | ate | lm_robust | 1.09 | 1.00 | 0.50 |
35 | iv_nonrandom | late | iv_robust | NA | 0.04 | NA |
37 | iv_nonrandom | late | lm_robust | NA | 1.00 | NA |
43 | iv_nonexcl | ate | iv_robust | 0.26 | 0.04 | 0.00 |
45 | iv_nonexcl | ate | lm_robust | 0.71 | 0.92 | 0.00 |
49 | iv_nonexcl | late | iv_robust | 0.26 | 0.04 | 0.00 |
51 | iv_nonexcl | late | lm_robust | 0.71 | 0.92 | 0.00 |
In the absence of the monotonicity assumption, our IV estimates are less biased for the LATE than the ATE estimand in the presence of heterogeneous treatment effects. Moreover, in these cases, the IV is best at estimating the weighted average of treatment effects for compliers and defiers (\(LATE_{het}\)).
Design Label | Estimand Label | Estimator Label | Bias | Power | Mean Estimand | |
---|---|---|---|---|---|---|
57 | iv_defiers | ate | iv_robust | -0.05 | 0.04 | 0.00 |
59 | iv_defiers | ate | lm_robust | -0.02 | 0.04 | 0.00 |
63 | iv_defiers | late | iv_robust | -0.05 | 0.04 | 0.00 |
65 | iv_defiers | late | lm_robust | -0.02 | 0.04 | 0.00 |
67 | iv_defiers | late_het | iv_robust | -0.05 | 0.04 | 0.00 |
69 | iv_defiers | late_het | lm_robust | -0.02 | 0.04 | 0.00 |
71 | iv_het | ate | iv_robust | 0.09 | 0.36 | 0.40 |
73 | iv_het | ate | lm_robust | -0.11 | 0.28 | 0.40 |
77 | iv_het | late | iv_robust | -0.01 | 0.36 | 0.50 |
79 | iv_het | late | lm_robust | -0.21 | 0.28 | 0.50 |
81 | iv_het | late_het | iv_robust | -0.18 | 0.36 | 0.67 |
83 | iv_het | late_het | lm_robust | -0.38 | 0.28 | 0.67 |