Chow–Liu tree

A first-order dependency tree representing the product on the left.

In probability theory and statistics Chow–Liu tree is an efficient method for constructing a second-order product approximation of a joint probability distribution, first described in a paper by (Chow Liu). The goals of such a decomposition, as with such Bayesian networks in general, may be either data compression or inference.

The Chow–Liu representation

The Chow–Liu method describes a joint probability distribution [math]\displaystyle{ P(X_{1},X_{2},\ldots,X_{n}) }[/math] as a product of second-order conditional and marginal distributions. For example, the six-dimensional distribution [math]\displaystyle{ P(X_{1},X_{2},X_{3},X_{4},X_{5},X_{6}) }[/math] might be approximated as

[math]\displaystyle{ P^{\prime }(X_{1},X_{2},X_{3},X_{4},X_{5},X_{6})=P(X_{6}|X_{5})P(X_{5}|X_{2})P(X_{4}|X_{2})P(X_{3}|X_{2})P(X_{2}|X_{1})P(X_{1}) }[/math]

where each new term in the product introduces just one new variable, and the product can be represented as a first-order dependency tree, as shown in the figure. The Chow–Liu algorithm (below) determines which conditional probabilities are to be used in the product approximation. In general, unless there are no third-order or higher-order interactions, the Chow–Liu approximation is indeed an approximation, and cannot capture the complete structure of the original distribution. (Pearl 1988) provides a modern analysis of the Chow–Liu tree as a Bayesian network.

The Chow–Liu algorithm

Chow and Liu show how to select second-order terms for the product approximation so that, among all such second-order approximations (first-order dependency trees), the constructed approximation [math]\displaystyle{ P^{\prime} }[/math] has the minimum Kullback–Leibler divergence to the actual distribution [math]\displaystyle{ P }[/math], and is thus the closest approximation in the classical information-theoretic sense. The Kullback–Leibler divergence between a second-order product approximation and the actual distribution is shown to be

[math]\displaystyle{ D(P\parallel P^{\prime })=-\sum I(X_{i};X_{j(i)})+\sum H(X_{i})-H(X_{1},X_{2},\ldots ,X_{n}) }[/math]

where [math]\displaystyle{ I(X_{i};X_{j(i)}) }[/math] is the mutual information between variable [math]\displaystyle{ X_{i} }[/math] and its parent [math]\displaystyle{ X_{j(i)} }[/math] and [math]\displaystyle{ H(X_{1},X_{2},\ldots ,X_{n}) }[/math] is the joint entropy of variable set [math]\displaystyle{ \{X_{1},X_{2},\ldots ,X_{n}\} }[/math]. Since the terms [math]\displaystyle{ \sum H(X_{i}) }[/math] and [math]\displaystyle{ H(X_{1},X_{2},\ldots ,X_{n}) }[/math] are independent of the dependency ordering in the tree, only the sum of the pairwise mutual informations, [math]\displaystyle{ \sum I(X_{i};X_{j(i)}) }[/math], determines the quality of the approximation. Thus, if every branch (edge) on the tree is given a weight corresponding to the mutual information between the variables at its vertices, then the tree which provides the optimal second-order approximation to the target distribution is just the maximum-weight tree. The equation above also highlights the role of the dependencies in the approximation: When no dependencies exist, and the first term in the equation is absent, we have only an approximation based on first-order marginals, and the distance between the approximation and the true distribution is due to the redundancies that are not accounted for when the variables are treated as independent. As we specify second-order dependencies, we begin to capture some of that structure and reduce the distance between the two distributions.

Chow and Liu provide a simple algorithm for constructing the optimal tree; at each stage of the procedure the algorithm simply adds the maximum mutual information pair to the tree. See the original paper, (Chow Liu), for full details. A more efficient tree construction algorithm for the common case of sparse data was outlined in (Meilă 1999).

Chow and Wagner proved in a later paper (Chow Wagner) that the learning of the Chow–Liu tree is consistent given samples (or observations) drawn i.i.d. from a tree-structured distribution. In other words, the probability of learning an incorrect tree decays to zero as the number of samples tends to infinity. The main idea in the proof is the continuity of the mutual information in the pairwise marginal distribution. More recently, the exponential rate of convergence of the error probability was provided.

Variations on Chow–Liu trees

The obvious problem which occurs when the actual distribution is not in fact a second-order dependency tree can still in some cases be addressed by fusing or aggregating together densely connected subsets of variables to obtain a "large-node" Chow–Liu tree (Huang King), or by extending the idea of greedy maximum branch weight selection to non-tree (multiple parent) structures (Williamson 2000). (Similar techniques of variable substitution and construction are common in the Bayes network literature, e.g., for dealing with loops. See (Pearl 1988).)

Generalizations of the Chow–Liu tree are the so-called t-cherry junction trees. It is proved that the t-cherry junction trees provide a better or at least as good approximation for a discrete multivariate probability distribution as the Chow–Liu tree gives. For the third order t-cherry junction tree see (Kovács Szántai), for the kth-order t-cherry junction tree see (Szántai Kovács). The second order t-cherry junction tree is in fact the Chow–Liu tree.

See also

• Bayesian network
• Knowledge representation

References

• Chow, C. K.; Liu, C.N. (1968), "Approximating discrete probability distributions with dependence trees", IEEE Transactions on Information Theory IT-14 (3): 462–467, doi:10.1109/tit.1968.1054142 .
• Huang, Kaizhu; King, Irwin; Lyu, Michael R. (2002), "Constructing a large node Chow–Liu tree based on frequent itemsets", in Wang, Lipo; Rajapakse, Jagath C.; Fukushima, Kunihiko et al., Proceedings of the 9th International Conference on Neural Information Processing ({ICONIP}'02), Singapore, pp. 498–502 .
• Pearl, Judea (1988), Probabilistic Reasoning in Intelligent Systems: Networks of Plausible Inference, San Mateo, CA: Morgan Kaufmann 
• Williamson, Jon (2000), "Approximating discrete probability distributions with Bayesian networks", Proceedings of the International Conference on Artificial Intelligence in Science and Technology, Tasmania, pp. 16–20 .
• Meilă, Marina (1999), "An Accelerated Chow and Liu Algorithm: Fitting Tree Distributions to High-Dimensional Sparse Data", Proceedings of the Sixteenth International Conference on Machine Learning, Morgan Kaufmann, pp. 249–257 .
• Chow, C. K.; Wagner, T. (1973), "Consistency of an estimate of tree-dependent probability distribution", IEEE Transactions on Information Theory IT-19 (3): 369–371, doi:10.1109/tit.1973.1055013 .
• Kovács, E.; Szántai, T. (2010), "On the Approximation of a Discrete Multivariate Probability Distribution Using the New Concept of t-Cherry Junction Tree", Coping with Uncertainty, Lecture Notes in Economics and Mathematical Systems, 633, Part 1, pp. 39–56, doi:10.1007/978-3-642-03735-1_3, ISBN 978-3-642-03734-4 .
• Szántai, T.; Kovács, E. (2010), "Hypergraphs as a mean of discovering the dependence structure of a discrete multivariate probability distribution", Annals of Operations Research .