# Difference between revisions of "stat946f11"

(→Basic Statistical Problems) |
|||

Line 869: | Line 869: | ||

* Clustering | * Clustering | ||

* Density Estimation | * Density Estimation | ||

− | + | ||

## Revision as of 10:38, 12 October 2011

## Contents

## Editor Sign Up

## Sign up for your presentation

## Assignments

## Introduction

### Notation

We will begin with short section about the notation used in these notes. Capital letters will be used to denote random variables and lower case letters denote observations for those random variables:

- [math]\{X_1,\ X_2,\ \dots,\ X_n\}[/math] random variables
- [math]\{x_1,\ x_2,\ \dots,\ x_n\}[/math] observations of the random variables

The joint *probability mass function* can be written as:

or as shorthand, we can write this as [math]p( x_1, x_2, \dots, x_n )[/math]. In these notes both types of notation will be used. We can also define a set of random variables [math]X_Q[/math] where [math]Q[/math] represents a set of subscripts.

### Example

Let [math]A = \{1,4\}[/math], so [math]X_A = \{X_1, X_4\}[/math]; [math]A[/math] is the set of indices for
the r.v. [math]X_A[/math].

Also let [math]B = \{2\},\ X_B = \{X_2\}[/math] so we can write

### Graphical Models

Graphs can be represented as a pair of vertices and edges: [math]G = (V, E).[/math]

Two branches of graphical representations of distributions are commonly used in graphical models; Bayesian networks and Markov networks. Both families encompass the properties of factorization and independence, but they differ in the factorization of the distribution that they induce.

- [math]V[/math] is the set of nodes (vertices).
- [math]E[/math] is the set of edges.

If the edges have a direction associated with them then we consider the graph to be directed as in Figure 1, otherwise the graph is undirected as in Figure 2.

We will use graphs in this course to represent the relationship between different random variables.

#### Directed graphical models (Bayesian networks)

In the case of directed graphs, the direction of the arrow indicates "causation". For example:

[math]A \longrightarrow B[/math]: [math]A\,\![/math] "causes" [math]B\,\![/math].

In this case we must assume that our directed graphs are *acyclic*. If our causation graph contains a cycle then it would mean that for example:

- [math]A[/math] causes [math]B[/math]
- [math]B[/math] causes [math]C[/math]
- [math]C[/math] causes [math]A[/math], again.

Clearly, this would confuse the order of the events. An example of a graph with a cycle can be seen in Figure 3. Such a graph could not be used to represent causation. The graph in Figure 4 does not have cycle and we can say that the node [math]X_1[/math] causes, or affects, [math]X_2[/math] and [math]X_3[/math] while they in turn cause [math]X_4[/math].

We will consider a 1-1 map between our graph's vertices and a set of random variables. Consider the following example that uses boolean random variables. It is important to note that the variables need not be boolean and can indeed be discrete over a range or even continuous.

Speaking about random variables, we can now refer to the relationship between random variables in terms of dependence. Therefore, the direction of the arrow indicates "conditional dependence". For example:

[math]A \longrightarrow B[/math]: [math]B\,\![/math] "is dependent on" [math]A\,\![/math].

#### Example

In this example we will consider the possible causes for wet grass.

The wet grass could be caused by rain, or a sprinkler. Rain can be caused by clouds. On the other hand one can not say that clouds cause the use of a sprinkler. However, the causation exists because the presence of clouds does affect whether or not a sprinkler will be used. If there are more clouds there is a smaller probability that one will rely on a sprinkler to water the grass. As we can see from this example the relationship between two variables can also act like a negative correlation. The corresponding graphical model is shown in Figure 5.

This directed graph shows the relation between the 4 random variables. If we have the joint probability [math]P(C,R,S,W)[/math], then we can answer many queries about this system.

This all seems very simple at first but then we must consider the fact that in the discrete case the joint probability function grows exponentially with the number of variables. If we consider the wet grass example once more we can see that we need to define [math]2^4 = 16[/math] different probabilities for this simple example. The table bellow that contains all of the probabilities and their corresponding boolean values for each random variable is called an *interaction table*.

**Example:**

Now consider an example where there are not 4 such random variables but 400. The interaction table would become too large to manage. In fact, it would require [math]2^{400}[/math] rows! The purpose of the graph is to help avoid this intractability by considering only the variables that are directly related. In the wet grass example Sprinkler (S) and Rain (R) are not directly related.

To solve the intractability problem we need to consider the way those relationships are represented in the graph. Let us define the following parameters. For each vertex [math]i \in V[/math],

- [math]\pi_i[/math]: is the set of parents of [math]i[/math]
- ex. [math]\pi_R = C[/math] \ (the parent of [math]R = C[/math])

