arxiv: v1 [math.st] 6 Jun 2016

Transcription

1 GENERALIZED PERMUTOHEDRA FROM PROBABILISTIC GRAPHICAL MODELS FATEMEH MOHAMMADI, CAROLINE UHLER, CHARLES WANG, AND JOSEPHINE YU arxiv:606.08v [math.st] 6 Jun 206 Abstract. A graphical model encodes conditional independence relations via the Markov properties. For an undirected graph these conditional independence relations can be represented by a simple polytope known as the graph associahedron, which can be constructed as a Minkowski sum of standard simplices. We show that there is an analogous polytope for conditional independence relations coming from a regular Gaussian model, and it can be defined using multiinformation or relative entropy. For directed acyclic graphical models we give a construction of this polytope as a Minkowski sum of matroid polytopes. Finally, we apply this geometric insight to construct a new ordering-based search algorithm for causal inference via directed acyclic graphical models.. Introduction A graphical model encodes conditional independence (CI) relations via the Markov properties. Our main goal is to understand the polyhedral geometry and combinatorics of the collection of CI relations encoded by a directed acyclic graph (DAG), a directed graph without directed cycles. It is natural, especially in view of causal inference, to associate to each conditional independence statement a collection of pairs of adjacent permutations of random variables that are compatible with that statement. Each of these pairs can be viewed as an edge of a permutohedron or a wall in the S n fan, which is the normal fan of the permutohedron. Removing these walls gives a coarsening of the fan and a natural question is whether this fan is the normal fan of a polytope. For undirected graphical models, the theory is well-understood. The coarsening of the S n fan corresponding to the CI relations encoded by an undirected graph is the normal fan of a polytope called a graph associahedron [MPS + 09]. These polytopes are Minkowski sums of standard simplices, and their facial structure has a nice description via tubings [CD06, PRW08]. In this paper we will show that the coarsened S n fan of any DAG is the normal fan of a polytope, which we call a DAG associahedron. We give two concrete constructions of DAG associahedra, one using multiinformation, or relative entropy, and another using matroids. We also show that in general DAG associahedra are not simple polytopes and cannot be realized as a Minkowski sum of standard simplices. Our main motivation for studying DAG associahedra is causal inference: Given a set of CI relations that are inferred from data, the goal is to estimate the underlying DAG model. A DAG is defined by an ordering of the nodes and an undirected graph. We show how our geometric insight on DAG associahedra can be applied to construct a new ordering-based search algorithm for causal inference. 2. Notation and background In this section, we discuss the relationship between CI relations, the S n fan, and generalized permutohedra. Please refer to the Appendix for a dictionary. Keywords: Graphical model, graphoid, permutohedron, causal inference, submodular function, matroid, entropy. MSC(200): 62H05 (primary); 52B2, 52B0 (secondary).

2 2 FATEMEH MOHAMMADI, CAROLINE UHLER, CHARLES WANG, AND JOSEPHINE YU Let [n] = {,..., n} and let P be a joint distribution on the random variables X i, i [n]. For I [n] we identify {X i : i I} with the index set I. For disjoint subsets I, J, K [n] we say that I is conditionally independent of J given K under P if the conditional probability P(A J, K) does not depend on J for any measurable set A in the sample space of X I. This statement is denoted by I P J K or simply I J K. If K =, we write I J. The set of CI relations arising from a distribution satisfies the following basic implications, known as the semigraphoid properties [Pea88]: (SG ) if I J L then J I L, (SG2 ) if I J L and U I, then U J L, (SG ) if I J L and U I, then I J (U L), (SG ) if I J L and I K J L, then I (J K) L. In this paper, CI relations can be considered as a formal construct and do not necessarily need any probabilistic interpretation. In addition, we will only work with relations in which I and J are both singletons, denoted by lower-case letters i, j. To simplify notation, we use concatenation to denote union among subsets and elements of [n], e.g. Lij means L {i, j}. Then a semigraphoid is a set of CI relations that satisfy the CI implications (SG) if i j L then j i L, (SG2) if i j L and i k jl, then i k L and i j kl, for i, j, k [n] distinct and L [n] \ {i, j, k}. For distributions with strictly positive densities such as the Gaussian distribution, the intersection axiom holds in addition to the semigraphoid axioms, namely (INT) if i j kl and i k jl, then i j L and i k L. The implications (SG), (SG2) and (INT) together are known as the graphoid properties. Note that these implications are not a complete list of CI implications that hold for distributions. In fact, Studený [Stu95] proved that there exists no finite such characterization. In [LM07], Lněničika and Matúš defined gaussoids as the graphoids satisfying the following additional axioms: (G) if i j L and i k L, then i j kl and i k jl, (G2) if i j L and i j kl, then i k L or j k L. The property (G) is the converse of the intersection axiom, and (G2) is called weak transitivity. The CI relations of any regular Gaussian distribution form a gaussoid, but not all gaussoids arise this way. The set of CI relations coming from probabilistic graphical models that we study in this paper can be faithfully represented by regular Gaussian distributions, so they are gaussoids. We will associate a geometric object to a collection of CI relations as follows: Consider the hyperplanes in R n defined by equations of the form x i = x j for all i < j n. The complement of these hyperplanes consists of points in R n with distinct coordinates, and they are partitioned into n! connected components corresponding to the permutations of [n] as follows: We identify a permutation (bijection) π on [n] with the linear order π() π(2) π(n). To every vector u R n with distinct coordinates, we associate a linear order on [n] by defining i j if and only if u i > u j. For example, the vector u = (25,, 6, 9) gives the linear order 2. We also denote this by the descent vector of the form ( 2). The closures of the n! cones and all their faces form a fan, which we will call the S n fan. It is also known as the A n fan or the braid arrangement fan. Each cone in the fan contains the line in direction (,,..., ) and is generated by a collection of 0/ vectors, every pair of which is nested (when a 0/ vector is identified with the set of coordinates that are ).

