Tutorials › Generalized Linear Models › The Identity Link and Gaussian Linear Regression

Gaussian GLMs and Weights · Tutorial 327 of 981

The Identity Link and Gaussian Linear Regression

Learn how the identity link makes the conditional mean a linear predictor, and how to interpret regression coefficients and fitted means in context.

Beginner 9 min read

What You'll Learn

  • Express the Gaussian identity-link model as a linear conditional mean.
  • Distinguish a conditional mean from an observed response and a fitted mean.
  • Interpret the intercept and slopes while holding other predictors fixed.
  • Calculate fitted means from supplied regression coefficients.
  • Explain how centering a predictor changes the intercept but not fitted means.
  • Recognize why identity-link predictions can fall outside a response’s restricted range.

From a Gaussian GLM to a Linear Mean

In a Gaussian generalized linear model, the link connects the conditional mean of a response to a linear predictor. With the identity link, that connection is especially direct: the conditional mean itself equals the linear predictor. This is the mean structure of ordinary Gaussian linear regression.

Let \(Y_i\) be the response for observation \(i\), and let \(x_i\) be a \(p\)-dimensional column vector of predictors, possibly including a constant term. The coefficient vector is \(\beta\). For the ordinary, unweighted Gaussian model, assume the responses are conditionally independent and

$$ Y_i \mid x_i \sim N(\mu_i,\phi), \qquad \phi>0, $$

where \(\phi\) is the common conditional variance. The identity link is \(g(\mu)=\mu\). The linear predictor is \(\eta_i=x_i^T\beta\), so the model specifies \(g(\mu_i)=\eta_i\).

Definition: A Gaussian identity-link model with predictor vector \(x_i\) and coefficient vector \(\beta\) has conditional mean \(\mu_i=E[Y_i\mid x_i]=x_i^T\beta\). In the ordinary model with common dispersion \(\phi\), it also has \(\operatorname{Var}(Y_i\mid x_i)=\phi\).

The link determines the scale on which the predictors combine. For the identity link, each term in \(x_i^T\beta\) contributes additively to the expected response. This statement concerns the conditional mean, not the observed response: individual observations can lie above or below that mean because of their random Gaussian errors.

What the Coefficients Mean

In the familiar multiple-regression form, suppose the predictor vector contains an intercept and \(q\) measured predictors:

$$ x_i=(1,x_{i1},\ldots,x_{iq})^T, \qquad E[Y_i\mid x_i]=\beta_0+\beta_1x_{i1}+\cdots+\beta_qx_{iq}. $$

The intercept \(\beta_0\) is the conditional mean when all predictors equal zero. A slope \(\beta_j\) describes the change in conditional mean associated with a one-unit increase in \(x_{ij}\), with the other predictors held fixed. This is a statement about the model’s mean function; by itself, it does not establish that changing a predictor causes a response change.

The intercept depends on how predictors are represented. If zero is far outside the range of observed values, the intercept may have little practical meaning, even though it is needed to specify the line or plane. Centering a predictor can make the intercept describe the mean at a more useful reference value. Indicator variables also give slopes a comparison interpretation: with a reference category coded as zero, an indicator coefficient is the difference in conditional means from that reference, when other predictors are held fixed.

Worked Example: Interpreting a Slope in a Simple Regression

Suppose a model for weekly study hours \(Y_i\), in hours, uses class attendance \(x_i\), measured as a percentage from 0 to 100: \(E[Y_i\mid x_i]=2+0.08x_i\). Interpret the intercept and slope, and calculate the conditional mean at 75 percent attendance.

Intercept. At \(x_i=0\), the model’s conditional mean is \(2+0.08(0)=2\) study hours per week.

Slope. A one-percentage-point increase in attendance is associated with a \(0.08\)-hour increase in the model’s conditional mean, since the predictor changes by one unit. For a ten-point increase, the mean changes by \(10(0.08)=0.8\) hours.

Conditional mean at 75 percent. Substitution gives \(2+0.08(75)=2+6=8\) hours per week. This is the expected response at that predictor value under the model, not a guarantee that every student attending 75 percent of classes studies exactly eight hours.

Fitted Means and Observed Responses

The coefficients in practice are usually estimated from data. Once coefficient estimates \(\hat{\beta}\) are available, the fitted mean for observation \(i\) is the model’s estimated conditional mean at that observation’s predictors:

$$ \hat{\mu}_i=x_i^T\hat{\beta}. $$

The fitted mean is not the observed response \(y_i\). Their difference, \(y_i-\hat{\mu}_i\), records how far that particular observation lies from its fitted mean. The identity link makes this plug-in calculation simple: multiply each predictor value by its corresponding coefficient estimate and add the terms, including the intercept.

Worked Example: Computing Fitted Means in Multiple Regression

Consider the estimated mean model \(\hat{\mu}_i=12+1.5x_{i1}-2x_{i2}\), where \(x_{i1}\) is a score in tens of points and \(x_{i2}\) is an indicator equal to 1 for a treatment group and 0 for a reference group. Compute fitted means for (a) score 6 in the reference group and (b) score 8 in the treatment group.

For (a), the score value in the model’s units is \(x_{i1}=6\), and the reference-group indicator is \(x_{i2}=0\). Thus \(12+1.5(6)-2(0)=12+9=21\). The fitted mean is 21 response units.

