Tutorials › Econometrics (Graduate) › Combining Multiple Instruments for Partial Identification

Instrumental Variables in Causal Analysis · Tutorial 866 of 952

Combining Multiple Instruments for Partial Identification

Translate each instrument’s information into restrictions on the same causal target, then intersect those restrictions or solve them jointly in a latent-type model.

Advanced 13 min read

What You'll Learn

  • State when instrument-specific restrictions can be combined for a common causal effect
  • Derive bounds from reduced-form contrasts with bounded direct effects
  • Intersect interval restrictions and handle negative first stages correctly
  • Represent multiple binary instruments with latent treatment-response types
  • Relax monotonicity by bounding the share of violating response types
  • Use a linear program to obtain sharp bounds under simultaneous restrictions

Why Multiple Instruments Can Sharpen Bounds

With several binary instruments, each instrument supplies information about how treatment and outcome change when that instrument is varied. If the assumptions and target are aligned, all those pieces of information must hold at once. Their joint implications can be more restrictive than the bounds from any one instrument alone.

The crucial qualification is that restrictions can be intersected only when they concern the same causal target under a coherent set of assumptions. Different instruments may shift treatment for different people, so their Wald ratios need not estimate one common effect when effects are heterogeneous. This tutorial develops two ways to combine information: a transparent common-effect sensitivity model for bounded exclusion violations, and a latent-type linear program that allows limited violations of monotonicity.

Definition: Given instruments \(Z_1,\ldots,Z_K\), a joint identified set for a target parameter is the set of parameter values that satisfy all observed-data restrictions and all maintained assumptions simultaneously. If instrument \(j\) implies a set \(\mathcal{T}_j\) for the same target, then the combined set is \(\bigcap_{j=1}^K \mathcal{T}_j\). In a latent-type model, the joint set is obtained by imposing all instrument-specific restrictions on the same latent population distribution.

Bounded Exclusion Violations for a Common Effect

Begin with a deliberately explicit sensitivity model. Let \(D\) be binary treatment and \(Y\) an outcome. For instrument \(j\), compare two assignment regimes that differ only in \(Z_j\), holding the other instruments at a specified reference setting. Define the population reduced-form contrasts \(\Delta D_j=E[D\mid Z_j=1]-E[D\mid Z_j=0]\) and \(\Delta Y_j=E[Y\mid Z_j=1]-E[Y\mid Z_j=0]\), where the expectations refer to those regimes.

Assume a common structural treatment effect \(\tau\) for the changes in treatment induced by these comparisons, and write \(\Delta Y_j=\tau\Delta D_j+\gamma_j\). Here \(\gamma_j\) is the instrument-specific direct-effect contribution to the outcome contrast. The sensitivity restriction is \(|\gamma_j|\leq\delta_j\), for a specified nonnegative bound \(\delta_j\). Instrument assignment must also be independent of the relevant potential outcomes and treatment responses, or the contrasts must be adjusted under a stated design that justifies the same relation.

Important qualification: This common-effect equation is a maintained sensitivity model, not a general consequence of instrumental-variable assumptions. With heterogeneous treatment effects, distinct instruments can shift treatment for distinct groups, and their reduced forms need not equal a common \(\tau\) times their first stages. In that setting, use a target-specific latent-type model or derive restrictions for the target actually shared across instruments.

The bounded direct effect gives a pair of inequalities for each instrument. The bounds from multiple instruments are imposed together, not chosen one at a time.

$$ |\Delta Y_j-\tau\Delta D_j|\leq\delta_j \quad\Longleftrightarrow\quad -\delta_j\leq \Delta Y_j-\tau\Delta D_j\leq\delta_j, \qquad j=1,\ldots,K. $$

When \(\Delta D_j>0\), rearranging gives \((\Delta Y_j-\delta_j)/\Delta D_j\leq\tau\leq (\Delta Y_j+\delta_j)/\Delta D_j\). When \(\Delta D_j<0\), dividing reverses the inequality directions; a safe method is to solve the two-sided inequality first and then order its endpoints. If \(\Delta D_j=0\), this instrument imposes no finite bound on \(\tau\) when \(|\Delta Y_j|\leq\delta_j\), and the restrictions are incompatible when \(|\Delta Y_j|>\delta_j\).

