Probability density function edit

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

${\displaystyle 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:

${\displaystyle 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 edit

The cumulative distribution function is given by:

${\displaystyle 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 edit

The ALD characteristic function is given by:

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

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

${\displaystyle \varphi ={\frac {1}{1/\varphi _{1}+1/\varphi _{2}-1}}}$ 

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

${\displaystyle {\frac {1}{\lambda ^{2}}}={\frac {1}{\lambda _{1}^{2}}}+{\frac {1}{\lambda _{2}^{2}}}}$ 

and the new skewness parameter κ obeys:

${\displaystyle {\frac {\kappa ^{2}-1}{\kappa \lambda }}={\frac {\kappa _{1}^{2}-1}{\kappa _{1}\lambda _{1}}}+{\frac {\kappa _{2}^{2}-1}{\kappa _{2}\lambda _{2}}}}$ 

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

${\displaystyle 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:

${\displaystyle 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)}$ 

${\displaystyle ={\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 ${\displaystyle E_{n}()}$  is the generalized exponential integral function ${\displaystyle E_{n}(x)=x^{n-1}\Gamma (1-n,x)}$ 

The first moment about zero is the mean:

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

The variance is:

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

and the skewness is:

${\displaystyle {\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}}}}$ 

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

${\displaystyle X=m-{\frac {1}{\lambda \,s\kappa ^{s}}}\log(1-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, λ, κ)

The differential entropy of the ALD is

${\displaystyle 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 ${\displaystyle (x-m)\,s\kappa ^{s}}$  where ${\displaystyle s=\operatorname {sgn}(x-m)}$ .

An alternative parametrization is made possible by the characteristic function:

${\displaystyle \varphi (t;\mu ,\sigma ,\beta )={\frac {e^{i\mu t}}{1-i\beta \sigma t+\sigma ^{2}t^{2}}}}$ 

where ${\displaystyle \mu }$  is a location parameter, ${\displaystyle \sigma }$  is a scale parameter, ${\displaystyle \beta }$  is an asymmetry parameter. This is specified in Section 2.6.1 and Section 3.1 of Lihn (2015). ^[3] Its probability density function is

${\displaystyle f(x;\mu ,\sigma ,\beta )={\frac {1}{2\sigma B_{0}}}{\begin{cases}\exp \left({\frac {x-\mu }{\sigma B^{-}}}\right)&{\text{if }}x<\mu \\[4pt]\exp(-{\frac {x-\mu }{\sigma B^{+}}})&{\text{if }}x\geq \mu \end{cases}}}$ 

where ${\displaystyle B_{0}={\sqrt {1+\beta ^{2}/4}}}$  and ${\displaystyle B^{\pm }=B_{0}\pm \beta /2}$ . It follows that ${\displaystyle B^{+}B^{-}=1,\P B^{+}-B^{-}=\beta }$ .

The n-th moment about ${\displaystyle \mu }$  is given by

${\displaystyle E[(x-\mu )^{n}]={\frac {\sigma ^{n}n!}{2B_{0}}}((B^{+})^{n+1}+(-1)^{n}(B^{-})^{n+1})}$ 

The mean about zero is:

${\displaystyle E[x]=\mu +\sigma \beta }$ 

The variance is:

${\displaystyle E[x^{2}]-E[x]^{2}=\sigma ^{2}(2+\beta ^{2})}$ 

The skewness is:

${\displaystyle {\frac {2\beta (3+\beta ^{2})}{(2+\beta ^{2})^{3/2}}}}$ 

The excess kurtosis is:

${\displaystyle {\frac {6(2+4\beta ^{2}+\beta ^{4})}{(2+\beta ^{2})^{2}}}}$ 

For small ${\displaystyle \beta }$ , the skewness is about ${\displaystyle 3\beta /{\sqrt {2}}}$  . Thus ${\displaystyle \beta }$  represents skewness in an almost direct way.

The Asymmetric Laplace distribution is commonly used with an alternative parameterization for performing quantile regression in a Bayesian inference context.^[4] Under this approach, the ${\displaystyle \kappa }$  parameter describing asymmetry is replaced with a ${\displaystyle p}$  parameter indicating the percentile or quantile desired. Using this parameterization, the likelihood of the Asymmetric Laplace Distribution is equivalent to the loss function employed in quantile regression. With this alternative parameterization, the probability density function is defined as:

${\displaystyle f(x;m,\lambda ,p)={\frac {p(1-p)}{\lambda }}{\begin{cases}\exp \left(-((p-1)/\lambda )(x-m)\right)&{\text{if }}x\leq m\\[4pt]\exp(-(p/\lambda )(x-m))&{\text{if }}x>m\end{cases}}}$ 

Where, m is a location parameter, λ > 0 is a scale parameter, and 0 < p < 1 is a percentile parameter.

The mean (${\displaystyle \mu }$ ) and variance (${\displaystyle \sigma ^{2}}$ ) are calculated as:

${\displaystyle \mu =m+{\frac {1-2p}{p(1-p)}}\lambda }$ 

${\displaystyle \sigma ^{2}={\frac {1-2p+2p^{2}}{p^{2}(1-p)^{2}}}\lambda ^{2}}$ 

The cumulative distribution function is given^[5] by:

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

The asymmetric Laplace distribution has applications in finance and neuroscience. For the example in finance, S.G. Kou developed a model for financial instrument prices incorporating an asymmetric Laplace distribution to address problems of skewness, kurtosis and the volatility smile that often occur when using a normal distribution for pricing these instruments.^[6] Another example is in neuroscience in which its convolution with normal distribution is considered as a model for brain stopping reaction times.^[7]