- [math]f_i(x_i, x_{\pi_i})[/math]: is the joint p.d.f. of [math]i[/math] and [math]\pi_i[/math] for which it is true that:
- [math]f_i[/math] is nonnegative for all [math]i[/math]
- [math]\displaystyle\sum_{x_i} f_i(x_i, x_{\pi_i}) = 1[/math]

**Claim**: There is a family of probability functions [math] P(X_V) = \prod_{i=1}^n f_i(x_i, x_{\pi_i})[/math] where this function is nonnegative, and

To show the power of this claim we can prove the equation (\ref{eqn:WetGrass}) for our wet grass example:

We want to show that

Consider factors [math]f(C)[/math], [math]f(R,C)[/math], [math]f(S,C)[/math]: they do not depend on [math]W[/math], so we can write this all as

since we had already set [math]\displaystyle \sum_{x_i} f_i(x_i, x_{\pi_i}) = 1[/math].

Let us consider another example with a different directed graph.

**Example:**

Consider the simple directed graph in Figure 6.

Assume that we would like to calculate the following: [math] p(x_3|x_2) [/math]. We know that we can write the joint probability as:

We can also make use of Bayes' Rule here:

We also need

Thus,

**Theorem 1.**

In our simple graph, the joint probability can be written as

Instead, had we used the chain rule we would have obtained a far more complex equation:

The *Markov Property*, or *Memoryless Property* is when the variable [math]X_i[/math] is only affected by [math]X_j[/math] and so the random variable [math]X_i[/math] given [math]X_j[/math] is independent of every other random variable. In our example the history of [math]x_4[/math] is completely determined by [math]x_3[/math].

By simply applying the Markov Property to the chain-rule formula we would also have obtained the same result.

Now let us consider the joint probability of the following six-node example found in Figure 7.

If we use Theorem 1 it can be seen that the joint probability density function for Figure 7 can be written as follows:

Once again, we can apply the Chain Rule and then the Markov Property and arrive at the same result.

### Independence

#### Marginal independence

We can say that [math]X_A[/math] is marginally independent of [math]X_B[/math] if:

#### Conditional independence

We can say that [math]X_A[/math] is conditionally independent of [math]X_B[/math] given [math]X_C[/math] if:

**Aside:** Before we move on further, we first define the following terms:

- I is defined as an ordering for the nodes in graph C.
- For each [math]i \in V[/math], [math]V_i[/math] is defined as a set of all nodes that appear earlier than i excluding its parents [math]\pi_i[/math].

Let us consider the example of the six node figure given above (Figure 7). We can define [math]I[/math] as follows:

We can then easily compute [math]V_i[/math] for say [math]i=3,6[/math].

while [math]\pi_i[/math] for [math] i=3,6[/math] will be.

We would be interested in finding the conditional independence between random variables in this graph. We know [math]X_i \perp X_{v_i} | X_{\pi_i}[/math] for each [math]i[/math]. In other words, given its parents the node is independent of all earlier nodes. So:

[math]X_1 \perp \phi | \phi[/math],

[math]X_2 \perp \phi | X_1[/math],

[math]X_3 \perp X_2 | X_1[/math],

[math]X_4 \perp \{X_1,X_3\} | X_2[/math],

[math]X_5 \perp \{X_1,X_2,X_4\} | X_3[/math],

[math]X_6 \perp \{X_1,X_3,X_4\} | \{X_2,X_5\}[/math]

To illustrate why this is true we can take a simple example. Show that:

Proof: first, we know [math]P(X_1,X_2,X_3,X_4,X_5,X_6) = P(X_1)P(X_2|X_1)P(X_3|X_1)P(X_4|X_2)P(X_5|X_3)P(X_6|X_5,X_2)\,\![/math]

then

The other conditional independences can be proven through a similar process.

#### Sampling

Even if using graphical models helps a lot facilitate obtaining the joint probability, exact inference is not always feasible. Exact inference is feasible in small to medium-sized networks only. Exact inference consumes such a long time in large networks. Therefore, we resort to approximate inference techniques which are much faster and usually give pretty good results.

In sampling, random samples are generated and values of interest are computed from samples, not original work.

As an input you have a Bayesian network with set of nodes [math]X\,\![/math]. The sample taken may include all variables (except evidence E) or a subset. Sample schemas dictate how to generate samples (tuples). Ideally samples are distributed according to [math]P(X|E)\,\![/math]

Some sampling algorithms:

- Forward Sampling
- Likelihood weighting
- Gibbs Sampling (MCMC)
- Blocking
- Rao-Blackwellised

- Importance Sampling

## Bayes Ball

The Bayes Ball algorithm can be used to determine if two random variables represented in a graph are independent. The algorithm can show that either two nodes in a graph are independent OR that they are not necessarily independent. The Bayes Ball algorithm can not show that two nodes are dependant. The algorithm will be discussed further in later parts of this section.

### Canonical Graphs

In order to understand the Bayes Ball algorithm we need to first introduce 3 canonical graphs.