Worked Example: A Bound From One Instrument

Suppose, hypothetically, that an instrument has \(\Delta D_1=0.20\) and \(\Delta Y_1=0.16\), in treatment-probability and outcome units respectively. Allow an absolute direct-effect contribution of at most \(\delta_1=0.02\) outcome units. Find the values of the common effect \(\tau\) consistent with this instrument.

The restriction is \(|0.16-0.20\tau|\leq0.02\). Equivalently, \(0.14\leq0.20\tau\leq0.18\). Dividing by the positive first stage gives \(0.70\leq\tau\leq0.90\).

Thus this instrument alone permits effects from 0.70 to 0.90 outcome units per treatment change under the common-effect model. The interval does not establish that the model is true; it reports the values compatible with the specified direct-effect tolerance and the hypothetical population contrasts.

Worked Example: Intersecting Two Instruments’ Restrictions

Suppose a second hypothetical instrument has \(\Delta D_2=0.10\), \(\Delta Y_2=0.09\), and \(\delta_2=0.015\), in the same units and for the same target \(\tau\). Combine its information with the first instrument’s restriction.

For instrument 2, \(|0.09-0.10\tau|\leq0.015\), so \(0.075\leq0.10\tau\leq0.105\). Hence its interval is \(0.75\leq\tau\leq1.05\). Instrument 1 permits \([0.70,0.90]\), while instrument 2 permits \([0.75,1.05]\). Both restrictions hold only on their intersection: \[ [0.70,0.90]\cap[0.75,1.05]=[0.75,0.90]. \]

The joint bound is therefore 0.75 to 0.90 outcome units per treatment change. It is narrower than either instrument’s interval, because the restrictions overlap but neither contains the other. If the intersection were empty, the maintained common-effect model and the stated direct-effect bounds would be mutually incompatible with the population contrasts.

Worked Example: A Negative First Stage

For a third hypothetical comparison, let \(\Delta D_3=-0.10\), \(\Delta Y_3=-0.08\), and \(\delta_3=0.01\). Find its allowed range for \(\tau\), taking care with the negative first stage.

Start without dividing: \(|-0.08-(-0.10)\tau|\leq0.01\), or \(|-0.08+0.10\tau|\leq0.01\). This is equivalent to \(-0.01\leq-0.08+0.10\tau\leq0.01\). Adding 0.08 throughout gives \(0.07\leq0.10\tau\leq0.09\), and therefore \(0.70\leq\tau\leq0.90\).

The interval is positive despite the negative first stage: both the treatment and outcome contrasts are negative. Dividing an inequality by \(-0.10\) without reversing its signs could incorrectly produce a reversed or negative interval.

Simultaneous Restrictions With Latent Response Types

A common-effect interval is useful when that model is substantively defensible, but it is not a substitute for modeling heterogeneous responses. A more general formulation treats the treatment response to all instrument settings as a latent type. For two binary instruments, write an assignment setting as \(z=(z_1,z_2)\in\{0,1\}^2\). A treatment-response pattern is \(r=(r_{00},r_{10},r_{01},r_{11})\), where \(r_z\in\{0,1\}\) is the treatment the unit would receive under setting \(z\).

Let \(Y(0)\) and \(Y(1)\) be binary potential outcomes under exclusion, and let \(\pi_{r,y_0,y_1}\) be the population share with response pattern \(r\), \(Y(0)=y_0\), and \(Y(1)=y_1\). These shares are nonnegative and sum to one. Under consistency, exclusion, and independence of the instrument vector from these latent potential outcomes and treatment responses, the observed conditional distribution in every assignment cell is a sum over the same latent shares.

$$ P(D=d,Y=y\mid Z=z) = \sum_{r,y_0,y_1} \pi_{r,y_0,y_1} \,\mathbf{1}\{r_z=d\} \,\mathbf{1}\{y_d=y\}, \qquad z\in\{0,1\}^2,\ d,y\in\{0,1\}. $$

