Asymmetric Laplace distribution

Asymmetric Laplace

Probability density function Asymmetric Laplace PDF with m = 0. Note that the κ =  2 and 1/2 curves are mirror images
Cumulative distribution function Asymmetric Laplace CDF with m = 0.
Parameters	$m$ location (real) $\lambda > 0$ scale (real) $\kappa > 0$ asymmetry (real)
Support	$x \in (-\infty; +\infty)\,$
PDF	(see article)
CDF	(see article)
Mean	$m+\frac{1-\kappa^2}{\lambda\kappa}$
Median	$m+(\kappa/\lambda)\log\left(\frac{1+\kappa^2}{2\kappa^2}\right)$
Variance	$\frac{1+\kappa^4}{\lambda^2\kappa^2}$
Skewness	$\frac{2 \left(1-\kappa ^6\right)}{\left(\kappa ^4+1\right)^{3/2}}$
Entropy	$\log\left(e\,\frac{1+\kappa^2}{\kappa\lambda}\right)$
CF	$\frac{ e^{i m t}}{(1+\frac{i t \kappa}{\lambda}) (1-\frac{i t}{\kappa\lambda})}$

In probability theory and statistics, the asymmetric Laplace distribution (ALD) is a continuous probability distribution which is a generalization of the Laplace distribution. Just as the Laplace distribution consists of two exponential distributions of equal scale back-to-back about x = m, the asymmetric Laplace consists of two exponential distributions of unequal scale back to back about x = m, adjusted to assure continuity and normalization. The difference of two variates exponentially distributed with different means and rate parameters will be distributed according to the ALD. When the two rate parameters are equal, the difference will be distributed according to the Laplace distribution.

Characterization

Probability density function

A random variable has an asymmetric Laplace(m, λ, κ) distribution if its probability density function is^[1][2]

$f(x;m,\lambda,\kappa)=\left(\frac{\lambda}{\kappa+1/\kappa}\right)\, e^{-(x-m)\lambda\, s \kappa^s}$

where s=sgn(x-m), or alternatively:

$f(x;m,\lambda,\kappa) = \frac{\lambda}{\kappa+1/\kappa} \begin{cases} \exp \left( (\lambda/\kappa)(x-m) \right) & \text{if }x < m \\[4pt] \exp ( -\lambda\kappa(x-m) ) & \text{if }x \geq m \end{cases}$

Here, m is a location parameter, λ > 0, is a scale parameter, and κ is an asymmetry parameter. When κ = 1, (x-m)s κ^s simplifies to |x-m| and the distribution simplifies to the Laplace distribution.

Cumulative distribution function

The cumulative distribution function is given by:

$F(x;m,\lambda,\kappa) = \begin{cases} \frac{\kappa^2}{1+\kappa^2}\exp ((\lambda/\kappa)(x-m)) & \text{if }x \leq m \\[4pt] 1-\frac{1}{1+\kappa^2} \exp (-\lambda\kappa(x-m)) & \text{if }x > m \end{cases}$

Characteristic function

The ALD characteristic function is given by:

$\varphi(t;m,\lambda,\kappa)=\frac{ e^{i m t}}{(1+\frac{i t \kappa}{\lambda}) (1-\frac{i t}{\kappa\lambda})}$

For m = 0, the ALD is a member of the family of geometric stable distributions with α = 2. It follows that if $\varphi_1$ and $\varphi_2$ are two distinct ALD characteristic functions with m = 0, then

$\varphi=\frac{1}{1/\varphi_1+1/\varphi_2-1}$

is also an ALD characteristic function with location parameter $m=0$. The new scale parameter λ obeys

$\frac{1}{\lambda^2}=\frac{1}{\lambda_1^2}+\frac{1}{\lambda_2^2}$

and the new skewness parameter κ obeys:

$\frac{\kappa^2-1}{\kappa\lambda}=\frac{\kappa_1^2-1}{\kappa_1\lambda_1}+\frac{\kappa_2^2-1}{\kappa_2\lambda_2}$

Moments, mean, variance, skewness

The n-th moment of the ALD about m is given by

$E[(x-m)^n]=\frac{n!}{\lambda^n(\kappa+1/\kappa)}\,(\kappa^{-(n+1)}-(-\kappa)^{n+1})$

From the binomial theorem, the n-th moment about zero (for m not zero) is then:

$E[x^n]=\frac{\lambda\,m^{n+1}}{\kappa+1/\kappa}\,\left( \sum_{i=0}^n\frac{n!}{(n-i)!}\,\frac{1}{(m\lambda\kappa)^{i+1}} -\sum_{i=0}^n\frac{n!}{(n-i)!}\,\frac{1}{(-m\lambda/\kappa)^{i+1}} \right)$

$= \frac{ \lambda \, m^{n+1}}{\kappa+1/\kappa} \left( e^{m \lambda\kappa} E_{-n}(m \lambda\kappa ) -e^{-m \lambda/\kappa } E_{-n}(-m \lambda /\kappa) \right)$

where $E_n( )$ is the generalized exponential integral function $E_n(x)=x^{n-1}\Gamma(1-n,x)$

The first moment about zero is the mean:

$\mu=E[x]=m-\frac{\kappa-1/\kappa}{\lambda}$

The variance is:

$\sigma^2=E[x^2]-\mu^2=\frac{1+\kappa^4}{\kappa^2\lambda^2}$

and the skewness is:

$\frac{E[x^3]-3\mu\sigma^2-\mu^3}{\sigma^3}=\frac{2 \left(1-\kappa ^6\right)}{\left(\kappa ^4+1\right)^{3/2}}$

Generating asymmetric Laplace variates

Asymmetric Laplace variates (X) may be generated from a random variate U drawn from the uniform distribution in the interval (-κ,1/κ) by:

$X=m-\frac{1}{\lambda\,s\kappa^s}\log(U\,s\kappa^S)$

where s=sgn(U).

They may also be generated as the difference of two exponential distributions. If X₁ is drawn from exponential distribution with mean and rate (m₁,λκ) and X₂ is drawn from an exponential distribution with mean and rate (m₂,λ/κ) then X₁ - X₂ is distributed according to the asymmetric Laplace distribution with parameters (m1-m2, λ, κ)

Entropy

The differential entropy of the ALD is

$H=-\int_{-\infty}^\infty f_{AL}(x)\log(f_{AL}(x)) dx = 1-\log\left(\frac{\lambda}{\kappa+1/\kappa}\right)$

The ALD has the maximum entropy of all distributions with a fixed value (1/λ) of $(x-m)\,s\kappa^s$ where $s=\sgn(x-m)$.

References

1. Kozubowski, Tomasz J.; Podgorski, Krzysztof (2000). "A Multivariate and Asymmetric Generalization of Laplace Distribution" (PDF). Computational Statistics 15: 531. Retrieved 2015-12-29.
2. Jammalamadaka, S. Rao; Kozubowski, Tomasz J. (2004). "New Families of Wrapped Distributions for Modeling Skew Circular Data" (PDF). Communications in Statistics – Theory and Methods 33 (9): 2059–2074. doi:10.1081/STA-200026570. Retrieved 2011-06-13.