Tutorials › Regression Analysis (Undergraduate) › How Specification Changes Alter the Target Association

Specification Checks and Residual Patterns · Tutorial 587 of 991

How Specification Changes Alter the Target Association

Learn to translate model terms into predicted outcome contrasts, and to identify which association a coefficient represents after a specification change.

Advanced 9 min read

What You'll Learn

  • Distinguish a model coefficient from the outcome contrast it represents.
  • Interpret a quadratic term using slopes and finite changes over specified predictor values.
  • Interpret an interaction by comparing group-specific slopes and matched predictions.
  • Explain what an adjustment-term coefficient means while holding other predictors fixed.
  • Compare specifications using the same predictor values and covariate settings.
  • Avoid treating a changed regression association as evidence of a changed causal effect.

Why the Same Coefficient Can Mean Something Different

A regression coefficient is not a description of an association in isolation. Its meaning depends on the model terms around it and on the comparison it implies. A coefficient of 2, for example, might describe a constant change in the fitted outcome per unit of a predictor, a slope only for a reference group, or a conditional slope holding another variable fixed. When a specification changes, a useful question is not only “Did the coefficient move?” but also “What comparison does each model now describe?”

The clearest way to answer that question is to compare matched predictions: predictions for two predictor settings that differ in the feature of interest but agree on other settings. Their difference states the association described by the model in context. This approach works whether the change involves adding curvature, an interaction, or an adjustment term.

Definition: A regression model’s target association is the outcome comparison the model is being used to describe, such as a change in predicted outcome across two values of a predictor. A matched prediction contrast compares predictions at two settings that differ in the predictor of interest and hold the specified comparison variables fixed. The contrast is model-dependent: changing the model can change the fitted predictions and the meaning of a coefficient.

Write \(\hat m(\mathbf{x})\) for a model’s fitted mean outcome at predictor settings \(\mathbf{x}\). For two settings that differ in predictor \(x\), a matched contrast is \( \hat m(x_1,\mathbf{z})-\hat m(x_0,\mathbf{z}) \), where \(\mathbf{z}\) denotes the other settings held fixed. It is measured in outcome units. A coefficient may equal this contrast, a contrast divided by a predictor change, or a slope at a particular setting; it need not represent one constant comparison in every model.

Curvature Changes a Constant Slope into a Setting-Specific One

In a linear mean model, the fitted mean changes at a constant rate as \(x\) changes, provided the other terms are held fixed. Adding a quadratic term allows that rate to vary with \(x\). Consider the two models

$$ \hat m_1(x)=\hat\beta_0+\hat\beta_1x \qquad\text{and}\qquad \hat m_2(x)=\hat\beta_0+\hat\beta_1x+\hat\beta_2x^2. $$

In the first model, \(\hat\beta_1\) is the predicted outcome change per one-unit increase in \(x\). In the quadratic model, it is the slope at \(x=0\), not a constant slope over the full range. At a general value \(x\), the slope is \(\hat\beta_1+2\hat\beta_2x\). For two values \(x_0\) and \(x_1\), the fitted change is

$$ \hat m_2(x_1)-\hat m_2(x_0) =(x_1-x_0)\left[\hat\beta_1+\hat\beta_2(x_1+x_0)\right]. $$

The change per unit between those two values is therefore \(\hat\beta_1+\hat\beta_2(x_1+x_0)\), which depends on the comparison’s endpoints. The derivative gives the local slope; a finite matched contrast gives the predicted change over a specific interval. These are related, but they are not interchangeable.

Worked Example: Comparing a Linear and a Quadratic Fit

Suppose a hypothetical analysis predicts a response score from practice time \(x\), measured in hours. A quadratic fit is \( \hat m_2(x)=12+2x-0.25x^2 \). The observed practice times include values from 4 to 9 hours, so the comparison below is within that range. What does the model predict when practice time increases from 6 to 8 hours?

