q-exponential distribution

From HandWiki
q-exponential distribution
Probability density function
Probability density plots of q-exponential distributions
Parameters q<2 shape (real)
λ>0 rate (real)
Support x[0,) for q1
x[0,1λ(1q)) for q<1
PDF (2q)λeqλx
CDF 1eqλx/q where q=12q
Mean 1λ(32q) for q<32
Otherwise undefined
Median qlnq(1/2)λ where q=12q
Mode 0
Variance q2(2q3)2(3q4)λ2 for q<43
Skewness 254q3q4q2 for q<54
Kurtosis 64q3+17q220q+6(q2)(4q5)(5q6) for q<65

The q-exponential distribution is a probability distribution arising from the maximization of the Tsallis entropy under appropriate constraints, including constraining the domain to be positive. It is one example of a Tsallis distribution. The q-exponential is a generalization of the exponential distribution in the same way that Tsallis entropy is a generalization of standard Boltzmann–Gibbs entropy or Shannon entropy.[1] The exponential distribution is recovered as q1.

Originally proposed by the statisticians George Box and David Cox in 1964,[2] and known as the reverse Box–Cox transformation for q=1λ, a particular case of power transform in statistics.

Characterization

Probability density function

The q-exponential distribution has the probability density function

(2q)λeq(λx)

where

eq(x)=[1+(1q)x]1/(1q)

is the q-exponential if q ≠ 1. When q = 1, eq(x) is just exp(x).

Derivation

In a similar procedure to how the exponential distribution can be derived (using the standard Boltzmann–Gibbs entropy or Shannon entropy and constraining the domain of the variable to be positive), the q-exponential distribution can be derived from a maximization of the Tsallis Entropy subject to the appropriate constraints.

Relationship to other distributions

The q-exponential is a special case of the generalized Pareto distribution where

μ=0,ξ=q12q,σ=1λ(2q).

The q-exponential is the generalization of the Lomax distribution (Pareto Type II), as it extends this distribution to the cases of finite support. The Lomax parameters are:

α=2qq1,λLomax=1λ(q1).

As the Lomax distribution is a shifted version of the Pareto distribution, the q-exponential is a shifted reparameterized generalization of the Pareto. When q > 1, the q-exponential is equivalent to the Pareto shifted to have support starting at zero. Specifically, if

XqExp(q,λ) and Y[Pareto(xm=1λ(q1),α=2qq1)xm],

then XY.

Generating random deviates

Random deviates can be drawn using inverse transform sampling. Given a variable U that is uniformly distributed on the interval (0,1), then

X=qlnq(U)λqExp(q,λ)

where lnq is the q-logarithm and q=12q.

Applications

Being a power transform, it is a usual technique in statistics for stabilizing the variance, making the data more normal distribution-like and improving the validity of measures of association such as the Pearson correlation between variables. It has been found to be an accurate model for train delays.[3] It is also found in atomic physics and quantum optics, for example processes of molecular condensate creation via transition through the Feshbach resonance.[4]

See also

Notes

  1. Tsallis, C. Nonadditive entropy and nonextensive statistical mechanics-an overview after 20 years. Braz. J. Phys. 2009, 39, 337–356
  2. Box, George E. P.; Cox, D. R. (1964). "An analysis of transformations". Journal of the Royal Statistical Society, Series B 26 (2): 211–252. 
  3. Keith Briggs and Christian Beck (2007). "Modelling train delays with q-exponential functions". Physica A 378 (2): 498–504. doi:10.1016/j.physa.2006.11.084. Bibcode2007PhyA..378..498B. 
  4. C. Sun; N. A. Sinitsyn (2016). "Landau-Zener extension of the Tavis-Cummings model: Structure of the solution". Phys. Rev. A 94 (3): 033808. doi:10.1103/PhysRevA.94.033808. Bibcode2016PhRvA..94c3808S. 

Further reading