3 GENERALIZED PERMUTOHEDRA FROM PROBABILISTIC GRAPHICAL MODELS To each CI relation i j K, where i, j [n] distinct and K [n] \ {i, j}, we associate pairs of adjacent permutations of the form () (a a k i j b b n k 2 ) and (a a k j i b b n k 2 ), where {a,..., a k } = K and {b,..., b n k 2 } = [n]\(k {i, j}). We will denote such a pair by (a a k i j b b n k 2 ). For each relation i j K there are K! (n K 2)! such pairs. A fan F in R n is said to be a coarsening of the S n fan if every cone in the S n fan is contained in a cone of F, or equivalently, if every cone of F is a union of some cones of the S n fan. In particular, maximal cones of F are unions of maximal cones of the S n fan, and we can think of constructing F from the S n fan by removing certain walls. This gives an equivalence relation on S n two permutations are equivalent if and only if their corresponding cones in the S n fan are contained in the same cone in F. Such an equivalence relation coming from a fan is called a convex rank test in [MPS + 09]. We identify a coarsening of the S n fan with the collection of walls that are removed. Each wall corresponds to an adjacent pair of permutations as in (), which gives a CI relation i j {a,..., a k }. It was shown in [MPS + 09, Theorem 6] that a set of walls form the missing walls in a fan that coarsens the S n fan if and only if the corresponding set of CI relations form a semigraphoid. In particular, if the wall associated to the pair () is not a wall in a coarsened S n fan F, then any pair obtained by permuting the a s and b s is also not a wall in F. A complete fan F in R n, the union of whose cones is all of R n, is called polytopal if it is the normal fan of a polytope. The S n fan itself is polytopal since it is the normal fan of a polytope called the (standard) permutohedron P n which is the convex hull of the n! permutations of [n] in R n. Two vertices of P n form an edge if and only if their descent vectors differ by an adjacent transposition as in (). Thus each CI relation corresponds to a certain set of edges of P n. A generalized permutohedron is a polytope whose normal fan is a coarsening of the S n fan, or equivalently, whose edges are in direction e i e j. See Figures and 2 for some examples. These polytopes are also called M-convex polytopes or base polyhedra [Mur0, (.)], and their projections along a coordinate direction give generalized polymatroids [Fuj05, Theorem.58]. We use the term generalized permutohedron to highlight the connection to permutations. Example 2. (Undirected graphical models and graph associahedra). Let G be an undirected graph with node set [n]. We associate a random variable X i to each node i of the graph. The 2 (0,, 0) (0,, ) (,, 0) (0, 0, ) (, 0, 0) (a) The S fan modulo the line (,, ). The maximal cones are labeled with permutations and the walls are labeled with CI relations. (, 0, ) (b) Permutohedron P with outer normals of its facets. Figure. The permutohedron P and its normal fan S.

4 FATEMEH MOHAMMADI, CAROLINE UHLER, CHARLES WANG, AND JOSEPHINE YU (a) The standard simplex conv{e, e 2, e }. (b) The matroid polytope of U 2,, conv{e + e 2, e + e, e 2 + e }. (c) The DAG associahedron of 2 with CI relation 2. This polytope is not a Minkowski sum of scaled standard simplices. Figure 2. Some generalized permutohedra. Compare with the fan in Figure. joint distribution P of the random vector X = (X,..., X n ) satisfies the undirected (global) Markov property with respect to G if I J K for all disjoint subsets I, J, K [n] such that K separates I and J in G, i.e. every path between nodes i I and j J passes through a node k K. If a distribution P satisfies exactly the CI relations corresponding to separations in the graph G, then P is called faithful or perfectly Markovian with respect to G. For any undirected graph there exist faithful regular Gaussian distributions; see [Lau96, Chapter ] for more details. Hence for any undirected graph G the corresponding CI relations defined by the Markov property satisfy the gaussoid axioms. The coarsened S n fan associated to the gaussoid of an undirected graph is the normal fan of a polytope, which can be realized as the Minkowski sum of standard simplices I = conv{e i : i I} where I runs over all sets of nodes that induce connected subgraphs of G [MPS + 09]. These polytopes are called graph associahedra and were studied in [Dev09, CD06, PRW08]. We now give a characterization of coarsened S n fans that are polytopal. We denote by 2 [n] the power set of [n], the set of all subsets of [n]. Definition 2.2. A function ω : 2 [n] R is called submodular if (2) ω(ki) + ω(kj) ω(kij) + ω(k) for all K [n] and i, j [n]\k. A semigraphoid on [n] is called submodular if there is a submodular function ω on 2 [n] with ω( ) = 0 such that equality in (2) is achieved if and only if the relation i j K is in the semigraphoid. It can be shown that a submodular function also satisfies ω(a) + ω(b) ω(a B) + ω(a B) for all A, B [n]. Note that a submodular function on 2 [n] is an L-convex function on the unit cube {0, } n [Mur0]. The following result follows from the conjugacy between L- and M- convex functions and also from [Mur0, Theorem.5] which was well known, but neither precise statement nor proof can be found in the literature [Mur0, Chapter Bibliographical Notes]. A part of it appeared in [MPS + 09, Proposition 2 and Theorem ]. Lemma 2.. A coarsening F of the S n fan is polytopal if and only if the corresponding semigraphoid is submodular. If ω : 2 [n] R is a submodular function, then the fan F corresponding to the semigraphoid defined by ω is the outer normal fan of the polytope defined by () x i ω(i) for each nonempty I [n], and x i = ω([n]). i [n] i I

