Exercise 20.3 - Message passing on a tree

Answers

The complete likelihood for this model is:

p ( X 1 , X 2 , X 3 , G 1 , G 2 , G 3 ) = p ( G 1 ) p ( G 2 | G 1 ) p ( G 3 | G 1 ) p ( X 1 | G 1 ) p ( X 2 | G 2 ) p ( X 3 | G 3 ) .

For question (a), we have the following variable elimination:

p ( X 2 = 50 ) = X 1 X 3 G 1 G 2 G 3 p ( X 1 , X 2 = 50 , X 3 , G 1 , G 2 ) = G 1 p ( G 1 ) G 2 p ( G 2 | G 1 ) p ( X 2 = 50 | G 2 ) G 3 p ( G 3 | G 1 ) X 3 p ( X 3 | G 3 ) X 1 p ( X 1 | G 1 ) = G 1 p ( G 1 ) G 2 p ( G 2 | G 1 ) p ( X 2 = 50 | G 2 ) =

On the other hand:

p ( X 2 = 50 , G 1 = 1 ) = X 1 X 3 G 2 G 3 p ( X 1 , X 2 = 50 , X 3 , G 1 = 1 , G 2 ) = p ( G 1 = 1 ) G 2 p ( G 2 | G 1 = 1 ) p ( X 2 = 50 | G 2 ) .

Hence:

p ( G 1 = 1 | X 2 = 50 ) = p ( G 1 = 1 , X 2 = 50 ) p ( X 2 = 50 ) = 0.45 + 0.05 × e 5 0.5 + 0.5 × e 5 0.9 .

For question (b), we have:

p ( X 2 = 50 , X 3 = 50 ) = G 1 p ( G 1 ) G 2 p ( G 2 | G 1 ) p ( X 2 = 50 | G 2 ) G 3 p ( G 3 | G 1 ) p ( X 3 = 50 | G 3 ) X 1 p ( X 1 | G 1 ) ,

similar for p ( G 1 = 1 , X 2 = 50 , X 3 = 50 ) , hence:

p ( G 1 = 1 | X 2 = 50 , X 3 = 50 ) = p ( G 1 = 1 , X 2 = 50 , X 3 = 50 ) p ( X 2 = 50 , X 3 = 50 ) = 0.405 + 0.09 × e 5 + 0.005 × e 10 0.41 + 0.18 × e 5 + 0.41 × e 10 0.99 .

Since both X 2 , X 3 equal fifty implies that both G 2 , G 3 are approximately unity, this confirms the belief that G 1 = 1 .

For question (c), the setting is symmetric to (b), hence:

p ( G 1 = 0 | X 2 = 60 , X 3 = 60 ) 0.01 .

For question (d), we follow the analysis as before:

p ( G 1 = 1 | X 2 = 50 , X 3 = 60 ) = p ( G 1 = 1 ) p ( X 2 = 50 , X 3 = 60 ) | G 1 ) p ( X 2 = 50 , X 3 = 60 ) .

Note that:

p ( X 2 = 50 , X 3 = 60 ) | G 1 ) = p ( X 2 = 50 , X 3 = 60 ) | G 10 )

by symmetry, hence:

p ( G 1 = 1 | X 2 = 50 , X 3 = 60 ) = 1 2 .

Since two contradictive evidence cancell each other.

This is a soft/fuzzy version of the ordinary logical reasoning, where we could have set p ( G 2 , 3 | G 1 ) to then identity matrices and the variance for generating X i to zero.

User profile picture
2021-03-24 13:42
Comments