Davis distribution: Difference between revisions

From formulasearchengine
Jump to navigation Jump to search
en>Helpful Pixie Bot
m ISBNs (Build KC)
 
en>BeyondNormality
Line 1: Line 1:
Got nothing to write about myself at all.<br>Finally a part of this site.<br>I just hope Im useful at all
In [[probability theory]] and [[statistics]], the '''split normal distribution''' also known as the '''two-piece normal distribution''' results from joining at the mode the corresponding halves of two [[normal distributions]] with the same mode but different variances.  The split normal distribution was introduced by Gibbons and Mylroie<ref name=Gibbons/> and by John.<ref name="John1982" />
 
{{Probability distribution
  | name      = Split-normal
  | type      = density
<!--  | pdf_image  = [[Image:Normal Distribution PDF.svg|350px|Probability density function for the normal distribution]]<br /><small>The red curve is the ''standard normal distribution''</small>
  | cdf_image  = [[Image:Normal Distribution CDF.svg|350px|Cumulative distribution function for the normal distribution]]<br /> -->
  | notation  = <math>\mathcal{SN}(\mu,\,\sigma_1,\sigma_2)</math>
  | parameters = <math>\mu \in \Re </math> — [[Mode (statistics)|mode]] ([[location parameter|location]], [[real number|real]])<br/><math>\sigma_1 > 0 </math> — left-hand-side [[standard deviation]] ([[scale parameter|scale]], [[real number|real]])<br/><math>\sigma_2 > 0 </math> — right-hand-side [[standard deviation]] ([[scale parameter|scale]], [[real number|real]])
  | support    = <math>x \in \Re </math>|
  | pdf        = <math> A \exp \left(- \frac {(x-\mu)^2}{2 \sigma_1^2}\right) \quad \text{if } x< \mu</math>
<br/><math> A \exp \left(- \frac {(x-\mu)^2}{2 \sigma_2^2}\right) \quad \text{otherwise,}</math>
<br/><math> \text{where} \quad A= \sqrt{2/\pi} (\sigma_1+\sigma_2)^{-1}</math>
  | mode      = <math> \mu </math>
  | mean      = <math> \mu+\sqrt{2 / \pi}(\sigma_2-\sigma_1)</math>
  | variance  = <math>(1-2 / \pi)(\sigma_2-\sigma_1)^2 + \sigma_1 \sigma_2</math>
  | skewness  = <math> \gamma_3 = \sqrt{\frac{2}{\pi}}(\sigma_2-\sigma_1)\left[\left(\frac{4}{\pi}-1\right)(\sigma_2-\sigma_1)^2 + \sigma_1 \sigma_2\right]</math>
  }}
 
== Definition ==
The split normal distribution arises from merging two opposite halves of two [[probability density function]]s (PDFs) of [[normal distribution]]s in their common [[Mode (statistics)|mode]].
 
The PDF of the split normal distribution is given by<ref name=Johnson1994/>
 
:<math>f(x;\mu,\sigma_1,\sigma_2)= A \exp (- \frac {(x-\mu)^2}{2 \sigma_1^2}) \quad \text{if } x< \mu
</math>
:<math> f(x;\mu,\sigma_1,\sigma_2)=  A \exp (- \frac {(x-\mu)^2}{2 \sigma_2^2}) \quad \text{otherwise}
</math>
where
:<math>\quad A= \sqrt{2/\pi} (\sigma_1+\sigma_2)^{-1}.</math>
 
=== Discussion ===
The split normal distribution results from merging two halves of normal distributions. In a general case the 'parent' normal distributions can have different variances which implies that the joined PDF would not be [[continuous function|continuous]]. To ensure that the resulting PDF [[Integral|integrates]] to 1, the [[normalizing constant]] '''A''' is used.
 
In a special case when <math>\sigma_1^2=\sigma_2^2=\sigma_{*}^2</math> the split normal distribution reduces to [[normal distribution]] with variance <math>\sigma_{*}^2</math>.
 
When σ<sub>2</sub>≠σ<sub>1</sub> the constant '''A''' it is different from the constant of normal distribution. However, when <math>\sigma_1^2=\sigma_2^2=\sigma_{*}^2</math> the constants are equal.
 
