Why Compare Adjustment Sets?
In multiple regression, an adjusted association describes the relationship between an exposure and an outcome after accounting for selected covariates. Choosing those covariates is not just a matter of adding every available variable. A variable may be a common cause of exposure and outcome, a mediator on the pathway between them, or a common effect of other variables. Those roles matter: adjusting for a common cause may help address confounding, while adjusting for a mediator can change the association of interest. Conditioning on a common effect can also introduce bias.
Sometimes subject-matter knowledge and the timing of measurements do not settle a variable’s role. In that situation, it is useful to identify a small number of plausible adjustment sets, state the assumptions behind each, and compare the exposure coefficient from each fitted model. The comparison documents how the adjusted association depends on the assumptions; it does not, by itself, tell us which assumptions are true.
What the Coefficient Means
For a continuous outcome \(Y\), an exposure \(X\), and covariates \(Z_1,\ldots,Z_k\), an additive multiple regression model can be written as:
Here, \(\beta_X\) is the model’s exposure coefficient. Under this model, it describes the difference in the fitted mean outcome associated with a one-unit difference in \(X\), comparing observations with the same included covariate values. It is an adjusted association. To interpret it as a causal effect requires further assumptions, including an appropriate adjustment set and adequate control of confounding; regression alone does not establish those assumptions.
For the coefficient interpretation to be appropriate, the model’s conditional mean must be represented adequately by its chosen form. In the equation above, that means a linear, additive relationship in the included predictors. The observations must also be analyzed on a consistent basis: comparisons of candidate sets should use the same outcome, exposure scale, and analysis sample. Otherwise, a coefficient difference could reflect a change in data or model specification unrelated to the covariate choice. Exact linear dependence among predictors must not prevent the coefficients from being separately estimated.
The coefficient is not required to stay the same across adjustment sets. It can change because the included variables account for different sources of association, because a variable lies on a pathway of interest, or because the model’s assumptions are imperfect. A changed estimate does not identify which explanation is correct.
Worked Example: Comparing Adjustment for a Possible Mediator
Suppose a hypothetical observational analysis examines the association between a treatment-dose indicator \(X\) (1 for the higher dose and 0 for the lower dose) and systolic blood pressure \(Y\), measured in mmHg. The analysis includes 240 people with complete observations. Investigators are unsure whether a follow-up measure of symptom improvement \(M\) should be adjusted for. They fit both models to the same 240 people, with the same outcome and exposure coding.
The hypothetical fitted coefficients are \(5.2\) mmHg for \(X\) without covariates, \(3.1\) mmHg for \(X\) after adjustment for baseline severity \(Z\), and \(1.4\) mmHg for \(X\) after adjustment for both \(Z\) and follow-up improvement \(M\).
The first adjusted model is \(Y_i=\beta_0+\beta_X X_i+\beta_Z Z_i+\varepsilon_i\). Its exposure coefficient is \(3.1\) mmHg: among observations with the same modeled baseline severity, the higher-dose group has a fitted mean systolic pressure \(3.1\) mmHg higher than the lower-dose group.
The second adjusted model also includes \(M\), and its exposure coefficient is \(1.4\) mmHg. Relative to the first adjusted estimate, the coefficient changes by \(1.4-3.1=-1.7\) mmHg. This arithmetic describes a change in the fitted adjusted association, not proof that the second model is better.
If follow-up improvement is caused partly by treatment and then affects blood pressure, it may be a mediator. Adjusting for it would then answer a different question from adjusting only for baseline severity. The study’s timing and subject-matter knowledge are needed to decide which interpretation is defensible. The two estimates should therefore be reported with their distinct assumptions, rather than choosing one solely because its coefficient is smaller.
Build Candidate Sets From Competing Assumptions
A useful comparison begins with the causal question and a diagram or verbal account of how the variables might relate. Earlier material on using subject-matter knowledge to select covariates and on distinguishing confounders, mediators, and prognostic variables provides the relevant concepts. Here the aim is to use them to organize a comparison when a variable’s role remains uncertain.
For each candidate set, write down what must be true for it to support the intended interpretation. For example, one set might adjust for a variable because it is assumed to precede and affect both exposure and outcome. A second might omit that variable because it is assumed to occur after exposure and lie on a pathway to the outcome. If the variable’s timing or causal role is unknown, neither choice should be presented as settled fact.
Keep the comparison focused. Each model should use the same outcome, exposure coding, and observations, and should differ only in the covariates whose roles are genuinely in question. This makes the comparison easier to interpret. It also avoids treating a broad collection of arbitrary models as if it were a principled analysis.
Worked Example: Baseline or Post-Exposure Body Mass Index?
Consider a hypothetical study of the association between a dietary program \(X\) (enrolled or not enrolled) and a later cholesterol measurement \(Y\), in mg/dL. Body mass index \(B\) was recorded, but its timing relative to program enrollment is unclear. The same 180 people are used for both fitted models. A model adjusting for age \(A\) gives an exposure coefficient of \(-6.4\) mg/dL. A model adjusting for age and \(B\) gives \(-4.1\) mg/dL.
In the first model, the coefficient is the fitted difference in mean cholesterol between enrolled and non-enrolled people with the same modeled age: \(-6.4\) mg/dL. In the second, it is the fitted difference for people with the same modeled age and body mass index: \(-4.1\) mg/dL. The coefficient changes by \(-4.1-(-6.4)=2.3\) mg/dL toward zero.
If \(B\) was measured before enrollment and was a common cause of enrollment and cholesterol, adjusting for it may be important to address confounding. If it was measured after enrollment and was affected by the program, it could instead be a mediator; adjusting for it would remove part of the pathway that contributes to the overall association. Because the timing is uncertain, the numerical comparison cannot resolve the role of \(B\). The defensible conclusion is conditional: the estimate is \(-6.4\) mg/dL under the age-only specification and \(-4.1\) mg/dL under the age-and-\(B\) specification, with the latter requiring a different assumption about what adjustment is appropriate.
Interpreting Differences Without Overclaiming
When the exposure coefficient changes across candidate sets, describe the change in the coefficient’s units and explain which covariate was added or removed. Then connect each model to its assumptions. Avoid saying that adjustment “revealed the true effect” or that the model with the most covariates is automatically least biased. The comparison is informative because it makes dependence on assumptions visible, not because regression can determine the causal structure from the coefficient estimates alone.
It can be tempting to select whichever model has the largest \(R^2\), the smallest p-value, or an exposure coefficient that matches expectations. These criteria do not establish that a set is causally appropriate. Adding predictors can improve prediction of the outcome while including a mediator or common effect that is unsuitable for the intended question. Likewise, an estimate’s direction or magnitude cannot prove that the model’s adjustment choices are correct.
Worked Example: Recording a Coefficient Comparison
A hypothetical analysis estimates the association between a one-hour-per-week increase in exercise \(X\) and resting heart rate \(Y\), in beats per minute. All three models use the same 210 observations. The candidate models and fitted exposure coefficients are:
| Model | Adjustment set | Fitted coefficient for exercise |
|---|---|---|
| A | None | \(-0.8\) beats per minute per exercise hour/week |
| B | Age and baseline health status | \(-1.3\) beats per minute per exercise hour/week |
| C | Age, baseline health status, and follow-up fitness | \(-0.6\) beats per minute per exercise hour/week |
Model B’s coefficient is \(-1.3\) beats per minute per additional exercise hour per week, conditional on age and baseline health status. Compared with model A, the difference is \(-1.3-(-0.8)=-0.5\) beats per minute per exercise hour per week. Model C’s coefficient differs from model B’s by \(-0.6-(-1.3)=0.7\) in the same units.
Suppose follow-up fitness may be influenced by exercise and may also be related to heart rate. Then model C’s adjustment set may condition on a post-exposure variable, while model B’s does not. If the target is an association that includes pathways through fitness, model B may align more closely with that target, subject to its other assumptions. The table alone cannot establish that follow-up fitness is a mediator or that model B is unbiased. A careful report would give the estimates, state the timing and role assumptions, and explain why the chosen set is the primary one.
A Practical Comparison and a Common Pitfall
A concise comparison can be organized in four steps:
Name the exposure, outcome, units, and whether the intended quantity includes or excludes pathways through possible mediators.
List the covariates in each model and state the variable-role assumptions that motivate including or excluding them.
Use the same observations, outcome, exposure scale, and model form so the coefficient comparison isolates the covariate choices.
Give each exposure coefficient in context, describe how it changes, and identify what cannot be resolved from the regression comparison alone.
Common mistake: “The model with the most covariates must control the most confounding, so I will report that model.” This reasoning treats all covariates as harmless adjustments. A post-exposure mediator or common effect can make an expanded set inappropriate for the target question. The correction is to justify each covariate’s inclusion using its timing and plausible causal role, then present alternative estimates when reasonable assumptions remain unsettled.
A clear report need not imply that all candidate sets are equally plausible. It can identify one set as primary because it best matches the available subject-matter knowledge and the target association, while showing how alternatives change the estimate. The limitation should be stated plainly: if the variable roles cannot be determined from the available information, the adjusted association remains dependent on those assumptions.
Check Your Understanding
Use the ideas in this tutorial to answer the following questions.
- Why should candidate adjustment models generally use the same analysis sample and exposure coding?
- A coefficient changes after adding a variable measured after exposure. What information is needed before deciding whether that adjustment is appropriate?
- In a model with an exposure and covariates, what does the exposure coefficient describe under the additive linear model?
- Why does a larger \(R^2\) not show that an adjustment set is more appropriate for a causal question?
- What should be documented when two plausible adjustment sets produce different exposure coefficients?