Exponential dispersion model explained

In probability and statistics, the class of exponential dispersion models (EDM), also called exponential dispersion family (EDF), is a set of probability distributions that represents a generalisation of the natural exponential family.^{[1] [2] [3]} Exponential dispersion models play an important role in statistical theory, in particular in generalized linear models because they have a special structure which enables deductions to be made about appropriate statistical inference.

Definition

Univariate case

There are two versions to formulate an exponential dispersion model.

Additive exponential dispersion model

In the univariate case, a real-valued random variable

X belongs to the additive exponential dispersion model with canonical parameter \theta and index parameter λ , X \sim ED ^{*(\theta,} λ ) , if its probability density function can be written asf _{X(x\mid\theta,} λ ) = h ^{*(λ,x)} \exp\left(\theta x - λ A(\theta)\right) .

Reproductive exponential dispersion model

The distribution of the transformed random variable

Y = (X)/(λ) is called reproductive exponential dispersion model, Y \sim ED (\mu, \sigma ^{2)} , and is given byf _{Y(y\mid\mu,} \sigma ^{2)} = h(\sigma ^{2,y)} \exp\left( (\theta y - A(\theta))/(\sigma ^{2}) \right) ,

with

\sigma ^{2} = (1)/(λ) and \mu = A'(\theta) , implying \theta = (A') ^{-1} (\mu) .The terminology dispersion model stems from interpreting \sigma ^{2} as dispersion parameter. For fixed parameter \sigma ^{2} , the ED (\mu, \sigma ^{2)} is a natural exponential family.

Multivariate case

In the multivariate case, the n-dimensional random variable

X has a probability density function of the following form^[1] f _{X} ( x |\boldsymbol{\theta}, λ ) = h( λ , x ) \exp\left( λ (\boldsymbol\theta ^{\top} x - A(\boldsymbol\theta))\right) ,

where the parameter

\boldsymbol\theta has the same dimension as X .

Properties

Cumulant-generating function

The cumulant-generating function of

Y\sim ED (\mu,\sigma ^{2)} is given byK(t;\mu,\sigma ^{2)} = log \operatorname{E}[e ^{tY} ] = (A(\theta + \sigma ^{2} t) - A(\theta))/(\sigma ^{2}) ,

with

\theta = (A') ^{-1} (\mu)

Mean and variance

Mean and variance of

Y\sim ED (\mu,\sigma ^{2)} are given by\operatorname{E}[Y] = \mu = A'(\theta) , \operatorname{Var}[Y] = \sigma ^{2} A''(\theta) = \sigma ^{2} V(\mu) ,

with unit variance function

V(\mu) = A''((A') ^{-1} (\mu)) .

Reproductive

If

Y _{1,\ldots,} Y _{n} are i.i.d. with Y _{i\sim ED \left(\mu, (\sigma ^{2})/(w _{i})} \right) , i.e. same mean \mu and different weights w _{i} , the weighted mean is again an ED with\sum ^{n} _{i =1} (w _{i} Y _{i})/(w _{\bullet}) \sim ED \left(\mu, (\sigma ^{2})/(w _{\bullet}) \right) ,

with

w _{\bullet} = \sum ^{n} _{i =1} w _{i} . Therefore Y _{i} are called reproductive.

Unit deviance

The probability density function of an

ED (\mu, \sigma ^{2)} can also be expressed in terms of the unit deviance d(y,\mu) asf _{Y(y\mid\mu,} \sigma ^{2)} = \tilde{h}(\sigma ^{2,y)} \exp\left( - (d(y,\mu))/(2 \sigma ^{2}) \right) ,

where the unit deviance takes the special form

d(y,\mu) = y f(\mu) + g(\mu) + h(y) or in terms of the unit variance function as d(y,\mu) = 2 \int ^{y} _{\mu} (y - t)/(V(t)) dt .

Examples

Many very common probability distributions belong to the class of EDMs, among them are: normal distribution, binomial distribution, Poisson distribution, negative binomial distribution, gamma distribution, inverse Gaussian distribution, and Tweedie distribution.

Notes and References

1. Jørgensen, B. (1987). Exponential dispersion models (with discussion). Journal of the Royal Statistical Society, Series B, 49 (2), 127 - 162.
2. Jørgensen, B. (1992). The theory of exponential dispersion models and analysis of deviance. Monografias de matemática, no. 51.
3. Marriott, P. (2005) "Local Mixtures and Exponential DispersionModels" pdf