5 GENERALIZED PERMUTOHEDRA FROM PROBABILISTIC GRAPHICAL MODELS 5 Remark 2.. If ω is a submodular function on 2 [n] with ω( ) = 0, then ω : 2 [n] R defined as ω (S) = ω([n]\s) ω([n]) is also submodular with ω ( ) = 0. The polytopes P and P, defined by ω and ω as in (), are related by P = P. Proof of Lemma 2.. We first explain how a halfspace description of a polytope gives a regular subdivision of the normal vectors of the halfspaces. See [DLRS0, 2.5] for details. Let P = {x R n : Ax + b 0}, where A is a k n matrix whose rows positively span R n and b R k is a column vector. Suppose the polytope P is non-empty and that all inequalities are tight but possibly redundant. Let C be the cone in R n R generated by be the rows of the concatenated matrix [A b]. The row (a i, b i ) of [A b] is called a lift of the vector a i. The dual cone C :={u : u T v 0 v C} is generated by {(x, ) : x P }. Thus the projection of faces of C onto R n form the inner normal fan of P. Since all inequalities are assumed to be tight, all lifted vectors (a i, b i ) lie on the boundary of C. All vectors in the dual cone C have positive last coordinates, so all proper faces of C are on the lower hull of C. In particular, if a vector (a, b) lies on the boundary of C, then (a, b + ε) does not lie on the boundary of C for any ε > 0. Now let F be a coarsening of the S n fan in R n. Every cone in F contains a line in direction (,,..., ) and is generated by this line together with some 0/ vectors. As shown above, the fan F is polytopal if and only if there exists a lift ω on the set of rays {e I I [n]} { e [n] } such that the faces of the cone C spanned by the lifted rays project precisely onto the cones of F and every lifted ray is on the boundary of C. Since all cones in F contain the line (,,..., ), it suffices to consider lifts ω such that ω(e [n] ) = ω( e [n] ). Such lifts can be identified with functions on 2 [n] with value 0 on. We wish to show submodularity of ω. For any I, J [n], the vectors e I, e I J, e I J lie in a common cone in the S n fan. Since F coarsens the S n fan, they also lie in a common cone in F. Similarly e J, e I J, e I J lie in a common cone of F. First, consider the case when e I and e J are lifted to the same proper face of C. Then this cone also contains e I J and e I J. Since we assumed that all lifted vectors lie on the boundary, hence a proper face, of C, and ω is linear on this face, we must have that ω(e I )+ω(e J ) = ω(e I J )+ω(e I J ). Now suppose that e I and e J are not lifted to the same proper face of C. Then ω is not linear on the vectors e I, e J, e I J, and e I J. We must then have that ω(e I ) + ω(e J ) > ω(e I J ) + ω(e I J ), because ω(e I ) + ω(e J ) < ω(e I J ) + ω(e I J ) would imply that (e I J + e I J, ω(e I J ) + ω(e I J )) > (e I + e J, ω(e I ) + ω(e J )), contradicting the fact that e I J and e I J are lifted to the same cone in the lower hull of C. For the converse, suppose ω is a submodular function on 2 [n] with ω( ) = 0 and consider the lift of e I to ω(i) for each I [n] and e [n] to ω([n]). Let F be the projection of the lower hull of the lifted cone C. The submodularity inequality ω(e I ) + ω(e J ) ω(e I J ) + ω(e I J ) ensures that whenever I and J are lifted to the same cone in the lower hull of C, then so are I J and I J. In other words, whenever a cone of F contains both, e I and e J, then it must also contain both, e I J and e I J, showing that F is a coarsening of the S n fan. Now suppose that the coarsened S n fan F is polytopal defined by a submodular function ω as above. The wall corresponding to the pair of adjacent permutations in () is not a cone in the fan F if and only if the two adjacent maximal cones are contained in the same cone of F. In particular, this happens if and only if e Ki and e Kj are in the same cone where K = {a,..., a k }. This is equivalent to the condition that (2) is attained at equality. Let P be the polytope defined by (). Its inner normal fan is obtained by lifting the rays e I to height ω(i) for nonempty I [n] and e [n] to height ω([n]). This is the negation of the fan F, which is obtained by lifting e I to ω(i) and e [n] to ω([n]). This shows that F is the outer normal fan of P.

6 6 FATEMEH MOHAMMADI, CAROLINE UHLER, CHARLES WANG, AND JOSEPHINE YU Example 2.5. Consider the submodular function ω on 2 [n] whose value is on all non-empty sets and 0 on the empty set. This is the rank function of the uniform rank one matriod on [n]. The generalized permutohedron defined by this submodular function is a standard simplex of dimension n whose outer normal vectors are e I for subsets I of size n. Any set of n 2 facet normals spans a wall in the normal fan, with pairs of the form (), where K =, corresponding to CI relations of the form i j. See Figures and 2. This characterization leads to the following questions for any given semigraphoid: Question A. Is a given semigraphoid submodular? And if so, can we construct a submodular function with the desired equalities as in Definition 2.2? In the following sections we will give a positive answer to these questions for semigraphoids coming from DAG models. Rank functions of matroids are submodular functions, so every matroid M on the ground set [n] gives a semigraphoid on [n] as follows: i j K rank(ki) + rank(kj) > rank(kij) + rank(k) Note that since a matroid rank function takes integer values and rank(aa) rank(a) + for any A [n] and a [n], we obtain () i j K rank(k) + = rank(ki) = rank(kj) = rank(kij). In this case, the coarsening of the S n fan is the outer normal fan of the matroid polytope, which is defined as the convex hull of the indicator functions of the bases of the matroid. For example, the standard simplex I = conv{e i : i I} is the matroid polytope of the rank one matroid in which each element of I forms a base. The intersection of two semigraphoids is again a semigraphoid. The corresponding operations on fans, polytopes, and submodular functions are common refinement, Minkowski sum, and sum, respectively. Question B. Which semigraphoids are submodular with respect to a sum of rank functions of matroids? Which fans are normal fans of Minkowski sums of matroid polytopes? For example, the Minkowski sum of all standard simplices I, for all non-empty subsets I [n], is affinely equivalent to the permutohedron P n, i.e. they have the same normal fan, which is the entire S n fan. This decomposition is not unique, however, e.g. P is a hexagon and can be decomposed as the Minkowski sum of either two triangles or three line segments, all of which are matroid polytopes. Example 2.6. Let G be the following DAG. 2 We will see in the next section that the Markov property on G defines a single CI relation, namely 2. Removing the corresponding wall in the S fan gives a fan with 5 maximal cones. Figure 2(c) depicts a polytope with this normal fan. It is straightforward to check that this fan is not the normal fan of a Minkowski sum of standard simplices, but it is the normal fan of the Minkowski sum of the simplex in Figure 2(b) with two additional line segments, which are all matroid polytopes.

