6

I am curious what the derivation for the expectation of the maximum of two jointly normal random variables $X$ and $Y$ with correlation coefficient $\rho$.

I could start with the following but the absolute value sign under expectation doesn't look like a walk in the park:

$\mathbb{E}\left[\text{max}(X,Y)\right] = \mathbb{E}\left[\frac{X+Y}{2}+\frac{|X-Y|}{2}\right] = \ ...$

kjetil b halvorsen
  • 63,378
  • 26
  • 142
  • 467
ambushed
  • 239
  • 1
  • 10
  • 4
    The distribution of $\max(X,Y)$ (for $X$ and $Y$ with equal variances and equal means) is given at http://stats.stackexchange.com/questions/139072/. From that you can compute the expectation (it looks like numerical methods have to be used). The general problem (for arbitrary variances and means) looks difficult: do you need a solution in that case? – whuber Apr 07 '16 at 20:58
  • What i had in mind was the general problem.. – ambushed Apr 07 '16 at 21:05
  • It looks like this question has been [asked](http://stats.stackexchange.com/questions/77110/mathbbe-and-variance-of-the-maximum-of-independent-mathcaln-mu-i-sig) before. Appologies. – ambushed Apr 07 '16 at 21:31
  • I'm not sure of that: I could not find an exact duplicate. The link in my comment was found after conducting three or four keyword searches of this site and inspecting several likely threads; it's the closest I could come. The link you found concerns the maximum of *independent* normal variables and its answer relies fundamentally on that assumption. – whuber Apr 07 '16 at 21:42
  • Indeed, you are right! – ambushed Apr 07 '16 at 21:50

2 Answers2

6

The solution is contained in the paper https://www.gwern.net/docs/conscientiousness/2008-nadarajah.pdf cited by @Lucas in his answer at Distribution of the maximum of two correlated normal variables

I will give the answer here, maybe I come back to add a proof ... Let $(X_1,X_2)$ be a bivariate random vector with a binormal distribution, with means $\mu_1, \mu_2$, standard deviations $\sigma_1, \sigma_2$ and correlation coefficient $\rho$. Then $X=\max(X_1, X_2)$ has probability density function $f(x) = f_1(-x)+f_2(-x)$ where $$ f_1(x)= \frac1{\sigma_1}\phi(\frac{x+\mu_1}{\sigma_1})\cdot \Phi\left( \frac{\rho(x+\mu_1)}{\sigma_1\sqrt{1-\rho^2}}-\frac{x+\mu_2}{\sigma_2\sqrt{1-\rho^2}} \right) \\ f_2(x)= \frac1{\sigma_2}\phi(\frac{x+\mu_2}{\sigma_2})\cdot \Phi\left( \frac{\rho(x+\mu_2)}{\sigma_2\sqrt{1-\rho^2}}-\frac{x+\mu_1}{\sigma_1\sqrt{1-\rho^2}} \right) $$ where $\phi, \Phi$ are the density and cumulative distribution function of the standard normal.

This paper also gives an exact expression for the expectation: $$ \DeclareMathOperator{\E}{\mathbb{E}} \E X = \mu_1 \Phi\left( \frac{\mu_1-\mu_2}{\theta} \right) + \mu_2 \Phi\left( \frac{\mu_2-\mu_1}{\theta} \right) + \theta \phi\left( \frac{\mu_1-\mu_2}{\theta} \right) $$ where $\theta = \sqrt{\sigma_1^2 +\sigma_2^2 - 2\rho\sigma_1\sigma_2}$. (the paper contains more, like the variance and moment generating functions).

kjetil b halvorsen
  • 63,378
  • 26
  • 142
  • 467
0

About the distribution |X-Y|, the Folded normal distribution may be helpful.

Sven Hohenstein
  • 6,285
  • 25
  • 30
  • 39
MyNT
  • 1