Tutorials › Generalized Linear Models › The Poisson Mean-Variance Relationship

Mean Variance and Dispersion · Tutorial 50 of 981

The Poisson Mean-Variance Relationship

See how the Poisson log-partition function yields equal mean and variance, and what “unit dispersion” means for this response family.

Beginner 9 min read

What You'll Learn

  • Write the Poisson distribution in canonical exponential-family form.
  • Identify its natural parameter and log-partition function.
  • Derive its mean and variance by differentiating the log-partition function.
  • Verify the moment calculations directly from Poisson probabilities.
  • Explain why the standard Poisson family has fixed unit dispersion.
  • Distinguish a distribution’s variance from the value of one observed count.

Why the Poisson Mean and Variance Matter

The Poisson distribution is a standard model for nonnegative integer counts. Its defining mean-variance relationship is unusually simple: the variance equals the mean. That equality is not an extra assumption added after choosing the distribution; it follows from the Poisson probability mass function and, in particular, from the derivatives of its log-partition function.

This connection matters in generalized linear models because the mean describes the response's expected count, while the variance describes its variability around that mean. Deriving both from the same function makes clear which part of the relationship is fixed by the Poisson family and why its dispersion is called one.

Canonical Form and Log-Partition Function

For a Poisson random variable \(Y\) with mean parameter \(\mu>0\), the probability mass function is

$$ P(Y=y)=\frac{e^{-\mu}\mu^y}{y!}, \qquad y=0,1,2,\ldots . $$

Set the natural parameter to \(\theta=\log\mu\), so that \(\mu=e^\theta\) and \(\theta\) ranges over all real numbers. Rewriting the mass function gives

$$ P(Y=y)=\exp\left\{y\theta-e^\theta-\log(y!)\right\}. $$

In the canonical exponential-family form \(\exp\{y\theta-b(\theta)+c(y)\}\), the log-partition function is \(b(\theta)=e^\theta\), and the data-only term is \(c(y)=-\log(y!)\). The support, the nonnegative integers, does not depend on \(\theta\).

Definition: The Poisson family in canonical form has natural parameter \(\theta=\log\mu\), log-partition function \(b(\theta)=e^\theta\), and \(c(y)=-\log(y!)\). Its mean parameter is \(\mu=e^\theta\), with \(\mu>0\).

The normalizing role of \(b\) can be checked by summing over the support:

$$ \sum_{y=0}^{\infty}\exp\{y\theta-e^\theta-\log(y!)\} =e^{-e^\theta}\sum_{y=0}^{\infty}\frac{(e^\theta)^y}{y!} =e^{-e^\theta}e^{e^\theta}=1. $$

Deriving the Mean and Variance

For a regular one-parameter exponential family with a fixed support and a differentiable log-partition function, differentiating \(b\) once gives the expectation, and differentiating twice gives the variance. Here the support is fixed and the relevant series can be differentiated term by term: on any bounded interval of \(\theta\), the series of derivatives is dominated by convergent exponential series. We can also verify the result directly from the Poisson mass function.

Theorem: If \(Y\) has the Poisson distribution with parameter \(\mu=e^\theta>0\), then $$ E[Y]=b'(\theta)=e^\theta=\mu,\qquad \operatorname{Var}(Y)=b''(\theta)=e^\theta=\mu. $$

Proof. Let \(Z(\theta)=\sum_{y=0}^{\infty}e^{y\theta}/y!\). For any bounded interval \(|\theta|\le M\), the \(k\)th derivative terms are bounded in absolute value by \(y^k e^{My}/y!\). The series of these bounds converges for each fixed nonnegative integer \(k\), by the ratio test. Thus the series may be differentiated term by term there. Since \(Z(\theta)=e^{e^\theta}\), its log is \(b(\theta)=\log Z(\theta)=e^\theta\).

$$ b'(\theta)=e^\theta,\qquad b''(\theta)=e^\theta. $$

To connect these derivatives to moments, differentiate \(Z\). The probability of \(Y=y\) is \(e^{y\theta}/(y!Z(\theta))\), so

$$ \frac{Z'(\theta)}{Z(\theta)} =\sum_{y=0}^{\infty}y\,P(Y=y) =E[Y]. $$

Differentiating the ratio once more yields

$$ \frac{d}{d\theta}\frac{Z'}{Z} =\frac{Z''}{Z}-\left(\frac{Z'}{Z}\right)^2 =E[Y^2]-E[Y]^2 =\operatorname{Var}(Y). $$

Because \(\log Z=b\), these identities say \(E[Y]=b'(\theta)\) and \(\operatorname{Var}(Y)=b''(\theta)\). Substituting \(b'(\theta)=b''(\theta)=e^\theta=\mu\) proves the result. \(\square\)

The variance conclusion can also be checked through the second factorial moment. For \(\mu>0\),

$$ E[Y(Y-1)] =\sum_{y=2}^{\infty}y(y-1)\frac{e^{-\mu}\mu^y}{y!} =\mu^2e^{-\mu}\sum_{j=0}^{\infty}\frac{\mu^j}{j!} =\mu^2. $$

Since \(Y^2=Y(Y-1)+Y\), it follows that \(E[Y^2]=\mu^2+\mu\). Therefore

$$ \operatorname{Var}(Y)=E[Y^2]-E[Y]^2=(\mu^2+\mu)-\mu^2=\mu. $$