7 GENERALIZED PERMUTOHEDRA FROM PROBABILISTIC GRAPHICAL MODELS 7. Bayesian networks Similarly to undirected graphs we can define probabilistic models on DAGs. models are also known as Bayesian networks. Such graphical Let G be a DAG with nodes [n]. If there is a directed edge from i to j in G, which we denote by i j in G or (i, j) G, the node i is called a parent of the node j. The set of all parent nodes of j is denoted by pa(j). We now review the concept of separation for DAGs. A path in G is an alternating sequence of nodes and edges, starting and ending at nodes, in which each edge is adjacent in the sequence to its two endpoints. The path may contain repeated edges and nodes. We do not assume that the direction of the edges is compatible with the ordering of the nodes in the path. Definition.. Let G be a DAG on [n] and let i, j [n] and K [n] \ {i, j}. A Bayes ball path from i to j given K in G is a path from i to j in G such that () if a b c or a b c or a b c is on the path, then b / K; (2) if a b c is on the path, then b K (where a and c need not be distinct). In this case the node b is called a collider along the path. See Figures and 5 for examples. Informally we think of a directed edge i j as pointing down from i to j. A Bayes ball rolls along edges of the DAG. It cannot roll through nodes that are in K, but it can bounce off them by going down, touching K, then going back up either along the same or a different edge. For subsets of nodes I, J, K [n], we say that I and J are directed separated or d-separated by K in G if there is no Bayes-ball path from any element of I to any element of J given K [VP90]. This led to the construction of the Bayes-Ball algorithm [Sha98], an algorithm for determining d-separation statements. Similarly as for undirected graphs, we can also associate a random vector with joint distribution P to the nodes of a DAG G. Then P satisfies the directed (global) Markov property with respect to G if I J K for all disjoint subsets I, J, K V such that K d-separates I and J in G. A faithful distribution to G, i.e. a distribution that satisfies exactly the CI relations corresponding to d-separation in G, can be realized by regular Gaussian distributions (see ). Hence, for any DAG G the CI relations of the form i j K, where i and j are d-separated given K in G, form a gaussoid, which we call a DAG gaussoid. It is important to note that while the set of separation statements uniquely determines an undirected graph, this is not the case for d-separation statements for DAGs. Two DAGs are called Markov equivalent if they satisfy the same d-separation statements. The Markov equivalence class is determined by the skeleton of a DAG and its V-structures triples of nodes (i, j, k) such that i k j and i, j are not adjacent [AMP97]. An essential graph [AMP97] (also called a completed partially directed acyclic graph or CPDAG in [Chi02] and a maximally oriented graph in [Mee95]) is a graph with undirected and directed edges that uniquely represents a Markov equivalence class of DAGs. It has the same skeleton as the DAGs in the Markov equivalence class and contains a directed edge i j if and only if each DAG in the Markov equivalence class contains the directed edge i j. The following is our main result and answers Questions A and B for DAG gaussoids. Theorem.2 (Main Theorem). Every DAG gaussoid is submodular. Equivalently, the associated coarsening of the S n fan is the normal fan of a polytope. Moreover, there is a realization of this polytope as a Minkowski sum of matroid polytopes. The equivalence of the first two statements was proven in Lemma 2. above. We call any such polytope resulting from a DAG gaussoid a DAG associahedron. DAG associahedra are uniquely

8 8 FATEMEH MOHAMMADI, CAROLINE UHLER, CHARLES WANG, AND JOSEPHINE YU 2 2 Figure. The DAG G (left) and its moral graph G (right) discussed in Example.. defined up to equivalence of normal fans, and they only depend on the DAG up to Markov equivalence. Remark.. Let G be a DAG. The normal fan of the DAG associahedron corresponding to G can be obtained by coarsening the normal fan of the graph associahedron corresponding to the moral graph of G the undirected graph with edges (i, j) if i j in G, j i in G, or i k j for some k in G; see Figure. In the next two sections, we will give two independent proofs for the submodularity of DAG gaussoids. In the first proof, in, we use multiinformation, or relative entropy, to give a formula for the submodular function and hence a realization of DAG associahedra. However, in general the constant terms of the inequalities in this construction are not rational. We will discuss some heuristic methods for finding exact combinatorial information from approximate inequalities. In the second proof, in 5, we give a realization of DAG associahedra as Minkowski sums of matroid polytopes, which are integral polytopes. The submodularity of a semigraphoid can be tested using linear programming [HMS + 08]. So our theorem states that the linear programs coming from DAG gaussoids are always feasible, and our proofs give an explicit construction of a feasible solution. We illustrate the concepts introduced so far by an example of a DAG model on nodes and describe the corresponding DAG associahedron. Example.. Consider the DAG G shown in Figure. An example of a Bayes ball path in G is the path from node to 2 given K = {}, since on the path 2 the node / K but K. The DAG gaussoid corresponding to G consists of the CI relations 2,, 2, {2, }, 2 {, }. The corresponding edges of the permutohedron are shown in green and blue in Figure (a). Since these CI relations form a semigraphoid, we obtain a coarsening of the S n fan by removing the edges (2 ), (2 ), ( 2), ( 2 ), (2 ), ( 2 ), ( 2) and ( 2). The resulting coarsening of the S n fan obtained by contracting the colored edges in the permutohedron is polytopal. The convex polytope corresponding to this DAG associahedron is shown in Figure (c). The moral graph G of G is shown in Figure (right). The gaussoid corresponding to G consists of the CI relations, 2, {2, }, 2 {, }. In general any DAG gaussoid contains the gaussoid of its moral graph. The edges corresponding to the CI relations for the moral graph are shown in green in Figure (a). By contracting the green edges in the permutohedron we obtain the graph associahedron corresponding to G shown in Figure (b). By further contracting also the blue edges, we obtain the DAG associahedron corresponding to G. As we will see in Proposition.6, the DAG associahedron in this example cannot be realized as a Minkowski sum of simplices. However, we will show in 5 that it can be realized as the following

