core28 min studyLesson 243 in Mathematics

Exponential families and completeness

Connect natural parameters, log normalizers, sufficient statistics, and a uniqueness property of unbiased estimation.

On this page

What you’ll understand

  • Express a Bernoulli family in natural exponential form.
  • Derive means and variances from a log normalizer under differentiation conditions.
  • Prove completeness for a binomial count using a polynomial identity.

Many familiar probability laws share an algebraic structure: the parameter couples linearly to a statistic inside an exponential, and a normalizer makes total probability one. This structure explains why sums of observations repeatedly appear in inference. A second property, completeness, answers a different question: can a nonzero function of a statistic have expectation zero for every parameter value? Combining sufficiency with completeness will later make an unbiased estimator unique. Neither property should be inferred merely because a density contains an exponential symbol.

Natural parameters expose the structure

A one-parameter natural exponential family has density h(x)exp⁡{ηT(x)−A(η)}h(x)\exp\{\eta T(x)-A(\eta)\} with respect to a fixed measure. Its natural parameter domain contains the values where the normalizing integral is finite. The data-only factor includes a fixed support restriction; a parameter-dependent support needs special care and may prevent this representation. For independent observations, the exponents add, producing the sufficient statistic equal to the sum of the individual statistics. This is a factorization argument, not a claim about every possible model.

For a Bernoulli observation with interior probability, choose the log odds η=log⁡[p/(1−p)]\eta=\log[p/(1-p)]. Solving gives p=eη/(1+eη)p=e^\eta/(1+e^\eta), and the log normalizer is A(η)=log⁡(1+eη)A(\eta)=\log(1+e^\eta). Every real natural parameter is allowed. The endpoint Bernoulli laws are limits of this representation rather than finite values of η\eta. Working with the open family is enough for the completeness argument below, and any almost-sure conclusion for its finite support also extends to the endpoint laws.

fη(x)=exp⁡{ηx−A(η)},x∈{0,1},A(η)=log⁡(1+eη)f_\eta(x)=\exp\{\eta x-A(\eta)\},\quad x\in\{0,1\},\qquad A(\eta)=\log(1+e^\eta)
The Bernoulli model is a natural exponential family with statistic x, counting measure, and a log normalizer that enforces total mass one.↔ Scroll wide equations horizontally.

The normalizer records moments

Differentiate the identity integrating the density to one. When differentiation can pass through the integral, the first derivative of the log normalizer equals the mean of the natural statistic. A second differentiation gives its variance. For a finite sample space these interchanges are automatic finite-sum manipulations. In continuous families, an open neighborhood of finite exponential moments supplies the needed domination. Do not take such interchanges for granted at the edge of a natural parameter domain.

In the Bernoulli family these derivatives yield pp and p(1−p)p(1-p). For nn independent observations, the log normalizer becomes nAnA and the count variance is np(1−p)np(1-p). Convexity of the log normalizer follows from the nonnegative variance. A zero variance can signal a redundant statistic direction; strict convexity needs a nondegeneracy condition. These identities link probability structure to the curvature used in estimation, but curvature is expressed in the chosen natural coordinate, not invariant as an unqualified second derivative.

A′(η)=Eη[T(X)],A′′(η)=Var⁡η(T(X));A′(η)=p,A′′(η)=p(1−p)A'(\eta)=E_\eta[T(X)],\qquad A''(\eta)=\operatorname{Var}_\eta(T(X));\qquad A'(\eta)=p,\quad A''(\eta)=p(1-p)
The first two identities require justified differentiation. The last two evaluate them for the single Bernoulli observation.↔ Scroll wide equations horizontally.

Completeness is an expectation uniqueness theorem

A statistic is complete for a family if every measurable function gg that is integrable under every model law and has expectation zero for every parameter is zero almost surely under every model law. Sufficiency concerns the conditional data distribution; completeness concerns functions of the statistic and expectations across the parameter space. They are different statements. A standard theorem gives completeness of the natural statistic in a full finite-dimensional exponential family when its natural parameter space contains a nonempty open subset of the full Euclidean parameter space, with integrability as stated. Curved subfamilies need not inherit this conclusion.

