The junction tree algorithm (also known as 'Clique Tree') is a method used in machine learning to extract marginalization in general graphs. In essence, it entails performing belief propagation on a modified graph called a junction tree. The graph is called a tree because it branches into different sections of data; nodes of variables are the branches.[1] The basic premise is to eliminate cycles by clustering them into single nodes. Multiple extensive classes of queries can be compiled at the same time into larger structures of data.[1] There are different algorithms to meet specific needs and for what needs to be calculated. Inference algorithms gather new developments in the data and calculate it based on the new information provided.[2]
Construct a junction tree from the triangulated graph (we will call the vertices of the junction tree "supernodes").
Propagate the probabilities along the junction tree (via belief propagation)
Note that this last step is inefficient for graphs of large treewidth. Computing the messages to pass between supernodes involves doing exact marginalization over the variables in both supernodes. Performing this algorithm for a graph with treewidth k will thus have at least one computation which takes time exponential in k. It is a message passing algorithm.[3] The Hugin algorithm takes fewer computations to find a solution compared to Shafer-Shenoy.
The Shafer-Shenoy algorithm is the sum product of a junction tree.[6] It is used because it runs programs and queries more efficiently than the Hugin algorithm. The algorithm makes calculations for conditionals for belief functions possible.[7]Joint distributions are needed to make local computations happen.[7]
The first step concerns only Bayesian networks, and is a procedure to turn a directed graph into an undirected one. We do this because it allows for the universal applicability of the algorithm, regardless of direction.
The second step is setting variables to their observed value. This is usually needed when we want to calculate conditional probabilities, so we fix the value of the random variables we condition on. Those variables are also said to be clamped to their particular value.
Example of a chordal graph
The third step is to ensure that graphs are made chordal if they aren't already chordal. This is the first essential step of the algorithm. It makes use of the following theorem:[8]
Theorem: For an undirected graph, G, the following properties are equivalent:
Graph G is triangulated.
The clique graph of G has a junction tree.
There is an elimination ordering for G that does not lead to any added edges.
Thus, by triangulating a graph, we make sure that the corresponding junction tree exists. A usual way to do this, is to decide an elimination order for its nodes, and then run the Variable elimination algorithm. The variable elimination algorithm states that the algorithm must be run each time there is a different query.[1] This will result to adding more edges to the initial graph, in such a way that the output will be a chordal graph.
All chordal graphs have a junction tree.[4] The next step is to construct the junction tree. To do so, we use the graph from the previous step, and form its corresponding clique graph.[9] Now the next theorem gives us a way to find a junction tree:[8]
Theorem: Given a triangulated graph, weight the edges of the clique graph by their cardinality, |A∩B|, of the intersection of the adjacent cliques A and B. Then any maximum-weight spanning tree of the clique graph is a junction tree.
So, to construct a junction tree we just have to extract a maximum weight spanning tree out of the clique graph. This can be efficiently done by, for example, modifying Kruskal's algorithm.
The last step is to apply belief propagation to the obtained junction tree.[10]
Usage: A junction tree graph is used to visualize the probabilities of the problem. The tree can become a binary tree to form the actual building of the tree.[11] A specific use could be found in auto encoders, which combine the graph and a passing network on a large scale automatically.[12]
^Ramirez, Julio C.; Munoz, Guillermina; Gutierrez, Ludivina (September 2009). "Fault Diagnosis in an Industrial Process Using Bayesian Networks: Application of the Junction Tree Algorithm". 2009 Electronics, Robotics and Automotive Mechanics Conference (CERMA). IEEE. doi:10.1109/cerma.2009.28.
Lepar, V., Shenoy, P. (1998). "A Comparison of Lauritzen-Spiegelhalter, Hugin, and Shenoy-Shafer Architectures for Computing Marginals of Probability Distributions." https://arxiv.org/ftp/arxiv/papers/1301/1301.7394.pdf
