# Posts / Expected value of the minimum of two standard half-normal distributions

Let $X_1, X_2$ be i.i.d. random variables that follow the standard normal distribution $\mathcal{N}(0,1^2)$. One day I wondered, what is the expected value of $Z = \min(|X_1|, |X_2|)$? It turned out to be a fun exercise. Let’s solve it together!

Let’s denote the absolute values of $X_1, X_2$ as $Y_1, Y_2$:

$$Y_1=|X_1|, \quad Y_2=|X_2|.$$

Thus, $Y_1$, $Y_2$ follow the standard half-normal distribution. The CDF of this distribution is well known:

$$F_Y(y) = \operatorname{erf}\bigg( \frac{y}{\sqrt{2}} \bigg), \quad \textrm{for}\,\,\, y\geq 0,$$

where $\operatorname{erf}$ is the error function. Let $\Phi$ be the CDF of the standard normal distribution:

$$F_X(x) = \Phi(x) = \frac{1}{2} \Bigg( 1 + \operatorname{erf} \bigg( \frac{x}{\sqrt{2}} \bigg) \Bigg).$$

Let also $\phi$ be the PDF of the standard normal distribution:

$$\phi(x) = \Phi'(x) = \frac{1}{\sqrt{2\pi}}e^{-x^2/2}$$

It’s easy to see that

$$F_Y(y) = 2\Phi(y)-1.$$

Suppose $U, V$ are independent random variables with CDFs $F_U$ and $F_V$, $W = \min(U,V)$. Let’s express the CDF $F_W$ of $W$ via $F_U$ and $F_W$:

$$\begin{split} F_W(w) & = \mathbb{P}(W \leq w) = 1 - \mathbb{P}(W > w) = 1 - \mathbb{P}(\min(U, V) > w) = 1 - \mathbb{P}(U > w, V > w) = \\ & = 1 - \mathbb{P}(U > w)\cdot \mathbb{P}(V > w) = 1 - (1-F_U(w))(1-F_V(w)) = \\ & = F_U(w) + F_V(w) - F_U(w)F_V(w). \end{split}$$

Now let’s apply this rule to $F_Z$:

$$F_Z(z) = 2F_Y(z) - F_Y^2(z) = 2(2\Phi(z) - 1) - (2\Phi(z) - 1)^2 = -4\Phi^2(z)+8\Phi(z)-3.$$

Now we can calculate the PDF of $Z$:

$$f_Z(z) = F_Z'(z) = (-4\Phi^2(z)+8\Phi(z)-3)' = -8\phi(z)\Phi(z)+8\phi(z) = 8\phi(z) (1-\Phi(z)).$$

Since $Y = |X|\geq 0$, the expected value of $Z$ is a definite integral from $0$ to $\infty$ of $xf_Z(x)dx$:

$$\mathbb{E}(Z) = \int_0^\infty xf_Z(x)dx = \int_0^\infty 8x \phi(z) (1-\Phi(z))dx$$

First of all, let’s calculate the indefinite integral of this expression:

$$\begin{split} \int 8x \phi(z) (1-\Phi(z))dx & = \int 8x \frac{1}{\sqrt{2\pi}}e^{-x^2/2} (1 - \frac{1}{2}(1+\operatorname{erf}(x/\sqrt{2}))dx =\\ & = \int -2\sqrt{\frac{2}{\pi}} x e^{-x^2/2} (\operatorname{erf}(x/\sqrt{2})-1)dx =\\ & = \int -2\sqrt{\frac{2}{\pi}} \bigg( 2\frac{e^{-x^2}}{\sqrt{2\pi}} + xe^{-x^2/2}(\operatorname{erf}(x/\sqrt{2})-1) - 2\frac{e^{-x^2}}{\sqrt{2\pi}} \bigg)dx =\\ & = \int -2\sqrt{\frac{2}{\pi}} \bigg( 2\frac{e^{-x^2}}{\sqrt{2\pi}} - \Big( -xe^{-x^2/2}(\operatorname{erf}(x/\sqrt{2})-1) + e^{-x^2/2}\frac{1}{\sqrt{2}} \frac{2}{\sqrt{\pi}} e^{-x^2/2} \Big) \bigg)dx =\\ & = \int -2\sqrt{\frac{2}{\pi}} \bigg( \Big(\frac{\operatorname{erf}(x)}{\sqrt{2}} \Big)' - \Big( e^{-x^2/2} (\operatorname{erf}(x/\sqrt{2})-1) \Big)' \bigg)dx =\\ & = \int \bigg( -2\sqrt{\frac{2}{\pi}} \Big( \frac{\operatorname{erf}(x)}{\sqrt{2}} - e^{-x^2/2}(\operatorname{erf}(x/\sqrt{2})-1) \Big) \bigg)'dx =\\ & = -2\sqrt{\frac{2}{\pi}} \Big( \frac{\operatorname{erf}(x)}{\sqrt{2}} - e^{-x^2/2}(\operatorname{erf}(x/\sqrt{2})-1) \Big). \end{split}$$

Now we are ready to finish the calculation of $\mathbb{E}(Z)$:

$$\begin{split} \mathbb{E}(Z) & = \int_0^\infty 8x \phi(z) (1-\Phi(z))dx =\\ & = -2\sqrt{\frac{2}{\pi}} \Big( \frac{\operatorname{erf}(x)}{\sqrt{2}} - e^{-x^2/2}(\operatorname{erf}(x/\sqrt{2})-1) \Big) \Bigg|_0^\infty =\\ & = -2\sqrt{\frac{2}{\pi}} \bigg( \frac{1}{\sqrt{2}} - 0 \bigg) - \Bigg( -2\sqrt{\frac{2}{\pi}} \bigg( 0 -e^0 (-1) \bigg) \Bigg) =\\ & = -\frac{2}{\sqrt{\pi}} + \frac{2}{\sqrt{\pi}} \sqrt{2} = \frac{2(\sqrt{2}-1)}{\sqrt{\pi}} \approx 0.467389954510218. \end{split}$$

Hooray, we have solved the initial problem:

$$\mathbb{E}(\min(X_1, X_2)) = \frac{2(\sqrt{2}-1)}{\sqrt{\pi}} \approx 0.467389954510218.$$

We can check the correctness of this result using the following R script that uses a Monte-Carlo simulation to evaluate the first three digits after the decimal point of the result:

set.seed(42)
round(mean(pmin(abs(rnorm(100000000)), abs(rnorm(100000000)))), 3)

## 0.467