Thursday, 7 August 2014

chapter-15 Concept of Inverse Probability

To understand the underlying concept of inverse probability, lets start with a very simple example.
Suppose that you have two bags \alpha \,\,{\rm{and}}\,\,\beta , one containing 10 red and 1 white balls, and the second containing 10 white and 1 red balls.
You play a game with your friend. The friend tosses a fair coin, without telling you the outcome and if he gets a Head he withdraws a ball from Bag-\alpha  while if he gets a Tail, he withdraws a ball from Bag–\beta , with you looking away all the time. After doing this once, he has a red ball in his hand. Which bag do think is the more likely one from which this ball was drawn? Intuition immediately tells us that it should be Bag–\alpha , since it has a large number of red balls. What we need to do now is quantify the inverse probability of the ball being drawn from Bag–\alpha  and Bag–\beta , the word ‘inverse’ being used since you are trying to find the probability of an event that has already taken place, using information from a subsequent event. We do this intuitively all the time: “India’s scorecard against Pakistan in yesterday’s match had a century.” “Oh, it would most likely have been Tendulkar!” So here, the speaker is expressing his conviction that the century-scorer must have been Tendulkar, since he has in his mind the information about the various players, and he thinks Tendulkar is the best.
Coming back to our bags and balls, let us draw a tree diagram highlighting the various possible actions your friend can take (the brackets show the probabilities of the corresponding paths)
Now comes the crucial part. Note that the total probability of selecting a red ball is the sum of the probabilities of the two darkened paths (one through Bag–\alpha, one through Bag–\beta ). This is
P\left( {{\rm{Red}}\,{\rm{Ball}}} \right) = \dfrac{1}{2} \times \dfrac{{10}}{{11}} + \dfrac{1}{2} \times \left( {\dfrac{1}{{11}}} \right) = \dfrac{1}{2}
Now, the probability of selecting a red ball through Bag –\alpha corresponds to the upper path only, and it equals
P\left( {{\rm{Red}}\,{\rm{Ball\, from \,Bag \,- }}\alpha } \right) = \dfrac{1}{2} \times \dfrac{{10}}{{11}} = \dfrac{5}{{11}}
Similarly,
P\left( {{\rm{Red}}\,{\rm{Ball\, from \,Bag \,- }}\beta } \right) = \dfrac{1}{2} \times \dfrac{1}{{11}} = \dfrac{1}{{22}}
Finally, it should now be intuitively obvious that
P\left\{ \begin{array}{l}  {\rm{Bag\, selected\, was\, Bag\, - }}\alpha {\rm{ }}\\  {\rm{given\, that \,ball\, is\, red}}  \end{array} \right\} = \dfrac{{P\left( {{\rm{Red\, ball\, from \,Bag\, - }}\alpha } \right)}}{{P\left( {{\rm{Red \,ball}}} \right)}}
 = \dfrac{{5/22}}{{1/2}}
 = \dfrac{{10}}{{11}}
while P\left\{ \begin{array}{l}  {\rm{Bag\, selected\, was\, Bag\, - }}\beta {\rm{ }}\\  {\rm{given \,that\, ball\, is\, red}}  \end{array} \right\} = \dfrac{{P\left( {{\rm{Red \,ball\, from \,Bag \,- }}\beta } \right)}}{{P\left( {{\rm{Red\, ball}}} \right)}}
 = \dfrac{{1/22}}{{1/2}}
 = \dfrac{1}{{11}}
