Half-normal distribution

Using the ${\displaystyle \sigma }$ parametrization of the normal distribution, the probability density function (PDF) of the half-normal is given by

${\displaystyle f_{Y}(y;\sigma )={\frac {\sqrt {2}}{\sigma {\sqrt {\pi }}}}\exp \left(-{\frac {y^{2}}{2\sigma ^{2}}}\right)\quad y\geq 0,}$

where ${\displaystyle E[Y]=\mu ={\frac {\sigma {\sqrt {2}}}{\sqrt {\pi }}}}$.

Alternatively using a scaled precision (inverse of the variance) parametrization (to avoid issues if ${\displaystyle \sigma }$ is near zero), obtained by setting ${\displaystyle \theta ={\frac {\sqrt {\pi }}{\sigma {\sqrt {2}}}}}$, the probability density function is given by

${\displaystyle f_{Y}(y;\theta )={\frac {2\theta }{\pi }}\exp \left(-{\frac {y^{2}\theta ^{2}}{\pi }}\right)\quad y\geq 0,}$

where ${\displaystyle E[Y]=\mu ={\frac {1}{\theta }}}$.

The cumulative distribution function (CDF) is given by

${\displaystyle F_{Y}(y;\sigma )=\int _{0}^{y}{\frac {1}{\sigma }}{\sqrt {\frac {2}{\pi }}}\,\exp \left(-{\frac {x^{2}}{2\sigma ^{2}}}\right)\,dx}$

Using the change-of-variables ${\displaystyle z=x/({\sqrt {2}}\sigma )}$, the CDF can be written as

${\displaystyle F_{Y}(y;\sigma )={\frac {2}{\sqrt {\pi }}}\,\int _{0}^{y/({\sqrt {2}}\sigma )}\exp \left(-z^{2}\right)dz=\operatorname {erf} \left({\frac {y}{{\sqrt {2}}\sigma }}\right),}$

where erf is the error function, a standard function in many mathematical software packages.

The quantile function (or inverse CDF) is written:

${\displaystyle Q(F;\sigma )=\sigma {\sqrt {2}}\operatorname {erf} ^{-1}(F)}$

where ${\displaystyle 0\leq F\leq 1}$ and ${\displaystyle \operatorname {erf} ^{-1}}$ is the inverse error function

The expectation is then given by

${\displaystyle E[Y]=\sigma {\sqrt {2/\pi }},}$

The variance is given by

${\displaystyle \operatorname {var} (Y)=\sigma ^{2}\left(1-{\frac {2}{\pi }}\right).}$

Since this is proportional to the variance σ² of X, σ can be seen as a scale parameter of the new distribution.

The differential entropy of the half-normal distribution is exactly one bit less the differential entropy of a zero-mean normal distribution with the same second moment about 0. This can be understood intuitively since the magnitude operator reduces information by one bit (if the probability distribution at its input is even). Alternatively, since a half-normal distribution is always positive, the one bit it would take to record whether a standard normal random variable were positive (say, a 1) or negative (say, a 0) is no longer necessary. Thus,

${\displaystyle h(Y)={\frac {1}{2}}\log _{2}\left({\frac {\pi e\sigma ^{2}}{2}}\right)={\frac {1}{2}}\log _{2}\left(2\pi e\sigma ^{2}\right)-1.}$

The half-normal distribution is commonly utilized as a prior probability distribution for variance parameters in Bayesian inference applications.[1][2]

Given numbers ${\displaystyle \{x_{i}\}_{i=1}^{n}}$ drawn from a half-normal distribution, the unknown parameter ${\displaystyle \sigma }$ of that distribution can be estimated by the method of maximum likelihood, giving

${\displaystyle {\hat {\sigma }}={\sqrt {{\frac {1}{n}}\sum _{i=1}^{n}x_{i}^{2}}}}$

The bias is equal to

${\displaystyle b\equiv \operatorname {E} {\bigg [}\;({\hat {\sigma }}_{\mathrm {mle} }-\sigma )\;{\bigg ]}=-{\frac {\sigma }{4n}}}$

which yields the bias-corrected maximum likelihood estimator

${\displaystyle {\hat {\sigma \,}}_{\text{mle}}^{*}={\hat {\sigma \,}}_{\text{mle}}-{\hat {b\,}}.}$

• The distribution is a special case of the folded normal distribution with μ = 0.
• It also coincides with a zero-mean normal distribution truncated from below at zero (see truncated normal distribution)
• If Y has a half-normal distribution, then (Y/σ)² has a chi square distribution with 1 degree of freedom, i.e. Y/σ has a chi distribution with 1 degree of freedom.
• The half-normal distribution is a special case of the generalized gamma distribution with d = 1, p = 2, a = ${\displaystyle {\sqrt {2}}\sigma }$.
• If Y has a half-normal distribution, Y ^-2 has a Levy distribution
• The Rayleigh distribution is a moment-tilted and scaled generalization of the half-normal distribution.

