Tutorials › Generalized Linear Models › Admissibility of Identity and Inverse Links

Choosing Continuous Families · Tutorial 641 of 981

Admissibility of Identity and Inverse Links

Check whether identity and reciprocal links keep fitted means in the allowable positive range across the covariates of interest.

Intermediate 10 min read

What You'll Learn

  • Define admissible fitted means for positive-response GLMs
  • Check an identity-link predictor over an entire covariate range
  • Use covariate-range vertices to verify positivity for affine predictors
  • Identify the zero and sign restrictions of the reciprocal link
  • Distinguish invalid mean predictions from response-family misspecification

Why Link Admissibility Matters

For a positive-response model, the conditional mean must itself be positive. A Gamma or inverse-Gaussian distribution, for example, cannot have a nonpositive mean parameter. The link function connects that mean to the linear predictor, so the link and the range of the predictor determine whether the model can produce allowable fitted means.

This is especially important for the identity and reciprocal links. Neither link gives positive means for every possible real-valued linear predictor. A fit can yield plausible estimates at the observed covariate values yet produce an invalid prediction elsewhere in the range where the model is intended to apply. The right check is therefore not merely whether the fitted values look reasonable in the data table, but whether the mean remains admissible over the relevant covariate range.

Definition: In a GLM, the conditional mean is \(\mu_i=E[Y_i\mid x_i]\), and the linear predictor is \(\eta_i=x_i^\mathsf{T}\beta\), possibly with an offset added. A link \(g\) specifies \(g(\mu_i)=\eta_i\). For a positive-response family, a fitted mean is admissible when it belongs to the family’s mean-parameter domain; in the settings considered here, this requires \(\mu_i>0\). A link fit is admissible over a covariate range if this condition holds at every point in that range.

The family and the link impose different requirements. The family determines which distributions and parameter values are available for the response; the link determines how the mean is related to covariates. A valid positive mean does not establish that the chosen family describes the data well. Conversely, a sensible family cannot make a negative fitted mean allowable.

Identity Link: Check the Predictor Itself

With the identity link, \(g(\mu)=\mu\), so the inverse link is also the identity function and \(\mu_i=\eta_i=x_i^\mathsf{T}\beta\). For a positive-response family, every relevant value of this linear predictor must be strictly greater than zero. A zero prediction is not enough: the mean parameter must be positive, not merely nonnegative.

For one covariate \(x\) in an interval \([a,b]\), an identity-link predictor has the form \(\eta(x)=\beta_0+\beta_1x\). Its minimum occurs at an endpoint. Thus, checking both endpoints suffices to establish positivity throughout the interval. With multiple predictors, checking only the observed rows may fail to check the full range or region of intended use.

Theorem: Let \(P\) be the convex hull of finitely many vertices \(v_1,\ldots,v_m\), and let \(\eta(x)=c+x^\mathsf{T}\beta\) be an affine predictor. Then \(\eta(x)>0\) for every \(x\in P\) if and only if \(\eta(v_j)>0\) for every vertex \(v_j\).

Proof. If \(\eta(x)>0\) for every \(x\in P\), then it is positive at each vertex because each vertex belongs to \(P\). For the converse, take any \(x\in P\). By the definition of a convex hull, there are weights \(\lambda_1,\ldots,\lambda_m\) with \(\lambda_j\geq 0\), \(\sum_{j=1}^m\lambda_j=1\), and \(x=\sum_{j=1}^m\lambda_jv_j\). Affineness gives \(\eta(x)=\sum_{j=1}^m\lambda_j\eta(v_j)\). If each vertex value is strictly positive, this weighted average is strictly positive. Therefore \(\eta(x)>0\) throughout \(P\). \(\square\)

Worked Example: Verifying an Identity Link on a Rectangle

Suppose a Gamma GLM uses \(\mu=\eta\), with predictor \(\eta(x_1,x_2)=1.1+0.15x_1-0.08x_2\). The intended range is the rectangle \(0\leq x_1\leq4\), \(0\leq x_2\leq3\). Check whether all fitted means in that range are positive.