First calculate the two fitted means: \( \hat m_2(6)=12+2(6)-0.25(6^2)=12+12-9=15 \) points, and \( \hat m_2(8)=12+2(8)-0.25(8^2)=12+16-16=12 \) points.

The matched contrast is \(12-15=-3\) points over 2 hours. Thus, for this model and this interval, the predicted score decreases by 3 points, or an average of \(-3/2=-1.5\) points per additional hour over the interval. The local slope at \(x=7\) is \(2+2(-0.25)(7)=2-3.5=-1.5\) points per hour. The equality here follows because 7 is the midpoint of 6 and 8 for a quadratic curve.

If a separate linear fit had estimated a slope of \(0.4\) points per hour, that coefficient would describe a constant fitted rate in the linear model, not the quadratic model’s rate around 6 to 8 hours. The two coefficients do not represent identical summaries. Comparing matched predictions over the same 6-to-8-hour interval makes the models’ implications directly understandable.

The coefficient \(\hat\beta_1\) in a quadratic model is also sensitive to where \(x=0\) is placed. If zero is far outside the data range, its “slope at zero” may be of little practical interest. Centering \(x\) at a meaningful value changes the interpretation of the linear coefficient, although it does not change the fitted curve when the model includes both the linear and quadratic terms.

Interactions Make Slopes Conditional on Another Predictor

An interaction allows the association between \(x\) and the outcome to vary with another predictor. With a binary group indicator \(G\), where \(G=0\) is the reference group and \(G=1\) is the comparison group, a model with an interaction is

$$ \hat m(x,G)=\hat\beta_0+\hat\beta_1x+\hat\beta_2G+\hat\beta_3xG. $$

For \(G=0\), the fitted slope in \(x\) is \(\hat\beta_1\). For \(G=1\), it is \(\hat\beta_1+\hat\beta_3\). The coefficient \(\hat\beta_3\) is the difference between those slopes. Meanwhile, \(\hat\beta_2\) is the group difference at \(x=0\); it is not necessarily the group difference at other values of \(x\).

Worked Example: Interpreting a Group-by-Predictor Interaction

In a hypothetical setting, a model predicts an outcome score from a dose \(x\) and a group indicator \(G\): \( \hat m(x,G)=40+3x+8G+2xG \). Dose is measured in units and is observed from 2 to 7 units in both groups. Interpret the slopes and compare the predicted changes from 4 to 6 units.

For the reference group, set \(G=0\): \( \hat m(x,0)=40+3x \). Its slope is 3 outcome points per dose unit.

For the other group, set \(G=1\): \( \hat m(x,1)=40+3x+8+2x=48+5x \). Its slope is 5 points per dose unit. The interaction coefficient, 2, is the difference in slopes: \(5-3=2\) points per dose unit.

For the reference group, the predictions at 4 and 6 units are \( \hat m(4,0)=40+3(4)=52 \) and \( \hat m(6,0)=40+3(6)=58 \). The matched change is \(58-52=6\) points. For the other group, they are \( \hat m(4,1)=48+5(4)=68 \) and \( \hat m(6,1)=48+5(6)=78 \), a change of \(78-68=10\) points.

The difference between these two changes is \(10-6=4\) points over the 2-unit increase, or 2 points per unit. This equals the interaction coefficient multiplied by the dose change: \(2(6-4)=4\). The coefficient 3 alone describes the reference-group slope, not a common slope for both groups.

With an interaction, a group comparison also depends on the chosen value of \(x\). In this example, the fitted group difference at dose \(x\) is \(8+2x\). It is 16 points at \(x=4\) and 20 points at \(x=6\). Naming the comparison value is therefore part of stating what the model says.

Adjustment Terms Define What Is Held Fixed

Adding an adjustment predictor changes the comparison represented by a coefficient. In a model such as \( \hat m(x,z)=\hat\beta_0+\hat\beta_1x+\hat\beta_2z \), the coefficient \(\hat\beta_1\) describes the fitted change per unit of \(x\) when \(z\) is held fixed. A model without \(z\) does not make that same conditional comparison. Its coefficient summarizes the fitted association between \(x\) and the outcome without explicitly holding \(z\) fixed.