9 GENERALIZED PERMUTOHEDRA FROM PROBABILISTIC GRAPHICAL MODELS 9 Minkowski sum of matroid polytopes: conv{e2, e, e2 } + conv{e2 + e + e2 + e + e2 }. As we will see in 5, the first three polytopes in the sum correspond to the three edges in the DAG, the next two correspond to the paths and 2, which have no colliders, and the last two correspond to the paths 2 (given ) and 2 (given ) respectively. We end this section with two observations about DAG associahedra. In the following example we show that unlike graph associahedra, DAG associahedra need not be simple. Example.5 (A non-simple DAG associahedron). Let G be the following DAG. 2 The corresponding DAG gaussoid consists only of the CI relation 2, which corresponds to a single edge 2 on the permutohedron P. Contracting this edge gives a vertex adjacent to edges on a -dimensional polyhedron, so the resulting polytope is not simple. In this case, the combinatorial operation of contracting the edge can be realized geometrically by pushing the two neighboring square facets towards each other until they meet at a point. Furthermore, as already mentioned in Example., unlike graph associahedra, DAG associahedra need not be Minkowski sums of standard simplices (MSS) in the sense of [MPS+ 09]. In fact, the following result shows that a DAG associahedron can be realized as a MSS if and only if the DAG gaussoid equals the gaussoid of its moral graph, or in other words, if and only if the DAG model coincides with an undirected graphical model. Proposition.6. The DAG associahedron associated to a DAG G is MSS if and only if G does not contain any V-structures, i.e. the DAG model coincides with an undirected graphical model. (a) Permutohedron P. The green edges correspond to CI relations in the moral graph G in Example.. The blue edges correspond to the additional CI relations in G. (b) The graph associahedron of the moral graph G obtained by contracting the green edges. (c) DAG associahedron for G obtained by contracting both, green and blue edges. Figure. The vertices are labeled by descent vectors of permutations, with s removed. The figures show how to combinatorics of the polytope changes as edges are contracted, but they are not drawn to be geometrically correct.

10 0 FATEMEH MOHAMMADI, CAROLINE UHLER, CHARLES WANG, AND JOSEPHINE YU We saw in Example 2.6 that a V-structure cannot be MSS. We can generalize this example to the following corollary of [MPS + 09, Proposition 20]: Lemma.7. If a semigraphoid arises from a Minkowski sum of standard simplices, then for any i, j [n] distinct and K K [n] \ {i, j}, we have i j K = i j K. Proof. For I [n], the standard simplex I is the matroid polytope of the rank one matroid on [n] whose loops are [n]\i. By () the semigraphoid corresponding to I contains i j K where i, j I and K I =. Taking a Minkowski sum of such simplices corresponds to taking the union of the associated conditional dependence statements. It follows that i j K = i j K for all K K [n] \ {i, j}. Using Lemma.7 we can now easily prove Proposition.6. Proof of Proposition.6. If G does not contain any V-structures, then the corresponding DAG gaussoid is equivalent to the gaussoid obtained from an undirected graph, namely the skeleton of G, so it is MSS. On the other hand, suppose that G contains a V-structure i l j. Let K [n] be the set of non-descendants of i and j in G, i.e. the set of k [n] such that there is no directed path from i to k or from j to k in G. Then the CI relation i j K is contained in the gaussoid corresponding to G. However, the CI relation i j K {l} is not in the gaussoid of G, since there is a Bayes ball path from i to j given K {l} in G. Hence by Lemma.7 above, the DAG associahedron corresponding to G is not MSS.. A construction of DAG associahedra from multiinformation The multiinformation of a probability measure P on [n] is a function m P : 2 [n] [0, ] defined by m P (S) = H(P Π i S P {i} ), where H denotes the relative entropy with respect to a product of one-dimensional marginals P {i}. Let P be a regular Gaussian measure on [n] with covariance matrix Σ. Let Γ be the correlation matrix of P a symmetric positive definite matrix obtained from Σ by simultaneously rescaling the rows and columns so that all the diagonal entries are equal to one. In other words, Γ = D /2 ΣD /2 where D = diag(σ). Then we have (5) i j K rank(γ Ki,Kj ) K where Γ A,B denotes the submatrix of Γ with rows and columns indexed by A and B respectively [Sul09]. By [Stu95, Corollary 2.6] the multiinformation m P (A) for A [n] is m P (A) = 2 log det(γ A,A). Since Γ is positive definite, all its principal minors det(γ A,A ) are non-zero. We define det(γ, ) to be. By [Stu95, Corollary 2.2] we have m P (A) = 0 for all A [n], A, and m P (ABC) + m P (C) m P (AC) + m P (BC) for all A, B, C [n] with equality if and only if A B C under P. We summarize this discussion in the following lemma.