The rectangle’s vertices are \((0,0)\), \((4,0)\), \((0,3)\), and \((4,3)\). Their predictor values are \(\eta(0,0)=1.1\), \(\eta(4,0)=1.1+0.60=1.70\), \(\eta(0,3)=1.1-0.24=0.86\), and \(\eta(4,3)=1.1+0.60-0.24=1.46\).

The minimum vertex value is \(0.86>0\). The predictor is affine and the range is a rectangle, which is the convex hull of these vertices. By the theorem, \(\eta\), and therefore the identity-link mean, remains positive throughout the rectangle. This establishes admissibility on the stated range, not goodness of fit or suitability of the Gamma family.

Worked Example: Finding Where an Identity-Link Prediction Fails

Consider a positive-response model with \(\eta(x)=2-0.25x\) and the identity link. The model is to be used for \(0\leq x\leq10\). Find the range over which its predicted means are admissible.

Because \(\mu(x)=\eta(x)=2-0.25x\), positivity requires \(2-0.25x>0\). Solving gives \(x<8\). At \(x=7\), the predicted mean is \(2-1.75=0.25\), which is positive. At \(x=8\), it is \(0\), which is inadmissible; at \(x=10\), it is \(2-2.5=-0.5\), also inadmissible.

Thus predictions are admissible only for \(0\leq x<8\) within the proposed range. Even if the observed data happened to have \(x\leq7\), the model would not provide valid positive means over the full intended range through \(10\). This is a link-and-predictor admissibility problem, not by itself evidence that the response family is wrong.

Reciprocal Link: Keep the Predictor Away from Zero

Here, “inverse link” means the reciprocal link, \(g(\mu)=1/\mu\), not simply the inverse function \(g^{-1}\) of an arbitrary link. Under this convention, the fitted mean is \(\mu=1/\eta\). For a positive mean, the predictor must be strictly positive. If \(\eta=0\), the expression is undefined; if \(\eta<0\), the resulting mean is negative.

Some texts or software use a signed reciprocal convention, such as \(g(\mu)=-1/\mu\). Under that convention, the predictor must instead be strictly negative to yield a positive mean. The sign condition must be matched to the stated link; the word “inverse” alone is not precise enough to identify it.

Theorem: Under the reciprocal link \(g(\mu)=1/\mu\), a fitted mean is positive and finite at \(x\) if and only if \(\eta(x)>0\). If \(\eta\) is continuous on a compact covariate range \(P\), then the fitted mean is positive and finite throughout \(P\) if and only if \(\eta(x)>0\) throughout \(P\).

Proof. If \(\eta(x)>0\), then its reciprocal exists, is finite, and is positive, so \(\mu(x)=1/\eta(x)>0\). If \(\eta(x)=0\), the reciprocal is undefined. If \(\eta(x)<0\), its reciprocal is negative. These cases prove the first claim. Applying the pointwise equivalence at every \(x\in P\) proves the second. In addition, if \(\eta\) is continuous and strictly positive on compact \(P\), it attains a minimum \(\delta\) there. Since it is positive everywhere, \(\delta>0\), and hence \(0<\mu(x)\leq1/\delta\) throughout \(P\). \(\square\)

Worked Example: Checking a Reciprocal-Link Model

A model uses \(\mu(x)=1/\eta(x)\), with \(\eta(x)=0.4+0.08x\), for \(0\leq x\leq10\). Determine whether its means are admissible across this range and calculate the endpoint predictions.

The predictor increases with \(x\), so its minimum occurs at \(x=0\), where \(\eta(0)=0.4>0\). Its maximum is \(\eta(10)=0.4+0.8=1.2\). Thus it is positive throughout the interval. The endpoint means are \(\mu(0)=1/0.4=2.5\) and \(\mu(10)=1/1.2=5/6\), approximately \(0.833\). Both are positive and finite, and the theorem ensures the same throughout the range.

