5
$\begingroup$

I have two normal distributions defined by their averages and standard deviations.

Sample 1: Mean=5.28; SD=0.91

Sample 2: Mean=8.45; SD=1.36

You can see how they look like in the next image:

enter image description here

How can I get the probability to obtain an individual from the overlapping area (green)? Is the probability the same as the area?

|cite|improve this question
$\endgroup$
  • 1
    $\begingroup$ What do you mean by the probability of obtaining an individual from the area? $\endgroup$ – Juho Kokkala Jun 18 '14 at 8:37
  • 2
    $\begingroup$ If you sample points from either normal distribution, you get points on the Perikymata-axis rather than on the 2-dimensional area. Furthermore, the green zone is infinitely wide, so all values sampled from either distribution are under the green zone, so in this sense the probability would be 1. $\endgroup$ – Juho Kokkala Jun 18 '14 at 8:47
  • $\begingroup$ @JuhoKokkala OP likely means integral over the green area. $\endgroup$ – Vladislavs Dovgalecs Nov 21 '17 at 18:40
  • $\begingroup$ @VladislavsDovgalecs What would the question "Is the probability the same as the area" then mean? $\endgroup$ – Juho Kokkala Nov 21 '17 at 19:21
  • 1
    $\begingroup$ @VladislavsDovgalecs I think you misunderstood my comment, I was not asking for an interpretation of the area but trying to parse the question. In any case, please note OP's comments to the accepted answer and a followup question posted here stats.stackexchange.com/questions/103821. $\endgroup$ – Juho Kokkala Nov 23 '17 at 6:14
21
$\begingroup$

It is not quite clear what you mean by probability to obtain an individual from the overlapping area. This solves for the area of the green zone in your diagram:

Let:

  • $X_1 \sim N(\mu_1,\sigma_1^2)$ with pdf $f_1(x_1)$ and cdf $ F_1(x_1)$ and

  • $X_2 \sim N(\mu_2,\sigma_2^2)$ with pdf $f_2(x_2)$ and cdf $ F_2(x_2)$,

where $\mu_1 < \mu_2$. In your example, the 'black variable' corresponds to $X_1$.

Let $c$ denote the point of intersection where the pdf's meet in the green zone of your plot Then, the area of your green intersection zone is simply:

$$P(X_1>c) + P(X_2<c) = 1 - F_1(c) + F_2(c) = 1-\frac{1}{2} \text{erf}\left(\frac{c-\mu _1}{\sqrt{2} \sigma _1}\right)+\frac{1}{2} \text{erf}\left(\frac{c-\mu _2}{\sqrt{2} \sigma _2}\right)$$

where erf(.) is the error function.

Point $c$ is the solution to $f_1(x) = f_2(x)$ within the green zone, which yields:

$$c = \frac{\mu _2 \sigma _1^2-\sigma _2 \left(\mu _1 \sigma _2+\sigma _1 \sqrt{\left(\mu _1-\mu _2\right){}^2+2 \left(\sigma _1^2-\sigma _2^2\right) \log \left(\frac{\sigma _1}{\sigma _2}\right)}\right)}{\sigma _1^2-\sigma _2^2}$$

For your example, with $ {\mu_1 = 5.28, \mu_2 = 8.45, \sigma_1 = 0.91, \sigma_2 = 1.36}$, this yields: $c = 6.70458...$,

and the area of the green section is: 0.158413 ...

|cite|improve this answer
$\endgroup$
  • 1
    $\begingroup$ @antecessor Depends on what you mean by 'overlapping among both'. In the question you talk about some probability, but it is unclear what probabilistic interpretation the area computed here would have. $\endgroup$ – Juho Kokkala Jun 18 '14 at 8:57
  • 1
    $\begingroup$ @JuhoKokkala Ok, I will write a new question to address the exact question. $\endgroup$ – antecessor Jun 18 '14 at 10:38
  • 1
    $\begingroup$ the resolution of the equation pdf1=pdf2 yields to two solutions, so we have two intersection points, why you delete the other point @wolfies $\endgroup$ – user61828 Dec 1 '14 at 11:51
  • 1
    $\begingroup$ Because, within the green zone, there is only one solution. $\endgroup$ – wolfies Dec 1 '14 at 12:32
  • 1
    $\begingroup$ @Dev These equations are fully explained in this post. Your request for a paper is a little bit like reading an answer that mentions 2+3+5=10 and asking for a reference to that fact. If you must have a reference, then, why not reference this very post? $\endgroup$ – whuber♦ Nov 14 '16 at 19:46
1
$\begingroup$

@abdelbasset, To improve @wolfies 's answer above, there are two intersection points c, let's call them c1 and c2. Here, c1=-1.2848 and c2=6.7046. @wolfies ignored c1 (it likely will be too far in the tails of the probability density functions to matter (especially if we round to just a few decimal points)).

The more correct way is to find both c1 and c2, and to find the area of overlap of both functions:

AREA OF OVERLAP = P(X1 > j1) + P(X2 < j1) - [P(X2 < j2)-P(X1 < j2)]
therefore, using cdf's
AREA OF OVERLAP = 1 - F1(j1,μ2,σ2)+F2(j1,μ2,σ2)-F2(j2,μ2,σ2)+F1(j2,μ1,σ1)

|cite|improve this answer
$\endgroup$
  • 1
    $\begingroup$ Welcome to the site, @luckapani. You cannot use the "Your Answer" field to comment or respond to comments. As a result, this would be deleted. However, there is the core of an answer here. Can you edit your post to make it more of an answer & not a comment? Since you're new here, you may want to take our tour, which contains information for new users. $\endgroup$ – gung♦ May 19 '15 at 16:12

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.