Distribution of the maximum of a multivariate normal random variable

For general $(\mu_X,\Sigma_X)$ The problem is quite difficult, even in 2D. Clearly $P(max(X_i)<x) = P(X_1<x \wedge X_2<x \cdots \wedge X_d <x)$, so to get the distribution function of the maximum one must integrate the joint density over that region... but that's not easy for a general gaussian, even in two dimensions. I'd bet there is no simple expression for the general case.

Some references:

http://www.springerlink.com/content/ca94xg2tdy7evdpb/

http://itc.ktu.lt/itc384/Aksom384.pdf

http://people.emich.edu/aross15/q/papers/bounds_Emax.pdf

Multivariate Skew-Normal Distributions and their Extremal Properties

Rolf Waeber February 8, 2008 Abstract In this thesis it is established that the distribution is a skew normal dist.

A paper by Nadarajah and Samuel Kotz gives the expression for the max of any bivariate normal F(x,y). IEEE TRANSACTIONS ON VERY LARGE SCALE INTEGRATION (VLSI) SYSTEMS, VOL. 16, NO. 2, FEBRUARY 2008 Exact Distribution of the Max/Min of Two Gaussian Random Variables Saralees Nadarajah and Samuel Kotz If F(x,y) is a standard normal (means=0 and variances=1, r>0) the dist of the maximum is a skew normal.