#### Markov Chain (also called serial connection)

In the following graph (Figure 8 X is independent of Z given Y.

We say that: [math]X[/math] [math]\perp[/math] [math]Z[/math] [math]|[/math] [math]Y[/math]

We can prove this independence:

Where

#### Hidden Cause (diverging connection)

In the Hidden Cause case we can say that X is independent of Z given Y. In this case Y is the hidden cause and if it is known then Z and X are considered independent.

We say that: [math]X[/math] [math]\perp[/math] [math]Z[/math] [math]|[/math] [math]Y[/math]

The proof of the independence:

The Hidden Cause case is best illustrated with an example:

In Figure 10 it can be seen that both "Shoe Size" and "Grey Hair" are dependant on the age of a person. The variables of "Shoe size" and "Grey hair" are dependent in some sense, if there is no "Age" in the picture. Without the age information we must conclude that those with a large shoe size also have a greater chance of having gray hair. However, when "Age" is observed, there is no dependence between "Shoe size" and "Grey hair" because we can deduce both based only on the "Age" variable.

#### Explaining-Away (converging connection)

Finally, we look at the third type of canonical graph:
*Explaining-Away Graphs*. This type of graph arises when a
phenomena has multiple explanations. Here, the conditional
independence statement is actually a statement of marginal
independence: [math]X \perp Z[/math].

In these types of scenarios, variables X and Z are independent. However, once the third variable Y is observed, X and Z become dependent (Fig. 11).

To clarify these concepts, suppose Bob and Mary are supposed to meet for a noontime lunch. Consider the following events:

If Mary is late, then she could have been kidnapped by aliens. Alternatively, Bob may have forgotten to adjust his watch for daylight savings time, making him early. Clearly, both of these events are independent. Now, consider the following probabilities:

We expect [math]P( late = 1 ) \lt P( aliens = 1 ~|~ late = 1 )[/math] since [math]P(
aliens = 1 ~|~ late = 1 )[/math] does not provide any information
regarding Bob's watch. Similarly, we expect [math]P( aliens = 1 ~|~
late = 1 ) \lt P( aliens = 1 ~|~ late = 1, watch = 0 )[/math]. Since
[math]P( aliens = 1 ~|~ late = 1 ) \neq P( aliens = 1 ~|~ late = 1, watch = 0 )[/math], *aliens* and
*watch* are not independent given *late*. To summarize,

- If we do not observe
*late*, then*aliens*[math]~\perp~ watch[/math] ([math]X~\perp~ Z[/math]) - If we do observe
*late*, then*aliens*[math] ~\cancel{\perp}~ watch ~|~ late[/math] ([math]X ~\cancel{\perp}~ Z ~|~ Y[/math])

### Bayes Ball Algorithm

**Goal:** We wish to determine whether a given conditional
statement such as [math]X_{A} ~\perp~ X_{B} ~|~ X_{C}[/math] is true given a directed graph.

The algorithm is as follows:

- Shade nodes, [math]~X_{C}~[/math], that are conditioned on.
- The initial position of the ball is [math]~X_{A}~[/math].
- If the ball cannot reach [math]~X_{B}~[/math], then the nodes [math]~X_{A}~[/math] and [math]~X_{B}~[/math] must be conditionally independent.
- If the ball can reach [math]~X_{B}~[/math], then the nodes [math]~X_{A}~[/math] and [math]~X_{B}~[/math] are not necessarily independent.

The biggest challenge in the *Bayes Ball Algorithm* is to
determine what happens to a ball going from node X to node Z as it
passes through node Y. The ball could continue its route to Z or
it could be blocked. It is important to note that the balls are
allowed to travel in any direction, independent of the direction
of the edges in the graph.

We use the canonical graphs previously studied to determine the route of a ball traveling through a graph. Using these three graphs we establish base rules which can be extended upon for more general graphs.

#### Markov Chain (serial connection)

A ball traveling from X to Z or from Z to X will be blocked at node Y if this node is shaded. Alternatively, if Y is unshaded, the ball will pass through.

In (Fig. 12(a)), X and Z are conditionally independent ( [math]X ~\perp~ Z ~|~ Y[/math] ) while in (Fig.12(b)) X and Z are not necessarily independent.

#### Hidden Cause (diverging connection)

A ball traveling through Y will be blocked at Y if it is shaded. If Y is unshaded, then the ball passes through.

(Fig. 13(a)) demonstrates that X and Z are conditionally independent when Y is shaded.

#### Explaining-Away (converging connection)

A ball traveling through Y is blocked when Y is unshaded. If Y is shaded, then the ball passes through. Hence, X and Z are conditionally independent when Y is unshaded.

### Bayes Ball Examples

#### Example 1

In this first example, we wish to identify the behavior of a ball going from X to Y in two-node graphs.

