Draft:Q-distribution

From Wikipedia, the free encyclopedia
Jump to navigation Jump to search

An asymmetric and heavy tailed probability distribution, Zhao's q-distribution [1] is a four parameter distribution for a real valued random variable. The distribution is proposed for modeling financial asset returns, which usually have heavy tails than Gaussian and also exhibit with asymmetry.

The four parameters of the q-distribution model explicitly represent the location, scaling parameters of positive and negative deviations and heaviness of tails. When one or two parameters reach their boundary, e.g. zero or infinity, the limiting distributions can be identical to one of the classical distributions. For instance,Gaussian, Laplace, logistic, and exponential are among the known special cases [2], which have numerous applications in many scientific fields.

Definition

[edit | edit source]

The q-distribution has a probability density function [1]

f(x)=1C(α,β,θ)[(1+exμαθ)(1+exμβθ)]θ

for x, where C(α,β,θ) is a compensating factor and it often needs to be computed numerically except for the symmetric cases when θ is an integer.

Density function of a symmetric q-distribution versus Gaussian distribution
Density of an asymmetric q-distribution versus Gaussian distribution

The q-distribution has four parameters:

  • μ - location parameter
  • α>0 - scaling parameter for positive domain or right tail
  • β>0 - scaling parameter for negative domain or left tail
  • θ>0 - shape parameter that controls kurtosis

Except the location parameter, all the other three parameters are positive.


Symmetric case

[edit | edit source]

A symmetric q-distribution arrives naturally when the scaling parameters α,β are identical. In this case, the probability density function is simplified to

f(x)=1αθB(θ,θ)[(1+exμαθ)(1+exμαθ)]θ

where the scaling parameters are equal and the compensator can be written as the Beta function. When the kurtosis parameter θ is an integer, the compensator can be easily calculated. For instance, when θ=1,B(1,1)=1, one gets the logistic distribution. The symmetric q-distribution can be directly derived from Beta distribution [2] though variable exchange with a logistic function.

The q2 distribution

[edit | edit source]

A very special case of the q-distribution is when θ=2, we call it a q2 distribution. It mimics the Gaussian distribution in many ways. The density function is simple and beautiful

f(x)=3α[(1+exμ2α)(1+exμ2α)]2

which is symmetric with exponential tails. Its MLE estimator is similar to Huber's [3] M-estimation, which is known to be robust and insensitive to outlier data points.

Other limiting cases

[edit | edit source]

A few classical distributions can be derived from q-distribution. These are the boundary cases when one or more of the parameters reaches zero or infinity.

f(x)12πσe(xμ)22σ2
f(x)1α[(1+e(xμ)α)(1+e(xμ)α)]1

When αβ, it defines an asymmetric logistic distribution[4], which is new and unseen from the literature.

f(x)1α+βe(xμ)+α(xμ)β

This is different from the definition from Kotz et al. [5] When α=β, it becomes the Laplace distribution.

f(x)1αe(xμ)+α

Properties

[edit | edit source]

The q-distribution has the following known properties:

  • Skewness ranges from 2 to 2.
  • Kurtosis ranges from 3 to 9.
  • Finite moments for any n>0
E|X|n<.
  • Finite expectation for λ<α1
E(eλX)<.

Computational method

[edit | edit source]

With a data sample {x1,x2,,xn}, the parameters of q-distribution can be estimated using maximum likelihood method.

In general, for asymmetric q-distribution, there is no closed formula for the cumulative distribution function (c.d.f.) and neither the compensator C(α,β,θ). However, expectation related to q-distribution can be computed effectively using Gaussian quadrature approach, because the distribution has exponential tails.

For any smooth function, the following integral can be effectively computed using Gauss–Laguerre quadrature[6] method:

0exG(x)=i=1kwiG(xi)

where {xi,wi} are the quadrature abscissas and weights.

Any expectation with respect to q-density can be written into the above format by variable exchanges [1] This makes it convenient to compute the compensator C(α,β,θ) as well as the c.d.f function and any other types of expected values for the related application.

Applications

[edit | edit source]

For many applications where Gaussian distribution is used, but is not perfect because the data set presents either asymmetry or long tails,  the q-distribution can jump in as a helper.

The four parameters of q-distribution have the flexibility to match moments up to the fourth order instead of the second order for Gaussian.  The distribution has potential applications in risk management, e.g. projection of tail risk, and many other scientific fields.


See also

[edit | edit source]

The following classical distributions are related to the q-distribution. Their properties are well studied.

References

[edit | edit source]
  1. ^ a b c Lua error in Module:Citation/CS1/Configuration at line 2172: attempt to index field '?' (a nil value).
  2. ^ a b Johnson N. L. and S. Lotz. (1970). Continuous univariate distributions I-II. John Wiley and Sons.
  3. ^ Lua error in Module:Citation/CS1/Configuration at line 2172: attempt to index field '?' (a nil value).
  4. ^ Lua error in Module:Citation/CS1/Configuration at line 2172: attempt to index field '?' (a nil value).
  5. ^ Lua error in Module:Citation/CS1/Configuration at line 2172: attempt to index field '?' (a nil value).
  6. ^ Lua error in Module:Citation/CS1/Configuration at line 2172: attempt to index field '?' (a nil value).