# Thread: Maximum of a set of independent normal random variables

1. ## Maximum of a set of independent normal random variables

Consider a set of random variables with unit varience but different means: . What is the probability that will have the maximum value of the set.

My initial guess was:

But if one considers a simple case where there are 3 random variables with mean 0 then obviously the probability of one of them being the maximum is 1/3.
But using the formula above will give a probablility of 0.5*0.5=0.25.
when

What am I doing wrong here? In the simple 3 variable example I suspect that the probabilities aren't independent or something, but I can't wrap my head around it.
Anyone have any thoughts?

2. Here is the thing I be thinkin

pr( X1 > X2 AND X1 > X3 ) = pr( X1 > X2 > X3 ONION X1 > X2 > X3)

Is the sum of the two disjoint events;

This accounts for two of the 3! equally likely outcomes of the eperiment

1/3

3. Originally Posted by fed1
pr( X1 > X2 AND X1 > X3 ) = pr( X1 > X2 > X3 ONION X1 > X2 > X3)
Onion? What's that mean?

4. Is it math operator or vegetable ?

5. LOL. Its a vegetable. Mightymergz, I'm a little confused about how you set up the problem though. If the distribution of each X has a different mean, wouldn't the probability of selecting the same number in each X be different?

Going back to your example, if you have 3 random variables, and lets say that their means are 0,1,2, then the probability of one of them being the max would not be 1/3.

6. They were working on the simplified example where each has a mean of 0 to test his result against intuition.

7. Indeed, it is a simplified version of the main problem. Also I've never heard of the ONION operator.. Onions are delicious to be sure, but they are less useful for calculating probabilities?

In the simplified example all events are equally proabable so there are only possible outcomes.
,
,
,
,
,

is two out of the 6 possibilities so the probability is 1/3. Unfortunately this method doesn't work if the random variables have different means. Is there another method or can I modify the above method to work when all choices aren't equally probable..

thanks for your help so far!

8. Can I ask what this is for? Just interested? Is it homework? Required for research?

Do we know anything else about the variables? What distribution they have maybe?

9. Hey Dason, ya, this problem is related to some research I am involved with for the summer. I am an undergrad co-op student on a summer research position

The problem I posted is just one way to see the problem, I suppose. The main problem is: consider a segment of a discrete signal from t=-T to T. The signal f(t) = A(t) + y(t) where A(t) is some known function and y(t) is a normal zero mean unit varience random variable.

The way I posed the problem is to just take A(t) as the mean of the random variables at each point and then calculate the statistics from there. We want to know with what probability will each point in the time series be the maximum.

Some things I would expect would be that the point with the highest probability of being the maximum would be the point with the highest mean/max(A(t)) (duh)
and normalization: ie the peak is somewhere:

But the initial method I tried in the OP definitely doesn't satisfy the normalization requirement. I suppose I could try just adding a normalization constant, but that doesn't seem right at all..

10. Ok...here's my reasoning...though I'd love for someone to chime in:

So, if the mu's were 0,1,2 and the sigma was 1, then we'd have

11. Hi Link, Thanks for your thoughts. I think I follow your train of thought but I numerically calculated the probabilities for 's equal to 0,1,2 and the probability for to be the maximum is closer to 0.73ish

12. really? And that's with a sigma of 1? I'm having a problem seeing where I might have went wrong. Though I relied on an online calculator to get me the p-values, I don't see them noticeably wrong.

13. I think I've got an answer. It's not necessarily pretty but I think I can do some work to simplify it. It matches intuition on that if all have the same mean and if there are 3 r.v.s then the probability is 1/3, if there are n then it's 1/n. But like I said when the means aren't the same it's not pretty so I'm going to see if it simplifies at all.

14. Originally Posted by mightymergz
Hi Link, Thanks for your thoughts. I think I follow your train of thought but I numerically calculated the probabilities for 's equal to 0,1,2 and the probability for to be the maximum is closer to 0.73ish
I did some simulations and I the probability around .728 so I think something went wrong Link. I think that P(X1 > X2) and P(X1 > X3) aren't necessarily independent because if you know that X1 > X2 then it's probably more likely that X1 > X3. But that's just a thought.

ya with var/stddev of 1. Someone suggested to me that aren't independent events. I.e. the probability of the second event depends on the result of the first one. i.e. the probability should be more like:

though I'm not sure how to calculate that second conditional probability.. In the case of the simplified example I know it should be 2/3 simply because the and 1/2 * 2/3 = 1/3 which is what we expect

@Dason
Ya I think it is something like that