Tutorials › Generalized Linear Models › Link-Scale and Response-Scale Effects

Marginal Mean Models · Tutorial 822 of 981

Link-Scale and Response-Scale Effects

Translate coefficients on identity, log, and logit scales into interpretable changes in the response mean, and see when those changes depend on the baseline.

Advanced 9 min read

What You'll Learn

  • Define a link-scale contrast and its corresponding response-scale contrast
  • Show why an identity-link coefficient gives a constant mean difference
  • Interpret a log-link coefficient as a mean ratio and calculate the resulting difference
  • Interpret a logit-link coefficient as an odds ratio and calculate probabilities
  • Derive how response-scale changes depend on the baseline mean
  • Avoid confusing a constant coefficient with a constant response-scale difference

Why the Scale of an Effect Matters

A coefficient in a marginal mean model describes a change on the model’s link scale. Its practical interpretation depends on the link: the same kind of coefficient may represent a difference, a ratio, or an odds ratio. Even when a coefficient is constant on the link scale, the corresponding difference between mean responses need not be constant.

Here, \(\mu\) denotes a marginal response mean, such as the probability of an outcome in the population or the average count for a measurement occasion. We compare two covariate values while holding the other parts of the comparison fixed. This tutorial concerns translating that mean contrast between scales; it does not by itself give a causal interpretation to the contrast.

Definition: Let \(g\) be a link function and let \(\mu_0\) and \(\mu_1\) be response means at two covariate values. The link-scale contrast is \(g(\mu_1)-g(\mu_0)\). The response-scale contrast is \(\mu_1-\mu_0\). Depending on the link, a useful response-scale summary may instead be a ratio of means or a ratio of odds.

For the comparisons below, suppose a covariate changes by \(\delta\), and the model’s linear predictor changes by \(\beta\delta\). Thus the link-scale contrast is \(\beta\delta\). For a one-unit increase, \(\delta=1\), it is \(\beta\). This statement concerns the particular comparison being made; it does not assume that every covariate or every model has the same form.

Identity Link: The Coefficient Is a Mean Difference

The identity link is \(g(\mu)=\mu\), so its inverse is \(g^{-1}(\eta)=\eta\). If the linear predictor changes by \(\beta\delta\), the response mean changes by precisely that amount:

$$ \mu_1-\mu_0 = g^{-1}\{g(\mu_0)+\beta\delta\}-\mu_0 = \mu_0+\beta\delta-\mu_0 = \beta\delta. $$

The difference does not depend on the baseline mean. For example, a coefficient of \(-2.4\) means a decrease of 2.4 response units per one-unit increase in the covariate, for comparisons to which that coefficient applies. The identity link does not guarantee that all fitted means are valid for every outcome: a model for a nonnegative response, for instance, may require attention to whether its predicted means are nonnegative.

Worked Example: A Constant Difference Under the Identity Link

Suppose an identity-link marginal mean model gives a coefficient \(\beta=-2.4\) for a covariate measured in units of one. Compare covariate values \(x=3\) and \(x=4\), with the baseline mean at \(x=3\) equal to 12.5.

The covariate change is \(\delta=4-3=1\), so the link-scale contrast is \(\beta\delta=-2.4(1)=-2.4\). Because the link is the identity, this is also the response-scale difference:

$$ \mu(4)-\mu(3)=-2.4, \qquad \mu(4)=12.5-2.4=10.1. $$

If instead the baseline mean were 7.0, the same one-unit increase would give a mean of \(7.0-2.4=4.6\). The differences are both \(-2.4\), although the two pairs of means differ.

Log Link: A Constant Ratio, Not a Constant Difference

For a positive mean, the log link is \(g(\mu)=\log(\mu)\), with inverse \(g^{-1}(\eta)=e^\eta\). A link-scale contrast of \(\beta\delta\) implies a ratio of response means:

$$ \log(\mu_1)-\log(\mu_0)=\beta\delta \quad\Longrightarrow\quad \frac{\mu_1}{\mu_0}=e^{\beta\delta}. $$