The four graphs in (Fig. 15 show different
scenarios. In (a), the ball is blocked at Y. In (b) the ball
passes through Y. In both of these cases, we use the rules of the
*Explaining Away Canonical Graph* (refer to
Fig. 14.) Finally, for the last two
graphs, we used the rules of the *Hidden Cause Canonical*
Graph* (Fig. 13). In (c), the ball passes through*
Y while in (d), the ball is blocked at Y.

#### Example 2

Suppose your home is equipped with an alarm system. There are two possible causes for the alarm to ring:

- Your house is being burglarized
- There is an earthquake

Hence, we define the following events:

The *burglary* and *earthquake* events are independent
if the alarm does not ring. However, if the alarm does ring, then
the *burglary* and the *earthquake* events are not
necessarily independent. Also, if the alarm rings then it is
possible for a police report to be issued.

We can use the *Bayes Ball Algorithm* to deduce conditional
independence properties from the graph. Firstly, consider figure
(16(a)) and assume we are trying to determine
whether there is conditional independence between the
*burglary* and *earthquake* events. In figure
(\ref{fig:AlarmExample1}(a)), a ball starting at the *burglary*
event is blocked at the *alarm* node.

Nonetheless, this does not prove that the *burglary* and
*earthquake* events are independent. Indeed,
(Fig. 16(b)) disproves this as we have found an
alternate path from *burglary* to *earthquake* passing
through *report*. It follows that [math]burglary
~\cancel{\amalg}~ earthquake ~|~ report[/math]

#### Example 3

Referring to figure (Fig. 17), we wish to determine whether the following conditional probabilities are true:

To determine if the conditional probability Eq.\ref{eq:c1} is true, we shade node [math]X_{2}[/math]. This blocks balls traveling from [math]X_{1}[/math] to [math]X_{3}[/math] and proves that Eq.\ref{eq:c1} is valid.

After shading nodes [math]X_{3}[/math] and [math]X_{4}[/math] and applying the *Bayes Balls Algorithm}, we find that the ball travelling from [math]X_{1}[/math] to [math]X_{5}[/math] is blocked at [math]X_{3}[/math]. Similarly, a ball going from [math]X_{5}[/math] to [math]X_{1}[/math] is blocked at [math]X_{4}[/math]. This proves that Eq.\ref{eq:c2* also holds.

#### Example 4

Consider figure (Fig. 18). Using the *Bayes Ball Algorithm* we wish to determine if each of the following
statements are valid:

To disprove Eq.\ref{eq:c3}, we must find a path from [math]X_{4}[/math] to [math]X_{1}[/math] and [math]X_{3}[/math] when [math]X_{2}[/math] is shaded (Refer to Fig. 19(a)). Since there is no route from [math]X_{4}[/math] to [math]X_{1}[/math] and [math]X_{3}[/math] we conclude that Eq.\ref{eq:c3} is true.

Similarly, we can show that there does not exist a path between [math]X_{1}[/math] and [math]X_{6}[/math] when [math]X_{2}[/math] and [math]X_{3}[/math] are shaded (Refer to Fig.19(b)). Hence, Eq.\ref{eq:c4} is true.

Finally, (Fig. 19(c)) shows that there is a route from [math]X_{2}[/math] to [math]X_{3}[/math] when [math]X_{1}[/math] and [math]X_{6}[/math] are shaded. This proves that the statement \ref{eq:c4} is false.

**Theorem 2.**

Define [math]p(x_{v}) = \prod_{i=1}^{n}{p(x_{i} ~|~ x_{\pi_{i}})}[/math] to be the factorization as a multiplication of some local probability of a directed graph.

Let [math]D_{1} = \{ p(x_{v}) = \prod_{i=1}^{n}{p(x_{i} ~|~ x_{\pi_{i}})}\}[/math]

Let [math]D_{2} = \{ p(x_{v}):[/math]satisfy all conditional independence statements associated with a graph [math]\}[/math].

Then [math]D_{1} = D_{2}[/math].

#### Example 5

Given the following Bayesian network (Fig.19 ): Determine whether the following statements are true or false?

a.) [math]x4\perp \{x1,x3\}[/math]

Ans. True

b.) [math]x1\perp x6\{x2,x3\}[/math]

Ans. True

c.) [math]x2\perp x3 \{x1,x6\}[/math]

Ans. False

## Undirected Graphical Model

Generally, the graphical model is divided into two major classes, directed graphs and undirected graphs. Directed graphs and its characteristics was described previously. In this section we discuss undirected graphical model which is also known as Markov random fields. We can define an undirected graphical model with a graph [math] G = (V, E)[/math] where [math] V [/math] is a set of vertices corresponding to a set of random variables and [math] E [/math] is a set of undirected edges as shown in (Fig.20)

#### Conditional independence