The Modified Half Normal Distribution

Notation	${\displaystyle {\text{MHN}}\left(\alpha ,\beta ,\gamma \right)}$
Parameters	${\displaystyle \alpha >0,\beta >0{\text{ and }}\gamma \in \mathbb {R} }$
Support	${\displaystyle x>0}$
PDF	${\displaystyle f(x)={\frac {2\beta ^{\frac {\alpha }{2}}x^{\alpha -1}\exp(-\beta x^{2}+\gamma x)}{\Psi {\left({\frac {\alpha }{2}},{\frac {\gamma }{\sqrt {\beta }}}\right)}}}}$, where ${\displaystyle \Psi (\alpha ,z)={}_{1}\Psi _{1}\left({\begin{matrix}\left(\alpha ,{\frac {1}{2}}\right)\\(1,0)\end{matrix}};z\right)}$ denotes the Fox-Wright Psi function.
CDF	${\displaystyle F_{_{\text{MHN}}}(x\mid \alpha ,\beta ,\gamma )={\frac {2\beta ^{\frac {\alpha }{2}}}{\Psi {\left({\frac {\alpha }{2}},{\frac {\gamma }{\sqrt {\beta }}}\right)}}}\sum _{i=0}^{\infty }{\frac {\gamma ^{i}}{2i!}}\beta ^{-{\frac {\alpha +i}{2}}}\gamma ({\frac {\alpha +i}{2}},\beta x^{2})}$, where ${\displaystyle \gamma (s,y)}$ denotes the lower incomplete gamma function.
Mean	${\displaystyle E(X)={\frac {\Psi \left({\frac {\alpha +1}{2}},{\frac {\gamma }{\sqrt {\beta }}}\right)}{\beta ^{\frac {1}{2}}\Psi \left({\frac {\alpha }{2}},{\frac {\gamma }{\sqrt {\beta }}}\right)}}}$ ${\displaystyle X_{\text{mode}}\leq E(X)\leq {\frac {\gamma +{\sqrt {\gamma ^{2}+8\alpha \beta }}}{4\beta }}{\text{ if }}\alpha >1.}$
Mode	${\displaystyle X_{\text{mode}}={\frac {\gamma +{\sqrt {\gamma ^{2}+8\beta (\alpha -1)}}}{4\beta }}{\text{ if }}\alpha >1}$.
Variance	${\displaystyle {\text{Var}}(X)={\frac {\Psi \left({\frac {\alpha +2}{2}},{\frac {\gamma }{\sqrt {\beta }}}\right)}{\beta \Psi \left({\frac {\alpha }{2}},{\frac {\gamma }{\sqrt {\beta }}}\right)}}-\left[{\frac {\Psi \left({\frac {\alpha +1}{2}},{\frac {\gamma }{\sqrt {\beta }}}\right)}{\beta ^{\frac {1}{2}}\Psi \left({\frac {\alpha }{2}},{\frac {\gamma }{\sqrt {\beta }}}\right)}}\right]^{2}}$.

The modified half-normal distribution (MHN)[3] is a three-parameter family of continuous probability distributions supported on the positive part of the real line. The truncated normal distribution, half-normal distribution, and square-root of the Gamma distribution are special cases of the MHN distribution.

The MHN distribution is used a probability model, additionally it appears in a number of Markov Chain Monte Carlo (MCMC) based Bayesian procedures including the Bayesian modeling of the Directional Data,[4][5] Bayesian Binary regression,[6] Bayesian Graphical model.[7] The MHN distribution occurs in the diverse areas of research [8][9][10] [11] signifying its relevance to the contemporary statistical modeling and associated computation. Additionally, the moments and its other moment based statistics (including variance, skewness) can be represented via the Fox-Wright Psi functions, denoted by ${\displaystyle \Psi (\cdot ,\cdot )}$. There exists a recursive relation between the three consecutive moments of the distribution.

• Half-t distribution
• Truncated normal distribution
• Folded normal distribution
• Rectified Gaussian distribution

• Leone, F. C.; Nelson, L. S.; Nottingham, R. B. (1961), "The folded normal distribution", Technometrics, 3 (4): 543–550, doi:10.2307/1266560, hdl:2027/mdp.39015095248541, JSTOR 1266560

• Half-Normal Distribution at MathWorld

(note that MathWorld uses the parameter ${\displaystyle \theta ={\frac {1}{\sigma }}{\sqrt {\pi /2}}}$