Stable distribution

Stable

Probability density function Symmetric α-stable distributions with unit scale factor Skewed centered stable distributions with unit scale factor
Cumulative distribution function CDFs for symmetric α-stable distributions CDFs for skewed centered stable distributions
Parameters	α ∈ (0, 2] — stability parameter β ∈ [−1, 1] — skewness parameter (note that skewness is undefined) c ∈ (0, ∞) — scale parameter μ ∈ (−∞, ∞) — location parameter
Support	x ∈ R, or x ∈ [μ, +∞) if α < 1 and β = 1, or x ∈ (-∞, μ] if α < 1 and β = −1
PDF	not analytically expressible, except for some parameter values
CDF	not analytically expressible, except for certain parameter values
Mean	μ when α > 1, otherwise undefined
Median	μ when β = 0, otherwise not analytically expressible
Mode	μ when β = 0, otherwise not analytically expressible
Variance	2c^2 when α = 2, otherwise infinite
Skewness	0 when α = 2, otherwise undefined
Ex. kurtosis	0 when α = 2, otherwise undefined
Entropy	not analytically expressible, except for certain parameter values
MGF	undefined
CF	$\exp\!\Big[\; it\mu - |c\,t|^\alpha\,(1-i \beta\,\mbox{sgn}(t)\Phi) \;\Big],$ where $\Phi = \begin{cases} \tan\tfrac{\pi\alpha}{2} & \text{if }\alpha \ne 1 \\ -\tfrac{2}{\pi}\log|t| & \text{if }\alpha = 1 \end{cases}$

In probability theory, a distribution or a random variable is said to be stable if a linear combination of two independent copies of a random sample has the same distribution, up to location and scale parameters. The stable distribution family is also sometimes referred to as the Lévy alpha-stable distribution, after Paul Lévy, the first mathematician to have studied it.^[1][2]

Of the four parameters defining the family, most attention has been focused on the stability parameter, α (see panel). Stable distributions have 0 < α ≤ 2, with the upper bound corresponding to the normal distribution, and α = 1 to the Cauchy distribution. The distributions have undefined variance for α < 2, and undefined mean for α ≤ 1. The importance of stable probability distributions is that they are "attractors" for properly normed sums of independent and identically-distributed (iid) random variables. The normal distribution defines a family of stable distributions. By the classical central limit theorem the properly normed sum of a set of random variables, each with finite variance, will tend towards a normal distribution as the number of variables increases. Without the finite variance assumption the limit may be a stable distribution. Mandelbrot referred to stable distributions that are non-normal as "stable Paretian distributions",^[3][4][5] after Vilfredo Pareto. Mandelbrot referred to "positive" stable distributions (meaning maximally skewed in the positive direction) with 1<α<2 as "Pareto-Lévy distributions".^[1] He also regarded this range as relevant for the description of stock and commodity prices.

Definition

A non-degenerate distribution is a stable distribution if it satisfies the following property:

Let X₁ and X₂ be independent copies of a random variable X. Then X is said to be stable if for any constants a > 0 and b > 0 the random variable aX₁ + bX₂ has the same distribution as cX + d for some constants c > 0 and d. The distribution is said to be strictly stable if this holds with d = 0.^[6]

Since the normal distribution, the Cauchy distribution, and the Lévy distribution all have the above property, it follows that they are special cases of stable distributions.

Such distributions form a four-parameter family of continuous probability distributions parametrized by location and scale parameters μ and c, respectively, and two shape parameters β and α, roughly corresponding to measures of asymmetry and concentration, respectively (see the figures).

Although the probability density function for a general stable distribution cannot be written analytically, the general characteristic function can be. Any probability distribution is given by the Fourier transform of its characteristic function φ(t) by:

$f(x)=\frac{1}{2\pi}\int_{-\infty}^\infty \varphi(t)e^{-ixt}\,dt$

A random variable X is called stable if its characteristic function can be written as^[6][7]

$\varphi (t;\alpha ,\beta ,c,\mu )=\exp \left[~it\mu \!-\!|ct|^{\alpha }\,(1\!-\!i\beta \,{\textrm {sgn}}(t)\Phi )~\right]$

where sgn(t) is just the sign of t and Φ is given by

$\Phi=\tan(\pi \alpha/2)\,$

for all α except α = 1 in which case:

$\Phi=-\frac{2}{\pi}\log|t|.\,$

μ ∈ R is a shift parameter, β ∈ [−1, 1], called the skewness parameter, is a measure of asymmetry. Notice that in this context the usual skewness is not well defined, as for α < 2 the distribution does not admit 2nd or higher moments, and the usual skewness definition is the 3rd central moment.