For directed graphs Bayes ball method was defined to determine the conditional independence properties of a given graph. We can also employ the Bayes ball algorithm to examine the conditional independency of undirected graphs. Here the Bayes ball rule is simpler and more intuitive. Considering (Fig.21) , a ball can be thrown either from x to z or from z to x if y is not observed. In other words, if y is not observed a ball thrown from x can reach z and vice versa. On the contrary, given a shaded y, the node can block the ball and make x and z conditionally independent. With this definition one can declare that in an undirected graph, a node is conditionally independent of non-neighbors given neighbors. Technically speaking, [math]X_A[/math] is independent of [math]X_C[/math] given [math]X_B[/math] if the set of nodes [math]X_B[/math] separates the nodes [math]X_A[/math] from the nodes [math]X_C[/math]. Hence, if every path from a node in [math]X_A[/math] to a node in [math]X_C[/math] includes at least one node in [math]X_B[/math], then we claim that [math] X_A \perp X_c | X_B [/math].

#### Question

Is it possible to convert undirected models to directed models or vice versa?

In order to answer this question, consider (Fig.22 ) which illustrates an undirected graph with four nodes - [math]X[/math], [math]Y[/math],[math]Z[/math] and [math]W[/math]. We can define two facts using Bayes ball method:

It is simple to see there is no directed graph satisfying both conditional independence properties. Recalling that directed graphs are acyclic, converting undirected graphs to directed graphs result in at least one node in which the arrows are inward-pointing(a v structure). Without loss of generality we can assume that node [math]Z[/math] has two inward-pointing arrows. By conditional independence semantics of directed graphs, we have [math] X \perp Y|W[/math], yet the [math]X \perp Y|\{W,Z\}[/math] property does not hold. On the other hand, (Fig.23 ) depicts a directed graph which is characterized by the singleton independence statement [math]X \perp Y [/math]. There is no undirected graph on three nodes which can be characterized by this singleton statement. Basically, if we consider the set of all distribution over [math]n[/math] random variables, a subset of which can be represented by directed graphical models while there is another subset which undirected graphs are able to model that. There is a narrow intersection region between these two subsets in which probabilistic graphical models may be represented by either directed or undirected graphs.

# Elimination Algorithm

## Elimination Algorithm on Directed Graphs

Given a graph G =(V,*E*), an evidence set E, and a query node F, we first choose an elimination ordering I such that F appears last in this ordering. The following figure shows the steps required to perform the elimination algorithm for probabilistic inference on directed graphs:

