The size distribution for Markov equivalence classes of acyclic digraph models

Bayesian networks, equivalently graphical Markov models determined by acyclic digraphs or ADGs (also called directed acyclic graphs or dags), have proved to be both effective and efficient for representing complex multivariate dependence structures in terms of local relations. However, model search and selection is potentially complicated by the many-to-one correspondence between ADGs and the statistical models that they represent. If the ADGs/models ratio is large, search procedures based on unique graphical representations of equivalence classes of ADGs could provide substantial computational efficiency. Hitherto, the value of the ADGs/models ratio has been calculated only for graphs with n = 5 or fewer vertices. In the present study, a computer program was written to enumerate the equivalence classes of ADG models and study the distributions of class sizes and number of edges for graphs up to n = 10 vertices. The ratio of ADGs to numbers of classes appears to approach an asymptote of about 3.7. Distributions of the classes according to number of edges and class size were produced which also appear to be approaching asymptotic limits. Imposing a bound on the maximum number of parents to any vertex causes little change if the bound is sufficiently large, with four being a possible minimum. The program also includes a new variation of orderly algorithm for generating undirected graphs. (C) 2002 Elsevier Science B.V. All rights reserved.