Jump to content

User:Stevelihn/sandbox/Lihn's Lambda Distribution

From Wikipedia, the free encyclopedia

Lihn's lambda distribution is a family of parametric continuous probability distributions on the real line. It is also called "version 1" of "generalized normal distribution", the exponential power distribution, or the generalized error distribution, with a different parametrization, but such nomenclature has not been standardized yet.

This distribution has been known as early as 1937.[1] This page serves as a reference to the particular parametrization and notation used in Lihn's study (from 2015 to 2018). Lihn's contribution is to reveal the richness of this distribution family as following.

  1. Connected this distribution family with the stable law via lambda decomposition;
  2. Established it as a continuous leptokurtic distribution in which the normal distribution is a special case;
  3. Applied this distribution family in the Hidden Markov Model, where the lepkurtotic nature of financial data can be better captured;
  4. Established it as a stationary solution of leptokurtic extension of the Ornstein–Uhlenbeck process;
  5. Derived an analytic solution for S&P 500 option that contains the feature of volatility smile;
  6. As a special case of a larger "stable lambda distribution" family, which aims to describe the daily returns distribution of financial data more suitably.

Definition

[edit]
Lihn's Lambda
Probability density function
Probability density plots of generalized normal distributions
Cumulative distribution function
Cumulative distribution function plots of generalized normal distributions
Parameters

- location (real)
- scale (positive, real)

- shape (positive, real)
Support
PDF
CDF
Mean
Median
Mode
Variance
Skewness 0
Excess kurtosis
MGF not analytically expressible, subject to tail truncation

Although Lihn studied both the symmetric and asymmetric distributions[2], only the symmetric formulation is mentioned here since the sequence of studies have showed that the symmetric formulation is far more important in applications.

The probability density function (PDF) is defined as[3]

where is the location parameter, is the scale parameter, and is the shape parameter.

The cumulative distribution function (CDF) is

where is the upper incomplete gamma function.

The parameter is called the order of the distribution, which is structured such that its integer values bear important meaning. It is connected to the stability index in the stable distribution via . Using the inverse of stability index makes many formulas cleaner for this distribution, especially when gamma functions are involved.

It is easy to see that it becomes a normal distribution when , and a Laplace distribution when . When , certain analytic solutions exist. Most importantly, is called "quartic lambda distribution", where there are many elegant analytic solutions.

This family allows for tails that are either heavier than normal (when ) or lighter than normal (when ). It is a useful way to parametrize a continuum of symmetric, platykurtic densities spanning from the normal () to the uniform density (), and a continuum of symmetric, leptokurtic densities spanning from the normal density to Laplace and beyond ().

The following table lists the variance and kurtosis of the most important integers :

Variance and Kurtosis
var kurtosis
1 3
2 6
3 12.257
4 25.2
The kurtosis is approximately in this limited range.


Properties

[edit]

Moments

[edit]

The odd moments are zero. The n-th moment is

Thus the variance is and the kurtosis is .

Moment generating function

[edit]

When all the moments are known, the moment generating function (MGF) can be calculated as

However, both the integral and the sum diverge when . Lihn attributed this divergence as the root cause of "moment explosion", which has been a major obstacle in many stochastic volatility models (See Section 2.4 of [2]). This regime is called "the local regime", where MGF solution is available as long as is reasonably small. The "tail" of the integral and the sum must be truncated. The procedures are described as following.

Truncation of MGF integral
[edit]

Assume in the following discussion. The truncation point in the right tail, that is, the upper limit of the integral, is where the integrand reaches its minimum: This leads to the solution, (See Sections 3.1 of [2])

Truncation of MGF summation
[edit]

The truncation point of the summation is where the summand reaches its minimum: This leads to the equation, (See Sections 2.5 of [2])

where is a digamma function. It has a large-n solution when ,

We note that , in other words, they are equal in a dimensionless setting when .

In Lihn's work of option pricing model, the risk-neutral drift is defined via MGF,

Therefore, it is essential that can be calculated properly in order to price an option. The meaning of can be illustrated by the trivial case of a normal distribution, where leads to . This is minus half of the variance, the well-known drift of a geometric Brownian motion.


Connection to the stable law

[edit]

Relation to symmetric alpha-stable distribution

[edit]

TBD

Relation to one-sided stable distribution

[edit]

TBD

Lambda Decomposition

[edit]

TBD


Applications

[edit]

This distribution can be used in modeling when the non-normality is prevailing. For instance, financial data is well known to violate normality almost all the time.

Leptokurtic extension of the Ornstein–Uhlenbeck process

[edit]

TBD

Leptokurtic extension of Hidden Markov Model

[edit]

TBD

Generalization - Stable Lambda Distribution

[edit]

TBD.


  1. ^ Bochner, Salomon (1937). "Stable laws of probability and completely monotone functions". Duke Mathematical Journal. 3 (4): 726--728.
  2. ^ a b c d Lihn, Stephen H. T. (2015-12-23). "The Special Elliptic Option Pricing Model and Volatility Smile". Rochester, NY. doi:10.2139/ssrn.2707810. {{cite journal}}: Cite journal requires |journal= (help)
  3. ^ Lihn, Stephen H. T. (2017-11-15). "A Theory of Asset Return and Volatility Under Stable Law and Stable Lambda Distribution". Rochester, NY. {{cite journal}}: Cite journal requires |journal= (help)