The reason this gives a stable distribution is that the characteristic function for the sum of two random variables equals the product of the two corresponding characteristic functions. Adding two random variables from a stable distribution gives something with the same values of α and β, but possibly different values of μ and c.

Not every function is the characteristic function of a legitimate probability distribution (that is, one whose cumulative distribution function is real and goes from 0 to 1 without decreasing), but the characteristic functions given above will be legitimate so long as the parameters are in their ranges. The value of the characteristic function at some value t is the complex conjugate of its value at −t as it should be so that the probability distribution function will be real.

In the simplest case β = 0, the characteristic function is just a stretched exponential function; the distribution is symmetric about μ and is referred to as a (Lévy) symmetric alpha-stable distribution, often abbreviated SαS.

When α < 1 and β = 1, the distribution is supported by [μ, ∞).

The parameter c > 0 is a scale factor which is a measure of the width of the distribution while α is the exponent or index of the distribution and specifies the asymptotic behavior of the distribution.

Parametrizations

The above definition is only one of the parametrizations in use for stable distributions; it is the most common but is not continuous in the parameters at α = 1.

A parametrization that is continuous is^[6]

${\displaystyle \varphi (t;\alpha ,\beta ,\gamma ,\delta )=\exp \left[~it\delta -|\gamma t|^{\alpha }\,(1-i\beta \,{\textrm {sgn}}(t)\Phi )~\right]}$

where Φ is given by

${\displaystyle \Phi =(|\gamma t|^{1-\alpha }-1)\tan(\pi \alpha /2)\,}$

for all α except α = 1 in which case:

${\displaystyle \Phi =-{\frac {2}{\pi }}\log |\gamma t|.}$

The ranges of α and β are the same as before, γ (like c) should be positive, and δ (like μ) should be real.

In either parametrization one can make a linear transformation of the random variable to get a random variable whose density is ${\displaystyle f(y;\alpha ,\beta ,1,0)}$. In the first parametrization, this is done by defining the new variable y as:

${\displaystyle y={\frac {x-\mu }{\gamma }}}$

when α ≠ 1, or

${\displaystyle y={\frac {x-\mu }{\gamma }}-\beta {\frac {2}{\pi }}\ln \gamma }$

when α = 1. For the second parametrization, we simply use

${\displaystyle y={\frac {x-\delta }{\gamma }}.}$

no matter what α is. In the first parametrization, if the mean exists (that is, α > 1) then it is equal to μ, whereas in the second parametrization when the mean exists it is equal to ${\displaystyle \delta -\beta \gamma \tan(\pi \alpha /2)}$.

The distribution

A stable distribution is therefore specified by the above four parameters. It can be shown that any non-degenerate stable distribution has a smooth (infinitely differentiable) density function.^[6] If $f(x;\alpha ,\beta ,c,\mu )$ denotes the density of X and Y is the sum of independent copies of X:

$Y = \sum_{i=1}^N k_i (X_i-\mu)\,$

then Y has the density $s^{{-1}}f(y/s;\alpha ,\beta ,c,0)$ with

$s=\left(\sum_{i=1}^N |k_i|^\alpha\right)^{1/\alpha}.\,$

The asymptotic behavior is described, for α< 2, by:^[6]

$f(x)\sim\frac{c^\alpha (1+\mbox{sgn}(x)\beta) \sin(\pi \alpha / 2)\Gamma(\alpha+1)/\pi}{|x|^{1+\alpha}}$

where Γ is the Gamma function (except that when α < 1 and β = ±1, the tail vanishes to the left or right, resp., of μ). This "heavy tail" behavior causes the variance of stable distributions to be infinite for all α < 2. This property is illustrated in the log-log plots below.

When α = 2, the distribution is Gaussian (see below), with tails asymptotic to exp(−x²/4c²)/(2c√π).

Properties

• All stable distributions are infinitely divisible.
• With the exception of the normal distribution (α = 2), stable distributions are leptokurtotic and heavy-tailed distributions.
• Closure under convolution

Stable distributions are closed under convolution for a fixed value of α. Since convolution is equivalent to multiplication of the Fourier-transformed function, it follows that the product of two stable characteristic functions with the same α will yield another such characteristic function. The product of two stable characteristic functions is given by:

$\exp\left[it\mu_1+it\mu_2 - |c_1 t|^\alpha - |c_2 t|^\alpha +i\beta_1|c_1 t|^\alpha\textrm{sgn}(t)\Phi +i\beta_2|c_2 t|^\alpha\,\textrm{sgn}(t)\Phi \right]$

