Natural exponential family

In probability and statistics, a natural exponential family (NEF) is a class of probability distributions that is a special case of an exponential family (EF).

Definition

Univariate case

The natural exponential families (NEF) are a subset of the exponential families. A NEF is an exponential family in which the natural parameter η and the natural statistic T(x) are both the identity. A distribution in an exponential family with parameter θ can be written with probability density function (PDF) $$f_X(x\mid \theta )=h(x)\exp (\eta (\theta )T(x)-A(\theta )),$$ where $h(x)$ and $A(\theta )$ are known functions. A distribution in a natural exponential family with parameter θ can thus be written with PDF $$f_X(x\mid \theta )=h(x)\exp (\theta x-A(\theta )).$$ [Note that slightly different notation is used by the originator of the NEF, Carl Morris. Morris uses ω instead of η and ψ instead of A.]

General multivariate case

Suppose that $𝐱\in 𝒳\subseteq {ℝ}^p$, then a natural exponential family of order p has density or mass function of the form: $$f_X(𝐱\mid \mathit{\theta })=h(𝐱)\exp ({\mathit{\theta }}^{\mathrm{T}}𝐱-A(\mathit{\theta })),$$ where in this case the parameter $\mathit{\theta }\in {ℝ}^p.$

Moment and cumulant generating functions

A member of a natural exponential family has moment generating function (MGF) of the form $$M_X(𝐭)=\exp (A(\mathit{\theta }+𝐭)-A(\mathit{\theta })).$$

The cumulant generating function is by definition the logarithm of the MGF, so it is $$K_X(𝐭)=A(\mathit{\theta }+𝐭)-A(\mathit{\theta }).$$

Kullback-Leibler divergence

The Kullback–Leibler divergence of two natural exponential families with parameters $\theta$ and $\lambda$ is

$D_{KL}(\theta ||\lambda )=(\theta -\lambda )A^{\prime }(\theta )-(A(\theta )-A(\lambda ))=\underset{\theta }{\overset{\lambda }{\int }}(\lambda -t)A^{″}(t)dt$

Examples

The five most important univariate cases are:

• normal distribution with known variance
• Poisson distribution
• gamma distribution with known shape parameter α (or k depending on notation set used)
• binomial distribution with known number of trials, n
• negative binomial distribution with known $r$

These five examples – Poisson, binomial, negative binomial, normal, and gamma – are a special subset of NEF, called NEF with quadratic variance function (NEF-QVF) because the variance can be written as a quadratic function of the mean. NEF-QVF are discussed below.

Distributions such as the exponential, Bernoulli, and geometric distributions are special cases of the above five distributions. For example, the Bernoulli distribution is a binomial distribution with n = 1 trial, the exponential distribution is a gamma distribution with shape parameter α = 1 (or k = 1 ), and the geometric distribution is a special case of the negative binomial distribution.

Some exponential family distributions are not NEF. The lognormal and Beta distribution are in the exponential family, but not the natural exponential family. The gamma distribution with two parameters is an exponential family but not a NEF and the chi-squared distribution is a special case of the gamma distribution with fixed scale parameter, and thus is also an exponential family but not a NEF (note that only a gamma distribution with fixed shape parameter is a NEF).

The inverse Gaussian distribution is a NEF with a cubic variance function.

The parameterization of most of the above distributions has been written differently from the parameterization commonly used in textbooks and the above linked pages. For example, the above parameterization differs from the parameterization in the linked article in the Poisson case. The two parameterizations are related by $\theta =\log (\lambda )$, where λ is the mean parameter, and so that the density may be written as $$f(k;\theta )=\frac{1}{k!}\exp (\theta k-\exp (\theta )),$$ for $\theta \in ℝ$, so $$h(k)=\frac{1}{k!},\text{ and }A(\theta )=\exp (\theta ).$$

This alternative parameterization can greatly simplify calculations in mathematical statistics. For example, in Bayesian inference, a posterior probability distribution is calculated as the product of two distributions. Normally this calculation requires writing out the probability distribution functions (PDF) and integrating; with the above parameterization, however, that calculation can be avoided. Instead, relationships between distributions can be abstracted due to the properties of the NEF described below.

An example of the multivariate case is the multinomial distribution with known number of trials.

Properties

The properties of the natural exponential family can be used to simplify calculations involving these distributions.

Multivariate case

In the multivariate case, the mean vector and covariance matrix areTemplate:Citation needed $$\mathrm{E}[X]=\mathrm{\nabla }A(\mathit{\theta })\text{ and }\mathrm{Cov}[X]=\mathrm{\nabla }{\mathrm{\nabla }}^{\mathrm{T}}A(\mathit{\theta }),$$ where$\mathrm{\nabla }$ is the gradient and $\mathrm{\nabla }{\mathrm{\nabla }}^{\mathrm{T}}$ is the Hessian matrix.

Natural exponential families with quadratic variance functions (NEF-QVF)

A special case of the natural exponential families are those with quadratic variance functions. Six NEFs have quadratic variance functions (QVF) in which the variance of the distribution can be written as a quadratic function of the mean. These are called NEF-QVF. The properties of these distributions were first described by Carl Morris.

$$\mathrm{Var}(X)=V(\mu )={\nu }_0+{\nu }_1\mu +{\nu }_2{\mu }^2.$$

The six NEF-QVFs

The six NEF-QVF are written here in increasing complexity of the relationship between variance and mean.

1. The normal distribution with fixed variance $X\sim N(\mu ,{\sigma }^2)$ is NEF-QVF because the variance is constant. The variance can be written $\mathrm{Var}(X)=V(\mu )={\sigma }^2$, so variance is a degree 0 function of the mean.
2. The Poisson distribution $X\sim \mathrm{Poisson}(\mu )$ is NEF-QVF because all Poisson distributions have variance equal to the mean $\mathrm{Var}(X)=V(\mu )=\mu$, so variance is a linear function of the mean.
3. The Gamma distribution $X\sim \mathrm{Gamma}(r,\lambda )$ is NEF-QVF because the mean of the Gamma distribution is $\mu =r\lambda$ and the variance of the Gamma distribution is $\mathrm{Var}(X)=V(\mu )={\mu }^2/r$, so the variance is a quadratic function of the mean.
4. The binomial distribution $X\sim \mathrm{Binomial}(n,p)$ is NEF-QVF because the mean is $\mu =np$ and the variance is $\mathrm{Var}(X)=np(1-p)$ which can be written in terms of the mean as

   $V(X)=-n{p}^2+np=-{\mu }^2/n+\mu .$

5. The negative binomial distribution $X\sim \mathrm{NegBin}(n,p)$ is NEF-QVF because the mean is $\mu =np/(1-p)$ and the variance is $V(\mu )={\mu }^2/n+\mu .$
6. The (not very famous) distribution generated by the generalizedTemplate:Clarify hyperbolic secant distribution (NEF-GHS) hasTemplate:Citation needed $V(\mu )={\mu }^2/n+n$ and $\mu >0.$

Properties of NEF-QVF

The properties of NEF-QVF can simplify calculations that use these distributions. Template:Ordered list Template:One source Template:Refimprove

See also

• Generalized linear model
• Pearson distribution
• Sheffer sequence
• Orthogonal polynomials

References

• Morris C. (1982) Natural exponential families with quadratic variance functions: statistical theory. Dept of mathematics, Institute of Statistics, University of Texas, Austin.