11 GENERALIZED PERMUTOHEDRA FROM PROBABILISTIC GRAPHICAL MODELS Lemma.. If P is a regular Gaussian distribution with correlation matrix Γ, then its semigraphoid is submodular, with submodular function given by A 2 log det(γ A,A). The submodularity of DAG gaussoids (i.e. the first part of Theorem.2) follows from the lemma above and the fact that any DAG gaussoid has a faithful regular Gaussian realization. See for example [DSS09,.], where the following construction is described. Let G be a DAG on the vertices [n]. Assume that the vertices are labeled so that if i j is an edge in G, then i < j. Let Λ be an upper-triangular matrix whose entries have the form if i = j, Λ i,j = l ij if i j is an edge in G, 0 otherwise, where l ij are real numbers. Let K = ΛΛ T and Σ = K. Then K is symmetric positive definite by construction, and so is Σ. For almost all choices of real numbers l ij (apart from an algebraic hypersurface), a Gaussian distribution P with covariance matrix Σ is faithful to the DAG gaussoid of G [URBY]. In fact, as explained in the following lemma, the inequalities for the desired generalized permutohedron can also be computed directly from minors of K = ΛΛ T instead of from the correlation matrix Γ = D /2 Λ T Λ D /2, where D = diag(λ T Λ ). This result simplifies computations considerably since we don t need to perform any matrix inversion on ΛΛ T. Lemma.2. Let K be a positive definite matrix and let ω be the submodular function on 2 [n] given by ω(a) = log det(k A,A ). Let P be the polytope defined as in (). Then P is the generalized permutohedron corresponding to the semigraphoid of a regular Gaussian distribution P with covariance matrix Σ = K. Proof. The polytope defined by the submodular function A log det(γ A,A ) is obtained from the polytope defined by the submodular function A log det(σ A,A ) by translation in each coordinate direction i by log Σ i,i. Thus these two polytopes have the same normal fans and encode the same semigraphoids. For A [n] and B = [n]\a, we have (Σ A,A ) = K A,A K A,B (K B,B ) K B,A, the Schur complement. Using the equality det(k) = det(k B,B ) det(k A,A K A,B (K B,B ) K B,A ), we obtain log det(σ A,A ) = log det(σ A,A ) = log det(k A,A K A,B (K B,B ) K B,A ) = log det(k B,B ) log det(k). Combining this with Remark 2., it follows that the polytopes given by A log det(σ A,A ) and by A log det(k A,A ) are negatives of each other. In other words, by using K instead of Σ we obtain the dual semigraphoid defined in [MPS + 09]. In particular, if a semigraphoid has a faithful regular Gaussian distribution, then so does its dual. Example. (Multiinformation of the -node DAG in Example.). We start by constructing Λ from G using edge weights (i.e. l ij = if i j is an edge in G). We then compute ΛΛ T : Λ = K = ΛΛ T =

12 2 FATEMEH MOHAMMADI, CAROLINE UHLER, CHARLES WANG, AND JOSEPHINE YU Taking the log of the principal minors, we arrive at the system of inequalities: x log 2 x 2 log 2 x log 2 x 0 x + x 2 log x + x log x + x log 2 x 2 + x log x 2 + x log 2 x + x 0 x + x 2 + x log x + x 2 + x log x + x + x 0 x 2 + x + x 0 x + x 2 + x + x = 0 ( ) 2 For instance, the submatrix K {,},{,} is, whose determinant is, giving the inequality x + x log. These inequalities give a realization of the DAG associahedron in Example. 2 that is geometrically different from but has the same normal fan as the realization obtained by matroid polytopes in 5. As seen in the example above, the constant terms in the inequalities will almost never be rational numbers, making it difficult to obtain exact combinatorial information such as the f-vector and the normal fan. We found that the following heuristic works well in practice to obtain exact combinatorial information from this polytope description: First, round off the real numbers to nearby rational numbers (e.g. using 52 bit precision). Then, use exact arithmetic to compute the vertices of the polytope defined by these approximate inequalities. This results in perturbations of the true vertices. Now, form an approximate slack matrix by evaluating each approximate inequality at each approximate vertex and replace the entries in the slack matrix by 0 or depending on whether the entry is approximately zero or not (e.g. by rounding off to 5 bit precision) to obtain an incidence matrix between the vertices and facets. By eliminating duplicate rows and columns from this matrix we obtain the incidence matrix of a new polytope. In our simulations, the incidence matrix obtained this way gives the correct number of vertices and facets of the DAG associahedron, but this does not immediately lead to a rational realization of the polytope. Our code is available on Github at 5. A construction of DAG associahedra as Minkowski sums of matroid polytopes In the following, we obtain a construction of DAG associahedra as Minkowski sums of matroid polytopes, resulting in a rational realization of these polytopes. Until now we viewed a semigraphoid as defined by CI relations. However, we can equivalently define a semigraphoid by its complimentary conditional dependence relations. Minkowski addition of generalized permutohedra translates to taking the union of the corresponding conditional dependence relations. Thus, for every relation i j K in the gaussoid defined by a DAG G, we need to find a matroid whose semigraphoid (as defined by its rank function; see ()) contains the given relation. We now describe how to construct these matroids: For any conditional dependence relation i j K in the gaussoid defined by a DAG G, there is a Bayes ball path from i to j given K. We partition the path into canyons and treks as follows: A trek along a path is a consecutive subpath that does not contain any colliders. A canyon along a path is a consecutive subpath that is palindromic, contains exactly one collider, and whose first edge is oriented away from the first vertex. If we think of the arrows as always pointing down, then a canyon is a path that first goes down and then backtracks up the same edges to the first vertex. A single collider can be considered a canyon by itself. See Figure 5 for an example.

13 GENERALIZED PERMUTOHEDRA FROM PROBABILISTIC GRAPHICAL MODELS Figure 5. In the DAG, is a Bayes ball path from to 8 given {, 7}. The treks and canyons along the path are overlined and underlined respectively. If there is a Bayes ball path from i to j given K, we claim that there is one in which no vertex is repeated except in the same canyon. We call such a Bayes ball path simple. Suppose there is a repeated vertex a. Then we can take the first edge into a and the last edge out of a. This is allowed except when a would become a collider on the new path and a is not in the conditioned set K. In this case there must be a descendant of a that is a collider, hence in K, so we can make a canyon between a and this collider. For example, in Figure 5 the Bayes ball path 8 7 from to given {8} has a repeated vertex, and simply removing the path between the two occurrences of would give, which is not a Bayes ball path given {8} since the collider is not in the conditioned set {8}. However, has a descendant, 8, which is a collider in the original Bayes Ball path, so we can create a canyon 8 and take the path 8 instead. Let α be a Bayes ball path from i to j given K. From the discussion above, we may assume that the nodes on the path can be partitioned into an alternating sequence of treks and canyons, starting and ending with treks. Let d be the number of canyons or colliders along the path. We define a matroid M α on the vertex set [n] of the DAG G, which can be represented by affine independence among points in R d as follows: The one-element circuits (loops) are precisely the vertices that are not on the path α. The two-element circuits are precisely the pairs of vertices that are in the same trek or the same canyon. For each consecutive sequence of trek canyon trek, the corresponding three points affinely span a line in R d. The d lines corresponding to d canyons are affinely independent. See Figure 6 for an example. Recall from 2 that every matroid gives a collection of conditional dependence relations of the form i j K, where i, j [n] and K [n] \ {i, j} satisfy the condition (), namely (6) rank(k) + = rank(ki) = rank(kij) = rank(kj) Lemma 5.. Let G be a DAG and let α be a simple Bayes ball path in G. Then the conditional dependence relations of the matroid M α are contained in the set of conditional dependence relations defined by the gaussoid corresponding to G. The active paths in [Sha98] can be obtained from simple Bayes ball paths by replacing each canyon with only the top of the canyon, e.g. for the Bayes ball path 8 we get an active path. We prefer to keep the canyons in the path because we need them for our matroid construction. 8