Thus, for a fixed covariate change, the mean ratio is the same at different baseline means. The difference is instead \(\mu_1-\mu_0=\mu_0\{e^{\beta\delta}-1\}\), so it is proportional to the baseline. If \(\beta>0\), the mean increases; if \(\beta<0\), it decreases. The ratio formula is exact, not an approximation.

Worked Example: Turning a Log-Link Coefficient Into a Difference

Suppose a log-link model has coefficient \(\beta=0.3\). At a baseline covariate value, the marginal mean is \(\mu_0=4\). Calculate the mean at a covariate value one unit higher, and distinguish the ratio from the difference.

For \(\delta=1\), the mean ratio is \(e^{0.3}\approx1.3498588\). Therefore \(\mu_1=4e^{0.3}\approx4(1.3498588)=5.3994352\). The response-scale difference is \(5.3994352-4=1.3994352\), approximately 1.3994 units.

At a different baseline mean of 10, the same coefficient and covariate change give \(\mu_1=10e^{0.3}\approx13.498588\), a difference of about 3.4986. The ratio remains \(1.3498588\), but the difference changes with the baseline.

Logit Link: A Constant Odds Ratio, Not a Constant Probability Difference

For a binary response with mean \(\mu=P(Y=1)\), the logit link is \(g(\mu)=\log\{\mu/(1-\mu)\}\), where \(0<\mu<1\). Its inverse is the logistic function:

$$ g^{-1}(\eta)=\frac{e^\eta}{1+e^\eta} =\frac{1}{1+e^{-\eta}}. $$

A link-scale contrast of \(\beta\delta\) is a log odds ratio. Exponentiating gives the odds ratio \(e^{\beta\delta}\). This does not make the probability difference constant: the inverse link is nonlinear, and a given change in log odds yields different probability changes at different baseline probabilities.

Worked Example: Odds Ratio and Probability Difference

Suppose the baseline linear predictor is \(-1\), and a one-unit covariate increase has logit coefficient \(\beta=\log 2\). Find the baseline probability, the probability after the increase, and their difference.

At baseline, the odds are \(e^{-1}=1/e\), so the probability is \(\mu_0=1/(1+e)\approx0.2689414\). The coefficient is \(\log 2\), so the odds are multiplied by \(e^{\log 2}=2\). The new odds are \(2/e\), and the new probability is \(\mu_1=(2/e)/(1+2/e)=2/(e+2)\approx0.4238831\).

Using the unrounded model probabilities, the difference is

$$ \mu_1-\mu_0 =\frac{2}{e+2}-\frac{1}{e+1} \approx0.1549. $$

The increase is about 15.49 percentage points. Subtracting the displayed probabilities after rounding to four decimal places would give \(0.4239-0.2689=0.1550\); that is not the correctly rounded difference from the unrounded model values. The odds ratio is exactly 2, but the probability difference is not a fixed consequence of that odds ratio alone: it also depends on the baseline probability.

Deriving the Dependence on the Baseline

The examples illustrate a general calculation. Let \(\eta_0=g(\mu_0)\) be the baseline linear predictor, and suppose the covariate change shifts it by \(c=\beta\delta\). The response mean after the change is \(g^{-1}(\eta_0+c)\). Consequently, the response-scale difference is

$$ \Delta\mu =g^{-1}(\eta_0+c)-g^{-1}(\eta_0). $$

This identity follows directly from the definition of the inverse link. Whether \(\Delta\mu\) depends on the baseline is determined by the shape of \(g^{-1}\). For the identity link it does not; for the log and logit links, the difference generally varies with the baseline.

Theorem: Suppose \(g\) is one-to-one on the relevant range of response means, \(\mu_0\) is in that range, and the link-scale contrast between two covariate values is \(c\). Then \(\mu_1=g^{-1}\{g(\mu_0)+c\}\). For the identity link, \(\mu_1-\mu_0=c\). For the log link, when \(\mu_0>0\), \(\mu_1/\mu_0=e^c\). For the logit link, when \(0<\mu_0<1\), the odds at \(\mu_1\) divided by the odds at \(\mu_0\) equal \(e^c\).

