Tutorials › Generalized Linear Models › Identity Links And Additive Mean Changes

Link Functions and Scale · Tutorial 110 of 981

Identity Links And Additive Mean Changes

Learn how identity-link coefficients translate into additive changes in conditional means, and how to check whether the resulting predictions are admissible.

Beginner 9 min read

What You'll Learn

  • Define the identity link and express the conditional mean through the linear predictor.
  • Calculate exact mean changes for continuous and indicator predictors.
  • Interpret coefficients in models with multiple predictors while holding other terms fixed.
  • Check whether identity-link predictions meet a response family’s mean constraints.
  • Identify how the modeled covariate range affects admissibility.

What an Identity Link Says About the Mean

A coefficient is easiest to interpret when we identify the scale on which it acts. In a generalized linear model, the linear predictor is built from the covariates, while the link connects that predictor to the conditional mean. With the identity link, the two are equal. Consequently, a change in the linear predictor is an equal-sized change in the mean—not a change that must first be transformed.

For observation \(i\), let \(\mathbf{x}_i\) be its design row and let \(\boldsymbol{\beta}\) be the coefficient vector. The linear predictor is \(\eta_i=\mathbf{x}_i^\mathsf{T}\boldsymbol{\beta}\). Writing the identity link as \(g(\mu)=\mu\), the model’s mean equation is \(\mu_i=\eta_i\). The distribution of the response still matters: it determines which values of \(\mu_i\) are admissible.

Definition: Under the identity link, \(g(\mu)=\mu\) and \(g^{-1}(\eta)=\eta\). Thus a generalized linear model with this link satisfies \(E[Y_i\mid\mathbf{x}_i]=\mu_i=\eta_i=\mathbf{x}_i^\mathsf{T}\boldsymbol{\beta}\), for covariate values where this mean is permitted by the response distribution.

The inverse link being the identity makes the link-scale and mean-scale values numerically equal. Therefore, in a model with a single untransformed predictor \(x\), such as \(\mu(x)=\beta_0+\beta_1x\), the coefficient \(\beta_1\) is the change in conditional mean for a one-unit increase in \(x\). More generally, a change of \(d\) units changes the mean by \(\beta_1d\), provided other predictors and model terms are held fixed.

Exact Additive Changes

The additive interpretation is an exact difference, not an approximation for small changes. In a model with several predictors, changing one design component while leaving the others unchanged gives a particularly direct comparison.

Result: Suppose \(\mu(\mathbf{x})=\mathbf{x}^\mathsf{T}\boldsymbol{\beta}\), and two design rows differ only in component \(j\), which changes from \(x_j\) to \(x_j+d\). Then the conditional mean changes by \(d\beta_j\).

