Tutorials › Causal Inference › Plausibly Exogenous Instruments and Conley Sensitivity Analysis

Instrumental Variable Extensions · Tutorial 566 of 972

Plausibly Exogenous Instruments and Conley Sensitivity Analysis

Translate a defensible range of direct instrument effects into a range of treatment effects, and assess how sensitive the conclusion is to the exclusion restriction.

Advanced 12 min read

What You'll Learn

  • State the additive model that permits a direct effect of the instrument on the outcome.
  • Relate the reduced-form effect, first stage, direct effect, and treatment effect.
  • Calculate the implied treatment-effect range from a proposed direct-effect range.
  • Handle endpoint ordering correctly when the first stage is negative.
  • Find the direct effect required to move an estimate past a chosen threshold.
  • Distinguish a sensitivity range from a confidence interval.

When Exclusion Is Plausible but Not Certain

Instrumental-variable analyses often rely on an exclusion restriction: the instrument affects the outcome only by changing treatment. In an encouragement study, for example, assignment to receive an offer may change treatment uptake, and the exclusion restriction says that assignment has no other route to the outcome. The design may make this assumption credible without making it certain.

The plausibly exogenous approach makes that uncertainty explicit. Rather than setting the instrument’s direct effect on the outcome to exactly zero, specify a defensible range for it and calculate the treatment effects compatible with each value in that range. The result is a sensitivity analysis: it shows how the treatment-effect estimate depends on the strength and direction of a possible exclusion violation.

Definition: In this tutorial, a plausibly exogenous instrument is an instrument for which independence from unmeasured outcome causes is assumed or supported by the design, while the exclusion restriction is relaxed. Its direct effect on the outcome is allowed to be nonzero, but is restricted to a substantively justified range.

The Model and the Sensitivity Relationship

Let \(Z\) denote a binary randomized encouragement, \(D\) a binary treatment, and \(Y\) a quantitative outcome. Define the first-stage effect \(\pi\) as the difference in mean treatment uptake between the encouraged and unencouraged groups. Define the reduced-form effect \(\rho\) as the corresponding difference in mean outcome. Both contrasts are oriented as \(Z=1\) minus \(Z=0\).

Use an additive structural mean model in which the average outcome contrast can be decomposed into a treatment component and a direct instrument component:

$$ \rho = \beta\pi+\gamma. $$

Here \(\beta\) is the treatment effect on the outcome scale, and \(\gamma\) is the direct effect of encouragement on the outcome scale, apart from any effect operating through treatment. With a binary encouragement, \(\gamma\) is the direct outcome contrast for assignment to encouragement rather than no encouragement, holding the treatment-mediated contribution apart. If \(Y\) is measured in points, then \(\rho\) and \(\gamma\) are in points; if \(D\) is a treatment indicator, \(\pi\) is a proportion and \(\beta\) is in outcome points per treatment.

This decomposition requires more than arithmetic. The encouragement must be independent of the relevant potential outcomes and treatment responses, as would be supported by valid random assignment. The first stage must be nonzero (\(\pi\ne 0\)). The additive model must be appropriate for the target effect: for example, it is satisfied by a constant treatment effect and a common additive direct effect. In a binary-instrument setting with heterogeneous treatment effects, the corresponding interpretation can be a complier-average treatment effect if monotonicity holds and the direct effect is common across compliance types. If direct effects vary across people or compliance types, a single \(\gamma\) may not capture the needed decomposition; the analyst must justify the chosen model rather than assume the formula automatically applies.

Under the exact exclusion restriction, \(\gamma=0\), and the equation gives the usual Wald ratio \(\beta=\rho/\pi\). Under plausible exogeneity, \(\gamma\) is not fixed at zero. Rearranging the model for any specified value of \(\gamma\) gives:

$$ \beta(\gamma)=\frac{\rho-\gamma}{\pi}. $$

Thus the sensitivity calculation is direct: subtract the posited direct effect from the reduced form, then divide by the first stage. The range of \(\beta(\gamma)\) as \(\gamma\) varies over a chosen interval is the implied treatment-effect sensitivity range.

Deriving the Effect Range

Suppose the analyst judges that the direct effect lies in \([\gamma_L,\gamma_U]\), where \(\gamma_L\le\gamma_U\). Because \(\beta(\gamma)\) is linear in \(\gamma\), its extrema over this interval occur at the endpoints. Its slope is \(-1/\pi\), so the endpoint order depends on the sign of the first stage.