The two coefficients can differ because the models use different information to describe the association. Such a change does not, by itself, show that either model is correct or that the association is causal. Interpretation also depends on the question, the predictor values that are meaningfully comparable, and whether the model’s mean structure is suitable for the data.

Worked Example: A Conditional Association After Adjustment

Suppose a hypothetical score is modeled using study sessions \(x\) and a baseline score \(z\), both measured in points or sessions as appropriate. The adjusted fitted model is \( \hat m(x,z)=10+2x+0.5z \). A simpler model without baseline score has an estimated slope of 3.2 points per session. Compare the fitted implications for 2 versus 5 sessions at a fixed baseline score of 60.

In the adjusted model, the prediction at \(x=2,z=60\) is \(10+2(2)+0.5(60)=10+4+30=44\) points. At \(x=5,z=60\), it is \(10+2(5)+0.5(60)=10+10+30=50\) points. The matched contrast is \(50-44=6\) points for three additional sessions, or 2 points per session. It is the same at any fixed \(z\) because this model has no \(x z\) interaction.

The simpler model’s slope of 3.2 implies a fitted increase of \(3.2(5-2)=9.6\) points over the same session difference, but it does not hold baseline score fixed. Its slope and the adjusted slope summarize different model comparisons. The numerical difference alone does not establish why they differ or which association is more relevant to a particular question.

These fitted associations are not automatically causal effects. For example, students who choose more sessions might differ in other ways, and the model equations alone do not establish that changing study time would change an individual’s score by the stated amount.

A Practical Rule for Comparing Specifications

When a specification changes, first state the target comparison in words, then calculate it from each model’s predictions using the same settings wherever possible. This makes the comparison explicit even when individual coefficients no longer align.

1
Name the outcome and comparison.
State which predictor changes, the two values being compared, and the outcome units.
2
Set the other predictors.
Specify the group, adjustment-variable values, or other settings held constant. If the model includes an interaction, these settings may change the contrast.
3
Calculate both predictions.
Substitute each predictor setting into the model and subtract the first fitted mean from the second.
4
Describe the model-specific result.
Report the predicted outcome change, its direction, and its units. Do not label it causal unless the study design and assumptions justify that claim.

For comparisons across fitted models, use the same outcome definition, observations, predictor units, and predictor settings when possible. If the models use different samples or units, a coefficient difference may reflect those changes as well as the changed terms. Matched predictions do not make two models equivalent; they make the comparison being drawn easier to see.

Common Mistake: Treating a Coefficient as a Universal Effect

  • Wrong reasoning: “The coefficient for \(x\) is 2, so increasing \(x\) by one unit always increases the outcome by 2.”
  • Correction: Check the model. With curvature, the slope depends on \(x\). With an interaction, it depends on the interacting predictor. With adjustment terms, it is conditional on the stated values of those predictors.
  • Wrong reasoning: “The coefficient changed after adjustment, so the adjusted model proves the true effect.”
  • Correction: A coefficient change shows that the fitted association depends on the specification. It does not alone identify which model is appropriate or establish causation.
Key takeaway: A regression coefficient gets its meaning from the full model. To explain a specification change, identify the comparison of interest, hold the relevant settings fixed, and use matched predictions to report the model’s implied outcome contrast in context.

Check Your Understanding

For each question, state what additional setting or calculation is needed to interpret the model’s association precisely.

  1. In a quadratic model for outcome against age, why is the coefficient of age not generally a constant slope across all ages?
  2. A model includes a group-by-dose interaction. Which slope does the dose coefficient represent, and how do you find the slope for the other group?
  3. A fitted model predicts outcome from exposure and baseline severity. What must be held fixed to interpret the exposure coefficient as a conditional association?
  4. For an interaction model, describe how you would compare the groups at a specified exposure value using matched predictions.
  5. Why does a coefficient changing after adding adjustment terms not, on its own, establish that one specification is correct or that the association is causal?