View source for Rayleigh distribution

{{Distinguish|Rayleigh mixture distribution}}

{{Probability distribution|
  name       =Rayleigh|
  type       =density|
  pdf_image =[[Image:Rayleigh distributionPDF.svg|325px|Plot of the Rayleigh PDF]]<br /><small></small>|
  cdf_image  =[[Image:Rayleigh distributionCDF.svg|325px|Plot of the Rayleigh CDF]]<br /><small></small>|
  parameters =scale: <math>\sigma>0</math>|
  support    =<math>x\in [0,\infty)</math>|
  pdf        =<math>\frac{x}{\sigma^2} e^{-x^2/\left(2\sigma^2\right)}</math>|
  cdf        =<math>1 - e^{-x^2/\left(2\sigma^2\right)}</math>|
  quantile   =<math>Q(F;\sigma)=\sigma \sqrt{-2\ln(1 - F)}</math>|
  mean       =<math>\sigma \sqrt{\frac{\pi}{2}}</math>|
  median     =<math>\sigma\sqrt{2\ln(2)}</math>|
  mode       =<math>\sigma</math>|
  variance   =<math>\frac{4 - \pi}{2} \sigma^2</math>|
  skewness   =<math>\frac{2\sqrt{\pi}(\pi - 3)}{(4-\pi)^{3/2}}</math>|
  kurtosis   =<math>-\frac{6\pi^2 - 24\pi +16}{(4-\pi)^2}</math>|
  entropy    =<math>1+\ln\left(\frac{\sigma}{\sqrt{2}}\right)+\frac{\gamma}{2}</math>|
  mgf        =<math>1+\sigma te^{\sigma^2t^2/2}\sqrt{\frac{\pi}{2}} \left(\operatorname{erf}\left(\frac{\sigma t}{\sqrt{2}}\right) + 1\right)</math>|
  char       =<math>1 - \sigma te^{-\sigma^2t^2/2}\sqrt{\frac{\pi}{2}} \left(\operatorname{erfi} \left(\frac{\sigma t}{\sqrt{2}}\right) - i\right)</math>|
}}

In [[probability theory]] and [[statistics]], the '''Rayleigh distribution''' is a [[continuous probability distribution]] for nonnegative-valued [[random variable]]s. It is essentially a [[chi distribution]] with two [[degrees of freedom]].