Here \(\mathbf{1}\{\cdot\}\) equals one when its condition holds and zero otherwise. Each observed-data equation is linear in the latent shares. The average treatment effect is also linear: \(E[Y(1)-Y(0)]=\sum_{r,y_0,y_1}\pi_{r,y_0,y_1}(y_1-y_0)\). Consequently, minimizing and maximizing this expression over the shares consistent with the observed cell probabilities gives lower and upper bounds on the average treatment effect.

Relaxing Monotonicity Inside the Joint Model

Monotonicity for instrument \(j\) means that changing \(Z_j\) from zero to one cannot reduce anyone’s treatment, with the other instrument setting held fixed. In the latent-type representation, this requires \(r_{(1,z_{-j})}\geq r_{(0,z_{-j})}\) for every setting \(z_{-j}\) of the other instrument. Rather than require this for every unit, allow a limited share \(\eta_j\) of types to violate it. Define \(V_j(r)=1\) if the response pattern violates this condition at any setting of the other instrument, and zero otherwise. The relaxation is the linear constraint \(\sum_{r,y_0,y_1}\pi_{r,y_0,y_1}V_j(r)\leq\eta_j\).

Both instruments’ restrictions are imposed on the same \(\pi\). Thus the combined analysis requires the observed-data equations for every instrument setting, nonnegative shares summing to one, and each specified limit on the share violating monotonicity. Exact monotonicity is the special case \(\eta_j=0\). Increasing \(\eta_j\) weakens that restriction and can widen the treatment-effect bounds; it cannot remove the observed-data equations.

1
Define the target and latent types.
Specify whether the target is the average treatment effect or another functional, and encode each unit’s treatment response to every instrument setting.
2
Translate assumptions into constraints.
State independence, exclusion, and any permitted share of monotonicity violations. Write the resulting observed-cell equations and inequalities for the common latent shares.
3
Optimize the target over the feasible set.
Find the smallest and largest target values among all nonnegative latent shares satisfying every restriction. These extrema are the joint bounds under the specified model.

This is a linear program: its objective and constraints are linear in the latent shares. The bounds are sharp relative to the stated model when every feasible share vector corresponds to a potential-outcome population satisfying the assumptions and observed probabilities. This approach avoids the invalid shortcut of calculating separate instrument-specific bounds and intersecting them when those bounds refer to different estimands or incompatible latent populations.

What the Intersection Does—and Does Not—Say

Intersecting restrictions can sharpen partial identification without selecting one instrument as the uniquely valid one. It also provides a model check: an empty feasible set indicates that at least one maintained restriction, the data inputs, or their compatibility must be reconsidered. It does not identify which assumption failed, and a nonempty set does not prove instrument validity.

Population restrictions and sampling uncertainty are distinct. In applications, contrasts and conditional cell probabilities are estimated from data. An interval formed by plugging in sample estimates is not automatically a confidence region for the partially identified target. Inference must account for uncertainty in the restrictions and for the fact that bounds can be sensitive to estimated endpoints.

Common pitfall: “Each instrument gives a Wald estimate, so I will average them or intersect intervals around them to bound the ATE.” This reasoning fails if the instruments move treatment for different subpopulations: their Wald estimands can be different LATEs, not noisy measurements of one ATE. First specify a shared target and the assumptions connecting each instrument to it. Only then combine the restrictions, either by a justified common-effect model or by imposing them jointly on a common latent-type distribution.

Check Your Understanding

Use the equations and assumptions above to answer the following questions.

  1. For a positive first stage, derive the interval for \(\tau\) from \(|\Delta Y-\tau\Delta D|\leq\delta\).
  2. Why can intervals from two instruments not be intersected automatically when treatment effects are heterogeneous?
  3. If \(\Delta D=0\), what determines whether an instrument’s bounded-exclusion restriction is compatible with the data?
  4. In the latent-type model, what does \(\eta_j=0\) impose, and what does a positive \(\eta_j\) permit?
  5. What does an empty joint feasible set imply, and what does it not identify?
Clinical Biostats