The reciprocal transformation also reverses the direction of the predictor’s effect. Differentiating gives \(d\mu/dx=-0.08/\eta(x)^2<0\), so the fitted mean decreases as \(x\) increases, even though the linear predictor increases.

A reciprocal predictor can fail even when its values appear moderate at some observed covariate values. For example, if \(\eta(x)=0.6-0.1x\) on \(0\leq x\leq8\), then \(\eta(0)=0.6\), but it reaches zero at \(x=6\) and is negative for \(x>6\). The mean is \(1/0.6=5/3\) at \(x=0\), undefined at \(x=6\), and negative at \(x=8\), where it equals \(-5\). The zero is a singularity, not a very large but valid prediction.

What an Invalid Prediction Does—and Does Not—Tell You

A nonpositive mean under either link indicates that the proposed mean structure cannot be used as specified across that covariate range for a positive-response family. The problem may be addressed by changing the link, restricting the prediction domain, or changing the mean structure, but the appropriate choice depends on the scientific purpose and data. A log link is one common way to enforce positivity because its inverse, \(\mu=\exp(\eta)\), is positive for every finite real \(\eta\); positivity alone, however, does not make a model otherwise appropriate.

It is important not to confuse an inadmissible prediction with a family mismatch. For example, negative means from an identity-link Gamma model do not show that Gamma variability or positive skewness is unsuitable. They show that the fitted predictor conflicts with the family’s positive mean domain. The reverse distinction matters too: if all means are positive, that does not establish that the family captures the response distribution or its variance pattern.

1
Specify the prediction domain.
State the covariate values or region where the fitted model is intended to be used; observed rows alone may not represent the full range of interest.
2
Compute the linear predictor.
Include the intercept and any offset. For an identity link, this is the mean itself; for the reciprocal link, it is the denominator of the mean.
3
Check the required sign everywhere.
For identity-link positive means require \(\eta>0\); for the reciprocal convention \(g(\mu)=1/\mu\), also require \(\eta>0\). On a convex polytope with an affine predictor, checking all vertices is sufficient.
4
Diagnose the scope of the failure.
A zero or negative predicted mean violates the parameter domain for a positive-response family. Investigate the link and mean structure before treating it as evidence against the family.

Common Mistakes

  • Checking only fitted values at observed rows. A model intended for a wider covariate range may predict zero or negative means between or beyond those rows. Define the range of use and check it explicitly.
  • Accepting a zero mean for a positive family. The condition is strict: \(\mu>0\). A zero identity-link prediction is not admissible for the positive mean domain considered here.
  • Assuming a reciprocal link is positive for any predictor. With \(g(\mu)=1/\mu\), \(\eta=0\) is undefined and \(\eta<0\) gives a negative mean.
  • Calling an inadmissible link prediction a family failure. It is first a conflict between the link-generated mean and the family’s parameter domain. A separate assessment is needed to judge the family’s fit.
Key takeaway: For positive-response models, admissibility is a domain check on predicted means across the covariate range of interest. The identity link requires a positive affine predictor; the reciprocal link \(g(\mu)=1/\mu\) requires a strictly positive denominator, with no zero crossing.

Check Your Understanding

Use the stated link convention and check the entire range, not just one selected prediction.

  1. For the identity link with \(\eta(x)=1.5-0.2x\) on \(0\leq x\leq6\), determine whether the predicted mean is admissible everywhere.
  2. An identity-link model has an affine predictor over a rectangular covariate region. Why can checking all vertices establish positivity throughout the region?
  3. Under \(g(\mu)=1/\mu\), what happens to the fitted mean when the predictor is zero? What if it is negative?
  4. A Gamma GLM yields a negative identity-link prediction at one covariate value. What conclusion follows immediately, and what conclusion about family fit does not follow?
  5. For a reciprocal-link model with \(\eta(x)=0.3+0.05x\) on \(0\leq x\leq4\), calculate the fitted means at both endpoints and state whether the predictions are admissible across the interval.
Clinical Biostats