Why IV Uses Several Equations
Instrumental variables (IV) analysis separates three relationships that are easy to conflate: how treatment relates to the outcome, how the instrument relates to treatment, and how the instrument relates to the outcome. Writing these relationships as equations makes their different roles explicit. In particular, the treatment is the exposure whose effect is of interest; the instrument is a variable used to generate or explain variation in that treatment.
The equations are a way to state a model, not a guarantee that its assumptions hold. An instrument’s presence in a first-stage equation does not by itself make it valid, and a treatment coefficient in an outcome equation is not automatically a causal effect. The equations help clarify what would need to be justified.
The Linear Outcome Equation
Let \(Y_i\) be an outcome for unit \(i\), \(D_i\) the treatment, \(Z_i\) an instrument, and \(X_i\) a vector of exogenous covariates. A common linear structural outcome equation is:
Here, \(\beta_0\) is an intercept, \(\beta\) is the coefficient on treatment, \(\gamma\) is a vector of covariate coefficients, and \(u_i\) collects influences on the outcome not represented by the included variables. The prime in \(X_i'\gamma\) denotes the sum of the covariate values multiplied by their corresponding coefficients.
The structural interpretation of \(\beta\) is clearest under a constant linear treatment-effect model. In that model, changing treatment from \(d\) to \(d+1\), holding the covariates and the unit’s other outcome determinants fixed, changes the outcome by \(\beta\) units. If \(D_i\) is binary, \(\beta\) represents the model’s treated-versus-untreated outcome contrast. If treatment is measured continuously, \(\beta\) is the modeled outcome change per one-unit increase in treatment.
The equation does not require treatment to be unrelated to the outcome error. In IV settings, \(D_i\) may be associated with \(u_i\), for example because unmeasured factors influence both treatment receipt and the outcome. That is one reason ordinary regression of \(Y_i\) on \(D_i\) and \(X_i\) may not recover the causal treatment coefficient.
The First-Stage Equation
The first stage describes how the instrument and covariates relate to treatment. For a single instrument, its linear form is:
The coefficient \(\pi_1\) describes the modeled change in treatment associated with a one-unit change in the instrument, conditional on the covariates. The error \(v_i\) collects treatment determinants not represented by the instrument and covariates. It is distinct from \(u_i\); the two errors may be related.
The first-stage equation can be used with either a binary or a continuous treatment. If treatment is binary, a linear equation is a linear-probability representation: when its conditional mean is correctly specified, its right-hand side gives \(P(D_i=1\mid Z_i,X_i)\). This is a model for the conditional probability, not a claim that each individual’s binary treatment value is itself continuous. If treatment is continuous, the equation instead describes its modeled conditional mean in the units of treatment.
An instrument is relevant when it predicts treatment after accounting for the included covariates; in this scalar equation, that means \(\pi_1\ne 0\). Relevance is a feature of the first stage, not a statement that the instrument is independent of unmeasured outcome causes or excluded from the outcome through every other pathway.
Worked Example: Binary Treatment and Binary Instrument
Worked Example: Binary Treatment and Binary Instrument
Suppose a hypothetical study considers enrollment in a counseling program, \(D_i\), as a binary treatment, and a randomized offer of a place, \(Z_i\), as a binary instrument. Let \(Y_i\) be a symptom score, with lower values indicating fewer symptoms, and let \(X_i\) be baseline symptom severity centered at its mean. Consider the following illustrative structural equations:
In the outcome equation, \(D_i\) is the treatment and \(-6\) is the modeled outcome contrast for enrollment versus non-enrollment, in symptom-score points, under the constant-effect structural interpretation. \(X_i\) adjusts for baseline severity. The instrument \(Z_i\) does not appear directly in this equation.
In the first stage, \(Z_i\) predicts enrollment. At \(X_i=0\) and \(Z_i=0\), the modeled treatment probability is \(0.20+0.15(0)+0.01(0)=0.20\), or 20%. At \(X_i=0\) and \(Z_i=1\), it is \(0.20+0.15(1)+0.01(0)=0.35\), or 35%. The model therefore represents a 15-percentage-point difference in treatment probability between the two instrument values at this covariate level.
These numbers are hypothetical model coefficients, not results from an actual study. The calculation identifies the roles of the variables; it does not establish that the offer affects symptoms only through enrollment or that the offer is independent of outcome determinants.
The Reduced-Form Equation
The reduced-form equation relates the outcome directly to the instrument and covariates, without including treatment as a regressor. Under the two equations above, it can be obtained algebraically by substituting the first-stage expression for \(D_i\) into the outcome equation. It is called “reduced form” because treatment has been removed from the right-hand side.
For the general equations, substitution gives:
Thus, the corresponding reduced-form equation is:
where \(\delta_0=\beta_0+\beta\pi_0\), \(\delta_1=\beta\pi_1\), \(\delta_2=\gamma+\beta\pi_2\), and \(\varepsilon_i=u_i+\beta v_i\). These equalities follow from substitution; they do not require the errors \(u_i\) and \(v_i\) to be independent.
If the equations represent conditional means, so that \(E(u_i\mid Z_i,X_i)=0\) and \(E(v_i\mid Z_i,X_i)=0\), then \(E(\varepsilon_i\mid Z_i,X_i)=0\) as well. In that case, the reduced-form equation also describes the conditional mean of the outcome given the instrument and covariates. Without the relevant mean conditions, the algebraic identity still holds, but the right-hand side need not equal that conditional mean.
Worked Example: Continuous Treatment and Instrument
Worked Example: Continuous Treatment and Instrument
Consider a hypothetical setting in which \(D_i\) is a medication dose in milligrams per day, \(Y_i\) is systolic blood pressure in millimeters of mercury (mmHg), \(Z_i\) is a continuous formulary-access score, and \(X_i\) is baseline disease severity. Suppose the illustrative equations are:
The outcome coefficient \(-0.8\) has units of mmHg per mg/day. Under the constant linear structural interpretation, a one mg/day increase in dose corresponds to a modeled decrease of \(0.8\) mmHg, holding \(X_i\) and the other outcome determinants fixed. The first-stage coefficient \(1.5\) has units of mg/day per score unit: a one-unit increase in \(Z_i\) corresponds to a modeled increase of \(1.5\) mg/day, conditional on \(X_i\).
For example, at \(Z_i=2\) and \(X_i=0\), the first-stage conditional mean is \(4+1.5(2)+0.5(0)=7\) mg/day. The outcome equation’s structural part at that dose and covariate value is \(142-0.8(7)+1.2(0)=136.4\) mmHg, before adding \(u_i\). This is an evaluation of the specified model, not evidence that formulary access is a valid instrument.
Worked Example: Deriving the Reduced Form Numerically
Worked Example: Deriving the Reduced Form Numerically
Suppose a hypothetical linear model has a binary treatment \(D_i\), an instrument \(Z_i\), and a covariate \(X_i\):
Substitute the second equation into the first, multiplying every term in the treatment equation by the treatment coefficient \(2\):
The reduced-form intercept is \(52\), the coefficient on \(Z_i\) is \(0.8\), the covariate coefficient is \(-2\), and the reduced-form error is \(u_i+2v_i\). For instance, the \(0.8\) is the product of the outcome treatment coefficient and the first-stage instrument coefficient: \(2(0.4)=0.8\). This calculation establishes an algebraic relationship among model coefficients; by itself, it does not determine whether any coefficient can be given a causal interpretation.
What the Equations Do—and Do Not—Establish
The roles of the variables can be checked by asking where they appear. Treatment \(D_i\) appears in the outcome equation and is predicted by the instrument in the first stage. Instrument \(Z_i\) appears in the first stage and in the reduced form, but is excluded from the structural outcome equation. Exogenous covariates \(X_i\) are included in both structural equations to account for measured differences relevant to treatment and outcome. The errors \(u_i\) and \(v_i\) represent different omitted influences and need not be unrelated.
The word “exogenous” is an assumption about a variable’s relationship to unmeasured determinants in the relevant model; it is not a label made true merely by including the variable. Likewise, excluding \(Z_i\) from the outcome equation encodes the exclusion restriction, but writing that restriction does not verify it. Instrument independence and relevance also require justification. These distinctions matter because a well-formed set of equations can still describe an invalid IV design.
A further modeling choice is linearity. For a binary treatment, the linear first stage is a convenient representation of a treatment probability when its conditional mean is linear in \(Z_i\) and \(X_i\). For a continuous treatment, linearity means that the specified conditional mean changes at a constant rate over the modeled range. In either case, the equation is an approximation or assumption about the relationship, not a consequence of the variables’ definitions.
Common Mistake: Treating the First Stage as a Causal Result
A frequent error is to reason, “The instrument predicts treatment, so it is valid and the treatment coefficient is causal.” The first clause concerns relevance only. It says nothing by itself about whether the instrument is associated with unmeasured outcome causes or whether it affects the outcome through a route other than treatment. The correction is to keep the claims separate: the first stage describes instrument–treatment association; the structural outcome equation states the proposed treatment–outcome relationship; and validity requires additional assumptions about the instrument and the outcome process.
Check Your Understanding
Use the equation roles and definitions above to answer these questions.
- In \(Y_i=\beta_0+\beta D_i+X_i'\gamma+u_i\), what does \(\beta\) represent under a constant linear treatment-effect model?
- In the first-stage equation, what does the condition \(\pi_1\ne 0\) say, and what does it not establish?
- Why can the errors \(u_i\) and \(v_i\) be related without invalidating the algebra used to derive the reduced form?
- For a binary treatment, what does a linear first-stage equation represent when it correctly specifies the conditional mean?
- Substitute \(D_i=2+0.3Z_i+X_i+v_i\) into \(Y_i=10+4D_i+2X_i+u_i\). Give the reduced-form coefficients and error term.