The sign of its third central moment is determined by the difference (σ<sub>2</sub>-σ<sub>1</sub>). If this difference is positive, the distribution is skewed to the right and if negative, then it is skewed to the left.
 
Other properties of the split normal density were discussed by Johnson et al.<ref name="Johnson1994" /> and Julio.<ref name="Julio2007" />
 
== Alternative formulations ==
The formulation discussed above originates from John.<ref name="John1982" /> The literature offers two mathematically equivalent alternative parameterizations . Britton, Fisher and Whitley<ref name="BFW1998" /> offer a parameterization if terms of mode, dispersion and normed skewness, denoted with <math>\mathcal{SN}(\mu,\, \sigma^2,\gamma)</math>. The parameter μ is the mode and has equivalent to the mode in John’s formulation. The parameter σ <sup>2</sup>>0 informs about the dispersion (scale) and should not be confused with variance. The third parameter, γ ∈ (-1,1), is the normalized  skew.
 
The second alternative parameterization is used in the [[Bank of England|Bank of England’s]] communication and is written in terms of mode, dispersion and unnormed skewness and is denoted with <math>\mathcal{SN}(\mu,\, \sigma^2,\xi)</math>. In this formulation  the parameter μ is the mode and is identical as in John’s <ref name="John1982" /> and Britton, Fisher and Whitley’s <ref name="BFW1998" /> formulation. The parameter σ <sup>2</sup> informs about the dispersion (scale) and is the same as in the Britton, Fisher and Whitley’s formulation. The parameter ξ equals the difference between the distribution’s mean and mode and can be viewed as unnormed measure of skewness.
 
The three parameterizations are mathematically equivalent, meaning that there is a strict relationship between the parameters and that it is possible to go from one parameterization to another. The following relationships hold:<ref name="BD2011" />
:<math>\begin{align}
\sigma^2 &= \sigma_1^2(1+\gamma)= \sigma_2^2(1-\gamma) \\
\gamma  &= \frac{\sigma_2-\sigma_1}{\sigma_2+\sigma_1} \\
\xi            &=\sqrt{2 / \pi}(\sigma_2-\sigma_1) \\
\gamma  &= \operatorname{sgn}(\xi) \sqrt{1-\left( \frac{\sqrt{1+2\beta}-1}{\beta} \right)^2}, \quad \text{where} \quad \beta = \frac{\pi\xi^2}{2\sigma^2}.
\end{align}</math>
 
== Multivariate Extensions ==
The multivariate generalization of the split normal distribution was proposed by Villani and Larsson.<ref name="VL2006"/> They assume that each of the [[Principal component analysis|principal components]] has univariate split normal distribution with a different set of parameters μ, σ<sub>2</sub> and σ<sub>1</sub>.
 
== Estimation of parameters ==
John<ref name="John1982" /> proposes to estimate the parameters using [[maximum likelihood]] method. He shows that the likelihood function can be expressed in an intensive form, in which the scale parameters σ<sub>1</sub>  and σ<sub>2</sub> are a function of the location parameter μ. The likelihood in its intensive form is:
:<math> L(\mu) = -\left[\sum_{x_i: x_i<\mu}  (x_i-\mu)^2 \right]^{1/3} - \left[\sum_{x_i: x_i>\mu}  (x_i-\mu)^2 \right]^{1/3}</math>
and has to be maximized numerically with respect to a single parameter μ only.
 
Given the maximum likelihood estimator <math>\hat{\mu}</math> the other parameters take values:
:<math> \hat{\sigma}_1^2 = \frac{-L(\mu)}{N} \left[\sum_{x_i: x_i<\mu}  (x_i-\mu)^2 \right]^{2/3},</math>
:<math> \hat{\sigma}_2^2 = \frac{-L(\mu)}{N} \left[\sum_{x_i: x_i>\mu}  (x_i-\mu)^2 \right]^{2/3},</math>
where '''N''' is the number of observations.
 
Villani and Larsson<ref name="VL2006" /> propose to use either  [[maximum likelihood]] method or [[Bayes estimator|bayesian estimation]] and provide some analytical results for either univariate and multivariate case.
 