For a fixed binomial sample size, completeness has an elementary proof. Set the expectation of g(K)g(K) to zero, divide by (1−p)n(1-p)^n, and write the odds as a positive variable zz. The result is a polynomial in zz that vanishes for every positive value. All coefficients must be zero, and the positive binomial coefficients force every value g(k)g(k) to be zero. Restricting the parameter to a few isolated probabilities would not justify this polynomial argument; the extent of the parameter space matters.

0=Ep[g(K)](1−p)n=∑k=0ng(k)(nk)zk,z=p1−p>00=\frac{E_p[g(K)]}{(1-p)^n}=\sum_{k=0}^n g(k)\binom nk z^k,\qquad z=\frac p{1-p}>0
A polynomial vanishing on the positive real axis is identically zero, establishing completeness of a binomial count for p in the open unit interval.↔ Scroll wide equations horizontally.

Check completeness for two trials

Let K∼Binomial⁡(2,p)K\sim\operatorname{Binomial}(2,p) for all 0<p<10<p<1. Suppose Ep[g(K)]=0E_p[g(K)]=0 for every such probability. Determine g(0),g(1),g(2)g(0),g(1),g(2).

  1. Expand the expectation as g(0)(1−p)2+2g(1)p(1−p)+g(2)p2=0g(0)(1-p)^2+2g(1)p(1-p)+g(2)p^2=0. This is an identity over the whole open probability interval, not an equation at one chosen probability.

  2. Divide by (1−p)2(1-p)^2 and substitute z=p/(1−p)z=p/(1-p). The identity becomes g(0)+2g(1)z+g(2)z2=0g(0)+2g(1)z+g(2)z^2=0 for every positive zz.

  3. A nonzero polynomial of degree at most two cannot have infinitely many roots. Therefore each coefficient is zero, giving g(0) = g(1) = g(2) = 0.

The result

The count is complete in this family; there is no nontrivial zero-expectation function of it across all interior probabilities.

Explore the model

Interactive lab

Interactive lab · probability

Expectations change across a family

Reset, set both p=0.25,n=2p=0.25,n=2, and use the exact bars to calculate the expectation of g(K)=K−1g(K)=K-1: predict −0.5-0.5. Reset, set both p=0.75,n=2p=0.75,n=2, and predict 0.50.5. The function is centered at p=0.5p=0.5 only, so it does not contradict completeness. Playback is outcome highlighting, not a proof over every probability.

Expectations change across a familyThese are exact binomial probabilities for independent trials with one fixed success probability. Playback tours the distribution; it does not resample outcomes. Current readings: Mean · np: 5; Variance · np(1−p): 2.5; Total probability: 1. All parameters are adjustable with the labeled controls below.03581000.070.150.220.3Successes kP(X = k)k = 0 · P = 9.77e-4
Exact probabilityHighlighted outcomeMean np
Mean · np
5
Variance · np(1−p)
2.5
Total probability
1
t = 0 s

These are exact binomial probabilities for independent trials with one fixed success probability. Playback tours the distribution; it does not resample outcomes. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forZero expectation at one parameter value is far weaker than zero expectation throughout a family.

Check your understanding

What is required in the definition of a complete statistic?

Try it yourself

For X∼Poisson⁡(λ)X\sim\operatorname{Poisson}(\lambda), λ>0\lambda>0, write the natural parameter and log normalizer. Use their first two derivatives to recover the mean and variance of XX.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Start from e−λλx/x!e^{-\lambda}\lambda^x/x! and put the parameter multiplying xx into one symbol.

Compare with the solution

Choose η=log⁡λ\eta=\log\lambda, so A(η)=eηA(\eta)=e^\eta and h(x)=1/x!h(x)=1/x! on nonnegative integers. Both A′(η)A'(\eta) and A′′(η)A''(\eta) equal eη=λe^\eta=\lambda. Thus the mean and variance are both λ\lambda. The natural domain is all real η\eta, where the exponential moments needed for the differentiations exist locally.

Take it with you

  • Natural exponential form separates a statistic from a log normalizer.
  • Normalizer derivatives encode moments under justified differentiation.
  • Completeness requires a zero-expectation identity across the entire model.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.