Kalpana Kalpana (Editor)

Skew normal distribution

Updated on
Edit
Like
Comment
Share on FacebookTweet on TwitterShare on LinkedInShare on Reddit
Skew normal distribution

Parameters
  
ξ {displaystyle xi ,} location (real) ω {displaystyle omega ,} scale (positive, real) α {displaystyle alpha ,} shape (real)

Support
  
x ∈ ( − ∞ ; + ∞ ) {displaystyle xin (-infty ;+infty )!}

PDF
  
1 ω π e − ( x − ξ ) 2 2 ω 2 ∫ − ∞ α ( x − ξ ω ) e − t 2 2   d t {displaystyle { rac {1}{omega pi }}e^{-{ rac {(x-xi )^{2}}{2omega ^{2}}}}int _{-infty }^{alpha left({ rac {x-xi }{omega }} ight)}e^{-{ rac {t^{2}}{2}}} dt}

CDF
  
Φ ( x − ξ ω ) − 2 T ( x − ξ ω , α ) {displaystyle Phi left({ rac {x-xi }{omega }} ight)-2Tleft({ rac {x-xi }{omega }},alpha ight)} T ( h , a ) {displaystyle T(h,a)} is Owen's T function

Mean
  
ξ + ω δ 2 π {displaystyle xi +omega delta {sqrt { rac {2}{pi }}}} where δ = α 1 + α 2 {displaystyle delta ={ rac {alpha }{sqrt {1+alpha ^{2}}}}}

Variance
  
ω 2 ( 1 − 2 δ 2 π ) {displaystyle omega ^{2}left(1-{ rac {2delta ^{2}}{pi }} ight)}

In probability theory and statistics, the skew normal distribution is a continuous probability distribution that generalises the normal distribution to allow for non-zero skewness.

Contents

Definition

Let ϕ ( x ) denote the standard normal probability density function

ϕ ( x ) = 1 2 π e − x 2 2

with the cumulative distribution function given by

Φ ( x ) = ∫ − ∞ x ϕ ( t )   d t = 1 2 [ 1 + erf ⁡ ( x 2 ) ] ,

where erf is the error function. Then the probability density function (pdf) of the skew-normal distribution with parameter α is given by

f ( x ) = 2 ϕ ( x ) Φ ( α x ) .

This distribution was first introduced by O'Hagan and Leonard (1976). A popular alternative parameterization is due to Mudholkar and Hutson (2000), which has a form of the c.d.f. that is easily inverted such that there is a closed form solution to the quantile function.

A stochastic process that underpins the distribution was described by Andel, Netuka and Zvara (1984). Both the distribution and its stochastic process underpinnings were consequences of the symmetry argument developed in Chan and Tong (1986), which applies to multivariate cases beyond normality, e.g. skew multivariate t distribution and others. The distribution is a particular case of a general class of distributions with probability density functions of the form f(x)=2 φ(x) Φ(x) where φ() is any PDF symmetric about zero and Φ() is any CDF whose PDF is symmetric about zero.

To add location and scale parameters to this, one makes the usual transform x → x − ξ ω . One can verify that the normal distribution is recovered when α = 0 , and that the absolute value of the skewness increases as the absolute value of α increases. The distribution is right skewed if α > 0 and is left skewed if α < 0 . The probability density function with location ξ , scale ω , and parameter α becomes

f ( x ) = 2 ω ϕ ( x − ξ ω ) Φ ( α ( x − ξ ω ) ) .

Note, however, that the skewness of the distribution is limited to the interval ( − 1 , 1 ) .

Estimation

Maximum likelihood estimates for ξ , ω , and α can be computed numerically, but no closed-form expression for the estimates is available unless α = 0 . If a closed-form expression is needed, the method of moments can be applied to estimate α from the sample skew, by inverting the skewness equation. This yields the estimate

| δ | = π 2 | γ ^ 1 | 2 3 | γ ^ 1 | 2 3 + ( ( 4 − π ) / 2 ) 2 3

where δ = α 1 + α 2 , and γ ^ 1 is the sample skew. The sign of δ is the same as the sign of γ ^ 1 . Consequently, α ^ = δ / 1 − δ 2 .

The maximum (theoretical) skewness is obtained by setting δ = 1 in the skewness equation, giving γ 1 ≈ 0.9952717 . However it is possible that the sample skewness is larger, and then α cannot be determined from these equations. When using the method of moments in an automatic fashion, for example to give starting values for maximum likelihood iteration, one should therefore let (for example) | γ ^ 1 | = min ( 0.99 , | ( 1 / n ) ∑ ( ( x i − x ¯ ) / s ) 3 | ) .

Concern has been expressed about the impact of skew normal methods on the reliability of inferences based upon them.

Differential equation

The differential equation leading to the pdf of the skew normal distribution is

ω 4 f ″ ( x ) + ( α 2 + 2 ) ω 2 ( x − ξ ) f ′ ( x ) + f ( x ) ( ( α 2 + 1 ) ( x − ξ ) 2 + ω 2 ) = 0 ,

with initial conditions

f ( 0 ) = exp ⁡ ( − ξ 2 2 ω 2 ) erfc ⁡ ( α ξ 2 ω ) 2 π ω  and f ′ ( 0 ) = exp ⁡ ( − ( α 2 + 1 ) ξ 2 2 ω 2 ) ( 2 α ω + 2 π ξ exp ⁡ ( α 2 ξ 2 2 ω 2 ) erfc ⁡ ( α ξ 2 ω ) ) 2 π ω 3 .

References

Skew normal distribution Wikipedia
(Text) CC BY-SA