For (b), \(x_{i1}=8\) and \(x_{i2}=1\). Thus \(12+1.5(8)-2(1)=12+12-2=22\). The fitted mean is 22 response units.

At a fixed score, the treatment-group indicator changes the fitted mean by \(-2\) response units relative to the reference group. For example, at score 8 the reference-group fitted mean would be \(12+1.5(8)=24\), while the treatment-group fitted mean is \(22\). This comparison follows from the chosen indicator coding.

Two Useful Consequences of the Identity Link

Theorem: In a Gaussian identity-link model with \(\eta_i=x_i^T\beta\), the conditional mean is \(E[Y_i\mid x_i]=x_i^T\beta\).

Proof. By definition of the link in a GLM, \(g(\mu_i)=\eta_i\), where \(\mu_i=E[Y_i\mid x_i]\). For the identity link \(g(\mu)=\mu\), so \(\mu_i=\eta_i\). Since \(\eta_i=x_i^T\beta\), it follows that \(E[Y_i\mid x_i]=x_i^T\beta\). \(\square\)

Thus the linear predictor is not merely a convenient quantity that is transformed to obtain the mean. With the identity link, it is the mean itself. A second consequence gives a precise interpretation for any comparison between two predictor settings.

Theorem: For any two predictor vectors \(x\) and \(x^\ast\) in the model, the difference in conditional means is \((x^\ast-x)^T\beta\).

Proof. The identity-link mean function at \(x\) is \(\mu(x)=x^T\beta\), and at \(x^\ast\) it is \(\mu(x^\ast)=(x^\ast)^T\beta\). Subtracting gives \(\mu(x^\ast)-\mu(x)=(x^\ast)^T\beta-x^T\beta=(x^\ast-x)^T\beta\). This establishes the stated comparison. \(\square\)

When only predictor \(j\) changes by one unit and the others remain fixed, \(x^\ast-x\) has a one in position \(j\) and zeros elsewhere. The theorem then gives a mean difference of \(\beta_j\). If several predictors change, their coefficient-weighted changes add. This is the algebraic basis for interpreting slopes and for calculating fitted-mean contrasts.

Worked Example: Centering a Predictor Changes the Intercept

Suppose \(E[Y\mid x]=10+3x\), where \(x\) is a measurement with reference value 5. Re-express the same mean using the centered predictor \(z=x-5\), and interpret the new intercept.

Since \(x=z+5\), substitution gives \(10+3x=10+3(z+5)=10+15+3z=25+3z\). Therefore the same mean function can be written \(E[Y\mid z]=25+3z\).

When \(z=0\), the original predictor is \(x=5\), so the new intercept 25 is the conditional mean at \(x=5\). The slope remains 3: increasing \(z\) by one unit is equivalent to increasing \(x\) by one unit.

For a check, at \(x=7\) the original expression gives \(10+3(7)=31\). The centered value is \(z=7-5=2\), and the re-expressed model gives \(25+3(2)=31\). The predicted mean is unchanged; only the reference point used to describe the intercept has changed.

Why the Identity Link Matters—and a Common Pitfall

The identity link preserves the response’s measurement scale. A slope can therefore be stated directly in response units per predictor unit, and the fitted mean is obtained by an ordinary linear combination. This makes the model a natural choice when additive changes in the expected response are scientifically meaningful and a Gaussian response model is appropriate.

A limitation follows from the same feature. A linear combination \(x_i^T\beta\) can be any real number, so the model does not constrain its mean predictions to be positive or to lie in a bounded interval. If the response is nonnegative or bounded, an identity-link fitted mean can fall outside its permissible range. That possibility should prompt a check of whether the model is suitable for the response and the predictor values being considered; it is not fixed merely by interpreting coefficients carefully.

  • Common mistake: treating the intercept as the average response. The intercept is the model’s conditional mean when all predictors are zero. It equals an overall average only under particular designs or parameterizations; the definition alone does not make it an average over the observed sample.
  • Correction: specify the predictor setting. Interpret the intercept at the stated zero values, or center predictors so that zero represents a meaningful reference setting.
  • Common mistake: describing a fitted mean as an observed outcome. A fitted mean is \(x_i^T\hat{\beta}\), an estimated conditional expectation. The observed response can differ from it.
  • Correction: retain the conditional language. Describe \(\hat{\mu}_i\) as the estimated mean response for observations with those predictor values, rather than as the response every such observation will have.
Key takeaway: Under the Gaussian identity link, the conditional mean equals \(x_i^T\beta\). Coefficients describe additive changes in that mean under specified predictor changes, and fitted means are obtained by substituting estimated coefficients into the same expression.

Check Your Understanding

Answer each question using the conditional-mean interpretation.

  1. In a model with an intercept and two predictors, what does the coefficient of the second predictor describe when the other predictor is held fixed?
  2. A fitted model is \(\hat{\mu}=4+0.5x\). Calculate and interpret the fitted mean at \(x=6\).
  3. Why is a fitted mean not the same quantity as an observed response?
  4. For the model \(E[Y\mid x]=7+2x\), what happens to the intercept if the predictor is recoded as \(z=x-3\)? Show the algebra.
  5. What limitation can arise when an identity-link model is used for a response that must be nonnegative?
Clinical Biostats