14 FATEMEH MOHAMMADI, CAROLINE UHLER, CHARLES WANG, AND JOSEPHINE YU 6,7 2, 5,8 (a) The matroid corresponding to the Bayes ball path , which goes from to 8 given {, 7}. 7 5,6,8 (b) The matroid corresponding to the Bayes ball path , which goes from to 8 given {, 7}. The element 2 is a loop in the matroid, i.e. {2} is dependent. Figure 6. Two matroids that are compatible with the DAG in Figure 5. Proof. A subset S of a matroid is called closed if adding any element to S increases the rank. The span or the closure of any subset of the matroid is the smallest closed subset containing it. The condition (6) can be translated as (7) span(k) span(ki) = span(kij) = span(kj). A subset S M α is closed if and only if it satisfies both of the following conditions: () If an element of a trek or a canyon is in S, then all the other nodes in the same trek or canyon are also in S. (2) For any consecutive sequence of trek -canyon-trek 2 along the path α, if two elements of {trek, trek 2, canyon} are in S then so is the third element. The span of a set can be computed by adding nodes along the path to satisfy these two conditions. Suppose the relation i j K comes from the matroid M α, i.e. (7) is satisfied. We wish to show that there is a Bayes ball path in G between i and j given K. Let us first consider the case when i and j are in the same trek or in the same canyon. Then K cannot contain any element from the same trek/canyon; otherwise both i and j would be in span(k), contradicting (7). Thus there is a Bayes ball path between i and j given K, along α. Now suppose that i and j are in different treks/canyons. Then K cannot contain any element in the treks/canyons containing i or j; otherwise i or j would be in span(k), contradicting (7). We claim that span(k) cannot contain any consecutive trek-canyon-trek sequence between i and j. Assume the contrary. If we compute span(ki) by adding to span(k) nodes along the path starting at i, then the process would terminate at the trek-canyon-trek sequence, before it reaches j. Thus j / span(ki), contradicting span(ki) = span(kj). It follows that K does not intersect any adjacent pair of a trek and a canyon or two treks separated by a canyon between i and j along α. Next we claim that K cannot skip an adjacent pair of a trek and a canyon between i and j along α. If K skips a canyon-trek sequence from i to j, then, as before, the computation of span(ki) stops before it reaches j; so again j / span(ki). Putting everything together we conclude that K intersects all the canyons and none of the treks that lie strictly between i and j along α. This means that the subpath between i and j along α forms a Bayes ball path given K in G, and the relation i j K is valid for the gaussoid of G. We now prove our main result.

15 GENERALIZED PERMUTOHEDRA FROM PROBABILISTIC GRAPHICAL MODELS 5 Proof of Theorem.2. For every conditional dependence relation i j K in the gaussoid corresponding to a DAG G, we find a simple Bayes ball path α from i to j given K and construct a matroid M α. By Lemma 5., all the conditional dependence relations coming from M α are among those in the gaussoid corresponding to G. By taking all the matroids of the form M α where α runs over all simple Bayes ball paths in G, we obtain exactly all the conditional dependence relations that are valid in the gaussoid of G. Taking the union of all these dependence relations translates into taking the Minkowski sum of all the corresponding matroid polytopes. Hence, DAG associahedra are obtained as Minkowski sums of matroid polytopes. 6. Relationship among families of semigraphoids Recall from Lemma 2. that submodular functions on 2 [n] correspond to polytopal coarsenings of the S n fan, and these are exactly the normal fans of generalized permutohedra. Graph associahedra and DAG associahedra are special classes of generalized permutohedra that are defined up to equivalence of normal fans. The former can be realized as Minkowski sums of standard simplices and the latter can be realized as Minkowski sums of matroid polytopes (MSMP). What additional classes of generalized permutohedra can be realized in this way? Since the standard simplices are matroid polytopes, MSS polytopes are also MSMP. Unfortunately, this question seems difficult to answer in general. For n =,, 5 respectively, the cone of submodular functions has 5, 7, and 7978 extreme rays of which only 5, 2, and 9 respectively correspond to (connected) matroid polytopes. It suffices to consider connected matroids because the direct sum of matroids corresponds to the Minkowski sum of the corresponding matroid polytopes. Thus the matroid polytope of a disconnected matroid, which is the direct sum of nontrivial matroids, is the Minkowski sum of the matroid polytopes of these direct summands. Although the structure of these extreme rays is unclear, it seems unlikely due to their sparsity that many submodular semigraphoids will arise in this way. Another interesting class of semigraphoids are gaussoids [LM07], an abstraction of regular Gaussian distributions in the language of CI relations; see 2. Since we have seen that probabilistic graphical models can be faithfully realized by regular Gaussian distributions, another natural question is whether all regular Gaussian models (also called representable gaussoids) or even all gaussoids are MSMP. Gaussoids and appear to be incompatible with the MSMP construction. We have computationally verified that for n 8 no submodular semigraphoid corresponding to a connected matroid on [n] is a gaussoid. Thus, none of the extreme matroidal rays of the submodular cone are gaussoids. Conversely, not all gaussoids, in fact not even all representable gaussoids, can be obtained via MSMP. For example, [DX0, table A.] lists all Gaussian CI models on four variables (up to equivalence) and examples 9, 20,, 50, 5 are not MSMP. On the other hand, the CI relations corresponding to graphical models in this list all correspond to generalized permutohedra arising as MSMP. This leads to the question: Do we have submodularity for semigraphoids of other graphical models such as for various mixed graphs or chain graphs? If so, are the polytopes realizable as MSMP? In Figure 7 we illustrate the relationship of all the different coarsenings of the S n fan discussed in this paper by a Venn diagram. We have seen that undirected graphical models give rise to MSSs, while DAG models can be realized by MSMPs. In Proposition.6 we showed that a DAG model is MSS if and only if it coincides with an undirected graphical model. As we have discussed above, gaussoids are incompatible with the MSMP construction. In fact, gaussoids are also incompatible