Since Φ is not a function of the μ, c or β variables it follows that these parameters for the convolved function are given by:

$\mu=\mu_1+\mu_2\,$

$|c|=(|c_1|^\alpha+|c_2|^\alpha)^{1/\alpha}\,$

${\displaystyle \beta ={\frac {\beta _{1}|c_{1}|^{\alpha }+\beta _{2}|c_{2}|^{\alpha }}{|c_{1}|^{\alpha }+|c_{2}|^{\alpha }}}}$

In each case, it can be shown that the resulting parameters lie within the required intervals for a stable distribution.

A generalized central limit theorem

Another important property of stable distributions is the role that they play in a generalized central limit theorem. The central limit theorem states that the sum of a number of independent and identically distributed (i.i.d.) random variables with finite variances will tend to a normal distribution as the number of variables grows.

A generalization due to Gnedenko and Kolmogorov states that the sum of a number of random variables with symmetric distributions having power-law tails (Paretian tails), decreasing as |x|^−α−1 where 0 < α < 2 (and therefore having infinite variance), will tend to a stable distribution $f(x;\alpha ,0,c,0)$ as the number of summands grows.^[8] If α>2 then the sum converges to a stable distribution with stability parameter equal to 2, i.e. a Gaussian distribution.^[9]

There are other possibilities as well. For example, if the characteristic function of the random variable is asymptotic to ${\displaystyle 1-a|t|^{\alpha }\ln(1/|t|)}$ at small t (positive or negative), then we may ask how the value of t, where the value of the characteristic function for the sum of n of these random variables equals some given value u, varies with n.

${\displaystyle \varphi _{\text{sum}}(t)=[\varphi (t)]^{n}=u}$

Assuming for the moment that t tends toward zero, we take the limit of the above as n → ∞:

${\displaystyle {\begin{aligned}\ln u&=\lim _{n\to \infty }[n\ln \varphi (t)]\\&=\lim _{n\to \infty }[-na|t|^{\alpha }\ln(1/|t|)]\end{aligned}}}$

${\displaystyle \ln(-\ln u)=\lim _{n\to \infty }[\ln(na)-\alpha \ln(1/|t|)+\ln \ln(1/|t|)]}$

This shows that ${\displaystyle \ln(1/|t|)}$ is asymptotic to ${\displaystyle (\ln n)/\alpha }$, so using the previous equation we have

${\displaystyle |t|\sim [(-\ln u)\alpha /(na\ln n)]^{1/\alpha }.}$