Note what a huge difference there is between the two probabilities, which was expected. Also expected is the fact that the two probabilities sum to 1.
This, then, is the essence of calculating inverse probabilities. We are given the information that an event E has occurred. This event E can occur through npaths {\rm{Path}}_1{\rm{Path}}_2\ldots{\rm{Path}}_n. We want to find the probability that Eoccurred through some particular path, say {\rm{Path}}_i, which is
P\left( {E\,{\rm{occurred \,through\, Pat}}{{\rm{h}}_i}} \right) = \dfrac{{P\left( {{\rm{Pat}}{{\rm{h}}_i}} \right)}}{{P\left( {{\rm{Pat}}{{\rm{h}}_1}} \right) + P\left( {{\rm{Pat}}{{\rm{h}}_2}} \right) + \ldots + P\left( {{\rm{Pat}}{{\rm{h}}_n}} \right)}}
Let us write this in standard terminology, which will give us the Baye’s theorem.
Suppose that the sample space consists of n mutually exclusive events {E_1},{E_2},\ldots ,{E_n}. Now, an event A occurs, which could have resulted from any of the events {E_i} (For example, think of A as obtaining a red ball in the previous example, while {E_1} and {E_2} are selecting Bag-\alpha and Bat-\beta respectively). We intend to find P({E_i}/A), i.e., the probability that {E_i} occurred given that A has occurred. There are now two ways to do the visualisation:
The left hand tree we have already explained. The right hand side shows that Ais an event that has occurred, which must have been a result of one of the {E_i}occurring (i.e. one of the {E_i}’s must occur for A to occur)
From the tree, evaluating P({E_i}/A) has already been explained:
P\left( {{E_i}/A} \right) = \dfrac{{P\left( {{\rm{Path}}\,\,{\rm{to}}\,A\,\,{\rm{through}}\,\,{E_i}} \right)}}{{\sum\limits_{j = 1}^n {P\left( {{\rm{Path}}\,\,{\rm{to}}\,A\,\,{\rm{through}}\,\,{E_j}} \right)} }}
 = \dfrac{{P\left( {{E_i}\,{\rm{occurs\, and }}\,{\bf{then}}\,\,A\,\,{\rm{occurs}}} \right)}}{{\sum\limits_{j = 1}^n {P\left( {{E_j}\,\,{\rm{occurs}}\,\,{\rm{and}}\,\,{\rm{then}}\,\,A\,\,{\rm{occurs}}} \right)} }}
 = \dfrac{{P\left( {{E_i}\,{\rm{occurs}}} \right) \times P\left\{ {A\,\,{\rm{occurs\, given\, that }}{E_i}\,{\rm{has \,occured}}} \right\}}}{{\sum\limits_{j = 1}^n {P\left( {{E_j}\,{\rm{occurs}}} \right) \times P\left\{ {A\,\,{\rm{occurs \,given\, that }}{E_j}\,{\rm{has \,occured}}} \right\}} }}
 = \dfrac{{P\left( {{E_i}} \right)P\left( {A/{E_i}} \right)}}{{\sum\limits_{j = 1}^n {P\left( {{E_j}} \right)P\left\{ {A\,/\,{E_j}} \right\}} }}
The same relation follows from the second figure:
P\left( {{E_i}/A} \right) = \dfrac{{P\left( {{\rm{Darkly\, shaded \, region\,}}} \right)}}{{P\left( {{\rm{Total \,shaded \, region}}} \right)}}
 = \dfrac{{P\left( {{\rm{A}} \cap {E_i}} \right)}}{{\sum\limits_{j{\rm{  = }}\,{\rm{1}}}^n {P\left( {A \cap {E_j}} \right)} }}
 = \dfrac{{P\left( {{E_i}} \right)P\left( {A/{E_i}} \right)}}{{\sum\limits_{j{\rm{  = }}\,{\rm{1}}}^n {P\left( {{E_j}} \right)P\left( {A/{E_j}} \right)} }}
Thus, the famous Baye’s theorem is
{P\left( {{E_i}/A} \right) = \dfrac{{P\left( {{E_i}} \right)P\left( {A/{E_i}} \right)}}{{\sum\limits_{j = 1}^n {P\left( {{E_j}} \right)P\left( {A/{E_j}} \right)} }}}
The name “inverse” stems from the fact that this relation gives us P({E_i}/A) in terms of P(A/{E_j}). The theorem is also known as a theorem on the probability of causes. The reason should be obvious.
In the examples following this, we’ll be using the tree diagram to calculate inverse probabilities. With sufficient practice, you’ll eventually not require drawing the tree diagram because by then you’ll be quite comfortable in using Baye’s theorem directly.

No comments:

https://www.youtube.com/TarunGehlot