Why the Adjustment Set Changes the Question
In longitudinal studies, a covariate measured after an earlier treatment can be both a consequence of that treatment and a predictor of a later treatment and outcome. Including it in a regression may help describe comparisons among people with the same measured covariate history. But it also changes what the treatment coefficient compares: the coefficient no longer includes outcome differences that operate through the covariate held fixed by the model.
This distinction matters because a regression coefficient is not automatically a total causal effect. Its statistical meaning is determined by the model and the variables held constant; its causal meaning additionally requires assumptions about how treatment and covariates were assigned. In particular, a coefficient after post-treatment adjustment generally answers a conditional question, not the question of what would happen if treatment changed and all its downstream consequences were allowed to unfold.
What a Treatment Coefficient Conditions On
Consider a regression with an additive conditional mean and no treatment-by-covariate interactions:
Here \(A\) is treatment, \(L\) is an included covariate, and \(Y\) is the outcome. Under this model, \(\beta_A\) is the difference in the model’s expected outcome for a one-unit increase in \(A\), holding \(L\) fixed. If \(A\) is binary, it is the modelled difference between \(A=1\) and \(A=0\) at the same value of \(L\). The units are outcome units per treatment unit. In a fitted regression, the estimated coefficient describes this adjusted association in the data, subject to the model specification.
The total effect instead compares outcomes under treatment strategies while allowing downstream variables to take the values they would take under each strategy. For a binary treatment at one time, its population average is \(E[Y(1)-Y(0)]\), where \(Y(a)\) denotes the outcome that would occur under treatment level \(a\). In a longitudinal setting, the strategies specify treatment at multiple times. A coefficient that holds a later covariate fixed does not generally make this total comparison.
Worked Example: A Treatment Effect Transmitted Through a Covariate
Worked Example: A Treatment Effect Transmitted Through a Covariate
Suppose a hypothetical study has a binary treatment \(A\), a later health measure \(L\), and a continuous outcome \(Y\), measured in symptom-score points. Assume the following structural equations describe the population:
Assume the errors have mean zero, are independent of treatment and of one another, and the equations correctly describe the relevant conditional means. Also assume treatment is assigned independently of the errors, and that the interventions \(A=0\) and \(A=1\) are well defined. These strong assumptions make the causal calculation transparent; they are not guaranteed merely by fitting a regression.
The regression of \(Y\) on \(A\) and \(L\) has conditional mean \(10+2A+4L\). Thus, at the same value of \(L\), its treatment contrast is:
That is the conditional treatment contrast. To calculate the total effect, allow \(L\) to respond to treatment. Increasing \(A\) from 0 to 1 increases \(L\) by 3 units. The change in \(Y\) is therefore the direct 2-point term plus the \(4\)-point-per-unit effect transmitted through the 3-unit change in \(L\):
Equivalently, substituting the equation for \(L\) into the equation for \(Y\) gives \(Y=18+14A+4U_L+U_Y\). The total effect is 14 points, whereas the coefficient on \(A\) in the regression adjusting for \(L\) is 2 points. They differ because the adjusted coefficient compares treatment levels at fixed \(L\), excluding the pathway through which treatment changes \(L\).
Why the Adjusted Coefficient Need Not Be the Total Effect
The example illustrates a general algebraic point. Suppose, for simplicity, that a treatment changes a later covariate by \(\alpha\) units per treatment unit, and that the covariate changes the outcome by \(\beta_L\) units per covariate unit. If the outcome equation also has a treatment term \(\beta_A A\), then substituting the covariate equation into the outcome equation adds the transmitted effect \(\alpha\beta_L\) to the treatment’s total effect.
The coefficient \(\beta_A\) is the contrast holding \(L\) fixed; the substituted coefficient \(\beta_A+\alpha\beta_L\) is the total effect in this simple linear model. This result relies on the equations being correctly specified, on the absence of relevant interactions, and on the stated error and treatment-assignment assumptions. It is not a universal rule for converting adjusted and unadjusted regression coefficients.
With time-varying treatment, the interpretation requires particular care. A later covariate may be affected by an earlier treatment while also predicting a subsequent treatment and the outcome. Adjusting for it conditions on a consequence of prior treatment, so the earlier treatment coefficient excludes pathways through that covariate. Further, the resulting conditional association need not be a clean causal direct effect: common causes of the covariate and outcome, or other consequences of treatment, can make such an interpretation invalid. A regression coefficient alone does not establish the assumptions needed for a causal reading.
Worked Example: A Two-Period Treatment History
Worked Example: A Two-Period Treatment History
Consider a hypothetical setting with an initial treatment \(A_0\), a later health measure \(L_1\), a subsequent treatment dose \(A_1\), and an outcome \(Y\) in functional-score points. Suppose the structural equations are:
Assume the equations are correct, the errors have mean zero, and the errors are mutually independent and independent of \(A_0\). Treat the displayed equations as a hypothetical structural model in which the interventions and dose changes are well defined. These assumptions let us compare the model’s conditional contrast with its total effect; they should not be inferred from an observed regression alone.
In the outcome equation, the coefficient of \(A_0\) is 3. It describes the change in predicted \(Y\), in functional-score points, for a one-unit increase in \(A_0\) while \(L_1\) and \(A_1\) are both held fixed. For example, at \(L_1=4\) and \(A_1=2\), the predicted outcome at \(A_0=0\) is \(40+2(4)+4(2)=56\). At \(A_0=1\) with those same covariate values, it is \(40+3+2(4)+4(2)=59\). The conditional contrast is \(59-56=3\) points.
The total effect allows both later variables to change when \(A_0\) changes. A one-unit increase in \(A_0\) changes \(L_1\) by 2 units. It also changes \(A_1\) directly by 0.5 dose units and indirectly through \(L_1\) by \(0.25\times2=0.5\) dose units, for a total change of 1 dose unit. Consequently:
For a direct check, set \(U_L=2\), \(U_A=0\), and \(U_Y=0\). When \(A_0=0\), \(L_1=4\), \(A_1=2\), and \(Y=56\). When \(A_0=1\), \(L_1=6\), \(A_1=3\), and \(Y=67\). The difference is \(67-56=11\) points. The coefficient 3 is the contrast holding later history fixed; the 11-point total effect includes the changes in both \(L_1\) and \(A_1\).
Interactions Change the Conditional Contrast
A coefficient is not even a single common conditional contrast when the regression includes a treatment-by-covariate interaction. In that case the treatment contrast depends on the covariate value. The coefficient on treatment alone gives the contrast at the covariate’s reference value, not a universal effect across all covariate levels.
Worked Example: A Treatment Interaction After Adjustment
Suppose a hypothetical outcome model, measured in recovery-score points, is:
Assume this conditional mean is correctly specified. The treatment coefficient 4 is the contrast at \(L=5\), because the interaction term is zero there. At \(L=7\), the contrast is instead:
If \(L\) is affected by treatment, neither conditional contrast automatically equals the total effect. For a concrete structural illustration, suppose the covariate takes \(L=4\) under \(A=0\) and \(L=6\) under \(A=1\), with no other sources of variation. The outcome under \(A=0,L=4\) is \(60+3(4-5)=57\). Under \(A=1,L=6\) it is \(60+4+3(6-5)+2(6-5)=69\). The total difference is \(69-57=12\) points. The model’s treatment coefficient of 4 is the conditional contrast at \(L=5\), while the total comparison uses different covariate values under the two treatment levels.
When Can a Causal Interpretation Be Defended?
A fitted coefficient always has a model-based statistical interpretation, but a causal interpretation requires more. At minimum, the treatment and covariate interventions must be well defined, the model must adequately describe the relevant conditional means, and the data must contain support for the treatment comparisons being made. Unmeasured common causes of treatment and outcome can undermine a causal interpretation. With a post-treatment covariate, additional assumptions are needed to justify conditioning on that covariate for the particular causal contrast of interest.
Even when a conditional coefficient does have a causal interpretation under a carefully specified model and assumptions, it is a conditional effect: it describes a comparison with the covariate fixed. That is a different target from a total effect that allows the covariate to respond to treatment. Naming the target before interpreting the coefficient helps prevent these two questions from being conflated.
Common Mistake
Incorrect reasoning: “The regression includes the later covariate, so the treatment coefficient controls confounding and is the treatment’s effect on the outcome.”
Correction: The coefficient compares treatment values conditional on the included covariate values. If the covariate was affected by earlier treatment, that comparison holds part of the treatment history fixed and generally leaves out pathways through it. Whether the coefficient has any causal interpretation also depends on the design, assumptions, and model; inclusion in a regression is not, by itself, proof of confounding control.
Check Your Understanding
Answer each question by identifying the contrast and the assumptions needed for a causal interpretation.
- In an additive regression of \(Y\) on binary \(A\) and a post-treatment covariate \(L\), what does the coefficient of \(A\) compare?
- Suppose treatment changes \(L\) by 2 units, and each unit of \(L\) changes the outcome by 3 units. In a no-interaction linear structural model, what pathway contribution is excluded when \(L\) is held fixed?
- Why does adjusting for a covariate affected by \(A_0\) change the question answered by the coefficient of \(A_0\)?
- In a model with an \(A\times L\) interaction, what does the treatment coefficient alone represent when \(L\) is centered at 5?
- Does a regression coefficient become a causal effect merely because the model includes measured covariates? Explain briefly.