This implies that the sum divided by ${\displaystyle [na(\ln n)/\alpha ]^{1/\alpha }}$ has a characteristic function whose value at some t' goes to u (with increasing n) when ${\displaystyle t'=(-\ln u)^{1/\alpha }}$. In other words, the characteristic function converges pointwise to ${\displaystyle \exp(-(t')^{\alpha })}$. and therefore by Lévy's continuity theorem the sum divided by ${\displaystyle [na(\ln n)/\alpha ]^{1/\alpha }}$ converges in distribution to the symmetric alpha-stable distribution with stability parameter α and scale parameter 1.

This can be applied to a random variable whose tails decrease as ${\displaystyle |x|^{-3}}$. This random variable has a mean but the variance is infinite. Let us take the following distribution:

${\displaystyle f(x)={\begin{cases}{\frac {1}{3}},&|x|<1,\\{\frac {1}{3}}x^{-3},&|x|>1\end{cases}}}$

We can write this as

${\displaystyle f(x)=\int _{1}^{\infty }{\frac {2}{w^{4}}}h(x/w)dw}$

where:

${\displaystyle h(x/w)={\begin{cases}{\frac {1}{2}},&|x/w|<1,\\0,&|x/w|>1\end{cases}}}$

We want to find the leading terms of the asymptotic expansion of the characteristic function. The characteristic function of the probability distribution ${\displaystyle {\frac {1}{w}}h\left({\frac {x}{w}}\right)}$ is

${\displaystyle {\frac {\sin(tw)}{tw}},}$

so the characteristic function for f(x) is

${\displaystyle \varphi (t)=\int _{1}^{\infty }{\frac {2\sin(tw)}{tw^{4}}}dw}$

and we have

${\displaystyle \varphi (t)-1=\int _{1}^{\infty }{\frac {2}{w^{3}}}\left[{\frac {\sin(tw)}{tw}}-1\right]dw.}$

The term in brackets is asymptotic to ${\displaystyle -{\frac {t^{2}w^{2}}{3!}}+{\frac {t^{4}w^{4}}{5!}}...}$, and we may break the integral into several parts. The first part is the integral of the leading term from w = 1 out to w = 1/|t|, the second term is the integral of the remainder, going from zero to 1/|t|, the third term subtracts out the integral of the same but from 0 to 1, and the last term is the integral from 1/|t| to infinity:

${\displaystyle {\begin{aligned}\varphi (t)-1=&\int _{1}^{1/|t|}{\frac {-t^{2}dw}{3w}}\\&+t^{2}\int _{0}^{1}{\frac {2}{|tw|^{3}}}\left[{\frac {\sin(tw)}{tw}}-1+{\frac {|tw|^{2}}{6}}\right]d|tw|\\&-\int _{0}^{1}{\frac {2}{w^{3}}}\left[{\frac {\sin(tw)}{tw}}-1+{\frac {t^{2}w^{2}}{6}}\right]dw\\&+t^{2}\int _{1}^{\infty }{\frac {2}{y^{3}}}\left({\frac {\sin y}{y}}-1\right)dy\end{aligned}}}$

The first term equals ${\displaystyle {\frac {t^{2}}{3}}\ln |t|}$, the second and fourth terms are constants times t² and the third term is of order t⁴. Therefore

${\displaystyle \varphi (t)\sim 1-{\frac {1}{3}}t^{2}\ln(1/|t|)}$

and according to what was said above (and the fact that the variance of f(x;2,0,1,0) is 2), the sum of n instances of this random variable, divided by ${\displaystyle {\sqrt {n(\ln n)/12}},}$ will converge in distribution to a Gaussian distribution with variance 1. But the variance at any particular n will still be infinite. Note that the width of the limiting distribution grows faster than in the case where the random variable has a finite variance (in which case the width grows as the square root of n). The average, obtained by dividing the sum by n, tends toward a Gaussian whose width approaches zero as n increases, in accordance with the Law of large numbers.

Special cases

Log-log plot of symmetric centered stable distribution PDF's showing the power law behavior for large x. The power law behavior is evidenced by the straight-line appearance of the PDF for large x, with the slope equal to −(α+1). (The only exception is for α = 2, in black, which is a normal distribution.)

Log-log plot of skewed centered stable distribution PDF's showing the power law behavior for large x. Again the slope of the linear portions is equal to −(α+1)

There is no general analytic solution for the form of p(x). There are, however three special cases which can be expressed in terms of elementary functions as can be seen by inspection of the characteristic function:^[6][7][10]

• For α = 2 the distribution reduces to a Gaussian distribution with variance σ² = 2c² and mean μ; the skewness parameter β has no effect.
• For α = 1 and β = 0 the distribution reduces to a Cauchy distribution with scale parameter c and shift parameter μ.
• For α = 1/2 and β = 1 the distribution reduces to a Lévy distribution with scale parameter c and shift parameter μ.

Note that the above three distributions are also connected, in the following way: A standard Cauchy random variable can be viewed as a mixture of Gaussian random variables (all with mean zero), with the variance being drawn from a standard Lévy distribution. And in fact this is a special case of a more general theorem ^[11] which allows any symmetric alpha-stable distribution to be viewed in this way (with the alpha parameter of the mixture distribution equal to twice the alpha parameter of the mixing distribution—and the beta parameter of the mixing distribution always equal to one).

A general closed form expression for stable PDF's with rational values of α is available in terms of Meijer G-functions.^[12] Fox H-Functions can also be used to express the stable probability density functions. For simple rational numbers, the closed form expression is often in terms of less complicated special functions. Several closed form expressions having rather simple expressions in terms of special functions are available. In the table below, PDF's expressible by elementary functions are indicated by an E and those that are expressible by special functions are indicated by an s.^[11]

	α
	1/3	1/2	2/3	1	4/3	3/2	2
β = 0	s	s	s	E	s	s	E
β = 1	s	E	s	s		s

Some of the special cases are known by particular names:

• For α = 1 and β = 1, the distribution is a Landau distribution which has a specific usage in physics under this name.
• For α = 3/2 and β = 0 the distribution reduces to a Holtsmark distribution with scale parameter c and shift parameter μ.

Also, in the limit as c approaches zero or as α approaches zero the distribution will approach a Dirac delta function δ(x − μ).

Series representation

The stable distribution can be restated as the real part of a simpler integral:^[13]

$f(x;\alpha,\beta,c,\mu)=\frac{1}{\pi}\Re\left[ \int_0^\infty e^{it(x-\mu)}e^{-(ct)^\alpha(1-i\beta\Phi)}\,dt\right].$

Expressing the second exponential as a Taylor series, we have:

$f(x;\alpha,\beta,c,\mu)=\frac{1}{\pi}\Re\left[ \int_0^\infty e^{it(x-\mu)}\sum_{n=0}^\infty\frac{(-qt^\alpha)^n}{n!}\,dt\right]$

where $q=c^\alpha(1-i\beta\Phi)$. Reversing the order of integration and summation, and carrying out the integration yields:

$f(x;\alpha,\beta,c,\mu)=\frac{1}{\pi}\Re\left[ \sum_{n=1}^\infty\frac{(-q)^n}{n!}\left(\frac{i}{x-\mu}\right)^{\alpha n+1}\Gamma(\alpha n+1)\right]$

which will be valid for x ≠ μ and will converge for appropriate values of the parameters. (Note that the n = 0 term which yields a delta function in x−μ has therefore been dropped.) Expressing the first exponential as a series will yield another series in positive powers of x−μ which is generally less useful.

Simulation of stable variables

Simulating sequences of stable random variables is not straightforward, since there are no analytic expressions for the inverse ${\displaystyle F^{-1}(x)}$ nor the CDF $F(x)$ itself.^[14][15] All standard approaches like the rejection or the inversion methods would require tedious computations. A much more elegant and efficient solution was proposed by Chambers, Mallows and Stuck (CMS),^[16] who noticed that a certain integral formula^[17] yielded the following algorithm:^[18]

• generate a random variable $U$ uniformly distributed on ${\displaystyle (-{\frac {\pi }{2}},{\frac {\pi }{2}})}$ and an independent exponential random variable $W$ with mean 1;
• for ${\displaystyle \alpha \neq 1}$ compute:

${\displaystyle X=(1+\zeta ^{2})^{\frac {1}{2\alpha }}{\frac {\sin\{\alpha (U+\xi )\}}{\{\cos(U)\}^{1/\alpha }}}\left[{\frac {\cos\{U-\alpha (U+\xi )\}}{W}}\right]^{\frac {1-\alpha }{\alpha }},}$

• for $\alpha =1$ compute:

${\displaystyle X={\frac {1}{\xi }}\left\{\left({\frac {\pi }{2}}+\beta U\right)\tan U-\beta \log \left({\frac {{\frac {\pi }{2}}W\cos U}{{\frac {\pi }{2}}+\beta U}}\right)\right\},}$

where

${\displaystyle \xi ={\begin{cases}{\frac {1}{\alpha }}\arctan(-\zeta ),&\alpha \neq 1,\\{\frac {\pi }{2}},&\alpha =1,\end{cases}}}$

and

${\displaystyle \zeta =-\beta \tan {\frac {\pi \alpha }{2}}}$.

This algorithm yields a random variable ${\displaystyle X\sim S_{\alpha }(\beta ,1,0)}$. For a detailed proof see.^[19]

Given the formulas for simulation of a standard stable random variable, we can easily simulate a stable random variable for all admissible values of the parameters $\alpha$, $c$, $\beta$ and $\mu$ using the following property. If ${\displaystyle X\sim S_{\alpha }(\beta ,1,0)}$ then

${\displaystyle Y={\begin{cases}cX+\mu ,&\alpha \neq 1,\\cX+{\frac {2}{\pi }}\beta c\log c+\mu ,&\alpha =1,\end{cases}}}$

is ${\displaystyle S_{\alpha }(\beta ,c,\mu )}$. It is interesting to note that for $\alpha =2$ (and $\beta =0$) the CMS method reduces to the well known Box-Muller transform for generating Gaussian random variables.^[20] Many other approaches have been proposed in the literature, including application of Bergström and LePage series expansions, see ^[21] and,^[22] respectively. However, the CMS method is regarded as the fastest and the most accurate.

Applications

Stable distributions owe their importance in both theory and practice to the generalization of the central limit theorem to random variables without second (and possibly first) order moments and the accompanying self-similarity of the stable family. It was the seeming departure from normality along with the demand for a self-similar model for financial data (i.e. the shape of the distribution for yearly asset price changes should resemble that of the constituent daily or monthly price changes) that led Benoît Mandelbrot to propose that cotton prices follow an alpha-stable distribution with α equal to 1.7.^[23] Lévy distributions are frequently found in analysis of critical behavior and financial data.^[7][24]

They are also found in spectroscopy as a general expression for a quasistatically pressure broadened spectral line.^[13]

The Lévy distribution of solar flare waiting time events (time between flare events) was demonstrated for CGRO BATSE hard x-ray solar flares in December 2001. Analysis of the Lévy statistical signature revealed that two different memory signatures were evident; one related to the solar cycle and the second whose origin appears to be associated with a localized or combination of localized solar active region effects.^[25]

Other analytic cases

A number of cases of analytically expressible stable distributions are known: In the following, the stable distribution is expressed as ${\displaystyle f(x;\alpha ,\beta ,c,\mu )}$. Besides the well-known distributions:

${\displaystyle f{\biggl (}x;1,0,1,0{\biggr )}}$ (the Cauchy Distribution)

${\displaystyle f{\biggl (}x;{\frac {1}{2}},-1,1,0{\biggr )}}$ (the Lévy distribution)

${\displaystyle f{\biggl (}x;2,0,1,0{\biggr )}}$ (the Normal distribution).

there are also:

${\displaystyle f{\biggl (}x;{\frac {1}{3}},0,1,0{\biggr )}=Re{\biggl (}{\frac {2\exp(-i\pi /4)}{3{\sqrt {3}}\pi }}x^{-3/2}S_{0,1/3}{\Bigl (}{\frac {2\exp(i\pi /4)}{3{\sqrt {3}}}}x^{-1/2}{\Bigr )}{\biggr )}}$

where ${\displaystyle S_{\mu ,\nu }(z)}$ is a Lommel function.^[26]

${\displaystyle f{\biggl (}x;{\frac {1}{2}},0,1,0{\biggr )}={\frac {\vert x\vert ^{-3/2}}{\sqrt {2\pi }}}{\Biggl (}\sin \left({\frac {1}{4\vert x\vert }}\right){\Biggl [}{\frac {1}{2}}-S{\biggl (}{\sqrt {\frac {1}{2\pi \vert x\vert }}}{\biggr )}{\Biggr ]}+\cos \left({\frac {1}{4\vert x\vert }}\right){\Biggl [}{\frac {1}{2}}-C{\biggl (}{\sqrt {\frac {1}{2\pi \vert x\vert }}}{\biggr )}{\Biggr ]}{\Biggr )}}$

where $S(x)$ and $C(x)$ are Fresnel Integrals. ^[27]

${\displaystyle f{\biggl (}x;{\frac {2}{3}},0,1,0{\biggr )}={\frac {1}{2{\sqrt {3\pi }}}}\vert x\vert ^{-1}\exp {\biggl (}{\frac {2}{27}}x^{-2}{\biggr )}W_{-1/2,1/6}{\biggl (}{\frac {4}{27}}x^{-2}{\biggr )}}$

where ${\displaystyle W_{k,\mu }(z)}$ is a Whittaker function. ^[28]

${\displaystyle f{\biggl (}x;{\frac {4}{3}},0,1,0{\biggr )}={\frac {3^{5/4}}{4{\sqrt {2\pi }}}}{\frac {\Gamma (7/12)\Gamma (11/12)}{\Gamma (6/12)\Gamma (8/12)}}\,_{2}F_{2}{\biggl (}{\frac {7}{12}},{\frac {11}{12}};{\frac {6}{12}},{\frac {8}{12}};{\frac {3^{3}x^{4}}{4^{4}}}{\biggr )}-{\frac {3^{11/4}x^{3}}{4^{3}{\sqrt {2\pi }}}}{\frac {\Gamma (13/12)\Gamma (17/12)}{\Gamma (18/12)\Gamma (15/12)}}\,_{2}F_{2}{\biggl (}{\frac {13}{12}},{\frac {17}{12}};{\frac {18}{12}},{\frac {15}{12}};{\frac {3^{3}x^{4}}{4^{4}}}{\biggr )}}$

where F is a hypergeometric function.^[26]

${\displaystyle f{\biggl (}x;{\frac {3}{2}},0,1,0{\biggr )}={\frac {1}{\pi }}\Gamma (5/3)\,_{2}F_{3}{\biggl (}{\frac {5}{12}},{\frac {11}{12}};{\frac {1}{3}},{\frac {1}{2}},{\frac {5}{6}};-{\frac {2^{2}x^{6}}{3^{6}}}{\biggr )}-{\frac {x^{2}}{3\pi }}\,_{3}F_{4}{\biggl (}{\frac {3}{4}},1,{\frac {5}{4}};{\frac {2}{3}},{\frac {5}{6}},{\frac {7}{6}},{\frac {4}{3}};-{\frac {2^{2}x^{6}}{3^{6}}}{\biggr )}+{\frac {7x^{4}}{3^{4}\pi ^{2}}}\Gamma (4/3)\,_{2}F_{3}{\biggl (}{\frac {13}{12}},{\frac {19}{12}};{\frac {7}{6}},{\frac {3}{2}},{\frac {5}{3}};-{\frac {2^{2}x^{6}}{3^{6}}}{\biggr )}}$

where F is a hypergeometric function. This is the Holtsmark distribution.^[26]

${\displaystyle f{\biggl (}x;{\frac {1}{3}},1,1,0{\biggr )}={\frac {1}{\pi }}{\frac {2{\sqrt {2}}}{3^{7/4}}}x^{-3/2}K_{1/3}{\Biggl (}{\frac {4{\sqrt {2}}}{3^{9/4}}}x^{-1/2}{\Biggr )}}$

where ${\displaystyle K_{v}(x)}$ is a modified Bessel function of the second kind.^[27]

${\displaystyle f{\biggl (}x;{\frac {2}{3}},1,1,0{\biggr )}={\frac {\sqrt {3}}{\sqrt {\pi }}}\vert x\vert ^{-1}\exp {\biggl (}-{\frac {16}{27}}x^{-2}{\biggr )}W_{1/2,1/6}{\biggl (}{\frac {32}{27}}x^{-2}{\biggr )}}$

where ${\displaystyle W_{k,\mu }(z)}$ is a Whittaker function.^[29]

${\displaystyle f{\biggl (}x;{\frac {3}{2}},1,1,0{\biggr )}=\left\{{\begin{array}{ll}{\frac {\sqrt {3}}{\sqrt {\pi }}}\vert x\vert ^{-1}\exp {\biggl (}{\frac {1}{27}}x^{3}{\biggr )}W_{1/2,1/6}{\biggl (}-{\frac {2}{27}}x^{3}{\biggr )}&x<0\\{\frac {1}{2{\sqrt {3\pi }}}}\vert x\vert ^{-1}\exp {\biggl (}{\frac {1}{27}}x^{3}{\biggr )}W_{-1/2,1/6}{\biggl (}{\frac {2}{27}}x^{3}{\biggr )}&x\geq 0\end{array}}\right.}$

where ${\displaystyle W_{k,\mu }(z)}$ is a Whittaker function.^[30]

See also

• Lévy flight
• Lévy process
• Fractional quantum mechanics
• Other "power law" distributions
  • Pareto distribution
  • Zeta distribution
  • Zipf distribution
  • Zipf–Mandelbrot distribution
• Stable and tempered stable distributions with volatility clustering – financial applications
• Multivariate stable distribution
• Discrete-stable distribution

Notes

• The STABLE program for Windows is available from John Nolan's stable webpage: http://academic2.american.edu/~jpnolan/stable/stable.html. It calculates the density (pdf), cumulative distribution function (cdf) and quantiles for a general stable distribution, and performs maximum likelihood estimation of stable parameters and some exploratory data analysis techniques for assessing the fit of a data set.
• Matlab codes for simulation of stable variables and estimation of stable parameters are available from RePEc: https://ideas.repec.org/e/pwe42.html#software

References

1. B. Mandelbrot, The Pareto-Lévy Law and the Distribution of Income, International Economic Review 1960 http://www.jstor.org/stable/2525289
2. Paul Lévy, Calcul des probabilités 1925
3. B.Mandelbrot, Stable Paretian Random Functions and the Multiplicative Variation of Income, Econometrica 1961 http://www.jstor.org/stable/pdfplus/1911802.pdf
4. B. Mandelbrot, The variation of certain Speculative Prices, The Journal of Business 1963
5. Eugene F. Fama, Mandelbrot and the Stable Paretian Hypothesis, The Journal of Business 1963
6. Nolan, John P. "Stable Distributions - Models for Heavy Tailed Data" (PDF). Retrieved 2009-02-21.
7. Voit, Johannes (2005). The Statistical Mechanics of Financial Markets - Springer. Springer. doi:10.1007/b137351.
8. B.V. Gnedenko, A.N. Kolmogorov. Limit distributions for sums of independent random variables, Cambridge, Addison-Wesley 1954 https://books.google.com/books/about/Limit_distributions_for_sums_of_independ.html?id=rYsZAQAAIAAJ&redir_esc=y
9. Vladimir V. Uchaikin, Vladimir M. Zolotarev, Chance and Stability: Stable Distributions and their Applications, De Gruyter 1999 https://books.google.com/books/about/Chance_and_Stability.html?id=Y0xiwAmkb_oC&redir_esc=y
10. Samorodnitsky, G.; Taqqu, M.S. (1994). Stable Non-Gaussian Random Processes: Stochastic Models with Infinite Variance. CRC Press. ISBN 9780412051715.
11. Lee, Wai Ha (2010). Continuous and discrete properties of stochastic processes. PhD thesis, University of Nottingham.
12. Zolotarev, V. (1995). "On Representation of Densities of Stable Laws by Special Functions". Theory of Probability & Its Applications. 39 (2): 354–362. doi:10.1137/1139025. ISSN 0040-585X.
13. Peach, G. (1981). "Theory of the pressure broadening and shift of spectral lines". Advances in Physics. 30 (3): 367–474. doi:10.1080/00018738100101467. ISSN 0001-8732.
14. Nolan, John P. (1997). "Numerical calculation of stable densities and distribution functions". Communications in Statistics. Stochastic Models. 13 (4): 759–774. doi:10.1080/15326349708807450. ISSN 0882-0287.
15. Matsui, Muneya; Takemura, Akimichi (2006). "Some Improvements in Numerical Evaluation of Symmetric Stable Density and Its Derivatives". Communications in Statistics - Theory and Methods. 35 (1): 149–172. doi:10.1080/03610920500439729. ISSN 0361-0926.
16. Chambers, J. M.; Mallows, C. L.; Stuck, B. W. (1976). "A Method for Simulating Stable Random Variables". Journal of the American Statistical Association. 71 (354): 340–344. doi:10.1080/01621459.1976.10480344. ISSN 0162-1459.
17. Zolotarev, V. M. (1986). One-Dimensional Stable Distributions. American Mathematical Society. ISBN 978-0-8218-4519-6.
18. Misiorek, Adam; Weron, Rafał (2012). Gentle, James E.; Härdle, Wolfgang Karl; Mori, Yuichi, eds. Heavy-Tailed Distributions in VaR Calculations. Springer Handbooks of Computational Statistics. Springer Berlin Heidelberg. pp. 1025–1059. doi:10.1007/978-3-642-21551-3_34. ISBN 978-3-642-21550-6.
19. Weron, Rafał (1996). "On the Chambers-Mallows-Stuck method for simulating skewed stable random variables". Statistics & Probability Letters. 28 (2): 165–171. doi:10.1016/0167-7152(95)00113-1.
20. Janicki, Aleksander; Weron, Aleksander (1994). Simulation and Chaotic Behavior of Alpha-stable Stochastic Processes. CRC Press. ISBN 9780824788827.
21. Mantegna, Rosario Nunzio (1994). "Fast, accurate algorithm for numerical simulation of L\'evy stable stochastic processes". Physical Review E. 49 (5): 4677–4683. doi:10.1103/PhysRevE.49.4677.
22. Janicki, Aleksander; Kokoszka, Piotr (1992). "Computer investigation of the Rate of Convergence of Lepage Type Series to α-Stable Random Variables". Statistics. 23 (4): 365–373. doi:10.1080/02331889208802383. ISSN 0233-1888.
23. Mandelbrot, B., New methods in statistical economics The Journal of Political Economy, 71 #5, 421-440 (1963).
24. Rachev, Svetlozar T.; Mittnik, Stefan (2000). Stable Paretian Models in Finance. Wiley. ISBN 978-0-471-95314-2.
25. Leddon, D., A statistical Study of Hard X-Ray Solar Flares
26. Garoni, T. M.; Frankel, N. E. (2002). "Lévy flights: Exact results and asymptotics beyond all orders". Journal of Mathematical Physic. 43 (5): 2670–2689. doi:10.1063/1.1467095.
27. Hopcraft, K. I.; Jakeman, E.; Tanner, R. M. J. (1999). "Lévy random walks with fluctuating step number and multiscale behavior". Physical Review E. 60 (5): 5327–5343. doi:10.1103/physreve.60.5327.
28. Uchaikin, V. V.; Zolotarev, V. M. (1999). "Chance And Stability - Stable Distributions And Their Applications". VSP. Utrecht, Netherlands.
29. Zlotarev, V. M. (1961). "Expression of the density of a stable distribution with exponent alpha greater than one by means of a frequency with exponent 1/alpha". Selected Translations in Mathematical Statistics and Probability (Translated from the Russian article: Dokl. Akad. Nauk SSSR. 98, 735-738 (1954)). 1: 163–167.
30. Zaliapin, I. V.; Kagan, Y. Y.; Schoenberg, F. P. (2005). "Approximating the Distribution of Pareto Sums". Pure and Applied Geophysics. 162 (6): 1187–1228. doi:10.1007/s00024-004-2666-3.