$$ \frac{d\beta(\gamma)}{d\gamma}=-\frac{1}{\pi}. $$

If \(\pi>0\), this slope is negative: increasing the assumed direct effect reduces the treatment effect needed to explain the same reduced form. The implied range is therefore \([\beta(\gamma_U),\beta(\gamma_L)]\). If \(\pi<0\), the slope is positive, so the range is \([\beta(\gamma_L),\beta(\gamma_U)]\). Equivalently, calculate both endpoint values and report the smaller first.

Key result: For fixed estimates \(\rho\) and \(\pi\ne0\), and a stipulated direct-effect range \([\gamma_L,\gamma_U]\), calculate \((\rho-\gamma_L)/\pi\) and \((\rho-\gamma_U)/\pi\). The smaller and larger of these values are the implied treatment-effect range. This is a sensitivity range conditional on the stipulated direct-effect range and the model; by itself, it is not a confidence interval.

Worked Example: A Positive First Stage

Suppose a hypothetical encouragement study estimates that mean treatment uptake is 0.62 in the encouraged group and 0.32 in the unencouraged group. Mean outcome is 18.4 points and 13.6 points, respectively. For illustration, assume random assignment, a nonzero first stage, and the additive model above. Subject-matter reasoning suggests that the direct effect of encouragement could plausibly be between \(-0.6\) and \(0.9\) outcome points.

The first stage and reduced form are:

$$ \pi=0.62-0.32=0.30,\qquad \rho=18.4-13.6=4.8\text{ points}. $$

The first stage is positive, so the implied treatment effect decreases as \(\gamma\) increases. At the lower endpoint:

$$ \beta(-0.6)=\frac{4.8-(-0.6)}{0.30} =\frac{5.4}{0.30}=18.0\text{ points}. $$

At the upper endpoint:

$$ \beta(0.9)=\frac{4.8-0.9}{0.30} =\frac{3.9}{0.30}=13.0\text{ points}. $$

Therefore, the implied treatment-effect range is 13.0 to 18.0 outcome points per treated person, under the stated model and direct-effect range. If the direct effect were set to zero, the estimate would be \(4.8/0.30=16.0\) points, which lies within this sensitivity range. The result does not establish that every value in the range is equally plausible; it shows the values implied by the specified endpoint assumptions.

Worked Example: A Negative First Stage

Now suppose a hypothetical instrument is coded so that \(Z=1\) reduces treatment uptake. The estimated first stage is \(\pi=-0.25\), and the reduced-form outcome contrast is \(\rho=-3.2\) points. Suppose a direct effect between \(-0.5\) and \(0.5\) points is considered plausible. Assume the same independence and additive-model conditions.

At the lower direct-effect endpoint:

$$ \beta(-0.5)=\frac{-3.2-(-0.5)}{-0.25} =\frac{-2.7}{-0.25}=10.8\text{ points}. $$

At the upper endpoint:

$$ \beta(0.5)=\frac{-3.2-0.5}{-0.25} =\frac{-3.7}{-0.25}=14.8\text{ points}. $$

The implied range is 10.8 to 14.8 points. Here \(-1/\pi\) is positive, so the treatment effect increases with the assumed direct effect. This example illustrates why endpoint values should be calculated explicitly rather than automatically assigning the lower treatment-effect value to \(\gamma_L\).

Choosing a Direct-Effect Range and Finding a Tipping Point

The calculations are only as credible as the range supplied for \(\gamma\). The range should be expressed in the outcome’s units and justified using substantive knowledge about how assignment could affect the outcome apart from treatment. Relevant considerations include the content of the encouragement, possible changes in attention or access, and evidence about the size of comparable pathways. A symmetric interval around zero is convenient, but it is not required: there may be a reason to expect a direct effect to be mostly positive, mostly negative, or asymmetric.

A useful complementary question is: how large would the direct effect have to be for the treatment effect to reach a particular threshold? Solving \(\beta=(\rho-\gamma)/\pi\) for \(\gamma\) gives:

$$ \gamma=\rho-\beta\pi. $$

This value is a tipping point for the chosen treatment-effect threshold. It does not tell us whether that direct effect is plausible; that judgment comes from the substantive context.

Worked Example: A Tipping-Point Calculation

Consider a hypothetical study with \(\rho=6.0\) outcome points and \(\pi=0.25\). Under exact exclusion, the treatment-effect estimate is \(6.0/0.25=24.0\) points. Suppose the analyst wants to know how large a positive direct effect would be needed to bring the implied treatment effect down to 10.0 points. Assume the additive model and independence conditions hold.

