Lets say I have a normal distribution $N(\mu, \sigma^2)$ from which I have drawn $n$ i.i.d. samples $x_1, \dots, x_n$.

Now, lets define a random variable $Y = max(x_1, \dots, x_n)$.

When $n=1$, the expected value of $Y$ is $\mu$. I would expect that as $n$ increases, the expected value of $Y$ should increase as well. Is it possible to determine the expected value of $Y$ for any value of $n$, in terms of $\mu$ and $\sigma$?


If we combine two of the answers here (Approximate order statistics for normal random variables), we have for the $r$th $\it{smallest}$ order statistic

$$E[r,n] \approx \mu + \sigma \ \Phi^{-1} \left( \frac{r-\frac{\pi}{8}}{n-\frac{\pi}{4}+1}\right) $$

For the largest value we want $r=n,$ so we have

$$E[Y] \approx \mu + \sigma \ \Phi^{-1} \left( \frac{n-\frac{\pi}{8}}{n-\frac{\pi}{4}+1}\right) $$


First note that\begin{align}Y_n=\max\{X_1,\ldots,X_n\}&=\max\{\sigma\epsilon_1+\mu,\ldots,\sigma\epsilon_n+\mu\}\\&=\sigma\max\{\epsilon_1,\ldots,\epsilon_n\}+\mu\\&=\sigma\xi_n+\mu\end{align} hence that $(\mu,\sigma)$ is also a location-scale parameter for the maximum. Asymptotically, the Normal distribution belongs to the domain of attraction of the Gumbel distribution, meaning that $$\sqrt{2\log(n)}(\xi_n-d_n)\stackrel{{\cal L}}{\longrightarrow} G_0$$with $G_0(x)=\exp\{-\exp(-x)\}$ the Gumbel pdf and $$d_n = \sqrt{2\log(n)}-\dfrac{\log\log n + \log(4\pi)}{2\sqrt{2\log(n)}}$$



I found this paper referenced in a thread on the maths stack exchange (Approximate order statistics for normal random variables), so I had a look. For the maximum, $r=n$.

"In a sample of size n the expected value of the rth largest order statistic is given by


where $\phi(x)=1/\sqrt(2\pi)exp(-\frac{1}{2}x^2)$ and $\Phi(x)=\int^x_{-\infty}\phi(z)dz.$"

  • Royston, J. P. (1982), 'Algorithm AS 177: Expected Normal Order Statistics (Exact and Approximate)', Journal of the Royal Statistical Society. Series C (Applied Statistics), 31(2):161-165.

So $Y$ is an order statistic. Let's label its density function $g_{(n)}(x)$, to indicate that it's the pdf of the variable in the nth position (i.e. its the pdf of the maximum in the sample). Let's also label the normal $N(\mu, \sigma^2)$ density function as $f(x)$. It's a standard result that $$g_{(n)}(x)=n[F(x)]^{n-1}f(x),$$ where $F(x)$ is the cumulative density function of $N(\mu, \sigma^2)$ (as a reference, I suggest Mathematical Statistics (7th ed.) by Wackerly, Mendenhall, and Scheaffer, p.333).

It is at this point that I'm unable to proceed - I don't know how to evaluate the expected value of $Y$, given that it has such a strange pdf. However, I'd advise you to search for "expected value of order statistic" - in particular, I found a thread on this topic on the maths stack exchange site:

EDIT: As pointed out by Khol, thread is for a uniform distribution, not a normal distribution. The uniform is apparently more straightforward to deal with. Apologies for the partial answer!


  • $\begingroup$ Linked thread is for a uniform random sample, which is distinctly simpler than a normal random sample. $\endgroup$ – khol May 2 '18 at 4:36

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.