This second derivation checks the same relationship from the probability mass function rather than relying only on the log-partition derivative rule.

Worked Calculations

Worked Example: A Poisson Count With Mean Four

Suppose \(Y\sim\operatorname{Poisson}(4)\). Find its canonical parameter, mean, variance, and standard deviation.

The canonical parameter is \(\theta=\log\mu=\log 4\), and the log-partition function is \(b(\theta)=e^\theta\). Therefore

$$ b'(\log 4)=e^{\log 4}=4,\qquad b''(\log 4)=e^{\log 4}=4. $$

Thus \(E[Y]=4\) and \(\operatorname{Var}(Y)=4\). The standard deviation is the square root of the variance, \(\sqrt{4}=2\). The mean is four counts, while the standard deviation is two counts; the variance and mean are numerically equal but describe different features.

Worked Example: A Small Mean and the Probability of Zero

Let \(Y\sim\operatorname{Poisson}(0.5)\). Determine its mean and variance, and calculate the probability of observing at least one count.

The log-partition function is \(b(\theta)=e^\theta\), where \(\theta=\log(0.5)\). Both derivatives equal \(e^\theta=0.5\), so \(E[Y]=0.5\) and \(\operatorname{Var}(Y)=0.5\). The probability of zero is

$$ P(Y=0)=\frac{e^{-0.5}(0.5)^0}{0!}=e^{-0.5}\approx0.6065. $$

Consequently,

$$ P(Y\geq1)=1-P(Y=0)=1-e^{-0.5}\approx1-0.6065=0.3935. $$

A mean below one does not require every observation to be zero or one. The Poisson model still has support on all nonnegative integers; the mean summarizes the distribution, rather than prescribing a possible observed value.

Worked Example: A Count Over a Longer Exposure

Suppose a count \(Y\) over an exposure of \(3.5\) units follows a Poisson model with rate \(0.8\) events per exposure unit. The mean count is \(\mu=(3.5)(0.8)=2.8\). Find the variance and the probability of observing exactly four events.

The natural parameter is \(\theta=\log(2.8)\). Differentiating \(b(\theta)=e^\theta\) gives \(E[Y]=\operatorname{Var}(Y)=e^\theta=2.8\). The probability of four events is

$$ P(Y=4)=\frac{e^{-2.8}(2.8)^4}{4!} =\frac{e^{-2.8}(61.4656)}{24} \approx0.1557. $$

Here the model's expected count and variance are both \(2.8\), even though the particular count \(4\) has probability about \(0.1557\). The equality of mean and variance is a statement about the distribution over repeated observations under the model, not an equality between the mean and every observed count.

Why the Poisson Family Has Unit Dispersion

In an exponential-dispersion representation, the response variance is written as \(\operatorname{Var}(Y)=\phi V(\mu)\), where \(\phi\) is a dispersion factor and \(V(\mu)\) is the variance function. For the standard Poisson family, the mass function has no additional free scale parameter: its canonical form is \(\exp\{y\theta-b(\theta)+c(y)\}\), with no nontrivial divisor of the canonical term. Thus the standard family fixes \(\phi=1\). Since the variance derived above is \(\mu\), its variance function is \(V(\mu)=\mu\), and the variance law is \(1\cdot\mu=\mu\).

Key result: In the standard Poisson family, the dispersion is fixed at \(\phi=1\), and the variance function is \(V(\mu)=\mu\). Therefore \(\operatorname{Var}(Y)=\phi V(\mu)=\mu\).

“Unit dispersion” does not mean the response variance is one. It means the multiplicative dispersion factor is one; the response variance still changes with the mean. For example, a Poisson response with mean \(4\) has variance \(4\), not variance \(1\). The value one belongs to the dispersion factor, not to every response variance.

This distinction is useful when interpreting a Poisson model: the family links the expected count to its variance. It also explains why checking the variance relative to the mean is meaningful when assessing whether the standard Poisson assumption is plausible. The present derivation establishes the model's implication; it does not assert that every real count dataset must follow it.

Common Mistake

A common error is to reason, “The Poisson dispersion is one, so the variance must be one.” This confuses the dispersion multiplier with the variance itself. The correct calculation is

$$ \operatorname{Var}(Y)=\phi V(\mu)=1\cdot\mu=\mu, $$

so a mean of \(0.5\) gives variance \(0.5\), and a mean of \(4\) gives variance \(4\). A related error is to say that an observed count \(y\) has variance \(\mu\). Variance describes the distribution of the random variable \(Y\), not the fixed value observed in one realization. Once \(Y=y\) has been observed, \(y\) is a number; the model variance remains a property of the distribution that generated it.

Check Your Understanding

Use the canonical parameterization and the Poisson mean-variance result to answer these questions.

  1. For a Poisson response with mean \(7\), what are \(\theta\), \(E[Y]\), and \(\operatorname{Var}(Y)\)?
  2. Write the Poisson log-partition function in terms of the natural parameter and differentiate it once and twice.
  3. In the exponential-dispersion variance law, why does unit dispersion not imply that the Poisson response variance equals one?
  4. If \(Y\sim\operatorname{Poisson}(1.2)\), what is \(E[Y^2]\)? Show which moment identity you use.
  5. Explain why an observed count of four does not contradict a model with mean and variance both equal to \(2.8\).
Clinical Biostats