Normal-inverse Gaussian distribution

Normal-inverse Gaussian (NIG)

parameters:	μ location (real) α tail heavyness (real) β asymmetry parameter (real) δ scale parameter (real) $\gamma = \sqrt{\alpha^2 - \beta^2}$
support:	$x \in (-\infty; +\infty)\!$
pdf:	$\frac{\alpha\delta K_1 \left(\alpha\sqrt{\delta^2 + (x - \mu)^2}\right)}{\pi \sqrt{\delta^2 + (x - \mu)^2}} \; e^{\delta \gamma + \beta (x - \mu)}$ Kj denotes a modified Bessel function of the second kind
mean:	μ + δβ / γ
variance:	δα2 / γ3
skewness:	$3 \beta /(\alpha \sqrt{\delta \gamma})$
ex.kurtosis:	3(1 + 4β2 / α2) / (δγ)
mgf:	$e^{\mu z + \delta (\gamma - \sqrt{\alpha^2 -(\beta +z)^2})}$
cf:	$e^{i\mu z + \delta (\gamma - \sqrt{\alpha^2 -(\beta +iz)^2})}$

Nig: The normal-inverse Gaussian distribution (NIG) is continuous probability distribution that is defined as the normal variance-mean mixture where the mixing density is the inverse Gaussian distribution. The tails of the distribution decrease more slowly than the normal distribution. It is therefore suitable to model phenomena where numerically large values are more probable than is the case for the normal distribution. Examples are returns from financial assets and turbulent wind speeds. The normal-inverse Gaussian distributions form a subclass of the generalised hyperbolic distributions.

The fact that there is a simple expression for the moment generating function implies that simple expressions for all moments are available. The class of normal-inverse Gaussian distributions is closed under convolution in the following sense. If X₁ and X₂ are independent random variables that are NIG-distributed with the same values of the parameters α and β, but possibly different values of the location and scale parameters, μ₁, δ₁ and μ₂, δ₂, respectively, then X₁ + X₂ is NIG-distributed with parameters α, β,μ₁ + μ₂ and δ₁ + δ₂.

The normal-inverse Gaussian distribution can also be seen as the marginal distribution of the normal-inverse Gaussian process which provides an alternative way of explicitly constructing it. Starting with a drifting Brownian motion (Wiener process), W^(γ)(t) = W(t) + γt, we can define the inverse Gaussian process A_t = inf{s > 0:W^(γ)(s) = δt}. Then given a second independent drifting Brownian motion, $W^{(\beta)}(t)=\tilde W(t)+\beta t$, the normal-inverse Gaussian process is the time-changed process X_t = W^(β)(A_t). The process X(t) at time 1 has the normal-inverse Gaussian distribution described above. The NIG process is a particular instance of the more general class of Lévy processes.

The parameters of the normal-inverse Gaussian distribution are often used to construct a heaviness and skewness plot called the NIG-triangle.