{{

Template:namespace detect

| type = style
| image =
| imageright =
| style =
| textstyle =
| text = This article **may require cleanup to meet Wikicoursenote's quality standards.** The specific problem is: *Please type the algorithms rather than having them as images from other sources.*. Please improve this article if you can. *(October 7 2011)*
| small =
| smallimage =
| smallimageright =
| smalltext =

**Example:**

For the graph in figure 21 [math]G =(V,''E'')[/math]. Consider once again that node [math]x_1[/math] is the query node and [math]x_6[/math] is the evidence node.

[math]I = \left\{6,5,4,3,2,1\right\}[/math] (1 should be the last node, ordering is crucial)

We must now create an active list. There are two rules that must be followed in order to create this list.

- For i[math]\in{V}[/math] put [math]p(x_i|x_{\pi_i})[/math] in active list.
- For i[math]\in[/math]{E} put [math]p(x_i|\overline{x_i})[/math] in active list.

Here, our active list is: [math] p(x_1), p(x_2|x_1), p(x_3|x_1), p(x_4|x_2), p(x_5|x_3),\underbrace{p(x_6|x_2, x_5)\delta{(\overline{x_6},x_6)}}_{\phi_6(x_2,x_5, x_6), \sum_{x6}{\phi_6}=m_{6}(x2,x5) }[/math]

We first eliminate node [math]X_6[/math]. We place [math]m_{6}(x_2,x_5)[/math] on the active list, having removed [math]X_6[/math]. We now eliminate [math]X_5[/math].

Likewise, we can also eliminate [math]X_4, X_3, X_2[/math](which yields the unnormalized conditional probability [math]p(x_1|\overline{x_6})[/math] and [math]X_1[/math]. Then it yields [math]m_1 = \sum_{x_1}{\phi_1(x_1)}[/math] which is the normalization factor, [math]p(\overline{x_6})[/math].

## Elimination Algorithm on Undirected Graphs

The first task is to find the maximal cliques and their associated potential functions.

maximal clique: [math]\left\{x_1, x_2\right\}[/math], [math]\left\{x_1, x_3\right\}[/math], [math]\left\{x_2, x_4\right\}[/math], [math]\left\{x_3, x_5\right\}[/math], [math]\left\{x_2,x_5,x_6\right\}[/math]

potential functions: [math]\varphi{(x_1,x_2)},\varphi{(x_1,x_3)},\varphi{(x_2,x_4)}, \varphi{(x_3,x_5)}[/math] and [math]\varphi{(x_2,x_3,x_6)}[/math]

[math] p(x_1|\overline{x_6})=p(x_1,\overline{x_6})/p(\overline{x_6})\cdots\cdots\cdots\cdots\cdots(*) [/math]

[math]p(x_1,x_6)=\frac{1}{Z}\sum_{x_2,x_3,x_4,x_5,x_6}\varphi{(x_1,x_2)}\varphi{(x_1,x_3)}\varphi{(x_2,x_4)}\varphi{(x_3,x_5)}\varphi{(x_2,x_3,x_6)}\delta{(x_6,\overline{x_6})} [/math]

The [math]\frac{1}{Z}[/math] looks crucial, but in fact it has no effect because for (*) both the numerator and the denominator have the [math]\frac{1}{Z}[/math] term. So in this case we can just cancel it.

The general rule for elimination in an undirected graph is that we can remove a node as long as we connect all of the parents of that node together. Effectively, we form a clique out of the parents of that node.
The algorithm used to eliminate nodes in an undirected graph is given in figure 23.

{{

Template:namespace detect

| type = style
| image =
| imageright =
| style =
| textstyle =
| text = This article **may require cleanup to meet Wikicoursenote's quality standards.** The specific problem is: *Please type the algorithms rather than having them as images from other sources.*. Please improve this article if you can. *(October 7 2011)*
| small =
| smallimage =
| smallimageright =
| smalltext =

**Example: **

For the graph G in figure 24

when we remove x1, G becomes as in figure 25

while if we remove x2, G becomes as in figure 26

An interesting thing to point out is that the order of the elimination matters a great deal. Consider the two results. If we remove one node the graph complexity is slightly reduced. But if we try to remove another node the complexity is significantly increased. The reason why we even care about the complexity of the graph is because the complexity of a graph denotes the number of calculations that are required to answer questions about that graph. If we had a huge graph with thousands of nodes the order of the node removal would be key in the complexity of the algorithm. Unfortunately, there is no efficient algorithm that can produce the optimal node removal order such that the elimination algorithm would run quickly.

## Moralization

So far we have shown how to use elimination to successively remove nodes from an undirected graph. We know that this is useful in the process of marginalization. We can now turn to the question of what will happen when we have a directed graph. It would be nice if we could somehow reduce the directed graph to an undirected form and then apply the previous elimination algorithm. This reduction is called moralization and the graph that is produced is called a moral graph.

To moralize a graph we first need to connect the parents of each node together. This makes sense intuitively because the parents of a node need to be considered together in the undirected graph and this is only done if they form a type of clique. By connecting them together we create this clique.

After the parents are connected together we can just drop the orientation on the edges in the directed graph. By removing the directions we force the graph to become undirected.

The previous elimination algorithm can now be applied to the new moral graph. We can do this by assuming that the probability functions in directed graph [math] P(x_i|\pi_{x_i}) [/math] are the same as the mass functions from the undirected graph. [math] \psi_{c_i}(c_{x_i}) [/math]

**Example:**

I = [math]\left\{x_6,x_5,x_4,x_3,x_2,x_1\right\}[/math]

When we moralize the directed graph in figure 27, we obtain the
undirected graph in figure 28.

# Elimination Algorithm on Trees

**Definition of a tree:**

A tree is an undirected graph in which any two vertices are connected by exactly one simple path. In other words, any connected graph without cycles is a tree.

If we have a directed graph then we must moralize it first. If the moral graph is a tree then the directed graph is also considered a tree.

## Belief Propagation Algorithm (Sum Product Algorithm)

One of the main disadvantages to the elimination algorithm is that the ordering of the nodes defines the number of calculations that are required to produce a result. The optimal ordering is difficult to calculate and without a decent ordering the algorithm may become very slow. In response to this we can introduce the sum product algorithm. It has one major advantage over the elimination algorithm: it is faster. The sum product algorithm has the same complexity when it has to compute the probability of one node as it does to compute the probability of all the nodes in the graph. Unfortunately, the sum product algorithm also has one disadvantage. Unlike the elimination algorithm it can not be used on any graph. The sum product algorithm works only on trees.

For undirected graphs if there is only one path between any two pair of nodes then that graph is a tree (Fig.29). If we have a directed graph then we must moralize it first. If the moral graph is a tree then the directed graph is also considered a tree (Fig.30).

For the undirected graph [math]G(v, \varepsilon)[/math] (Fig.30) we can write the joint probability distribution function in the following way.

We know that in general we can not convert a directed graph into an undirected graph. There is however an exception to this rule when it comes to trees. In the case of a directed tree there is an algorithm that allows us to convert it to an undirected tree with the same properties.

Take the above example (Fig.30) of a directed tree. We can write the joint probability distribution function as:

If we want to convert this graph to the undirected form shown in (Fig. \ref{fig:UnDirTree}) then we can use the following set of rules. \begin{thinlist}

- If [math]\gamma[/math] is the root then: [math] \psi(x_\gamma) = P(x_\gamma) [/math].
- If [math]\gamma[/math] is NOT the root then: [math] \psi(x_\gamma) = 1 [/math].
- If [math]\left\lbrace i \right\rbrace[/math] = [math]\pi_j[/math] then: [math] \psi(x_i, x_j) = P(x_j | x_i) [/math].

\end{thinlist} So now we can rewrite the above equation for (Fig.30) as:

## Elimination Algorithm on a Tree

We will derive the Sum-Product algorithm from the point of view of the Eliminate algorithm. To marginalize [math]x_1[/math] in Fig.31,

where,

which is essentially (potential of the node)[math]\times[/math](potential of the edge)[math]\times[/math](message from the child).

The term "[math]m_{ji}(x_i)[/math]" represents the intermediate factor between the eliminated variable, *j*, and the remaining neighbor of the variable, *i*. Thus, in the above case, we will use [math]m_{53}(x_3)[/math] to denote [math]m_5(x_3)[/math], [math]m_{42}(x_2)[/math] to denote
[math]m_4(x_2)[/math], and [math]m_{32}(x_2)[/math] to denote [math]m_3(x_2)[/math]. We refer to the
intermediate factor [math]m_{ji}(x_i)[/math] as a "message" that *j*
sends to *i*. (Fig. \ref{fig:TreeStdEx})

## Elimination To Sum Product Algorithm

The Sum-Product algorithm allows us to compute all marginals in the tree by passing messages inward from the leaves of the tree to an (arbitrary) root, and then passing it outward from the root to the leaves, again using the above equation at each step. The net effect is that a single message will flow in both directions along each edge. (See Fig.32) Once all such messages have been computed using the above equation, we can compute desired marginals.

As shown in Fig.32, to compute the marginal of [math]X_1[/math] using elimination, we eliminate [math]X_5[/math], which involves computing a message [math]m_{53}(x_3)[/math], then eliminate [math]X_4[/math] and [math]X_3[/math] which involves messages [math]m_{32}(x_2)[/math] and [math]m_{42}(x_2)[/math]. We subsequently eliminate [math]X_2[/math], which creates a message [math]m_{21}(x_1)[/math].

Suppose that we want to compute the marginal of [math]X_2[/math]. As shown in Fig.33, we first eliminate [math]X_5[/math], which creates [math]m_{53}(x_3)[/math], and then eliminate [math]X_3[/math], [math]X_4[/math], and [math]X_1[/math], passing messages [math]m_{32}(x_2)[/math], [math]m_{42}(x_2)[/math] and [math]m_{12}(x_2)[/math] to [math]X_2[/math].

Since the messages can be "reused", marginals over all possible elimination orderings can be computed by computing all possible messages which is small in numbers compared to the number of possible elimination orderings.

The Sum-Product algorithm is not only based on the above equation, but also *Message-Passing Protocol*.
**Message-Passing Protocol** tells us that a node can
send a message to a neighboring node when (and only when) it has
received messages from all of its other neighbors.

### For Directed Graph

Previously we stated that:

Using the above equation (\ref{eqn:Marginal}), we find the marginal of [math]\bar{x}_E[/math].

Now we denote:

Since the sets, *F* and *E*, add up to [math]\mathcal{V}[/math],
[math]p(x_v)[/math] is equal to [math]p(x_F,x_E)[/math]. Thus we can substitute the
equation (\ref{eqn:Dir8}) into (\ref{eqn:Marginal}) and (\ref{eqn:Dir7}), and they become:

We are interested in finding the conditional probability. We substitute previous results, (\ref{eqn:Dir9}) and (\ref{eqn:Dir10}) into the conditional probability equation.

[math]p^E(x_v)[/math] is an unnormalized version of conditional probability, [math]p(x_F|\bar{x}_E)[/math].

### For Undirected Graphs

We denote [math]\psi^E[/math] to be:

## Max-Product

We would like to find the Maximum probability that can be achieved by some set of random variables given a set of configurations. The algorithm is similar to the sum product except we replace the sum with max.

[math]p(x_F|\bar{x}_E)[/math]

**Example:**
Consider the graph in Figure \ref{fig:MaxProdEx}.

## Maximum configuration

We would also like to find the value of the [math]x_i[/math]s which produces the largest value for the given expression. To do this we replace the max from the previous section with argmax.

[math]m_{53}(x_5)= argmax_{x_5}\psi{(x_5)}\psi{(x_5,x_3)}[/math]

[math]\log{m^{max}_{ji}(x_i)}=\max_{x_j}{\log{\psi^{E}{(x_j)}}}+\log{\psi{(x_i,x_j)}}+\sum_{k\in{N(j)\backslash{i}}}\log{m^{max}_{kj}{(x_j)}}[/math]

In many cases we want to use the log of this expression because the numbers tend to be very high. Also, it is important to note that this also works in the continuous case where we replace the summation sign with an integral.

# Basic Statistical Problems

In statistics there are a number of different 'standard' problems that always appear in one form or another. They are as follows: \begin{thinlist}

- Regression
- Classification
- Clustering
- Density Estimation

## Regression

In regression we have a set of data points [math] (x_i, y_i) [/math] for [math] i = 1...n [/math] and we would like to determine the way that the variables x and y are related. In certain cases such as (Fig. \ref{img:regression.eps}) we try to fit a line (or other type of function) through the points in such a way that it describes the relationship between the two variables.

Once the relationship has been determined we can give a functional value to the following expression. In this way we can determine the value (or distribution) of y if we have the value for x. [math]P(y|x)=\frac{P(y,x)}{P(x)} = \frac{P(y,x)}{\int_{y}{P(y,x)dy}}[/math]

## Classification

In classification we also have a set of points [math] (x_i, y_i) [/math] for [math] i = 1...n [/math] but we would like to use the x and y values to determine if a certain point belongs in group A or in group B. Consider the example in (Fig. \ref{img:Classification.eps}) where two sets of points have been divided into the set + and the set - by a line. The purpose of classification is to find this line and then place any new points into one group or the other.

We would like to obtain the probability distribution of to following equation where c is the class and x and y are the data points. In simple terms we would like to find the probability that this point is in class c when we know that the values of X and Y are x and y.

## Clustering

Clustering is somewhat like classification only that we do not know the groups before we gather and examine the data. We would like to find the probability distribution of the following equation without knowing the value of y.

We can use graphs to represent the three types of statistical problems that have been introduced so far. The first graph (Fig. \ref{fig:RegClass} can be used to represent either the Regression or the Classification problem because both the X and the Y variables are known. The second graph (Fig. \ref{fig:Clustering}) we see that the value of the Y variable is unknown and so we can tell that this graph represents the Clustering situation.

**Classification example: Naive Bayes classifier**

First define a set of boolean random variables [math]X_i[/math] and [math]Y[/math] for [math]i = 1...n[/math].

Then we will say that a certain pattern of Xs can either be classified as a 1 or a 0. The result of this classification will be represented by the variable Y. The graphical representation is shown in (Fig. \ref{img:classifi.eps}). One important thing to note here is that the two diagrams represent the same graph. The one on the right uses plate notation to simplify the representation of the graph for variables that are indexed. Such plate notation will also be used later in these notes.

\begin{tabular}{ccc}
[math] \stackrel{x}{\underbrace{\lt 01110\gt }_{n}}[/math] & [math] \rightarrow [/math] & [math]\stackrel{Y}{1}[/math]

[math] \lt 01110\gt [/math] & [math] \rightarrow [/math] & [math]0[/math]
\end{tabular}

We are interested in finding the following:

The classification is very intuitive in this case. We will calculate the probability that we are in class 1 and we will calculate the probability that we are in class 0. The higher probability will decide the class. For example if we have a higher probability of being in class 1 then we will place this set of Xs in class 1.

\begin{tabular}{ ccc }
[math] \widehat{y}=1 [/math] & [math] \Leftrightarrow [/math] & [math] P(y=1|x_1.....x_n) \gt P(y=0|x_1.....x_n) [/math]

[math] \widehat{y}=1 [/math] & [math] \Leftrightarrow [/math] & [math] \frac{P(y=1|x_1.....x_n)}{P(y=0|x_1.....x_n)} \gt 1 [/math]

& [math]\Leftrightarrow[/math] & [math] \log{\frac{P(y=1)}{P(y=0)}} + \sum_{i=1..n}{\log{\frac{P(x_i|y=1)}{P(x_i|y=0)}}}\gt 0 [/math]
\end{tabular}

Now if we define the following:

[math]P(y=1) =p[/math]

[math]P(x_i|y=1)=P_{i1}[/math]

[math]P(x_i|y=0)=P_{i0}[/math]

We can continue with the above simplification and we arrive at the solution:

\begin{tabular}{ ccc }
[math]\widehat{y}=1[/math] & [math]\Leftrightarrow[/math] & [math]x_i\log{\frac{P_{i1}}{P_{i0}}}+ (1-x_i)\log{\frac{(1-P_{i1})}{(1-P_{i0})}} \gt 0[/math]

& [math]\Leftrightarrow[/math] & [math] =x_i\underbrace{\log{\frac{P_{i1}(1-P_{i0})}{P_{i0}(1-P_{i1})}}}_{slope} + \underbrace{ \log{\frac{(1-P_{i1})}{(1-P_{i0})}} }_{intercept} [/math]
\end{tabular}