Tutorials › Generalized Linear Models › Normal Responses And Constant Variance Functions

Mean Variance and Dispersion · Tutorial 51 of 981

Normal Responses And Constant Variance Functions

Write a normal response in exponential-family form and distinguish its constant variance function from the dispersion-scaled response variance.

Beginner 9 min read

What You'll Learn

  • Express a normal density in exponential-dispersion form.
  • Identify its mean parameter, dispersion parameter, and log-partition function.
  • Derive the normal response mean and variance from its moment-generating function.
  • Calculate the variance function and explain why it equals one.
  • Interpret the dispersion parameter as the response variance.
  • Recognize common errors involving variance, standard deviation, and dispersion.

Why the Normal Variance Function Is Constant

For a normal response, the variance does not change as the mean changes, provided the dispersion parameter is held fixed. The exponential-dispersion representation makes that fact especially transparent: the variance function is \(V(\mu)=1\), while the response variance is \(\phi V(\mu)=\phi\). The distinction matters because the variance function describes how variance depends on the mean, whereas dispersion supplies the overall scale.

We derive both quantities from the normal density. In the course’s convention, \(N(\mu,\sigma^2)\) uses its second argument for the variance. We write the normal dispersion parameter as \(\phi>0\), so that a response with mean \(\mu\) and dispersion \(\phi\) has distribution \(N(\mu,\phi)\).

The Normal Exponential-Dispersion Form

The density of \(Y\sim N(\mu,\phi)\), for \(y\in\mathbb R\), is

$$ f(y;\mu,\phi)=\frac{1}{\sqrt{2\pi\phi}} \exp\!\left(-\frac{(y-\mu)^2}{2\phi}\right), \qquad \mu\in\mathbb R,\quad \phi>0. $$

Expanding the square in the exponent separates the terms involving the response \(y\) and the parameter \(\mu\):

$$ -\frac{(y-\mu)^2}{2\phi} = \frac{y\mu-\mu^2/2}{\phi}-\frac{y^2}{2\phi}. $$

Thus the density can be written in the exponential-dispersion form \(f(y;\theta,\phi)=\exp\{[y\theta-b(\theta)]/\phi+c(y,\phi)\}\). For the normal response, \(\theta=\mu\), \(b(\theta)=\theta^2/2\), and \(c(y,\phi)=-y^2/(2\phi)-\tfrac12\log(2\pi\phi)\).

