Understanding Conditional Probability in Standard Normal Distributions

Click For Summary
SUMMARY

This discussion focuses on calculating the conditional probability P(Z>0|Y>0) where Z is defined as Z=√qY+√(1-q)X-Φ⁻¹(p) with X and Y being standard normal random variables. The user attempted two numerical approaches in Mathematica, both yielding results that contradict their expectations, with P(Z>0|Y>0) not exceeding 0.5. The user also referenced Bayes' Theorem and the joint normal bivariate distribution, indicating confusion about the correct application of these concepts in their calculations.

PREREQUISITES
  • Understanding of conditional probability and Bayes' Theorem
  • Familiarity with standard normal distributions and random variables
  • Proficiency in using Mathematica for numerical calculations
  • Knowledge of joint normal bivariate distributions
NEXT STEPS
  • Study the application of Bayes' Theorem in continuous probability distributions
  • Learn about joint normal distributions and their properties
  • Explore numerical integration techniques in Mathematica for probability calculations
  • Investigate the implications of covariance in bivariate normal distributions
USEFUL FOR

Statisticians, data scientists, and anyone involved in probability theory or statistical modeling, particularly those working with normal distributions and conditional probabilities.

mariank
Messages
1
Reaction score
0
Dear all,

I am in need of an advice regarding a conditional probability... Although it is very simple, I haven't found anywhere anything similar so I am not sure I get it right...I tried numerically calculating this conditional probability using 2 approaches (in Mathematica) and I get 2 completely different results, which are both very different from my intuition... I would be very happy if you could help me figure it out.

Let X, Y be 2 standard normal r.v. and Z=[tex]\sqrt{q}Y+\sqrt{1-q}X-\Phi^{-1}(p).[/tex]
s.t he unconditional probability of P(Z>0)=1-p

What I would like to compute is P(Z>0|Y>0).

here is what I have tried:

1. P(Z>0|Y>0)=P(Z>0, Y>0)/P(Y>0)=[tex]\int_{0}^\infty \int_{0}^\infty \phi_{zy}(z,y) dz dy/0.5[/tex] where [tex]\phi_{zy}(z,y)[/tex] is the joint normal birvariate distribution of z and y, where the covariance between z and y is[tex]\sqrt{q}.[/tex]
Weird enough, computing this numerically turns out that P(Z>0|Y>0) is at most 0.5(even when I put p=0.01, that means P(Z>0)=0.99) (which is very weird, right?)

2. The second thing I was advised to do is to write P(Z>0|Y>0)=[tex]\int_0^\infty P(Z>0|Y=y) \phi(y) dy[/tex][tex](P(Z>0|Y=y)=P(X>\phi^{-1}(p)-y\sqrt{q}/{\sqrt{1-q}}))[/tex]...these results are even smaller than the previous one...

So what is that I am doing wrong?

I would very much appreciate your help. Thanks a lot,folks!
 
Last edited:
Physics news on Phys.org
mariank said:
Dear all,here is what I have tried:

1. P(Z>0|Y>0)=P(Z>0, Y>0)/P(Y>0)

So what is that I am doing wrong?

I don't understand this expression. The correct expression for Bayes Theorem is:

P(Z>0|Y>0)=P(Y>0|P(Z>0)P(Z>0)/P(Y>0)
 
Last edited:
For the continuous case, I think you want:

[tex]f_{Z,Y}(z,y)=f_{Z|Y}(z|y)f_{Y}(y); Y>0, Z>0[/tex]
 
Last edited:

Similar threads

  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 5 ·
Replies
5
Views
3K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 4 ·
Replies
4
Views
4K
  • · Replies 4 ·
Replies
4
Views
1K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 6 ·
Replies
6
Views
2K