16 6 FATEMEH MOHAMMADI, CAROLINE UHLER, CHARLES WANG, AND JOSEPHINE YU gaussoids representable gaussoids DAG models undirected graphical models MSMP MSS submodular semigraphoids semigraphoids Figure 7. Venn diagram representing the relationship of all the different coarsenings of the S n fan discussed in this paper. with the MSS construction. For example, it is easy to check that the standard simplex in Figure 2(a) is not a gaussoid. While every representable gaussoid is a submodular gaussoid as shown in Lemma., this is not the case for gaussoids. The semigraphoid studied in [HMS + 08, Section ] is a gaussoid that is not submodular. 7. Causal inference In this section, we describe how DAG associahedra can be used to perform causal inference. The main problem in causal inference is the following: We obtain data from an unobserved DAG G. From this data we infer a set of CI relations C. Under the faithfulness assumption, which we will assume throughout this section, C coincides with the gaussoid of G. The goal is to learn G from C. This problem is ill-defined since d-separation does not uniquely identify a DAG. So instead the problem is to learn G up to Markov equivalence, or in other words, to learn from C the essential graph, which is a partially directed graph with the same skeleton as G where an edge is directed if and only if it is directed the same way in every DAG in the Markov equivalence class. A popular algorithm for learning the Markov equivalence class of a DAG is Greedy Equivalence Search (GES) [Mee97, Chi02b], a greedy algorithm that searches through the space of DAGs by maximizing a scoring criterion such as the Bayesian Information Criterion (BIC). Under the faithfulness assumption GES is known to be consistent, i.e. it learns the correct essential graph [Mee97,Chi02b]. To reduce computation time, Teyssier and Koller [TK05] suggested to replace the greedy search in DAG space by a greedy search in the space of all orderings; a scoring criterion such as BIC is optimized by performing a walk on the edges of the permutohedron. Although no consistency guarantees were given for this greedy algorithm, simulations suggest that the greedy ordering-based search has a similar performance and lower computational costs as compared to GES [TK05]. In the following, we use our geometric insight on DAG associahedra to develop a new greedy orderingbased search with consistency guarantees. Let F be a coarsening of the S n fan. Each cone in F is defined by inequalities of the form x i x j and can be labeled a poset on [n]. Then we get a map from permutations of [n] to the set of partial orders on [n], derived from the map sending a maximal S n cone to the maximal cone F containing it. The preimage permutation (total order) is a linear extension of its image partial order. Hence the maximal cones of the coarsened S n fan or the vertices of the generalized permutohedron if the fan is polytopal can be labeled by posets so that every permutation is a linear extension of

17 GENERALIZED PERMUTOHEDRA FROM PROBABILISTIC GRAPHICAL MODELS 7 Algorithm Greedy SP algorithm on the permutohedron Input: A set of CI relations C on n random variables and a starting permutation π S n Output: An essential graph G. () Set t := 0 and π (0) := π. (2) Set t := t +. Randomly select a permutation π (t) that differs from π (t ) in a single adjacent transposition such that G π (t) is at least as sparse as G π (t ). () Iterate (2) until convergence to the sparsest Markov equivalence class and output the corresponding essential graph. exactly one of the posets. If two permutations π and τ are mapped to the same partial order, then we denote this by π τ. A semigraphoid C on [n] also gives a map from S n to the set of DAGs on nodes [n] as described in [RU]: To every permutation π we associate a DAG G π with (8) (π i, π j ) G π i < j and π i π j {π,..., π max(i,j) } \ {π i, π j }. In other words, the edge directions in the graph must be compatible with the ordering π = (π π 2 π n ), and the existence of an edge means that the two nodes are not independent given all the nodes that come before them in the ordering. We call π a topological ordering of G if any edge (i, j) in G implies that i j in π. Note that if the semigraphoid comes from a DAG G and π is a topological ordering of G, then G = G π. In [RU], it was proposed to use the number of edges of G π as a scoring criterion. It was shown that an algorithm that outputs the Markov equivalence class of G π with the fewest number of edges is consistent, i.e. it outputs the correct Markov equivalence class, under strictly weaker conditions than faithfulness. A permutation π giving a sparsest DAG is called a sparsest permutation. However the sparsest permutation (SP) algorithm is problematic from a computational point of view since it requires searching over all permutations. Instead, similarly as suggested in [TK05], we can perform a greedy search by traversing the edges of the permutohedron, using the number of edges of G π as a scoring function (see Algorithm ). Theorem 7.. The Greedy SP algorithm is consistent under the faithfulness assumption. We omit this proof since it is identical to the proof of Theorem 7.5, the consistency of the computationally more efficient Algorithm 2. Algorithm requires searching through neighboring permutations even when they give rise to the same DAG. For example, the neighboring permutations π = ( 2 ) and τ = (2 ) in Example. give rise to the same DAG G π = G τ = G shown in Figure (left). We will show that we can reduce the search space and hence computation time by performing the greedy search on the smaller DAG associahedron instead of the full permutohedron and keep the same consistency guarantees. The difficulty is that this needs to be done without having access to the DAG G on which the DAG associahedron is based. In order to do this, we give a description of the vertices and edges of a DAG associahedron in terms of the DAGs G π that are associated to its vertices. Theorem 7.2. For any fixed graphoid and two permutations π and τ, we have π τ G π = G τ. Moreover, the equivalence class of π consists of all topological orderings of G π. Proof. Suppose π τ. We may assume that π = (a a k i j b... b n k 2 ) and τ = (a a k j i b... b n k 2 ),