== Applications ==
The split normal distribution has been used mainly in econometrics and time series. A remarkable area of application is the construction of the [[Fan chart (time series)|fan chart]], a representation of the [[inflation]] forecast distribution reported by [[inflation targeting]] central banks around the globe.<ref name="Julio2007"/><ref name="BoEInfRep"/>
 
== References ==
{{Reflist|refs=
<ref name="BoEInfRep">Bank of England, ''[http://www.bankofengland.co.uk/publications/inflationreport/irfanch.htm Inflation Report]''</ref>
 
<ref name=Gibbons>{{cite journal
|last1=Gibbons|first1=J.F.
|last2=Mylroie|first2=S.
|year=1973
|title=Estimation of impurity profiles in ion-implanted amorphous targets using joined half-Gaussian distributions
|journal= Applied Physics Letters
|volume= 22 |pages= 568&ndash;569}}
</ref>
 
<ref name=BD2011>{{Cite conference
| last = Banerjee
| first = N.
| coauthors = A. Das
| title = Fan Chart: Methodology and its Application to Inflation Forecasting in India
| series = Reserve Bank of India Working Paper Series
| year = 2011
}}
</ref>
 
<ref name=BFW1998>{{Cite journal
| volume = February 1998
| pages =  30–37
| last = Britton
| first = E.
| coauthors = P. Fisher, Whitley, J.
| title = The inflation report projections: understanding the fan chart
| journal = Quarterly Bulletin
| year = 1998
}}
</ref>
 
<ref name=VL2006>{{Cite journal
| issn = 0361-0926
| volume = 35
| issue = 6
| pages = 1123–1140
| last = Villani
| first = Mattias
| coauthors = Rolf Larsson
| title = The Multivariate Split Normal Distribution and Asymmetric Principal Components Analysis
| journal = Communications in Statistics - Theory and Methods
| year = 2006
}}
</ref>
 
<ref name="John1982">{{Cite journal
| journal = Communications in Statistics - Theory and Methods
| volume = 11
| issue = 8
| pages = 879–885
| last = John |first=S.
| title = The three-parameter two-piece normal family of distributions and its fitting
| year = 1982
| url = http://www.informaworld.com/10.1080/03610928208828279
}}</ref>
 
<ref name="Johnson1994">{{Cite book
| publisher = John Wiley & Sons
| last =  Johnson, N.L., Kotz, S. and Balakrishnan, N.
| title = Continuous Univariate Distributions, Volume 1
|page=173 |isbn=0-471-58495-9
| year = 1994
}}</ref>
 
<ref name="Julio2007">{{Cite conference
| publisher = Banco de la República
| last = Juan Manuel Julio
| title = The Fan Chart: The Technical Details Of The New Implementation
| accessdate = 2010-09-11
| year = 2007
| url = http://ideas.repec.org/p/col/000094/004294.html
| postscript =, [http://www.banrep.gov.co/docum/ftp/borra468.pdf direct link]
}}</ref>
 
}}
{{ProbDistributions|continuous-infinite}}
 
[[Category:Continuous distributions]]
[[Category:Normal distribution]]
[[Category:Probability distributions]]

Revision as of 04:36, 10 December 2013

In probability theory and statistics, the split normal distribution also known as the two-piece normal distribution results from joining at the mode the corresponding halves of two normal distributions with the same mode but different variances. The split normal distribution was introduced by Gibbons and Mylroie[1] and by John.[2]

Template:Probability distribution

Definition

The split normal distribution arises from merging two opposite halves of two probability density functions (PDFs) of normal distributions in their common mode.

The PDF of the split normal distribution is given by[3]

f(x;μ,σ1,σ2)=Aexp((xμ)22σ12)if x<μ
f(x;μ,σ1,σ2)=Aexp((xμ)22σ22)otherwise

where

A=2/π(σ1+σ2)1.

Discussion

The split normal distribution results from merging two halves of normal distributions. In a general case the 'parent' normal distributions can have different variances which implies that the joined PDF would not be continuous. To ensure that the resulting PDF integrates to 1, the normalizing constant A is used.

In a special case when σ12=σ22=σ*2 the split normal distribution reduces to normal distribution with variance σ*2.

When σ2≠σ1 the constant A it is different from the constant of normal distribution. However, when σ12=σ22=σ*2 the constants are equal.

The sign of its third central moment is determined by the difference (σ21). If this difference is positive, the distribution is skewed to the right and if negative, then it is skewed to the left.

Other properties of the split normal density were discussed by Johnson et al.[3] and Julio.[4]

Alternative formulations

The formulation discussed above originates from John.[2] The literature offers two mathematically equivalent alternative parameterizations . Britton, Fisher and Whitley[5] offer a parameterization if terms of mode, dispersion and normed skewness, denoted with 𝒮𝒩(μ,σ2,γ). The parameter μ is the mode and has equivalent to the mode in John’s formulation. The parameter σ 2>0 informs about the dispersion (scale) and should not be confused with variance. The third parameter, γ ∈ (-1,1), is the normalized skew.

The second alternative parameterization is used in the Bank of England’s communication and is written in terms of mode, dispersion and unnormed skewness and is denoted with 𝒮𝒩(μ,σ2,ξ). In this formulation the parameter μ is the mode and is identical as in John’s [2] and Britton, Fisher and Whitley’s [5] formulation. The parameter σ 2 informs about the dispersion (scale) and is the same as in the Britton, Fisher and Whitley’s formulation. The parameter ξ equals the difference between the distribution’s mean and mode and can be viewed as unnormed measure of skewness.

The three parameterizations are mathematically equivalent, meaning that there is a strict relationship between the parameters and that it is possible to go from one parameterization to another. The following relationships hold:[6]

σ2=σ12(1+γ)=σ22(1γ)γ=σ2σ1σ2+σ1ξ=2/π(σ2σ1)γ=sgn(ξ)1(1+2β1β)2,whereβ=πξ22σ2.

Multivariate Extensions

The multivariate generalization of the split normal distribution was proposed by Villani and Larsson.[7] They assume that each of the principal components has univariate split normal distribution with a different set of parameters μ, σ2 and σ1.

Estimation of parameters

John[2] proposes to estimate the parameters using maximum likelihood method. He shows that the likelihood function can be expressed in an intensive form, in which the scale parameters σ1 and σ2 are a function of the location parameter μ. The likelihood in its intensive form is:

L(μ)=[xi:xi<μ(xiμ)2]1/3[xi:xi>μ(xiμ)2]1/3

and has to be maximized numerically with respect to a single parameter μ only.

Given the maximum likelihood estimator μ^ the other parameters take values:

σ^12=L(μ)N[xi:xi<μ(xiμ)2]2/3,
σ^22=L(μ)N[xi:xi>μ(xiμ)2]2/3,

where N is the number of observations.

Villani and Larsson[7] propose to use either maximum likelihood method or bayesian estimation and provide some analytical results for either univariate and multivariate case.

Applications

The split normal distribution has been used mainly in econometrics and time series. A remarkable area of application is the construction of the fan chart, a representation of the inflation forecast distribution reported by inflation targeting central banks around the globe.[4][8]

References

43 year old Petroleum Engineer Harry from Deep River, usually spends time with hobbies and interests like renting movies, property developers in singapore new condominium and vehicle racing. Constantly enjoys going to destinations like Camino Real de Tierra Adentro. 55 yrs old Metal Polisher Records from Gypsumville, has interests which include owning an antique car, summoners war hack and spelunkering. Gets immense motivation from life by going to places such as Villa Adriana (Tivoli).

my web site - summoners war hack no survey ios

  1. Cite error: Invalid <ref> tag; no text was provided for refs named Gibbons
  2. 2.0 2.1 2.2 2.3 Cite error: Invalid <ref> tag; no text was provided for refs named John1982
  3. 3.0 3.1 Cite error: Invalid <ref> tag; no text was provided for refs named Johnson1994
  4. 4.0 4.1 Cite error: Invalid <ref> tag; no text was provided for refs named Julio2007
  5. 5.0 5.1 Cite error: Invalid <ref> tag; no text was provided for refs named BFW1998
  6. Cite error: Invalid <ref> tag; no text was provided for refs named BD2011
  7. 7.0 7.1 Cite error: Invalid <ref> tag; no text was provided for refs named VL2006
  8. Cite error: Invalid <ref> tag; no text was provided for refs named BoEInfRep