(17 intermediate revisions by the same user not shown)
Line 1: Line 1:
 
[[Category:ECE302Spring2013Boutin]] [[Category:ECE]] [[Category:ECE302]] [[Category:probability]] [[Category:problem solving]] [[Category:conditional probability]] [[Category:independence]]
 
[[Category:ECE302Spring2013Boutin]] [[Category:ECE]] [[Category:ECE302]] [[Category:probability]] [[Category:problem solving]] [[Category:conditional probability]] [[Category:independence]]
  
 +
<center>
 
==Pixel = 'Pic'ture 'El'ement==
 
==Pixel = 'Pic'ture 'El'ement==
 
[[Image:Smiley_3_1_6.PNG]]
 
[[Image:Smiley_3_1_6.PNG]]
 +
</center>
  
 
Above is an RGB image with a pixel color distribution of roughly the following:
 
Above is an RGB image with a pixel color distribution of roughly the following:
Line 10: Line 12:
  
 
Therefore, if we were to randomly sample a pixel color from the image ('im1'), we would have the following probabilities of sampling each color:
 
Therefore, if we were to randomly sample a pixel color from the image ('im1'), we would have the following probabilities of sampling each color:
  P(R|im1) = .3
+
  <math style='inline'>P(R|im_1) = .3</math>
  P(G|im1) = .1
+
  <math style='inline'>P(G|im_1) = .1</math>
  P(B|im1) = .6
+
  <math style='inline'>P(B|im_1) = .6</math>
  
 
Given the following probabilities of two more RGB images with the same number of pixels as im1, solve '''(a)''', '''(b)''', and '''(c)'''.  (Note: we assume for this problem that pixels are sampled with replacement.)
 
Given the following probabilities of two more RGB images with the same number of pixels as im1, solve '''(a)''', '''(b)''', and '''(c)'''.  (Note: we assume for this problem that pixels are sampled with replacement.)
  
  P(R|im1) = .3 P(R|im2) = .8 P(R|im3) = .4
+
  <math style='inline'>P(R|im_1) = .3</math>    <math style='inline'>P(R|im_2) = .8</math>    <math style='inline'>P(R|im_3) = .4</math>
  P(G|im1) = .1 P(G|im2) = .1 P(G|im3) = .2     } where R = red pixel sampled, etc, and im1 = pixel sampled is in image 1, etc.
+
  <math style='inline'>P(G|im_1) = .1</math>    <math style='inline'>P(G|im_2) = .1</math>    <math style='inline'>P(G|im_3) = .2</math>    } where R = red pixel sampled, etc, and im1 = pixel sampled is in image 1, etc.
  P(B|im1) = .6 P(B|im2) = .1 P(B|im3) = .4
+
  <math style='inline'>P(B|im_1) = .6</math>    <math style='inline'>P(B|im_2) = .1</math>    <math style='inline'>P(B|im_3) = .4</math>
  
  (a) Assuming that we have an equal probability of sampling a pixel from each image, what is the probability that a sampled red
+
  (a) Assuming that we have an equal probability of sampling a pixel from each image (ie <math style='inline'>P(im_1) = P(im_2) = P(im_3) = 1/3</math>),
    pixel came from image 3?  
+
    what is the probability that a sampled red pixel came from image 3?  
  
 
  (b) Now, assume that after the first pixel is sampled, our remaining pixel samples must come from the same image as the first
 
  (b) Now, assume that after the first pixel is sampled, our remaining pixel samples must come from the same image as the first
 
     pixel.  Given that we first sampled a red pixel, then a green and blue, what is the probability that we have been sampling
 
     pixel.  Given that we first sampled a red pixel, then a green and blue, what is the probability that we have been sampling
     from image 1? (Note that because pixels are not removed from the image after being sampled, the order of pixels sampled
+
     from image 1?
    after the first pixel does not matter.)
+
  
 
  (c) Given the same {R,G,B} scenario as (b), and where pixels after the first must come from the same image, what is the
 
  (c) Given the same {R,G,B} scenario as (b), and where pixels after the first must come from the same image, what is the
Line 35: Line 36:
 
==Solution==
 
==Solution==
  
  '''(a)''' We have been asked to find P(im3|R):
+
  '''(a)'''
<math style='inline'>P(im3|R) = \frac{P(im3 \cap R)}{P(R)}</math> from conditional probability
+
We have been asked to find P(im3|R):
<math style='inline'>= \frac{P(im3)P(R|im3)}{P(R)}</math> from Bayes' rule
+
 
 +
<math style='inline'>P(im_3|R) = \frac{P(im_3 \cap R)}{P(R)}</math> from conditional probability
 +
<math style='inline'>= \frac{P(im_3)P(R|im_3)}{P(R)}</math> from Bayes' rule
  
 
Because the numerators of the previous equation are both given, we are only left to find P(R).  Since the sample space is partitioned into disjoint subsets im1, im2, and im3 (see graphic below), we can use the law of total probability:
 