Proof. Let \(\mathbf{x}\) denote the first design row and \(\mathbf{x}'\) the second. By assumption, \(x'_j=x_j+d\) and \(x'_k=x_k\) for every \(k\ne j\). Subtract the first conditional mean from the second:

$$ \mu(\mathbf{x}')-\mu(\mathbf{x}) =(\mathbf{x}'-\mathbf{x})^\mathsf{T}\boldsymbol{\beta} =d\beta_j. $$

This difference does not depend on the starting value \(x_j\). It does depend on what counts as changing only component \(j\): if changing a raw covariate also changes other design terms, the full design-row difference must be used instead. The claimed result is about a single design component, not automatically every possible change in the original covariates. \(\square\)

Worked Example: A Continuous Predictor

Suppose the modeled mean of a continuous outcome is \(\mu(x)=12+1.5x\), where \(x\) is a dose measured in units. Find the conditional mean at \(x=2\) and \(x=5\), and calculate the change.

At \(x=2\), the mean is \(\mu(2)=12+1.5(2)=15\). At \(x=5\), it is \(\mu(5)=12+1.5(5)=19.5\). The predictor increased by \(5-2=3\) units, so the model implies a mean change of \(1.5(3)=4.5\). Direct subtraction confirms it: \(19.5-15=4.5\). The coefficient \(1.5\) is the additive mean change per one-unit increase in dose.

For multiple predictors, a coefficient comparison is conditional on the other design components. If a predictor is centered, rescaled, or represented by more than one design column, interpretation must follow that representation. The additive change comes from multiplying the change in the design row by the coefficient vector.

Worked Example: A Continuous Predictor and a Treatment Indicator

Suppose a model for mean systolic blood pressure is \(\mu(a,T)=30+2.4(a-50)-6T\), where \(a\) is age in years and \(T=1\) denotes treatment while \(T=0\) denotes control. Compare the modeled means for a 60-year-old in each group, then compare a 50-year-old and a 60-year-old within the treatment group.

For a 60-year-old in control, \(a-50=10\) and \(T=0\), so \(\mu(60,0)=30+2.4(10)-6(0)=54\). For a 60-year-old receiving treatment, \(\mu(60,1)=30+2.4(10)-6(1)=48\). At this age, the treatment indicator’s change from 0 to 1 corresponds to an additive difference of \(-6\) in the modeled mean.

Within the treatment group, the 50-year-old has mean \(\mu(50,1)=30+2.4(0)-6=24\), whereas the 60-year-old has mean 48. The age increase is 10 years, and \(2.4(10)=24\), matching the difference \(48-24=24\). These are model-implied conditional mean comparisons; interpreting them causally would require additional assumptions about the study and treatment assignment.

Mean Changes Must Still Be Admissible

An additive mean equation does not automatically respect the response distribution. The identity inverse link returns any real-valued linear predictor unchanged, but a response family may allow only a restricted set of means. For example, a Bernoulli response with success probability \(p\) requires \(0\leq p\leq1\), and a Poisson response with mean \(\mu\) requires \(\mu\geq0\). A Gaussian response with unrestricted real support has no corresponding positivity or upper-bound constraint on its mean.

This distinction is important: the linear predictor may be calculated for any design row, but that calculation is a valid conditional mean only if it is in the mean space allowed by the response distribution. A coefficient can have the additive interpretation and yet produce an inadmissible prediction at some covariate values.

Result: Let \(\mathcal{C}\) be the set of covariate values at which a model is intended to apply, and let \(\mathcal{M}\) be the set of admissible means for its response distribution. Under the identity link, the model defines an admissible mean throughout \(\mathcal{C}\) if and only if \(\mathbf{x}^\mathsf{T}\boldsymbol{\beta}\in\mathcal{M}\) for every \(\mathbf{x}\in\mathcal{C}\).

Proof. Under the identity link, the proposed conditional mean at \(\mathbf{x}\) is exactly \(\mathbf{x}^\mathsf{T}\boldsymbol{\beta}\). By definition, a value is an admissible mean for the specified response distribution exactly when it belongs to \(\mathcal{M}\). Thus, for each \(\mathbf{x}\in\mathcal{C}\), the proposed mean is admissible if and only if \(\mathbf{x}^\mathsf{T}\boldsymbol{\beta}\in\mathcal{M}\). Requiring this for every covariate value in \(\mathcal{C}\) gives the stated condition. \(\square\)

In the special case of one predictor, \(\mu(x)=\beta_0+\beta_1x\), and a covariate range \([a,b]\), checking the two endpoints is enough when \(\mathcal{M}\) is an interval. For any \(x\in[a,b]\), write \(x=(1-t)a+tb\) for some \(t\in[0,1]\). Then \(\mu(x)=(1-t)\mu(a)+t\mu(b)\). If both endpoint means belong to an interval \(\mathcal{M}\), every such weighted average belongs to \(\mathcal{M}\) as well. This argument applies, for instance, to the Bernoulli mean interval \([0,1]\).

Worked Example: An Identity-Link Bernoulli Model

Consider a Bernoulli response and the proposed identity-link mean \(p(x)=0.15+0.22x\), intended for covariate values \(0\leq x\leq4\). Check whether the model gives valid Bernoulli means throughout that range.

At the left endpoint, \(p(0)=0.15+0.22(0)=0.15\), which is in \([0,1]\). At the right endpoint, \(p(4)=0.15+0.22(4)=0.15+0.88=1.03\), which is greater than 1. Thus the model is not admissible on the stated range. In particular, at \(x=4\), the linear predictor gives \(1.03\), but a Bernoulli distribution with success probability \(1.03\) does not exist.

The coefficient \(0.22\) still describes the proposed additive change in the mean for a one-unit increase in \(x\), wherever the model gives an admissible mean. That interpretation does not make the prediction at \(x=4\) a valid probability. Restricting the covariate range or imposing coefficient constraints could change the admissibility issue; simply calling the output a fitted mean cannot.

Checking the Scope of an Additive Interpretation

The identity link is especially transparent, but its mean-scale consequences should be checked across the covariate region relevant to the analysis. A model may produce valid means for observed design rows and invalid means for new rows. Conversely, a model that is admissible over a narrow observed range may fail under extrapolation. The check concerns the mean equation and the intended covariate domain, not just the values encountered in the fitting data.

For response families whose means are unrestricted real numbers, this particular constraint check imposes no additional bounds. For restricted mean spaces, verify the fitted or proposed means against the family’s requirements. For a single linear predictor and interval covariate range, endpoint checks can suffice as shown above; with multiple predictors, the relevant check must cover the intended design region.

Common Mistake

  • Wrong reasoning: “The identity link has an inverse that accepts every real number, so every identity-link prediction is a valid mean.”
  • Correction: The inverse link returns a real number, but the response distribution may permit only a subset of real numbers as means. For Bernoulli outcomes, a predicted mean must be between 0 and 1. Check the response family’s mean space over the covariate values where the model will be used.
  • Wrong reasoning: “The coefficient is the additive effect of changing the original predictor, no matter how the model is written.”
  • Correction: The exact additive change is the coefficient vector multiplied by the change in the design row. A coefficient for a transformed or interaction column must be interpreted according to the change in that column; changing a raw predictor may alter several columns at once.
Key takeaway: With the identity link, the conditional mean equals the linear predictor, so changing a design component by \(d\) changes the mean by \(d\) times its coefficient. This is an exact additive statement, subject to the response distribution’s mean constraints.

Check Your Understanding

Use the mean equation and the response distribution together when answering.

  1. In \(\mu(x)=4-0.8x\), what is the modeled mean change when \(x\) increases by 2.5 units?
  2. A design component changes by 1, while all other design components stay fixed. What conditional mean change does its coefficient imply?
  3. For a Bernoulli response, why is a predicted identity-link mean of 1.04 inadmissible?
  4. If an affine mean function is admissible at both endpoints of an interval covariate range, why is it admissible throughout the interval when the allowed mean space is an interval?
  5. Why does an additive coefficient interpretation alone not establish that predictions are valid over an extrapolated covariate range?
Clinical Biostats