Proof. By the definition of link-scale contrast, \(g(\mu_1)-g(\mu_0)=c\), so \(g(\mu_1)=g(\mu_0)+c\). Since \(g\) is one-to-one, applying its inverse gives \(\mu_1=g^{-1}\{g(\mu_0)+c\}\). Under the identity link, this is \(\mu_1=\mu_0+c\), hence the difference is \(c\). Under the log link, \(\log(\mu_1)=\log(\mu_0)+c\); exponentiation gives \(\mu_1=\mu_0e^c\), and division by the positive \(\mu_0\) gives the stated ratio. Under the logit link, write \(o_j=\mu_j/(1-\mu_j)\) for \(j=0,1\). The link equation gives \(\log(o_1)-\log(o_0)=c\). Exponentiating yields \(o_1/o_0=e^c\). Each conclusion follows from the specified link and the stated domain conditions. \(\square\)

A local change gives another view of why the baseline matters. If \(g^{-1}\) is differentiable, the derivative of the mean with respect to the linear predictor is \((g^{-1})'(\eta)\). Along a covariate change in which the linear predictor changes at rate \(\beta\), the derivative of the mean with respect to that covariate is \(\beta(g^{-1})'(\eta)\).

Result: If \(g^{-1}\) is differentiable and the linear predictor is \(\eta(x)=a+\beta x\) over the comparison of interest, then \(\frac{d\mu(x)}{dx}=\beta(g^{-1})'\{\eta(x)\}\). This is a local response-scale change per unit of \(x\); it need not equal the finite difference over a full unit.

Derivation. Since \(\mu(x)=g^{-1}\{\eta(x)\}=g^{-1}(a+\beta x)\), the chain rule gives \(\frac{d\mu(x)}{dx}=(g^{-1})'(a+\beta x)\frac{d(a+\beta x)}{dx}=\beta(g^{-1})'\{\eta(x)\}\). For the identity inverse link, the derivative is 1. For the log inverse link, it is \(e^\eta=\mu\), so the local change is \(\beta\mu\). For the logistic inverse link, it is \(\mu(1-\mu)\), so the local change is \(\beta\mu(1-\mu)\). These last two derivatives depend on the current mean. \(\square\)

Interpreting Effects Carefully

For marginal mean models, the means being compared are marginal means, not automatically subject-specific conditional means. The scale translation does not change which mean the model represents. Nor does a coefficient, by itself, establish that changing a covariate would cause the response to change: causal claims require additional design and assumptions.

When reporting an effect, name both the scale and the comparison. For example, “the mean is multiplied by \(e^\beta\) per unit” is a log-link interpretation; “the odds are multiplied by \(e^\beta\) per unit” is a logit-link interpretation. If a response-scale difference is important, calculate it at a stated baseline mean or probability rather than treating it as constant.

Common Mistake

Wrong reasoning: “The coefficient is \(\beta\), so the response mean increases by \(\beta\) everywhere.” This is valid for a coefficient contrast under the identity link, but not in general. Under a log link, \(\beta\) is a log mean ratio for a one-unit increase; under a logit link, it is a log odds ratio. Neither generally equals the change in the mean or probability.

Correction: Identify the link, transform the coefficient to its link-appropriate comparison, and use the baseline mean when calculating a response-scale difference. Report the covariate change and units as well: a coefficient per one unit cannot be interpreted as a per-ten-unit effect without multiplying the link-scale contrast by ten.

Key takeaway: Identity-link coefficients give constant mean differences. Log-link coefficients give constant mean ratios, and logit-link coefficients give constant odds ratios. For the latter two links, response-scale differences depend on the baseline mean.

Check Your Understanding

Use the link and baseline specified in each question; keep ratios distinct from differences.

  1. Under an identity link, a coefficient is \(1.7\) per unit. What is the mean difference for a covariate increase of three units?
  2. Under a log link, a coefficient is \(-0.2\) per unit. What mean ratio corresponds to a two-unit increase?
  3. For the log-link comparison in question 2, explain why the mean difference cannot be determined without a baseline mean.
  4. Under a logit link, what odds ratio corresponds to a coefficient of \(\log(1.5)\) for a one-unit increase?
  5. Why can the same logit coefficient yield different probability differences at different baseline probabilities?
Clinical Biostats