Because the numerators of the previous equation are both given, we are only left to find P(R).  Since the sample space is partitioned into disjoint subsets im1, im2, and im3 (see graphic below), we can use the law of total probability:
<math style='inline'>P(R) = \sum_{i=1}^{3} P(im_i)P(R|im_i) = \frac{1}{3}(.3 + .8 + .4) = .5</math>
 
  
<math style='inline'>\therefore P(im3|R) = \frac{(\frac{1}{3})(.4)}{.5}</math>
+
<math>P(R) = \sum_{i=1}^{3} P(im_i)P(R|im_i) = \frac{1}{3}(.3 + .8 + .4) = .5</math>
  
'''(b)''' [[Image:Ec1_partition.PNG | Visualization of Our Sample Space]]
+
<math>\therefore P(im_3|R) = \frac{(\frac{1}{3})(.4)}{.5}  = \frac{.4}{1.5}</math>
  
The key to working (b) and (c) is realizing that there is independence of subsequent pixel samplings within an image (intra-image, "locally"),but dependence of subsequent pixel samplings between images (inter-image, "globally").
+
<center>
 +
[[Image:Ec1_partition.PNG | Visualization of Our Sample Space]]
 +
</center>
 +
'''(b)'''
  
Here we define can define a new event comprised of other events {R,G,B}.  First, we need to determine the probability that the first (red) pixel was sampled from image 1.  Using the same method as in '''(a)''':
+
We know that once an image has been sampled from, future color order does not matter since pixels are sampled with replacement.  However, one might think that does matter for the first pixel (when the image number is not yet known).  However, after performing the following calculations, it surprisingly turns out that order does not matter for this.  (I don't know why.)
 +
 
 +
<math>P(im_1|\{R,G,B\}) = \frac{P(im_1|R)P(\{G,B\}|im_1)}{P(\{R,G,B\})}</math>
 +
 
 +
The calculation of <math style='inline'>P(\{R,G,B\})</math> is defined similarly to the numeratorThat is:
 +
 
 +
<math>P(\{R,G,B\}) = \sum_{i=1}^{3} P(im_i|R)P(\{G,B\}|im_i)</math>
 +
 
 +
Since <math style='inline'>P(\{G,B\}|im_i)</math> decomposes into <math style='inline'>P(G|im_i)P(B|im_i)</math> due to independence, the only terms we are missing are P(im1|R) and P(im2|R).  However, we calculated this for image 3 in '''(a)'''.  Using the same method:
 +
 
 +
<math>P(im1|R) = \frac{(\frac{1}{3})(.3)}{.5} = \frac{.3}{1.5}</math>
 +
 
 +
<math>P(im2|R) = \frac{(\frac{1}{3})(.8)}{.5} = \frac{.8}{1.5}</math>
 +
 
 +
<math>\therefore P(im_1|\{R,G,B\}) = \frac{P(im_1|R)P(\{G,B\}|im_1)}{P(\{R,G,B\})} = \frac{(\frac{.3}{1.5})(.1)(.6)}{(\frac{1}{1.5})((.3)(.1)(.6) + (.8)(.1)(.1) + (.4)(.2)(.4))} = \frac{(.3)(.1)(.6)}{((.3)(.1)(.6) + (.8)(.1)(.1) + (.4)(.2)(.4))} \approx .31</math>
 +
 
 +
'''(c)'''
 +
 
 +
Here we're looking for:
 +
 
 +
<math>P(R|\{R,G,B\}) = P(R|im_1)P(im_1|\{R,G,B\}) + P(R|im_2)P(im_2|\{R,G,B\}) + P(R|im_3)P(im_3|\{R,G,B\})</math>
 +
 
 +
<math>= \frac{(.3)(.3)(.1)(.6)}{((.3)(.1)(.6) + (.8)(.1)(.1) + (.4)(.2)(.4))} + \frac{(.8)(.8)(.1)(.1)}{((.3)(.1)(.6) + (.8)(.1)(.1) + (.4)(.2)(.4))} + \frac{(.4)(.4)(.2)(.4)}{((.3)(.1)(.6) + (.8)(.1)(.1) + (.4)(.2)(.4))}</math>
 +
 
 +
<math>\approx (.3)(.31) + (.8)(.14) + (.4)(.55) = .425</math>
 +
 
 +
----
  
<math style='inline'>P(im1|R) = \frac{(\frac{1}{3})(.3)}{.5} = \frac{1}{5}</math>
+
==Notes:==
  
 +
If this problem is 'too long/difficult', just choose '''(a)'''.
  
 +
Although I'm confident in these answers, I don't understand why order does not matter for the calculation of '''(b)'''.
 +
(Perhaps my notation of the conditional probabilities is inaccurate, e.g. I'm not comfortable with events like {R,G,B}.)
  
 
----
 
----
  
 
[[Bonus_point_1_ECE302_Spring2012_Boutin|Back to first bonus point opportunity, ECE302 Spring 2013]]
 
[[Bonus_point_1_ECE302_Spring2012_Boutin|Back to first bonus point opportunity, ECE302 Spring 2013]]

Latest revision as of 18:36, 27 January 2013


Pixel = 'Pic'ture 'El'ement

Smiley 3 1 6.PNG

Above is an RGB image with a pixel color distribution of roughly the following:

red = 30%
green = 10%
blue = 60%

Therefore, if we were to randomly sample a pixel color from the image ('im1'), we would have the following probabilities of sampling each color:

$ P(R|im_1) = .3 $
$ P(G|im_1) = .1 $
$ P(B|im_1) = .6 $

Given the following probabilities of two more RGB images with the same number of pixels as im1, solve (a), (b), and (c). (Note: we assume for this problem that pixels are sampled with replacement.)

$ P(R|im_1) = .3 $     $ P(R|im_2) = .8 $     $ P(R|im_3) = .4 $
$ P(G|im_1) = .1 $     $ P(G|im_2) = .1 $     $ P(G|im_3) = .2 $    } where R = red pixel sampled, etc, and im1 = pixel sampled is in image 1, etc.
$ P(B|im_1) = .6 $     $ P(B|im_2) = .1 $     $ P(B|im_3) = .4 $
(a) Assuming that we have an equal probability of sampling a pixel from each image (ie $ P(im_1) = P(im_2) = P(im_3) = 1/3 $),
    what is the probability that a sampled red pixel came from image 3? 
(b) Now, assume that after the first pixel is sampled, our remaining pixel samples must come from the same image as the first
    pixel.  Given that we first sampled a red pixel, then a green and blue, what is the probability that we have been sampling
    from image 1?
(c) Given the same {R,G,B} scenario as (b), and where pixels after the first must come from the same image, what is the
    probability that our next (fourth) pixel sampled is red?

Solution

(a)

We have been asked to find P(im3|R):

$ P(im_3|R) = \frac{P(im_3 \cap R)}{P(R)} $ from conditional probability $ = \frac{P(im_3)P(R|im_3)}{P(R)} $ from Bayes' rule

Because the numerators of the previous equation are both given, we are only left to find P(R). Since the sample space is partitioned into disjoint subsets im1, im2, and im3 (see graphic below), we can use the law of total probability:

$ P(R) = \sum_{i=1}^{3} P(im_i)P(R|im_i) = \frac{1}{3}(.3 + .8 + .4) = .5 $

$ \therefore P(im_3|R) = \frac{(\frac{1}{3})(.4)}{.5} = \frac{.4}{1.5} $

Visualization of Our Sample Space

(b)

We know that once an image has been sampled from, future color order does not matter since pixels are sampled with replacement. However, one might think that does matter for the first pixel (when the image number is not yet known). However, after performing the following calculations, it surprisingly turns out that order does not matter for this. (I don't know why.)

$ P(im_1|\{R,G,B\}) = \frac{P(im_1|R)P(\{G,B\}|im_1)}{P(\{R,G,B\})} $

The calculation of $ P(\{R,G,B\}) $ is defined similarly to the numerator. That is:

$ P(\{R,G,B\}) = \sum_{i=1}^{3} P(im_i|R)P(\{G,B\}|im_i) $

Since $ P(\{G,B\}|im_i) $ decomposes into $ P(G|im_i)P(B|im_i) $ due to independence, the only terms we are missing are P(im1|R) and P(im2|R). However, we calculated this for image 3 in (a). Using the same method:

$ P(im1|R) = \frac{(\frac{1}{3})(.3)}{.5} = \frac{.3}{1.5} $

$ P(im2|R) = \frac{(\frac{1}{3})(.8)}{.5} = \frac{.8}{1.5} $

$ \therefore P(im_1|\{R,G,B\}) = \frac{P(im_1|R)P(\{G,B\}|im_1)}{P(\{R,G,B\})} = \frac{(\frac{.3}{1.5})(.1)(.6)}{(\frac{1}{1.5})((.3)(.1)(.6) + (.8)(.1)(.1) + (.4)(.2)(.4))} = \frac{(.3)(.1)(.6)}{((.3)(.1)(.6) + (.8)(.1)(.1) + (.4)(.2)(.4))} \approx .31 $

(c)

Here we're looking for:

$ P(R|\{R,G,B\}) = P(R|im_1)P(im_1|\{R,G,B\}) + P(R|im_2)P(im_2|\{R,G,B\}) + P(R|im_3)P(im_3|\{R,G,B\}) $

$ = \frac{(.3)(.3)(.1)(.6)}{((.3)(.1)(.6) + (.8)(.1)(.1) + (.4)(.2)(.4))} + \frac{(.8)(.8)(.1)(.1)}{((.3)(.1)(.6) + (.8)(.1)(.1) + (.4)(.2)(.4))} + \frac{(.4)(.4)(.2)(.4)}{((.3)(.1)(.6) + (.8)(.1)(.1) + (.4)(.2)(.4))} $

$ \approx (.3)(.31) + (.8)(.14) + (.4)(.55) = .425 $


Notes:

If this problem is 'too long/difficult', just choose (a).
Although I'm confident in these answers, I don't understand why order does not matter for the calculation of (b).
(Perhaps my notation of the conditional probabilities is inaccurate, e.g. I'm not comfortable with events like {R,G,B}.)

Back to first bonus point opportunity, ECE302 Spring 2013

Alumni Liaison

Correspondence Chess Grandmaster and Purdue Alumni

Prof. Dan Fleetwood