A Rayleigh distribution is often observed when the overall magnitude of a vector is related to its directional [[Euclidean vector#Vector components|components]]. One example where the Rayleigh distribution naturally arises is when [[wind]] velocity is analyzed in [[Plane (geometry)|two dimensions]].
Assuming that each component is [[uncorrelated]], [[Normal distribution|normally distributed]] with equal [[variance]], and zero [[mean]], then the overall wind speed ([[Euclidean vector|vector]] magnitude) will be characterized by a Rayleigh distribution.   
A second example of the distribution arises in the case of random complex numbers whose real and imaginary components are independently and identically distributed [[normal distribution|Gaussian]] with equal variance and zero mean. In that case, the absolute value of the complex number is Rayleigh-distributed.

The distribution is named after [[John Strutt, 3rd Baron Rayleigh|Lord Rayleigh]] ({{IPAc-en|ˈ|r|eɪ|l|i}}).<ref>"The Wave Theory of Light", ''Encyclopedic Britannica'' 1888; "The Problem of the Random Walk", ''Nature'' 1905 vol.72 p.318</ref>

==Definition==

The [[probability density function]] of the Rayleigh  distribution is<ref Name=PP>Papoulis, Athanasios; Pillai, S. (2001) ''Probability, Random Variables and Stochastic Processes''. {{isbn|0073660116}}, {{isbn|9780073660110}} {{Page needed|date=April 2013}}</ref>

:<math>f(x;\sigma) = \frac{x}{\sigma^2} e^{-x^2/(2\sigma^2)}, \quad x \geq 0,</math>

where <math>\sigma</math> is the [[scale parameter]] of the distribution. The [[cumulative distribution function]] is<ref Name=PP/>

:<math>F(x;\sigma) = 1 - e^{-x^2/(2\sigma^2)}</math>

for <math>x \in [0,\infty).</math>

==Relation to random vector length==

Consider the two-dimensional vector <math> Y = (U,V) </math> which has components that are normally distributed, centered at zero, and independent.  Then <math>U</math> and <math>V</math> have density functions

:<math>f_U(x; \sigma) = f_V(x;\sigma) = \frac{e^{-x^2/(2\sigma^2)}}{\sqrt{2\pi\sigma^2}}.</math>

Let <math>X</math> be the length of <math>Y</math>. That is, <math>X = \sqrt{U^2 + V^2}.</math>  Then <math>X</math> has cumulative distribution function

:<math>F_X(x; \sigma) = \iint_{D_x} f_U(u;\sigma) f_V(v;\sigma) \,dA,</math>

where <math>D_x</math> is the disk

:<math>D_x = \left\{(u,v) : \sqrt{u^2 + v^2} < x\right\}.</math>

Writing the [[multiple integral|double integral]] in [[polar coordinate system|polar coordinates]], it becomes

:<math>F_X(x; \sigma) = \frac{1}{2\pi\sigma^2} \int_0^{2\pi} \int_0^x r e^{-r^2/(2\sigma^2)} \,dr\,d\theta = \frac 1 {\sigma^2} \int_0^x r e^{-r^2/(2\sigma^2)} \,dr.
</math>

Finally, the probability density function for <math>X</math> is the derivative of its cumulative distribution function, which by the [[fundamental theorem of calculus]] is

:<math>f_X(x;\sigma) = \frac d {dx} F_X(x;\sigma) = \frac x {\sigma^2} e^{-x^2/(2\sigma^2)},</math>

which is the Rayleigh distribution. It is straightforward to generalize to vectors of dimension other than 2. 
There are also generalizations when the components have [[Hoyt distribution|unequal variance]] or correlations, or when the vector ''Y'' follows a [[multivariate t-distribution|bivariate Student ''t''-distribution]].<ref>{{cite journal|last=Röver|first=C.|title=Student-t based filter for robust signal detection|journal=Physical Review D|volume=84|issue=12|year=2011|pages=122004|doi=10.1103/physrevd.84.122004|arxiv=1109.0442|bibcode=2011PhRvD..84l2004R}}</ref>

==Properties==
The [[moment (mathematics)|raw  moments]] are given by:
: <math>\mu_j = \sigma^j2^{j/2}\,\Gamma\left(1 + \frac j 2\right),</math> 
where <math>\Gamma(z)</math> is the [[gamma function]].

The [[mean]] of a  Rayleigh  random variable  is thus <!--(<math>k=1, \Gamma\left(\tfrac32\right) = \tfrac12 \sqrt{\pi}\,</math>)-->:
:<math>\mu(X) = \sigma \sqrt{\frac{\pi}{2}}\ \approx 1.253\ \sigma.</math>

The [[variance]]   of a  Rayleigh  random variable  is :

:<math>\operatorname{var}(X) = \mu_2-\mu_1^2 =  \left(2-\frac{\pi}{2}\right) \sigma^2 \approx 0.429\ \sigma^2</math>

The [[mode (statistics)|mode]] is <math>\sigma,</math> and the maximum pdf is

:<math> f_{\max} = f(\sigma;\sigma) = \frac{1}{\sigma} e^{-1/2} \approx \frac{0.606}{\sigma}.</math>

The [[skewness]] is given by:

:<math>\gamma_1 = \frac{2\sqrt{\pi}(\pi - 3)}{(4 - \pi)^{3/2}} \approx 0.631</math>

The excess [[kurtosis]] is given by:

:<math>\gamma_2 = -\frac{6\pi^2 - 24\pi + 16}{(4 - \pi)^2} \approx 0.245</math>

The [[characteristic function (probability theory)|characteristic function]] is given by:

:<math>\varphi(t) = 1 - \sigma te^{-\frac{1}{2}\sigma^2t^2}\sqrt{\frac{\pi}{2}} \left[\operatorname{erfi}\left(\frac{\sigma t}{\sqrt{2}}\right) - i\right]</math>

where <math>\operatorname{erfi}(z)</math> is the imaginary [[error function]]. The [[moment generating function]] is given by

:<math>
  M(t) = 1 + \sigma t\,e^{\frac{1}{2}\sigma^2t^2}\sqrt{\frac{\pi}{2}}
           \left[\operatorname{erf}\left(\frac{\sigma t}{\sqrt{2}}\right) + 1\right]</math>

where <math>\operatorname{erf}(z)</math> is the [[error function]].

===Differential entropy===

The [[differential entropy]] is given by{{Citation needed|date=April 2013}}

:<math>H = 1 + \ln\left(\frac \sigma {\sqrt{2}}\right) + \frac \gamma 2 </math>

where <math>\gamma</math> is the [[Euler–Mascheroni constant]].

== Parameter estimation ==

Given a sample of ''N'' [[independent and identically distributed]] Rayleigh random variables <math>x_i</math> with parameter <math>\sigma</math>,

: <math>\widehat{\sigma}^2\approx \!\,\frac{1}{2N}\sum_{i=1}^N x_i^2</math> is the [[maximum likelihood estimation|maximum likelihood]] estimate and also is [[bias of an estimator|unbiased]].

:<math>\widehat{\sigma}\approx \sqrt{\frac 1 {2N} \sum_{i=1}^N x_i^2}</math> is a biased estimator that can be corrected via the formula

:<math>\sigma = \widehat{\sigma} \frac {\Gamma(N)\sqrt{N}} {\Gamma(N + \frac 1 2)} = \widehat{\sigma} \frac {4^N N!(N-1)!\sqrt{N}} {(2N)!\sqrt{\pi}}</math><ref>[https://archive.org/details/jresv68Dn9p1005 Siddiqui, M. M. (1964) "Statistical inference for Rayleigh distributions", ''The Journal of Research of the National Bureau of Standards, Sec. D: Radio Science'', Vol. 68D, No. 9, p. 1007]</ref>

=== Confidence intervals ===
To find the (1&nbsp;&minus;&nbsp;''α'') confidence interval, first find the bounds <math>[a,b]</math> where:
:&nbsp; <math>P(\chi_{2N}^2 \leq a) = \alpha/2, \quad P(\chi_{2N}^2 \leq b) = 1 - \alpha/2</math>
then the scale parameter will fall within the bounds
:&nbsp; <math>\frac{{N}\overline{x^2}}{b} \leq {\widehat\sigma}^2 \leq \frac{{N}\overline{x^2}}{a}</math><ref>[http://nvlpubs.nist.gov/nistpubs/jres/66D/jresv66Dn2p167_A1b.pdf Siddiqui, M. M. (1961) "Some Problems Connected With Rayleigh Distributions", ''The Journal of Research of the National Bureau of Standards; Sec. D: Radio Propagation'', Vol. 66D, No. 2, p. 169]</ref>

== Generating random variates ==

Given a random variate ''U'' drawn from the [[uniform distribution (continuous)|uniform distribution]] in the interval <nowiki>(0,&nbsp;1)</nowiki>, then the variate

:<math>X=\sigma\sqrt{-2 \ln U}\,</math>

has a Rayleigh distribution with parameter <math>\sigma</math>. This is obtained by applying the [[inverse transform sampling]]-method.

==Related distributions==

*<math>R \sim \mathrm{Rayleigh}(\sigma)</math> is Rayleigh distributed if <math>R = \sqrt{X^2 + Y^2}</math>, where <math>X \sim N(0, \sigma^2)</math> and <math>Y \sim N(0, \sigma^2)</math> are independent [[Normal distribution|normal random variables]].<ref>[https://web.archive.org/web/20131105232146/http://home.kpn.nl/jhhogema1966/skeetn/ballist/sgs/sgs.htm#_Toc96439743 Hogema, Jeroen (2005) "Shot group statistics"]</ref> (This gives motivation to the use of the symbol "sigma" in the above parametrization of the Rayleigh density.)
*The [[chi distribution]] with ''v''&nbsp;=&nbsp;2 is equivalent to the Rayleigh Distribution with&nbsp;''&sigma;''&nbsp;=&nbsp;1.  I.e., if <math>R \sim \mathrm{Rayleigh} (1)</math>, then <math>R^2</math> has a [[chi-squared distribution]] with parameter <math>N</math>, degrees of freedom, equal to two (''N''&nbsp;=&nbsp;2)
::  <math>[Q=R^2] \sim \chi^2(N)\ .</math>

*If <math>R \sim \mathrm{Rayleigh}(\sigma)</math>, then <math>\sum_{i=1}^N R_i^2</math> has a [[gamma distribution]] with parameters <math>N</math> and <math>2\sigma^2</math>
:: <math>\left[Y=\sum_{i=1}^N R_i^2\right] \sim \Gamma(N,2\sigma^2) .</math>

* The [[Rice distribution]] is a generalization of the Rayleigh distribution: <math> \mathrm{Rayleigh}(\sigma) = \mathrm{Rice}(0,\sigma) </math>.
* The [[Weibull distribution]] with the "shape parameter" ''k''=2 yields a Rayleigh distribution.  Then the Rayleigh distribution parameter <math>\sigma</math> is related to the Weibull scale parameter according to <math>\lambda = \sigma \sqrt{2} .</math>
* The [[Maxwell–Boltzmann distribution]] describes the magnitude of a normal vector in three dimensions.
* If <math>X</math> has an [[exponential distribution]] <math>X \sim \mathrm{Exponential}(\lambda)</math>, then <math>Y=\sqrt{X} \sim \mathrm{Rayleigh}(1/\sqrt{2\lambda}) .</math>

== Applications ==

An application of the estimation of σ can be found in [[magnetic resonance imaging]] (MRI). As MRI images are recorded as [[complex numbers|complex]] images but most often viewed as magnitude images, the background data is Rayleigh distributed. Hence, the above formula can be used to estimate the noise variance in an MRI image from background data.<ref>{{cite journal | last1 = Sijbers | first1 = J. | last2 = den Dekker | first2 = A. J. | last3 = Raman | first3 = E. | last4 = Van Dyck | first4 = D. | year = 1999 | title = Parameter estimation from magnitude MR images | url = | journal = International Journal of Imaging Systems and Technology | volume = 10 | issue = 2| pages = 109–114 | doi=10.1002/(sici)1098-1098(1999)10:2<109::aid-ima2>3.0.co;2-r| citeseerx = 10.1.1.18.1228 }}</ref>
<ref>{{cite journal | last1 = den Dekker | first1 = A. J. | last2 = Sijbers | first2 = J. | year = 2014 | title = Data distributions in magnetic resonance images: a review | journal = Physica Medica | volume = 30 | issue = 7| pages = 725–741 | doi=10.1016/j.ejmp.2014.05.002| pmid = 25059432 }}</ref> 

The Rayleigh distribution was also employed in the field of [[Nutrition]] for linking [[Diet (nutrition)|dietary]] [[nutrient]] levels and [[human]] and [[Animal husbandry|animal]] responses. In this way, the [[parameter]] σ may be used to calculate nutrient response relationship.<ref>{{Cite journal|last=Ahmadi|first=Hamed|date=2017-11-21|title=A mathematical function for the description of nutrient-response curve|journal=PLOS ONE|volume=12|issue=11|pages=e0187292|doi=10.1371/journal.pone.0187292|pmid=29161271|issn=1932-6203|bibcode=2017PLoSO..1287292A|pmc=5697816}}</ref>

==See also==
*[[Rayleigh fading]]
*[[Rayleigh mixture distribution]]
*[[Circular error probable]]

{{More footnotes|date=April 2013}}

== References ==

{{reflist}}

{{ProbDistributions|continuous-semi-infinite}}

{{DEFAULTSORT:Rayleigh Distribution}}
[[Category:Continuous distributions]]
[[Category:Exponential family distributions]]