Substitute the target \(\beta=10.0\) and \(\pi=0.25\) into the tipping-point equation:

$$ \gamma=6.0-(10.0)(0.25) =6.0-2.5=3.5\text{ points}. $$

A direct effect of 3.5 points would therefore yield an implied treatment effect of 10.0 points:

$$ \beta(3.5)=\frac{6.0-3.5}{0.25} =\frac{2.5}{0.25}=10.0\text{ points}. $$

For comparison, if the plausible direct-effect range were \([0,2]\) points, the implied treatment-effect endpoints would be \((6.0-2)/0.25=16.0\) and \((6.0-0)/0.25=24.0\) points. The target of 10.0 points is outside that sensitivity range; reaching it would require a direct effect larger than the stipulated maximum. This comparison is informative only if the maximum of 2 points has a defensible substantive basis.

Interpreting the Sensitivity Analysis

A concise report should give the observed reduced form and first stage, the assumed direct-effect range and its units, and the resulting treatment-effect range. It should also state the assumptions that connect these quantities. For example: “With a 0.30 increase in treatment uptake and a 4.8-point reduced-form increase, allowing a direct encouragement effect from \(-0.6\) to \(0.9\) points implies treatment effects from 13.0 to 18.0 points under the additive model.” This keeps the conclusion tied to both the data and the sensitivity choice.

The range calculated from point estimates does not account for sampling variability. A confidence interval for a treatment effect must also account for uncertainty in the reduced form and first stage; dividing uncertain estimates can be especially problematic when the first stage is weak. Therefore, do not label the point-estimate sensitivity range a confidence interval or infer that it has a particular coverage probability. If inferential intervals are reported, their construction must reflect uncertainty in the estimates and the sensitivity parameter, using methods appropriate to the strength of the instrument.

The method matters because exact exclusion can be scientifically demanding, while silently ignoring a possible direct pathway can make an instrumental-variable estimate look more certain than it is. Making \(\gamma\) explicit reveals the direction and scale of the assumption: for a positive first stage, a more positive direct effect leaves less of the reduced-form association to attribute to treatment; a negative direct effect leaves more. Reporting how the conclusion moves across a justified range is more transparent than presenting only the exact-exclusion estimate.

Common Mistakes

  • Calling the sensitivity range a confidence interval. The endpoint calculations above vary \(\gamma\) while holding the estimated \(\rho\) and \(\pi\) fixed. They do not, by themselves, quantify sampling uncertainty.
  • Reversing the endpoints without checking the first stage. With \(\pi>0\), increasing \(\gamma\) decreases \(\beta\); with \(\pi<0\), it increases \(\beta\). Calculate both values and order them numerically.
  • Using incompatible units. The direct effect \(\gamma\) must be measured on the same outcome scale as \(\rho\). Dividing their difference by the treatment uptake contrast gives an outcome effect per treatment, not another outcome contrast per encouragement.
  • Treating a chosen range as evidence. A range is an assumption informed by the design and substantive knowledge, not something guaranteed by the observed reduced form. Explain why its endpoints are plausible.
  • Assuming the formula solves any form of exclusion violation. The decomposition requires a suitable additive model and a defensible interpretation of \(\gamma\). Heterogeneous or type-specific direct effects may not be represented by one common value.
Key takeaway: Plausible exogeneity replaces the claim that the instrument’s direct effect is exactly zero with an explicit range for that effect. Under a justified additive decomposition, each value \(\gamma\) implies \(\beta(\gamma)=(\rho-\gamma)/\pi\). The resulting sensitivity range makes the dependence of the treatment-effect conclusion on exclusion visible.

Check Your Understanding

Use the model and conventions in this tutorial to answer the following questions.

  1. If \(\rho=5.0\), \(\pi=0.20\), and \(\gamma\) ranges from \(-0.5\) to \(0.5\), what treatment-effect range is implied? Show both endpoint calculations.
  2. Why does the order of the treatment-effect endpoints depend on whether the first stage is positive or negative?
  3. In a study where the outcome is measured in symptom-score points, what units should be used for \(\gamma\), and what are the units of \(\beta\) when \(D\) is binary?
  4. For \(\rho=4.0\), \(\pi=0.25\), what direct effect would make the implied treatment effect equal to 8.0 points?
  5. Explain why a range obtained by varying \(\gamma\) over its stipulated values is not, by itself, a confidence interval.
Clinical Biostats