Definition: For a regular exponential-dispersion family written as \(f(y;\theta,\phi)=\exp\{[y\theta-b(\theta)]/\phi+c(y,\phi)\}\), the mean parameter is \(\mu=b'(\theta)\). When this relationship can be inverted, the variance function is \(V(\mu)=b''(\theta(\mu))\). The response variance is \(\phi V(\mu)\). For a normal response, \(\theta=\mu\), \(b(\theta)=\theta^2/2\), and \(\phi>0\).

The support is the whole real line, independent of \(\mu\) and \(\phi\), and \(b(\theta)\) is defined for every real \(\theta\). Differentiating the log-partition function gives \(b'(\theta)=\theta\), so the mean parameter is indeed \(\mu=\theta\). The same simple quadratic form will let us obtain the variance directly.

Worked Example: Identifying the Normal Family Components

Suppose \(Y\sim N(5,9)\). Identify the mean and dispersion and write the density in exponential-dispersion form.

The second argument of \(N(5,9)\) is the variance, so \(\mu=5\) and \(\phi=9\). Substituting into the normal density gives

$$ f(y;5,9)=\frac{1}{\sqrt{18\pi}} \exp\!\left(-\frac{(y-5)^2}{18}\right). $$

Here \(\theta=5\) and \(b(5)=5^2/2=12.5\). The exponent in the exponential-dispersion representation is

$$ \frac{5y-12.5}{9}-\frac{y^2}{18}-\frac12\log(18\pi). $$

Combining the first two terms gives \(- (y-5)^2/18\), since \((5y-12.5)/9-y^2/18=-(y^2-10y+25)/18\). This agrees with the original density. In particular, the dispersion is \(9\), not the standard deviation \(3\).

Deriving the Mean and Variance

A direct route from the exponential-dispersion form to moments is the moment-generating function. For the normal representation, multiplying the density by \(e^{ty}\) adds \(ty\) to the exponent. This is equivalent to replacing \(\theta\) with \(\theta+\phi t\), so long as the resulting parameter is in the family’s domain. Here that domain is all of \(\mathbb R\), and therefore the calculation is valid for every real \(t\).

Result: If \(Y\) has the normal exponential-dispersion density above, then its moment-generating function is \(M_Y(t)=\exp(\mu t+\phi t^2/2)\) for every real \(t\). Consequently, \(E[Y]=\mu\) and \(\operatorname{Var}(Y)=\phi\).

Derivation. Using the density’s exponential-dispersion form and the fact that it integrates to one at every real value of \(\theta\),

$$ \begin{aligned} M_Y(t) &=E[e^{tY}]\\ &=\int_{\mathbb R}\exp\!\left\{ \frac{y\theta-b(\theta)}{\phi}+c(y,\phi)+ty \right\}\,dy\\ &=\exp\!\left\{\frac{b(\theta+\phi t)-b(\theta)}{\phi}\right\} \int_{\mathbb R}\exp\!\left\{ \frac{y(\theta+\phi t)-b(\theta+\phi t)}{\phi} +c(y,\phi) \right\}\,dy\\ &=\exp\!\left\{\frac{b(\theta+\phi t)-b(\theta)}{\phi}\right\}. \end{aligned} $$

The integral in the third line is one: it is the density normalized at parameter value \(\theta+\phi t\). For the normal response, \(b(\theta)=\theta^2/2\) and \(\theta=\mu\). Expanding the difference in the last line gives

$$ \frac{b(\mu+\phi t)-b(\mu)}{\phi} = \frac{(\mu+\phi t)^2-\mu^2}{2\phi} =\mu t+\frac{\phi t^2}{2}. $$

Hence \(M_Y(t)=\exp(\mu t+\phi t^2/2)\). This moment-generating function is finite in a neighborhood of zero, so it can be differentiated there to obtain the first two moments. Writing \(K_Y(t)=\log M_Y(t)\), we have \(K_Y'(0)=E[Y]\) and \(K_Y''(0)=\operatorname{Var}(Y)\). In this case,

$$ K_Y(t)=\mu t+\frac{\phi t^2}{2}, \qquad K_Y'(0)=\mu, \qquad K_Y''(0)=\phi. $$

This establishes the mean and variance directly from the exponential-family form. The quadratic term in \(b\) is what makes the variance independent of the mean.

Deriving the Constant Variance Function

The variance function records the part of the response variance that varies with the mean parameter; it does not include the dispersion multiplier. For the normal family, the mean-parameter relationship and its inverse are especially simple.

Result: For the normal exponential-dispersion family, \(b'(\theta)=\theta\), so \(\theta(\mu)=\mu\). Therefore \(V(\mu)=b''(\theta(\mu))=1\) for every \(\mu\in\mathbb R\), and the response variance is \(\operatorname{Var}(Y)=\phi V(\mu)=\phi\).

Derivation. From \(b(\theta)=\theta^2/2\), differentiation gives \(b'(\theta)=\theta\) and \(b''(\theta)=1\). The mean parameter is \(\mu=b'(\theta)=\theta\), and the inverse is \(\theta(\mu)=\mu\), which exists for every real mean. Substituting into the definition of the variance function gives

$$ V(\mu)=b''(\theta(\mu))=b''(\mu)=1. $$

The exponential-dispersion mean-variance law then gives \(\operatorname{Var}(Y)=\phi V(\mu)=\phi\cdot1=\phi\). Thus the variance function is constant across all possible means, while the response variance equals the dispersion parameter. This derivation also agrees with the moment-generating-function calculation above. \(\square\)

Worked Example: Calculating the Variance at Different Means

Consider two normal responses with the same dispersion, \(Y_1\sim N(2,4)\) and \(Y_2\sim N(10,4)\). Calculate their variance functions and response variances.

For \(Y_1\), the mean is \(\mu_1=2\) and \(\phi=4\). Since \(V(\mu)=1\), \(\operatorname{Var}(Y_1)=\phi V(\mu_1)=4\cdot1=4\).

For \(Y_2\), the mean is \(\mu_2=10\) and \(\phi=4\). Again, \(\operatorname{Var}(Y_2)=\phi V(\mu_2)=4\cdot1=4\).

The means differ by \(10-2=8\), but both responses have variance \(4\) because their dispersions match. Their standard deviations are both \(\sqrt{4}=2\). The value of the mean changes the location of each distribution, not its variance function.

Dispersion Sets the Variance Scale

Because \(V(\mu)=1\) regardless of \(\mu\), changing the mean alone does not change the variance in this family. Changing \(\phi\), on the other hand, changes the response variance directly. For instance, \(N(3,1)\) has variance \(1\), whereas \(N(3,16)\) has variance \(16\), despite having the same mean. Their standard deviations are \(1\) and \(4\), respectively.

The parameter \(\phi\) is a variance scale, not a standard-deviation scale: the standard deviation is \(\sqrt{\phi}\). It is also useful to keep the two layers of the variance law visible:

$$ \underbrace{\operatorname{Var}(Y)}_{\text{response variance}} = \underbrace{\phi}_{\text{dispersion}} \underbrace{V(\mu)}_{\text{variance function}} = \phi\cdot1. $$

Worked Example: Changing Measurement Units

Suppose \(Y\sim N(12,9)\) measures a quantity in one unit, and define \(Z=10Y\), so that \(Z\) is measured in units ten times as small. Find the mean and variance of \(Z\), then identify its variance function and dispersion.

For a constant \(c\), \(E[cY]=cE[Y]\) and \(\operatorname{Var}(cY)=c^2\operatorname{Var}(Y)\). With \(c=10\), this gives \(E[Z]=10(12)=120\) and \(\operatorname{Var}(Z)=10^2(9)=900\).

A linear rescaling of a normal variable is still normal, so \(Z\sim N(120,900)\). In the normal exponential-dispersion representation its mean is \(120\), its dispersion is \(900\), and its variance function remains \(V(120)=1\). The response variance is therefore \(900\cdot1=900\), as calculated directly. The variance function remains one; the dispersion changes with the scale of the measurement.

Common Mistakes and Why the Distinction Matters

  • Calling the variance function \(\phi\). The normal variance function is \(V(\mu)=1\); the response variance is \(\phi V(\mu)=\phi\). These quantities have different roles.
  • Using the standard deviation as the dispersion. If a normal response has standard deviation \(3\), its variance and dispersion are \(9\), not \(3\).
  • Assuming that a constant variance function means every normal model has the same variance. The variance is constant with respect to \(\mu\) when \(\phi\) is held fixed. Different dispersions give different response variances.
  • Thinking \(V(\mu)=1\) means \(\operatorname{Var}(Y)=1\). That conclusion omits the dispersion multiplier. The variance equals one only when \(\phi=1\).

This distinction is useful whenever a model’s mean and variance are described separately. The mean parameter specifies the response’s location; the variance function describes how the mean relates to variance within the family; and dispersion supplies the scale multiplying that function. For a normal response the relationship is particularly simple: the mean can vary freely, the variance function stays at one, and the dispersion is the response variance.

Check Your Understanding

Use the exponential-dispersion representation and the distinction between variance function and dispersion to answer these questions.

  1. Write the normal density \(N(\mu,\phi)\) in exponential-dispersion form and identify \(b(\theta)\), \(\theta\), and \(c(y,\phi)\).
  2. For \(Y\sim N(-2,25)\), state the mean, variance, standard deviation, and variance function value.
  3. Explain why the variance function remains constant when the mean changes but the dispersion is fixed.
  4. If a normal response has \(V(\mu)=1\) and variance \(7\), what is its dispersion? What is its standard deviation?
  5. Two normal responses have means \(1\) and \(8\), and both have dispersion \(3\). What are their variances, and why?
Clinical Biostats