## CONDITIONAL INDEPENDENCE IN STATIONARY DIFFUSIONS

TOBIAS BOEGE, MATHIAS DRTON, BENJAMIN HOLLERING, SARAH LUMPP, PRATIK MISRA, AND DANIELA SCHKODA

Abstract. Stationary distributions of multivariate diffusion processes have recently been proposed as probabilistic models of causal systems in statistics and machine learning. Motivated by these developments, we study stationary multivariate diffusion processes with a sparsely structured drift. Our main result gives a characterization of the conditional independence relations that hold in a stationary distribution. The result draws on a graphical representation of the drift structure and pertains to conditional independence relations that hold generally as a consequence of the drift's sparsity pattern.

## 1. Introduction

Consider an R n -valued stationary stochastic process { X ( ) t } t ≥ 0 that arises as the solution of a stochastic differential equation (SDE)

$$\mathrm d \mathbb { X } ( t ) = f ( \mathbb { X } ( t ) ) \, \mathrm d t + D \, \mathrm d \mathbb { W } ( t ),$$

understood in the Itˆ sense [1]. o Here, f : R n → R n is a Lipschitz-continuous map that specifies the drift, D ∈ R n × n is a diagonal and invertible matrix, and W is a standard n -dimensional Brownian motion. Assuming its unique existence, the stationary (marginal) distribution of the process X depends on D only through the diffusion matrix C = DD T , and we denote the distribution by P f,C . Our interest lies in examining the conditional independence (CI) relations that hold in P f,C when the drift is sparse, i.e., when each component function f i only depends on some of its input coordinates. Let X be an R n -valued random vector distributed according to P f,C . Let I, J, K ⊆ [ n ] := { 1 , . . . , n } be disjoint index sets, and denote subvectors of X as X I = ( X i ) i ∈ I . Then we ask which sparsity restrictions on f imply that X I and X J are conditionally independent given X K ; in symbols, X I ⊥ ⊥ X J | X K or briefer I ⊥ ⊥ J | K . The main results of this paper answer this question via a characterization based on a directed graph that encodes the sparsity assumptions about the map f .

Before proceeding, we firm up the framework. If f is Lipschitz-continuous, the SDE in (1.1), initialized at a random vector X (0) that is independent of the driving Brownian motion, has a unique but not necessarily stationary solution. Existence and uniqueness of a stationary solution may be guaranteed in the framework provided by [11]. To this end, we assume that f has component functions f i that are locally m -integrable for m&gt;n , and that f satisfies a Lyapunov condition. Referring to the entries of the diagonal diffusion matrix C = ( c ij ), define the adjoint Fokker-Planck operator

$$\dot { \mathcal { L } } = \sum _ { i = 1 } ^ { n } \frac { 1 } { 2 } c _ { i i } \partial _ { i } ^ { 2 } + f _ { i } \partial _ { i }.$$

Key words and phrases. Conditional independence, graphical model, Lyapunov equation, Markov process, Ornstein-Uhlenbeck process.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00583v1-figure-60574d12-c92e-41bd-a49f-55540ad21e82.jpg)

2

3

1

4

ˆ

G

- (a) Directed graph on four nodes including all self-loops.

(b) Corresponding bidirected trek graph.

Figure 1. Directed graph G and its trek graph ˆ . G

The Lyapunov condition then requires existence of a non-negative function U ∈ C 2 ( R n ), the Lyapunov function, fulfilling

$$\lim _ { \| x \| \to \infty } U ( x ) = \infty, \text{ and } \lim _ { \| x \| \to \infty } \sup _ { \mathcal { L } U } \mathcal { L } U ( x ) < 0.$$

In the sequel, we use graphs to represent the sparsity structure of a drift map f . Specifically, we consider directed graphs G = ( V, E ) whose vertex set V = [ n ] indexes the coordinate processes of X . For all graphs considered in the paper, we tacitly assume that the edge set includes all self-loops i → i . Let Pa G ( ) = i { j ∈ [ n ] : ( j, i ) ∈ E } denote the set of parents of node i ∈ V in graph G . By our convention on self-loops, i ∈ Pa ( ). G i

Definition 1.1. We call a drift map f : R n → R n G -compatible if it is Lipschitz-continuous, satisfies the above local integrability and Lyapunov conditions, and has its sparsity pattern given by G , that is, each coordinate function f i : R n → R only depends the values x Pa G ( ) i of input coordinates indexed by Pa G ( ), i i ∈ [ n ].

̸

We now introduce our main theorem. If there is a directed path from node i to node j in the directed graph G , then i is an ancestor of j . Write An G ( ) for the set of all ancestors of i i . Note that i ∈ An ( ). Given G i G , we construct an auxiliary graph ˆ = ( G V, E ˆ ) that we call the trek graph . It is obtained by including a bidirected edge i ↔ j in ˆ if An E G ( ) i ∩ An ( ) = G j ∅ .

Theorem 1.2. Let G = ( V, E ) be a directed graph on V = [ n ] , and let ˆ G be its trek graph. Let I, J, K ⊆ [ n ] be three disjoint index sets. Then the following two statements are equivalent:

- (i) The set I is separated from J by V \ ( I ∪ J ∪ K ) in the trek graph ˆ G .
- (ii) The conditional independence X I ⊥ ⊥ X J | X K holds in any random vector X that is distributed according to a stationary distribution P f,C given by a G -compatible drift map f : R n → R n and a positive definite diagonal matrix C ∈ R n × n .

Example 1.3. Consider the graph G depicted in Figure 1a, and let X be a random vector following a stationary distribution given by a G -compatible drift map and a positive definite diagonal diffusion matrix. The graph G has the ancestor sets,

$$\ A n _ { \mathcal { G } } ( 1 ) = \{ 1 \}, \ \ A n _ { \mathcal { G } } ( 2 ) = \{ 1, 2 \}, \ \ A n _ { \mathcal { G } } ( 3 ) = \{ 1, 2, 3, 4 \}, \ \ A n _ { \mathcal { G } } ( 4 ) = \{ 4 \}.$$

As only the pairs of nodes (1 , 4) and (2 , 4) do not have a common ancestor, the trek graph ˆ G is missing precisely the edges 1 ↔ 4 and 2 ↔ 4; see Figure 1b. For subsets of nodes I , J , and S , write I ⊥ ˆ G J | S if I and J are separated by S in ˆ (i.e., every path from a node G i ∈ I to a node j ∈ J intersects S ). Then Theorem 1.2 implies the following marginal and conditional independencies:

$$\text{es} \colon \\ \{ 1, 2 \} \perp _ { \hat { \varphi } } \{ 4 \} \, | \, \{ 3 \} \quad \Longrightarrow \quad X _ { \{ 1, 2 \} } \oplus X _ { 4 }, \\ \{ 1 \} \perp _ { \hat { \varphi } } \{ 4 \} \, | \, \{ 2, 3 \} \quad \Longrightarrow \quad X _ { 1 } \oplus X _ { 4 }, \\ \{ 2 \} \perp _ { \hat { \varphi } } \{ 4 \} \, | \, \{ 1, 3 \} \quad \Longrightarrow \quad X _ { 2 } \oplus X _ { 4 }, \\ \{ 1 \} \perp _ { \hat { \varphi } } \{ 4 \} \, | \, \{ 3 \} \quad \Longrightarrow \quad X _ { 1 } \oplus X _ { 4 } \, | \, X _ { 2 }, \\ \{ 2 \} \perp _ { \hat { \varphi } } \{ 4 \} \, | \, \{ 3 \} \quad \Longrightarrow \quad X _ { 2 } \oplus X _ { 4 } \, | \, X _ { 1 }. \\ \text{t listed separation implies the other separations graphically, and } t$$

Here, the first listed separation implies the other separations graphically, and the first listed marginal independence implies the other independencies probabilistically by the semi-graphoid axioms of conditional independence [20, § 2.2].

Theorem 1.2 further states that no marginal or conditional independencies other than the ones listed above hold for all considered stationary distributions. For instance, Theorem 1.2 implies existence of a G -compatible drift map f and a positive definite diagonal diffusion matrix C such that

̸

$$X _ { 1 } \not \perp X _ { 4 } \, | \, X _ { 3 }$$

for X ∼ P f,C . We demonstrate how to construct such distributions in Section 4. △

The proof of Theorem 1.2 will be developed in the subsequent sections. The fact that separation in the trek graph implies conditional independence will be derived with the help of the Fokker-Planck equation that characterizes the density of P f,C (Section 2). In constructing counterexamples that exhibit conditional dependence when separation fails, we strive hard to structure them as simple as possible. Specifically, our counterexamples take f to be linear, which gives an Ornstein-Uhlenbeck process for which we let the diffusion matrix be a multiple of the identity (Section 4).

Remark 1.4 (Applied motivation) . Stationary diffusion processes have emerged as a new framework for statistical analysis of causal systems. Indeed, in many applications it is difficult if not impossible to observe a system more than once; e.g., in biology, measurement may destroy the observed collection of cells. The available data then constitute a sample of independent observations, for which the new paradigm retains an implicit temporal perspective by considering one-time snapshots of independent stationary dynamic processes. This approach was pioneered using stationary Ornstein-Uhlenbeck processes, resulting in Gaussian models termed Lyapunov models [9, 24]. Research progresses on issues like identifiability of drift parameters, statistical estimation, and new applications [3,4,19,25]. Moreover, extensions to non-linear diffusions were proposed in [14]. Theorem 1.2 advances this line of research via a characterization of conditional independence that can be used to disprove parameter identifiability [4]. It also reveals that even for general drift maps, the trek graph ˆ may be G estimated by learning a marginal independence graph [2,24].

Remark 1.5 (Relation to structural equation models) . Working with stationary distributions of sparse multivariate diffusion processes has particular appeal when modeling causal systems with feedback loops, which present mathematical and interpretational challenges for the

classical structural equation models (SEMs) [6,17]. For the latter, the well-known d-separation criterion yields a graphical characterization of conditional independence relations [15]. In our Theorem 1.2, a marginal independence X I ⊥ ⊥ X J (with conditioning set K = ∅ ) arises from the absence of common ancestors of nodes in I and nodes in J , which is precisely the same phenomenon as in d-separation. However, as will be evident from our later proof, Theorem 1.2 shows that for stationary diffusions the marginal independencies are the only source of independence relations, which radically differs from the case of d-separation.

The remainder of the paper is organized as follows. Section 2 treats the characterization of marginal independence and describes how to infer the conditional independence relations implied by these marginal independencies by means of a connected set Markov property for the trek graph. In Section 3, we set up notation, formally introduce the Lyapunov model and review some facts on Gaussian distributions. In Section 4, we prove the main result by showing that all conditional independence relations that hold in the Lyapunov model originate from marginal independencies. To this end, we reduce the problem to graphs obtained from 1-clique sums of treks. We then use a combination of algebraic and analytic methods to establish the result in this simpler case.

## 2. Marginal Independence and Implied Conditional Independencies

This section proves one direction in Theorem 1.2, namely, that V \ ( I ∪ J ∪ K ) separating I and J in the trek graph ˆ implies the conditional independence G X I ⊥ ⊥ X J | X K . We start with the case K = ∅ , i.e., marginal independence statements. Then, we show that these statements imply all claimed conditional independence statements. Before giving the proofs, we introduce graphical concepts, as used in this as well as later sections of the paper.

2.1. Graphs. For a directed graph G = ( V, E ) with edge set E ⊆ V × V , we write i → j if ( i, j ) ∈ E . We denote the set of nodes or edges of a specific graph G by V ( G ) or E ( G ), respectively. The set Pa G ( ) = i { j ∈ V : j → } i collects the parents of a node i ∈ V , and Ch ( ) = G i { j ∈ V : i → } j is its set of children. A directed path is a sequence of nodes ( i 1 , . . . , i m ) such that i l → i l +1 ∈ E for 1 ≤ l &lt; m , and m -1 is the length of this path. A path of length 0 consists of a single node. As defined in the introduction, An G ( ) is the set i of ancestors of i , i.e., the nodes j from which there is a directed path to i . This definition is conjunctively applied to sets, so An ( G I ) = ⋃ i ∈ I An ( ). G i The set of descendants of a node i is denoted De G ( ) and consists of all nodes i j such that there is a directed path from i to j . Typically, it is clear which graph G is concerned, and we drop the subscript, writing An( ) i ≡ An ( ) and, similarly, Pa( ) and Ch( ). G i i i The respective sets where the node itself is excluded are denoted in lowercase, e.g., pa( ) = Pa( ) i i \ { } i and ch( ) = Ch( ) i i \ { } i .

Combining two directed paths yields a trek: An ( ℓ, r )-trek between nodes i ℓ and j r is a sequence of nodes ( i , . . . , i ℓ 1 , t, j 1 , . . . , j r ) such that both ( t, i 1 , . . . , i ℓ ) and ( t, j 1 , . . . , j r ) are directed paths. Node t is the top node of the trek. If there is a trek between two nodes i and j , then the top node is a common element of An( ) and An( i j ). Thus, the trek graph ˆ from G Theorem 1.2 has two nodes adjacent precisely if there is a trek between them.

A subset of nodes I ⊆ V is ancestral if An( I ) = I . A trivial example of an ancestral set is a source node, meaning a node without parents. If a directed path starts and ends at the same node, then it forms a directed cycle. If the graph G does not contain any directed

̸

cycles, then it is termed a directed acyclic graph (DAG). The vertex set of a DAG may be ordered topologically as ( i 1 , . . . , i n ) such that for l = k , i l ∈ An( i k ) implies that l &lt; k .

- 2.2. Marginal Independence. The first result treats the implication ( ) i ⇒ ( ii ) in Theorem 1.2 when the conditioning set K is empty, in which case the independence is marginal.

Theorem 2.1. Let G = ( V, E ) be a directed graph on V = [ n ] . If f is G -compatible, then X ∼ P f,C satisfies X I ⊥ ⊥ X J for all sets of nodes I, J ⊆ V without trek between them.

This theorem generalizes an observation concerning Ornstein-Uhlenbeck processes made in [24], and we prove it using the Fokker-Planck equation. To this end, we draw on the next lemma, which follows from combining Theorem A and Proposition 2.1 in [11]. Here, C ∞ 0 ( R n ) is the set of smooth functions with compact support.

Lemma 2.2. If the Lyapunov condition is fulfilled, then (1.1) has a unique stationary solution admitting a density p . This density is the unique density fulfilling the Fokker-Planck equation

$$\int _ { \mathbb { R } ^ { n } } \mathcal { L } \phi ( x ) p ( x ) \, \mathrm d x = 0 \text{ for all } \phi \in C _ { 0 } ^ { \infty } ( \mathbb { R } ^ { n } ).$$

Proof of Theorem 2.1. By Lemma 2.2, the stationary distribution P f,C has a density p . We begin by showing that for any ancestral set A ⊆ V , the marginal density p A of X A satisfies the Fokker-Planck equation for the SDE restricted to A . Let A c = V \ A be the set of all remaining nodes. Fix a function ϕ ∈ C ∞ 0 ( R n ) with ϕ x ( ) only depending on x A . Then the Fokker-Planck equation for the overall density p yields that

$$0 = \int _ { \mathbb { R } ^ { n } } \mathcal { L } \phi ( x ) p ( x ) \, \mathrm d x = \int _ { \mathbb { R } ^ { n } } \left ( \sum _ { i = 1 } ^ { n } \frac { 1 } { 2 } c _ { i i } \partial _ { i } ^ { 2 } \phi ( x ) + f _ { i } ( x ) \partial _ { i } \phi ( x ) \right ) p ( x ) \, \mathrm d x \\ = \int _ { \mathbb { R } ^ { n } } \left ( \sum _ { i \in A } \frac { 1 } { 2 } c _ { i i } \partial _ { i } ^ { 2 } \phi ( x ) + f _ { i } ( x ) \partial _ { i } \phi ( x ) \right ) p ( x ) \, \mathrm d x \\ \text{the same as a $\mathcal{a}(\mathcal{a}(\mathcal{a})\mathcal{e}0$ for $i$} \, \mathcal{a}$ as $not$ only random$} \, on $r$. \quad \text{New the ancestrality of $A$} \, \text{immilies}$$

because ∂ ϕ x i ( ) ≡ 0 for i / A ∈ as ϕ only depends on x A . Now, the ancestrality of A implies that f i for i ∈ A only depends on x A . Rewriting the above integral, we obtain, by Fubini,

$$f _ { i } \text{ for } i \in A & \text{ only depends on } x _ { A }. \text{ Rewriting the above integral, we obtain, by Full } \\ 0 & = \int _ { \mathbb { R } ^ { A } } \int _ { \mathbb { R } ^ { A ^ { c } } } \left ( \sum _ { i \in A } \frac { 1 } { 2 } c _ { i i } \partial _ { i } ^ { 2 } \phi ( x ) + f _ { i } ( x ) \partial _ { i } \phi ( x ) \right ) p ( x ) \, d x _ { A ^ { c } } \, d x _ { A } \\ & = \int _ { \mathbb { R } ^ { A } } \left ( \sum _ { i \in A } \frac { 1 } { 2 } c _ { i i } \partial _ { i } ^ { 2 } \phi ( x ) + f _ { i } ( x ) \partial _ { i } \phi ( x ) \right ) \int _ { \mathbb { R } ^ { A ^ { c } } } p ( x ) \, d x _ { A ^ { c } } \, d x _ { A } \\ & = \int _ { \mathbb { R } ^ { A } } \left ( \sum _ { i \in A } \frac { 1 } { 2 } c _ { i i } \partial _ { i } ^ { 2 } \phi ( x _ { A } ) + f _ { i } ( x _ { A } ) \partial _ { i } \phi ( x _ { A } ) \right ) p _ { A } ( x _ { A } ) \, d x _ { A }, \\ \text{ is the Fokker-Planck equation for a restricted SDE of the form }$$

which is the Fokker-Planck equation for a restricted SDE of the form

$$\mathrm d \mathbb { X } _ { A } ( t ) = f _ { A } ( \mathbb { X } _ { A } ( t ) ) \, \mathrm d t + D _ { A, A } \, \mathrm d \mathbb { W } ( t )$$

with D A,A D T A,A = C A,A . Applying this result to A = An( ) I ∪ An( J ) gives a Fokker-Planck equation for the marginal density of X An( ) I ∪ An( J ) . To show the independence statement, we show that the product of marginal densities p An( ) I p An( J ) fulfills the same Fokker-Planck

equation. Let ˜ = An( I I ), ˜ = An( J J ), and ϕ ∈ C ∞ 0 ( R A ). Since there is no trek between I and J , the sets ˜ and I ˜ are disjoint and partition J A = ˜ I ∪ ˜ . J Therefore,

$$\int _ { \mathbb { R } ^ { A } } \left ( \sum _ { i \in A } \frac { 1 } { 2 } c _ { i i } \partial _ { i } ^ { 2 } \phi ( x _ { A } ) + f _ { i } ( x _ { A } ) \partial _ { i } \phi ( x _ { A } ) \right ) ( p _ { i } p _ { j } ) ( x _ { A } ) \, d x _ { A } \\ = & \int _ { \mathbb { R } ^ { A } } \left ( \sum _ { i \in I } \frac { 1 } { 2 } c _ { i i } \partial _ { i } ^ { 2 } \phi ( x _ { A } ) + f _ { i } ( x _ { A } ) \partial _ { i } \phi ( x _ { A } ) \right ) ( p _ { i } p _ { j } ) ( x _ { A } ) \, d x _ { A } \\ + & \int _ { \mathbb { R } ^ { A } } \left ( \sum _ { i \in \bar { J } } \frac { 1 } { 2 } c _ { i i } \partial _ { i } ^ { 2 } \phi ( x _ { A } ) + f _ { i } ( x _ { A } ) \partial _ { i } \phi ( x _ { A } ) \right ) ( p _ { i } p _ { j } ) ( x _ { A } ) \, d x _ { A }. \\ \text{e rewrite the first integral as}$$

We rewrite the first integral as

$$\text{We rewrite the first integral as} \\ & \quad \int _ { \mathbb { R } ^ { J } } \left ( \sum _ { i \in \bar { I } } \frac { 1 } { 2 } c _ { i i } \partial _ { i } ^ { 2 } \phi ( x _ { A } ) + f _ { i } ( x _ { A } ) \partial _ { i } \phi ( x _ { A } ) \right ) ( p _ { \bar { I } } p _ { \bar { J } } ) ( x _ { A } ) \, \mathrm d x _ { A } \\ = & \int _ { \mathbb { R } ^ { J } } \int _ { \mathbb { R } ^ { \bar { I } } } \left ( \sum _ { i \in \bar { I } } \frac { 1 } { 2 } c _ { i i } \partial _ { i } ^ { 2 } \phi ( x _ { A } ) + f _ { i } ( x _ { A } ) \partial _ { i } \phi ( x _ { A } ) \right ) ( p _ { \bar { I } } p _ { \bar { J } } ) ( x _ { A } ) \, \mathrm d x _ { \bar { I } } \, \mathrm d x _ { \bar { J } } \\ = & \int _ { \mathbb { R } ^ { J } } \int _ { \mathbb { R } ^ { \bar { I } } } \left ( \sum _ { i \in \bar { I } } \frac { 1 } { 2 } c _ { i i } \partial _ { i } ^ { 2 } \phi ( x _ { A } ) + f _ { i } ( x _ { A } ) \partial _ { i } \phi ( x _ { A } ) \right ) p _ { \bar { I } } ( x _ { \bar { I } } ) \, \mathrm d x _ { \bar { I } } p _ { \bar { J } } ( x _ { \bar { J } } ) \, \mathrm d x _ { \bar { J } }. \\ \text{The Fokker-Plank equation for $\tilde{I}$, with test function $x_{\bar{I}}\mapsto \phi(x_{\bar{I}},x_{\bar{J}})$ for fixed $x_{\bar{J}}$, yields}$$

$$\mathbb { K } \overset { \mathbb { K } } { \mathbb { K } }$$

↦

The Fokker-Planck equation for ˜ , I with test function x ˜ I → ϕ x , x ( ˜ I ˜ J ) for fixed x ˜ J , yields that the inner integral vanishes. Hence, the overall integral vanishes as well. Reversing the roles of I and J shows that also the second summand in (2.1) is zero. In conclusion, p An( ) I p An( J ) fulfills the Fokker-Planck equation for p An( ) I ∪ An( J ) , i.e., p An( ) I p An( J ) is a density of X An( ) I ∪ An( J ) . Consequently, X An( ) I ⊥ ⊥ X An( J ) and the decomposition property of conditional independence implies our claim that X I ⊥ ⊥ X J because I ⊆ An( ) and I J ⊆ An( J ). □

2.3. Implied Conditional Independencies. Let ˆ = ( G V, E ˆ ) be the trek graph associated to the given directed graph G . By definition, ( i, j ) ∈ ˆ if and only if there is a trek between E nodes i and j . Observe that ( i, j ) ∈ ˆ if and only if ( E j, i ) ∈ ˆ , E and, reflecting the structure of treks, we draw the edge bidirected as i ↔ j . Theorem 2.1 clarifies that the trek graph encodes marginal independencies, in particular,

$$i \leftrightarrow j \notin \hat { E } \implies X _ { i } \, \mathbb { \perp } \, X _ { j }$$

in random vectors X ∼ P f,C with f a G -compatible drift map. A more refined interpretation of the trek graph allows us to derive all conditional independencies that are implied by the marginal statements from Theorem 2.1. Recall the notation I ⊥ ˆ G J | S from Example 1.3, which indicates that the subsets of nodes I and J are separated by subset S in ˆ . G

Theorem 2.3. Let G = ( V, E ) be a directed graph on V = [ n ] , and consider a random vector X ∼ P f,C , where f : R n → R n is a G -compatible drift map and C ∈ R n × n is positive definite

Figure 2. A bidirected graph ˆ . G

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00583v1-figure-24cf5155-eb71-4c72-b44b-7cd9d74cbe8e.jpg)

and diagonal. Let I, J, K ⊆ [ n ] be three disjoint index sets, and let ˆ G be the trek graph of G . Then

$$I \perp _ { \hat { g } } J \, | \, V \, \langle I \cup J \cup K ) \ \implies \ X _ { I } \perp X _ { J } \, | \, X _ { K }.$$

Proof. In the literature on probabilistic graphical models, a random vector X that satisfies (2.2) is said to satisfy the pairwise Markov property relative to the bidirected graph ˆ [7]. G However, Theorem 2.1 shows that a richer set of marginal independencies hold. Call a subset A ⊆ V connected in ˆ if any two vertices G i, j ∈ A are connected by a path in ˆ all of whose G vertices are in A (i.e., A induces a connected subgraph of ˆ ). G Let Sp ˆ G ( A ) be the set of all nodes that are in A or are adjacent to an element of A via an edge of ˆ . G Then Theorem 2.1 implies that X satisfies the connected set Markov property relative to ˆ , G which requires that

$$X _ { A } \, \mathbb { \perp } \, X _ { V \vee \text{Sp} _ { \tilde { \mathcal { G } } } ( A ) } \ \forall \, A \subset V, \text{ non-empty and connected in } \hat { \mathcal { G } }.$$

By general properties of conditional independence, the connected set Markov property is equivalent to the global Markov property for bidirected graphs [7]. The global Markov property states that

$$I \perp _ { \hat { g } } J \, | \, V \, \langle I \cup J \cup K ) \ \implies \ X _ { I } \perp X _ { J } \, | \, X _ { K },$$

as was the claim of our theorem.

□

We illustrate the theorem in another example.

Example 2.4. Let ˆ be the bidirected graph from Figure 2. G The connected set Markov property requires

$$1 \perp \{ 3, 4 \}, \ \{ 1, 2 \} \perp 4.$$

As 1 is separated from 4 by 3 in ˆ , G the global Markov property additionally states 1 ⊥ ⊥ 4 | 2. The other conditional (and not marginal) independencies required by the global Markov property are

$$1 \, \perp \, 3 \, | \, 4, \quad 1 \, \perp \, 4 \, | \, 3, \quad 2 \, \perp \, 4 \, | \, 1.$$

In the next sections, we show that Theorem 2.3 already captures all conditional independencies that hold for general G -compatible drift. To this end, we consider Ornstein-Uhlenbeck processes and the graphical Lyapunov models obtained from their stationary distributions.

## 3. Ornstein-Uhlenbeck Processes and Lyapunov Models

This section briefly reviews Ornstein-Uhlenbeck processes and the associated graphical continuous Lyapunov models. We also recall basic facts regarding Gaussian conditional independence that will be used in Section 4.

- 3.1. Ornstein-Uhlenbeck Processes. The Ornstein-Uhlenbeck process is the only nontrivial time-homogeneous Markov process that is both Gaussian and stationary [5]. In the n -dimensional case, it is defined by the SDE

$$\mathrm d \mathbb { X } ( t ) = M ( \mathbb { X } ( t ) - \mu ) \, \mathrm d t + D \, \mathrm d \mathbb { W } ( t )$$

with an invertible drift matrix M ∈ R n × n encoding the interactions among the coordinates of X ( ). t The affine linear drift term introduces a restoring force towards the mean vector µ . Again we only consider the case with D and the diffusion matrix C = DD T diagonal.

Evidently, (3.1) is a special case of (1.1) with f being an affine function given by M and µ . In this case, the Lyapunov condition on f is satisfied if M is a stable matrix, that is, the real parts of the eigenvalues of M are negative. For M stable, X ( ) admits a unique stationary t distribution that is Gaussian with mean µ [12, 18]. The positive definite covariance matrix of the distribution is the unique solution Σ ≡ Σ( M,C ) of the continuous Lyapunov equation

$$M \Sigma + \Sigma M ^ { T } = - C.$$

We note that the solution Σ is invariant under joint rescaling of M and C .

↦

Given our interest in conditional independence structure, we restrict ourselves, without loss of generality, to stationary distributions of Ornstein-Uhlenbeck processes with mean µ = 0. For a directed graph G , the class of Gaussian distributions of interest naturally corresponds to the set of covariance matrices Σ that can be achieved using stable matrices M for which the linear map x → Mx is G -compatible. In other words, we consider stable matrices M whose support (i.e., non-zero entries) is contained in the edge set E ( G ).

Definition 3.1 ([4, 24]) . Let G = ( V, E ) be a directed graph on V = [ n ]. Let R E be the set of matrices M = ( m ij ) ∈ R n × n such that m ji = 0 if i → j / E ∈ and denote by PD n the cone of n × n positive definite matrices. Then for any fixed C ∈ PD , the (graphical continuous) n Lyapunov model of G is the family of all multivariate Gaussian distributions N (0 Σ) with , covariance matrix Σ in the set

$$\mathcal { M } _ { \mathcal { G }, C } = \{ \S$$

In the sequel, we identify the model with its set of covariance matrices and refer to M G ,C as the Lyapunov model. We denote the entries of a covariance matrix Σ by σ ij .

Convention 3.2. We denote the n -dimensional identity matrix by I n or equivalently by I V for a set V with | V | = n . We will typically fix C = 2 I n and write M G = M G , 2 I n . This assumption guarantees that I n ∈ M G ,C which is a desirable property for any Gaussian model. While the drift matrix M is typically thought of as the weighted adjacency matrix of G , we will usually suppress the self-loops in subsequent drawings of graphs for simplicity.

3.2. Conditional Independence in Multivariate Gaussians. Consider a symmetric matrix Σ with rows and columns indexed by a finite set V . For I, J ⊆ V , let Σ I,J denote the submatrix with rows indexed by I and columns by J . Submatrices of the form Σ K := Σ K,K are principal . Especially in subscripts, we will often abbreviate a union I ∪ J of index sets to IJ . A single index i is used interchangeably with its singleton set { } i in these conventions. So, iK := { } ∪ i K and K \ i := K \ { } i .

The starting point of our subsequent study of conditional independence is the following collection of results. Proofs can be found, for instance, in § 4.1 of [22].

Theorem 3.3. Let Σ ∈ PD n . The conditional independence (CI) statement I ⊥ ⊥ J | K holds in the Gaussian random vector X ∼ N ( µ, Σ) if and only if submatrix Σ IK,JK has rank | K | .

Since Gaussian distributions satisfy the compositional graphoid axioms, the CI statement I ⊥ ⊥ J | K holds if and only if i ⊥ ⊥ j | K holds for all i ∈ I and j ∈ J [13]. Thus, to characterize the set of CI statements that a Gaussian distribution satisfies, it suffices to characterize all statements of the form i ⊥ ⊥ j | K . This is summarized with the following corollary, which combines this fact with the previous theorem.

Corollary 3.4. For a Gaussian random vector X ∼ N ( µ, Σ) with Σ ∈ PD n , the following statements are equivalent:

- (i) I ⊥ ⊥ J | K ,
- (ii) i ⊥ ⊥ j | K for all i ∈ I and j ∈ J , and
- (iii) det(Σ iK,jK ) ≡ | Σ iK,jK | = 0 for all i ∈ I and j ∈ J .

This characterization of conditional independence for multivariate Gaussians will be used extensively in Section 4.

## 4. Conditional Independence in Lyapunov Models

Our focus is now on the Lyapunov model of a directed graph G = ( V, E ), which is given by the set of covariance matrices

- M G = { Σ ∈ PD n : M Σ+Σ M T +2 I n = 0 for a stable matrix M ∈ R E } .

Theorem 2.1 applies to M G and shows that the absence of a trek in G leads to a marginal independence, and Theorem 2.3 yields the implied additional conditional independence (CI) statements for M G . The aim of this section is to prove that the resulting collection of CI statements is complete , by showing that for any other CI statement, there is a distribution Σ ∈ M G which violates it. More precisely, our main result is as follows.

Theorem 4.1. Let G be a directed graph with trek graph ˆ G . Then the Lyapunov model M G is Markov-perfect to the marginal independence graph ˆ G . In other words, i ⊥ ⊥ j | K holds for all Σ ∈ M G if and only if V \ ( { i, j } ∪ K ) separates i and j in ˆ G .

The proof of this theorem will be completed at the end of this section. As noted above, the direction to be proven is the fact that a CI statement holds for all Σ ∈ M G only if we have graphical separation.

̸

4.1. Outline of the Proof. We will prove the contrapositive, i.e., we will show that if V \ ( { i, j } ∪ K ) does not separate i and j in ˆ , G then there exists a matrix Σ ∈ M G such that i ⊥ ⊥ j | K . To do this, we first clarify that we may reduce the problem to studies of suitable subgraphs, and we introduce a trek rule which gives a polynomial parameterization of a submodel of M G in terms of the entries of M (Section 4.2). We then show how to apply this restricted trek rule to the case where G is a single trek in order to produce covariance matrices Σ in M G such that the endpoints of the trek are conditionally dependent given all other nodes on the trek Via a trick of creating a limiting scenario of perfect correlation , we extend this result to more general conditioning sets comprised of arbitrary subsets of the other nodes on the trek (Section 4.3). Finally, to complete the proof, we extend the previous results from a single trek to zig-zag graphs, which correspond to paths that are connecting

for the global Markov property described in the proof of Theorem 2.3 (Section 4.4). We note that these zig-zag graphs also appear in the description of the d-separation criterion for structural equation models; cf. [16, Section 4.3].

- 4.2. Subgraphs and a Restricted Trek Rule. The following lemma clarifies how to pass to subgraphs for the construction of counterexamples to CI statements.

̸

Lemma 4.2. Let G = ( V, E ) be a directed graph, and let G ′ = ( V , E ′ ′ ) be a subgraph of G . If i ⊥ ⊥ j | K ′ on M G ′ , then i ⊥ ⊥ j | K on M G for all K ⊆ V with K ∩ V ′ = K ′ .

̸

̸

̸

Proof. The hypothesis that i ⊥ ⊥ j | K ′ means that there exists Σ ′ ∈ M G ′ for which | Σ ′ iK ,jK ′ ′ | = 0. Let M ′ be the corresponding parameter matrix, i.e., M ′ Σ +Σ ′ ′ M ′ T +2 I V ′ = 0. Denote V ′′ = V \ V ′ . Then the matrix equation with blocks indexed by V ′ and V ′′

$$\begin{array} {$$

̸

shows that Σ ∈ M G . Since K ∩ V ′ = K ′ , it follows that | Σ iK,jK | = Σ | ′ iK ,jK ′ ′ | = 0. □

The previous lemma will be used throughout this section to reduce the problem of determining whether a conditional dependence holds to a minimal subgraph which exhibits it. The next proposition is then our major tool to study the resulting minimal subgraphs.

Proposition 4.3 (Restricted trek rule) . Let G = ( V, E ) be a DAG, and consider the slice of the model M G obtained by restricting the drift matrix M = ( m ij ) to have diagonal entries m ii = -1 / 2 ζ for all i ∈ V . Then the following trek rule holds for the entries of the covariance Σ = ( σ ij ) associated to M :

$$\sigma _ { i j } = \sum _ { \substack { T \colon ($$

where the trek monomial m T associated to a trek T is given by m T = ∏ e ∈ T m e , i.e., the product over all the edges e in that trek.

Proof. The Lyapunov equation M Σ+Σ M T = -C is linear in Σ, and the ij -th element of the matrix equation reads

$$( m _ { i i } + m _ { j j } ) \sigma _ { i j$$

Note that the coefficient m ii + m jj of σ ij is always -1 / ζ . Solving for σ ij , we obtain

$$\sigma _ { i j } = \zeta \cdot \left ( \stackrel { \ddot { \sim } } { C } _ { i j } + \sum _ { l \in \text{pa} ( j ) } m _ { j l } \sigma _ { l i } + \sum _ { k \in \text{pa} ( i ) } m _ { i k } \sigma _ { k j } \right ).$$

These relations motivate a proof of (4.1) by induction.

We fix a topological ordering of the DAG G and extend it to a lexicographic ordering on the entries of Σ via σ ij &lt; σ kl if and only if j &lt; l or j = l and i &lt; k . We perform induction with respect to this ordering, noting that the right-hand side of (4.2) contains only Σ-entries which are strictly smaller than σ ij on the left-hand side. The smallest Σ-entries

Figure 3. The diamond graph discussed in Example 4.4.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00583v1-figure-cf521890-6695-4422-b48e-83062f32c747.jpg)

are σ ii where i is a source node. In this case σ ii = ζC ii = 2 ζ follows immediately from the Lyapunov equation for our choice of C = 2 I V . This verifies the claim as there is only the empty trek on i . Theorem 2.1 or equation (4.2) can be used to see that if i and k are distinct source nodes (having in particular no treks between them), then σ ki = 0. Now let j be a node all of whose parents are source nodes; the edges into such nodes j form the next smallest class of Σ-entries in the lexicographic ordering. The Lyapunov equation simplifies to σ ij = ζ ∑ k ∈ pa( ) j m σ jk ki = ζm σ ji ii = 2 ζ 2 m ji which corresponds to the unique trek between i and j of length ℓ + r = 1. These two cases prove the initial induction hypothesis.

Fix any two (not necessarily distinct) vertices i, j and a non-empty ( ℓ, r )-trek T between them. This trek decomposes into a unique ( ℓ, r -1)-trek T 1 from i to l , for l ∈ pa( ), and the j edge l → j . It also decomposes uniquely into an ( ℓ -1 , r )-trek T 2 from k to j , for k ∈ pa( ), i together with the edge k → i . By (4.2), the coefficient of m T in σ ij is thus ζ times the sum of the coefficients of m T 1 in σ il and m T 2 in σ kj . Using the induction hypothesis (since σ il and σ kj are both less than σ ij ) gives

$$\text{coefficient of } m _ { T } = \zeta \cdot \left ( 2 \zeta ^ { \ell + r } \binom { \ell + r - 1 } { \ell } + 2 \zeta ^ { \ell + r } \binom { \ell + r - 1 } { \ell - 1 } \right ) = 2 \zeta ^ { \ell + r + 1 } \binom { \ell + r } { \ell },$$

̸

as claimed. This immediately proves the induction step when i = . In case j i = , there is j also the empty trek from i to i to account for and its trek term arises as required from (4.2) via C ii = 2. □

Example 4.4. Consider the graph in Figure 3. Using the trek rule formulation for σ ij from Proposition 4.3, we get

$$\begin{array} { r c l }. 3, \, \text{we get} \\ \sigma _ { 1 1 } & = & 2 \zeta, \\ \sigma _ { 2 2 } & = & 4 \zeta ^ { 3 } m _ { 2 1 } ^ { 2 } + 2 \zeta, \\ \sigma _ { 3 3 } & = & 4 \zeta ^ { 3 } m _ { 3 1 } ^ { 2 } + 2 \zeta, \\ & \vdots \\ \sigma _ { 3 4 } & = & 6 \zeta ^ { 4 } m _ { 2 1 } m _ { 3 1 } m _ { 4 2 } + 6 \zeta ^ { 4 } m _ { 3 1 } ^ { 2 } m _ { 4 3 } + 2 \zeta ^ { 2 } m _ { 4 3 }. \\ \end{array}$$

For instance, observe that there are two treks from the vertex 2 to itself, the empty trek (corresponding to the monomial 2 ζ for σ 22 ) and the trek 2 ← 1 → 2 (4 ζ 3 m 2 21 ). Similarly, there are three treks between the vertices 3 and 4, which are the edge 3 → 4 (2 ζ 2 m 43 ), the trek 2 ← 1 → 3 → 4 (6 ζ 4 m m m 21 31 42 ), and 3 ← 1 → 3 → 4 (6 ζ 4 m m 2 31 43 ). For the last trek 3 ← 1 → 3 → 4, we have ℓ = 1 and r = 2, which gives us the coefficient 6 ζ 4 . △

Figure 4. A trek T which exhibits no conditional independence statements and a zig-zag graph Z which exhibits some marginal independencies.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00583v1-figure-a22067ef-8da5-4f77-a8ac-51c7f0b3c049.jpg)

This trek rule is reminiscent of the original trek rule for linear structural equation models [21, 23]. However, now each trek monomial is weighted by a binomial coefficient. While the trek rule we present here only parameterizes a proper subset of M G and is less general than that of [24], it provides an explicit rational formula for the entries of Σ in terms of the parameters M which makes it significantly easier to construct explicit covariance matrices in M G with the desired dependence properties. We remark that the very recent preprint of [10] gives a generalization of Proposition 4.3.

- 4.3. Conditional Independence on Treks. We now turn to the case where the considered graph is a single trek, for which we show that the associated Lyapunov model does not satisfy any CI relations. Before proceeding, we note that if T is a trek from node 1 to n with top node t , then a drift matrix M ∈ R E ( T ) has zero-pattern

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00583v1-figure-bfd296e5-b795-457b-97c5-8c04ce1c075c.jpg)

with diagonal entries d , . . . , d 1 n and the symbol ∗ indicating arbitrary values. Such matrices are stable if and only if the diagonal entries are negative.

̸

Proposition 4.5. Let T be an ( ℓ, r ) trek between i and j . Let K = V ( T ) \ { i, j } be the remaining vertices. Then i ⊥ ⊥ j | K in M T .

̸

The proof of this proposition and its several supporting lemmas can be found in Appendix A.1. The following example provides an explicit family of covariance matrices in the Lyapunov model for a small trek which exhibit the conditional dependence stated in Proposition 4.5. The proof of the previous proposition generalizes this family of covariance matrices to larger treks and shows that the conditional dependence still holds by inductively computing certain entries in the adjoint matrix of Σ iK,jK . The Schur complement expansion of | Σ iK,jK | with respect to the principal minor Σ K,K can then be used to show that | Σ iK,jK | = 0 by inductively computing the lowest degree term which appears in | Σ iK,jK | when it is viewed as a polynomial in ζ .

Example 4.6. Let T be the trek as shown in Figure 4. Setting every edge parameter m ji to 1 for all i → j ∈ E ( T ), we get the following covariance matrix from the trek rule which lies in the model:

$$\Sigma = \begin{pmatrix} 2 \zeta & 2 \zeta ^ { 2 } & 2 \zeta ^ { 2 } & 2 \zeta ^ { 3 } \\ 2 \zeta ^ { 2 } & 4 \zeta ^ { 3 } + 2 \zeta & 4 \zeta ^ { 3 } & 6 \zeta ^ { 4 } \\ 2 \zeta ^ { 2 } & 4 \zeta ^ { 3 } & 4 \zeta ^ { 3 } + 2 \zeta & 6 \zeta ^ { 4 } + 2 \zeta ^ { 2 } \\ 2 \zeta ^ { 3 } & 6 \zeta ^ { 4 } & 6 \zeta ^ { 4 } + 2 \zeta ^ { 2 } & 1 2 \zeta ^ { 5 } + 4 \zeta ^ { 3 } + 2 \zeta \end{pmatrix}.$$

̸

̸

Observe that the submatrix Σ { 1 2 3 , , } { , 2 3 4 , , } has determinant -8 ζ 6 , which is non-zero for ζ = 0, implying that 2 ⊥ ⊥ 4 | { 1 3 , } in the model M T . △

We now introduce a technique to extend the previous proposition to all conditioning sets K ⊆ V ( T ) \ { i, j } . In Lyapunov models with fixed diffusion matrix C = 2 I V , we cannot simply pass to marginal distributions as, in general, removing a row and corresponding column in a covariance matrix Σ ∈ M T does not yield a distribution in the Lyapunov model over the shortened trek in which the removed node is skipped. Instead, we provide a construction in which the random variables indexed by nodes in V ( T ) \ K become perfectly correlated with random variables indexed by elements of K . This construction then allows us to argue that the counterexample from Proposition 4.5 can be 'lifted' to counterexamples on longer treks between i and j that contain additional nodes not in the conditioning set K .

To introduce the construction, we recall the following notation. For a node k for which pa( k ) is a singleton set, we identify pa( k ) with its single element, i.e., the unique node j with j → k . Similarly, we identify ch( k ) with the unique node j such that k → j . Moreover, for k ∈ N , we define a map that shifts indices as

↦

$$\psi _ { - k } \, \colon \mathbb { N } \, \left \{ k \right \} \to \mathbb { N }, \quad x \mapsto \begin{cases} x & \text{if $x < k$;} \\ x - 1 & \text{if $x > k$}. \end{cases}$$

Now, consider a trek T between i and j with top node t and of length | V ( T ) | ≥ 4. If | ch( ) t | = 2, let c l and c r be the child nodes on the left and on the right side of the trek, respectively. Let k ∈ V ( T ) \ { i, j } be a node on the trek, and consider the subset of nodes V -k = V ( T ) \ { k } obtained by removing k . Define a shorter trek

$$\mathcal { T } _ { - k } = ( V _ { - k }, E _ { - k } )$$

̸

by constructing the edge set E -k as follows. If k = , replace the 2-path pa( t k ) → k → ch( k ) by the single edge pa( k ) → ch( k ). If k = , the trek t T has two non-trivial sides as k / ∈ { i, j } . In this case, if the right side of T has length 1, we replace c l ← k → c r by c l → c r . If k = t and the right side consists of more than one edge, we replace c l ← k → c r by c l ← c r . The construction is depicted in Figure 5. In all cases, the resulting graph is shorter by one node, as the k -th node is removed. Note that a similar definition is possible by adding the reverse edge c l → c r in the case where k = t and both sides of the trek consist of more than one edge. The distinction between the second and third case is necessary to keep the trek structure when removing a top node from a trek where one side consists of exactly one edge. Otherwise, the resulting graph could be a directed path, e.g., T -k could be i ← c l ← c r in Figure 5b. In the case that the nodes V ( T ) = [ n ] are numbered consecutively from left to right, we implicitly apply ψ -k to the elements of V -k and the start- and endpoint of each edge in E -k so that we have V -k = [ n -1] and E -k ⊆ [ n -1] × [ n -1].

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00583v1-figure-387b2075-8be6-4e53-a624-2e4901943dc3.jpg)

̸

- (a) Case t = k with ch( k ) =
- { } i and pa( k ) = { } t .

(c) Case t = k with ride side longer than length 1.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00583v1-figure-9f3d0b83-c5d1-4fd4-9183-97838f6cedf9.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00583v1-figure-310ce08d-6742-498f-b6dd-652bb8768793.jpg)

Figure 5. Construction of T -k depending on the position of the node k in T . The node k and its in- and outgoing edges (blue in T ) are replaced by a single edge oriented accordingly (blue in T -k ). The dashed edge indicates the two nodes that will be perfectly correlated in the following construction of a distribution in M T given a distribution in M T -k .

Given a matrix Σ, we write

$$\Sigma _ { - k, - k } = \Sigma _ { [ n ] \langle k \}, [ n ] \langle k \} }$$

for the submatrix obtained by removing the k th row and column of Σ. Given a set of matrices M⊆ R n × n , we denote its topological closure by M .

Lemma 4.7. Let T be a trek between i and j containing more than 3 nodes. Let t be the top node of T , and let k / ∈ { i, t, j } be a further node on T . Then for every matrix Σ ∗ ∈ M T -k , there exists a matrix Σ ∈ M T such that

$$\Sigma _ { - k, - k } = \Sigma _ { - \text{ pa} ( k ), - \text{ pa} ( k ) } = \Sigma ^ { * }.$$

Proof. Let n ≥ 4 be the number of nodes on the trek T . Without loss of generality, assume that the nodes on T are numbered consecutively from left to right, i.e., i = 1, j = n , and k ∈ { 2 , . . . , n - } \ { } 1 t .

From Σ ∗ ∈ M T -k , we are to construct a suitable larger covariance matrix Σ ∈ PD . To do n this, we make the new node at index k a copy of its parent on the trek by duplicating the respective row and column from Σ . In the resulting distribution, ∗ k and pa( k ) are perfectly correlated. Note that pa( k ) = { k +1 } for k &lt; t , and pa( k ) = { k - } 1 for k &gt; t . Thus, in the first case, the k th and k +1 st row and column of Σ are duplicates of each other, while in the second case, the k th and k -1 st row and column of Σ are duplicates of each other. The proof in both cases is analogous as the drift matrix of a single trek is defined in a symmetric way around its top node as seen in (4.3). We only give the proof for the second case, so assume k &gt; t . Furthermore, let U := { 1 , . . . , k - } 1 , W := { k +1 , . . . , n } , and write ψ ≡ ψ -k for the map from (4.4). Then, ψ U ( ) = U and ψ W ( ) = { k, . . . , n - } 1 . We note that the case of a directed path, i.e., a trek with i = t = 1, is included here as k &gt; 1 by assumption.

Given Σ ∗ ∈ M T -k , there exists a stable matrix M ∗ ∈ R ( n -1) × ( n -1) such that the Lyapunov equation

$$M ^ { * } \Sigma ^ { * } + \Sigma ^ { * } ( M ^ { * } ) ^ { T } + 2 I _ { n - 1 } = 0$$

is fulfilled. In the following, we show that our construction of inserting a duplicate of the row and column of Σ ∗ at index pa( k ) in T indeed yields a matrix Σ ∈ M T .

To achieve this, we construct a sequence of drift matrices M ( m ) as follows. First, partition

$$\Sigma ^ { * } = \left ( \frac { \Sigma _ { U } ^ { * } \left | \Sigma _ { U, \psi ( W ) } ^ { * } \right ) } { \Sigma _ { \psi ( W ), U } ^ { * } \right | \Sigma _ { \psi ( W ) } ^ { * } } \right )$$

as well as

$$M ^ { * } = \left ( \frac { M _ { U } ^ { * } \, \left | \, 0 } { M _ { \psi ( W ), U } ^ { * } \, \left | \, M _ { \psi ( W ) } ^ { * } } \right ), \quad \text{where} \quad M _ { \psi ( W ), U } ^ { * } = \left ( \begin{matrix} 0 & \cdots & 0 & m _ { k, k - 1 } ^ { * } \\ & & 0 \\ & 0 & \vdots \\ & & & 0 \end{matrix} \right ).$$

Considering (4.5) block-wise results in the following four equations:

- (a) M ∗ U Σ ∗ U +Σ ∗ U M ∗ U T +2 I k -1 = 0,
- (b) M ∗ U Σ ∗ U,ψ W ( ) +Σ ∗ U,ψ W ( ) M ∗ ψ W ( ) T +Σ ∗ U M ∗ ψ W ,U ( ) T = 0,
- (c) M ∗ ψ W ( ) Σ ∗ ψ W ( ) +Σ ∗ ψ W ( ) M ∗ ψ W ( ) T +2 I n -k + M ∗ ψ W ,U ( ) Σ ∗ U,ψ W ( ) +Σ ∗ ψ W ,U ( ) M ∗ ψ W ,U ( ) T = 0,
- (d) M ∗ ψ W ( ) Σ ∗ ψ W ,U ( ) +Σ ∗ ψ W ,U ( ) M ∗ U T + M ∗ ψ W ,U ( ) Σ ∗ U = 0,

where equations (b) and (d) are transposes of each other. Given the stable matrix M ∗ , these equations are uniquely fulfilled by the corresponding submatrices of Σ . ∗

Now, let m ∈ N &gt; 0 and define

$$M ^ { ( m ) } = \begin{pmatrix} M _ { U } ^ { * } & 0 & 0 \\ \frac { 0 - m } { 0 } & - m & 0 \\ 0 & m _ { k, k - 1 } ^ { * } & M _ { \psi ( W ) } ^ { * } \end{pmatrix} \in \mathbb { R } ^ { n \times n }.$$

Further, let

$$\Sigma ^ { ( m ) } = \begin{pmatrix} \Sigma ^ { ( m ) } _ { U } & \Sigma ^ { ( m ) } _ { U, k } & \Sigma ^ { ( m ) } _ { U, W } \\ \Sigma ^ { ( m ) } _ { k, U } & \Sigma ^ { ( m ) } _ { k, k } & \Sigma ^ { ( m ) } _ { k, W } \\ \Sigma ^ { ( m ) } _ { W, U } & \Sigma ^ { ( m ) } _ { W, k } & \Sigma ^ { ( m ) } _ { W } \\ \tau \in \end{pmatrix} \in \mathbb { R } ^ { n \times n }$$

be the solution of the resulting Lyapunov equation

$$M ^ { ( m ) } \Sigma ^ { ( m ) } + \Sigma ^ { ( m ) } ( M ^ { ( m ) } ) ^ { T } + 2 I _ { n } = 0$$

induced by M ( m ) for every m ∈ N &gt; 0 .

Due to the symmetry of the Lyapunov equation, any unique solution of (4.6) is symmetric. Since M ( m ) is stable, the solution Σ ( m ) is unique and therefore symmetric, so we only have to solve six matrix equations that arise by block-matrix multiplication, denoted by letters in brackets as follows:

$$\text{ets as follows:} \\ M ^ { ( m ) } \Sigma ^ { ( m ) } + \Sigma ^ { ( m ) } ( M ^ { ( m ) } ) ^ { T } + 2 I _ { n } = 0 \equiv \begin{pmatrix} ( A ) & ( B ) & ( C ) \\ & ( D ) & ( E ) \\ & & ( F ) \end{pmatrix}.$$

The upper left block Σ ( m ) U of the partitioned matrix Σ ( m ) is constant as it is the unique solution of equation (a). Proceeding recursively, the other blocks of Σ ( m ) can be shown to converge for m →∞ as well, with limits given by the relevant parts of Σ . ∗ The details of these calculations can be found in Appendix A.2. Combining the results, we have as desired.

$$\text{se calculations can be found in Appendix A.2. Combining the results, we have } \\ \lim _ { m \to \infty } \Sigma ^ { ( m ) } = \begin{pmatrix} \Sigma _ { U } ^ { * } & \Sigma _ { U, \beta ( k ) } ^ { * } & \Sigma _ { U, \psi ( W ) } ^ { * } \\ \frac { \Sigma _ { \text{pa} ( k ), U } ^ { * } & \Sigma _ { \text{pa} ( k ), \beta ( k ) } ^ { * } & \Sigma _ { \text{pa} ( k ), \psi ( W ) } ^ { * } \\ \Sigma _ { \text{\psi} ( W ), U } ^ { * } & \Sigma _ { \text{\psi} ( W ), \beta ( k ) } ^ { * } & \Sigma _ { \text{\psi} ( W ) } ^ { * } \end{pmatrix} \in \overline { M } _ { \mathcal { T } }, \\ \text{desired.} & \quad \Box \\ \text{rollary } 4.8. \ I. e t \, \mathcal { T } \, \text{ be a tree with endnoints } i \, \text{ and } i \, \text{ } \ I. e t \, \emptyset \neq K \, \text{ } V ( \mathcal { T } ) \, \vee \{ i \, i \} \text{ and } \text{ let}$$

Corollary 4.8. Let T be a trek with endpoints i and j . Let ∅ ̸ = K ⊊ V ( T ) \ { i, j } and let k / ∈ { i, j } ∪ K be a node on the trek. If there is

̸

$$\Sigma ^ { * } \in \mathcal { M } _ { \mathcal { T } _ { - k } } \ s u c h \ t h a t \left | \Sigma ^ { * } _ { \psi _ { - k } ( i K ), \psi _ { - k } ( j K ) } \right | \neq 0,$$

then there is

̸

$$\Sigma \in \mathcal { M } _ { \mathcal { T } } \ s u c h \ t h a t \, | \Sigma _ { i K, j K } | \neq 0.$$

̸

Proof. By continuity of the determinant, it suffices to construct a matrix Σ ̂ ∈ M T such that | ̂ Σ iK,jK | = 0. Considering the closure of the model is crucial as it allows us to apply Lemma 4.7 for the construction of Σ. Then as ̂ k / ∈ { i, j } ∪ K , we find that

̸

$$| \widehat { \Sigma } _ { i K, j K } | & = | \Sigma _ { \substack { \psi _ { - k } ( i K ), \psi _ { - k } ( j K ) } } ^ { * } | \neq 0. & \Box$$

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00583v1-figure-f54b198f-5fc4-4672-aca1-4c838ea248fd.jpg)

- (a) Iterative construction of p -trek zig-zag graph obtained via a 1-clique sum of a trek and p -1-trek zig-zag.

(b) The resulting p -trek zig-zag graph.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00583v1-figure-d5b40419-70af-43c9-a071-99cef03a8115.jpg)

Figure 6. A zig-zag graph Z with endpoints i and j where each edge in the above graph represents a directed path. It is obtained by taking a 1-clique sum along the nodes k α of the p treks T α with top node t α and endpoints s α and s α +1 . The last step of this iterative process is pictured in (a) while the result is pictured in (b).

- 4.4. Conditional Independence in Zig-Zag Graphs. We now introduce the family of graphs that will be the key element needed to complete the proof of Theorem 4.1.

Definition 4.9. Let T 1 , . . . , T p be treks such that for i ∈ [ p -1], there exists a node k i on the right branch of T i which is the endpoint of the left branch of T i +1 . A zig-zag graph Z (Figure 6) is a directed acyclic graph obtained by performing an iterated 1-clique sum of these treks along the nodes k i , i.e.

$$\mathcal { Z } = \mathcal { T } _ { 1 } \oplus _ { k _ { 1 } } \mathcal { T } _ { 2 } \oplus _ { k _ { 2 } } \cdots \mathcal { T } _ { p - 1 } \oplus _ { k _ { p - 1 } } \mathcal { T } _ { p }.$$

̸

Proposition 4.10. Let Z be a zig-zag graph with endpoints i and j , and let K ⊆ V ( Z \{ ) i, j } be a set such that for each collider k ∈ Z it holds that s ∈ K where s is the unique sink in De( k ) . Then i ⊥ ⊥ j | K in the Lyapunov model associated to Z .

̸

Proof. To see this we use Proposition 4.5 inductively on the number p of treks that Z is obtained from. If Z is simply a trek then the result already holds by Proposition 4.5. So suppose that Z p is a zig-zag graph with tails obtained from combining p treks and that for any set K ′ ⊆ V ( Z p ) of the required form it holds that k ⊥ ⊥ j | K ′ where k and j are the endpoints of Z p . Let T be a trek from i to s containing a node labelled k which is not on

̸

̸

the same branch as i and let Z p +1 be the 1-clique sum of T and Z p at k . Recall that this is exactly how every zig-zag graph is produced as pictured in Figure 6. Note that if k and i are on the same branch then checking that i ⊥ ⊥ j | K on Z p +1 naturally reduces to checking i ⊥ ⊥ j | K on a smaller zig-zag Z p which immediately follows from the induction hypothesis. To complete the induction we must show that i ⊥ ⊥ j | K where K is of the form K = T ∪ L ∪ K ′ such that T ⊆ V ( T ) \ { } ∪ i De( k ), s ∈ L ⊆ De( k ) and K ′ ⊆ Z p is a set such that k ⊥ ⊥ j | K ′ . Lastly, let l ∈ L be the node closest to k . Then observe that after rearranging the rows and columns of Σ iK,jK it has the block form

$$\Sigma _ { i K, j K } = \begin{array} {$$

Now observe that if we specialize this matrix by setting m ab = 0 and m la = 0 for all a, b ∈ L \ l then the above block matrix immediately simplifies to

$$\begin{array} {$$

̸

Thus we see that for any such matrix obtained from the submodel corresponding to this specialization, | Σ iK,jK | = Σ | lK ,jK ′ ′ | | Σ iT,lT | | Σ L l,L \ \ l | . Now observe that | Σ lK ,jK ′ ′ | = 0 by our induction hypothesis since this corresponds to conditional dependence of a zig-zag with endpoints l and j which is composed of only p treks. Further observe that | Σ iT,lT | = 0 by Corollary 4.8. The last factor is a principal minor and thus is nonzero for any Σ ∈ M Z p +1 . So any Σ which comes from the above submodel satisfies | Σ iK,jK | = 0 and thus i ⊥ ⊥ j | K for the Lyapunov model M Z . □

̸

̸

̸

̸

Example 4.11. Consider the zig-zag graph Z pictured in Figure 4. Then 3 ⊥ ⊥ 5 | { 1 4 , } since for any Σ ∈ M Z , it holds that

$$\Sigma _ { \{ 3, 1, 4 \}, \{ 5, 1, 4$$

As was shown in the previous proof, this matrix is block upper-triangular with blocks σ 45 and Σ { 3 1 , } { , 4 1 , } . Thus | Σ { 3 1 4 , , } { , 5 1 4 , , } | = σ 45 | Σ { 3 1 , } { , 4 1 , } | where the second term in the product is guaranteed to be nonzero for almost all choices of parameters by Proposition A.6. △

We now have all the necessary tools needed to complete the proof of Theorem 4.1.

Proof of Theorem 4.1. If V \ ( { i, j } ∪ K ) separates i and j in the marginal independence graph ˆ then Section 2.3 immediately implies that G i ⊥ ⊥ j | K holds for all Σ ∈ M G . It remains to show that if i ⊥ ⊥ j | K then V \ ( { i, j } ∪ K ) separates i in j in ˆ . G We instead show the contrapositive: if i and j are connected given V \ ( { i, j }∪ K ) in ˆ , G then there exists

̸

̸

Figure 7. The four possible intersections of a trek from i over t 1 to k and a trek from l over t 2 to j and how to shorten the walk from i to j (indicated in blue). Each arrow represents a directed path.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00583v1-figure-e271df4d-e020-48cc-972c-10a2501aad2e.jpg)

̸

̸

̸

Σ ∈ M G such that i ⊥ ⊥ j | K . We will prove this by constructing a specific Σ using the trek rule described in Proposition 4.3 and then showing that | Σ iK,jK | = 0 which by Corollary 3.4 implies that the conditional dependence i ⊥ ⊥ j | K holds.

If i and j are connected in ˆ given G V \ ( { i, j }∪ K ) then there is path in ˆ from G i to j which uses only the vertices in { i, j }∪ K . Without loss of generality, this path is simple, i.e., it does not repeat a vertex. Let i = k , k 0 1 , . . . , k s , k s +1 = j be the sequence of vertices on this path where { k , . . . , k 1 s } ⊆ K . Each edge k p -k p +1 in ˆ corresponds to a trek in G G . Even though the path in ˆ is simple, these treks can still intersect in G G . Suppose that the edges k p -k p +1 and k q -k q +1 with p &lt; q correspond to treks T 1 and T 2 in G which intersect. There are four cases depending on where this intersection occurs which are depicted in Figure 7. In the first three cases, there is a trek from k p to k q +1 through the intersection vertex - corresponding to a trek from i to j through m in Figure 7. This allows us to shorten the path in ˆ by G omitting k p + q , . . . , k q . The resulting path in ˆ still has all its vertices contained in G K . The fourth case Figure 7d goes from i over m to l , back up to m and then to j , omitting k from the path in ˆ . G The walk in G traced out by this is a zig-zag with a tail m → l whose endpoint l is still in K . For every pair of treks, multiple of these cases can occur simultaneously and each case can occur multiple times. One can freely choose any such intersection and resolve it according to Figure 7 which results in a shortened path in ˆ . G Iterating this procedure

produces a simple path in ˆ from G i to j whose inner vertices are all contained in K . The union of the corresponding treks in G is a zig-zag with tails Z ⊆ G as in Figure 6.

̸

̸

By construction the subgraph Z has the following property: for each collider m on Z , there is tail starting at m and completely contained in Z whose end vertex is in K . Hence, Proposition 4.10 applies and shows that i ⊥ ⊥ j | K ′ on M Z where K ′ is the set of vertices on Z which are in K . Since Z is a subgraph of G , Lemma 4.2 implies that i ⊥ ⊥ j | K on M G . □

## Acknowledgements

Tobias Boege was partially supported by the Academy of Finland grant number 323416 and by the Wallenberg Autonomous Systems and Software Program (WASP) funded by the Knut and Alice Wallenberg Foundation. Benjamin Hollering was supported by the Alexander von Humboldt Foundation. Mathias Drton, Pratik Misra, and Daniela Schkoda acknowledge support from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement No 883818). Sarah Lumpp acknowledges support by the Munich Data Science Institute (MDSI) at Technical University of Munich as part of the focus topic CAUSE. Daniela Schkoda and Sarah Lumpp were also supported by the DAAD programme Konrad Zuse Schools of Excellence in Artificial Intelligence, sponsored by the Federal Ministry of Education and Research.

Part of this research was performed while some of the authors were visiting the Institute for Mathematical and Statistical Innovation (IMSI), which is supported by the National Science Foundation (Grant No. DMS-1929348). Benjamin Hollering and Tobias Boege also thank the Combinatorial Coworkspace (funded by Combinatorial Synergies SPP 2458).

We thank Shiva Darshan for helpful discussions regarding the proof of Theorem 2.1.

## References

- [1] Rabi Bhattacharya and Edward C. Waymire. Continuous parameter Markov processes and stochastic differential equations , volume 299 of Graduate Texts in Mathematics . Springer, Cham, 2023.
- [2] Danai Deligeorgaki, Alex Markham, Pratik Misra, and Liam Solus. Combinatorial and algebraic perspectives on the marginal independence structure of Bayesian networks. Algebraic Statistics , 14(2):233-286, 2023.
- [3] Philipp Dettling, Mathias Drton, and Mladen Kolar. On the lasso for graphical continuous Lyapunov models. 3rd Conference on Causal Learning and Reasoning (CLeaR), Proceedings of Machine Learning Research , 236:514-550, 2024.
- [4] Philipp Dettling, Roser Homs, Carlos Am´ndola, Mathias Drton, and Niels Richard Hansen. Identifiability e in continuous Lyapunov models. SIAM Journal on Matrix Analysis and Applications , 44(4):1799-1821, 2023.
- [5] Joseph Leonard Doob. The Brownian movement and stochastic equations. Annals of Mathematics , 43(2):351-369, 1942.
- [6] Mathias Drton. Algebraic problems in structural equation modeling. In The 50th anniversary of Gr¨bner o bases , volume 77 of Adv. Stud. Pure Math. , pages 35-86. Math. Soc. Japan, Tokyo, 2018.
- [7] Mathias Drton and Thomas S. Richardson. Binary models for marginal independence. Journal of the Royal Statistical Society Series B: Statistical Methodology , 70(2):287-309, 02 2008.
- [8] Alan Edelman and Gilbert Strang. Pascal matrices. The American Mathematical Monthly , 111(3):189-197, 2004.
- [9] Katherine Fitch. Learning directed graphical models from Gaussian data, 2019. arXiv:1906.08050.
- [10] Niels Richard Hansen. A trek rule for the Lyapunov equation, 2024. arXiv:2407.21223.
- [11] Wen Huang, Min Ji, Zhenxin Liu, and Yingfei Yi. Steady states of Fokker-Planck equations: I. Existence. Journal of Dynamics and Differential Equations , 27(3-4):721-742, 2015.

| [12]   | Martin Jacobsen. A brief account of the theory of homogeneous Gaussian diffusions in finite dimensions. In Proceedings of the Third Finnish-Soviet symposium on probability theory and mathematical statistics. Turku, Finland, August 13-16, 1991 , pages 86-94. Utrecht: VSP; Moscow: TVP, 1993.   |
|--------|------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [13]   | Steffen Lauritzen and Kayvan Sadeghi. Unifying Markov properties for graphical models. Ann. Statist. 46(5):2251-2278, 2018.                                                                                                                                                                          |
| [14]   | Lars Lorch, Andreas Krause, and Bernhard Sch¨lkopf. o Causal modeling with stationary diffusions. In Proc. International Conference on Artificial Intelligence and Statistics (AISTATS) , pages 1927-1935. PMLR, April 2024.                                                                         |
| [15]   | Marloes Maathuis, Mathias Drton, Steffen Lauritzen, and Martin Wainwright, editors. Handbook of graphical models . Chapman & Hall/CRC Handbooks of Modern Statistical Methods. CRC Press, 2019.                                                                                                      |
| [16]   | Frantiˇek s Mat´ˇ. us Conditional independence structures examined via minors. Annals of Mathematics and Artificial Intelligence , 21(1):99-30, 1997.                                                                                                                                                |
| [17]   | Jonas Peters, Dominik Janzing, and Bernhard Sch¨lkopf. o Elements of causal inference . Adaptive Computation and Machine Learning. MIT Press, Cambridge, MA, 2017.                                                                                                                                   |
| [18]   | H. Risken. The Fokker-Planck equation , volume 18 of Springer Series in Synergetics . Springer-Verlag, Berlin, second edition, 1989. Methods of solution and applications.                                                                                                                           |
| [19]   | Martin Rohbeck, Brian Clarke, Katharina Mikulik, Alexandra Pettet, Oliver Stegle, and Kai Ueltzh¨ffer. o Bicycle: Intervention-based causal discovery with cycles. 3rd Conference on Causal Learning and Reasoning, Proceedings of Machine Learning Research , 236:209-242, 2024.                    |
| [20]   | Milan Studen´. y Probabilistic Conditional Independence Structures . Information Science and Statistics. Springer, 2005.                                                                                                                                                                             |
| [21]   | Seth Sullivant. Algebraic geometry of Gaussian Bayesian networks. Advances in Applied Mathematics 40(4):482-513, 2008.                                                                                                                                                                               |
| [22]   | Seth Sullivant. Algebraic statistics , volume 194 of Graduate Studies in Mathematics . American Mathe- matical Society, Providence, RI, 2018.                                                                                                                                                        |
| [23]   | Seth Sullivant, Kelli Talaska, and Jan Draisma. Trek separation for Gaussian graphical models. The Annals of Statistics , 38(3):1665-1685, 2010.                                                                                                                                                     |
| [24]   | Gherardo Varando and Niels Richard Hansen. Graphical continuous Lyapunov models. Proceedings of the 36th Conference on Uncertainty in Artificial Intelligence (UAI), Proceedings of Machine Learning Research , 124:989-998, 2020.                                                                   |
| [25]   | Lingfei Wang, Nikolaos Trasanidis, Ting Wu, Guanlan Dong, Michael Hu, Daniel E. Bauer, and Luca Pinello. Dictys: Dynamic gene regulatory network dissects developmental continuum with single-cell                                                                                                   |

## Appendix A. Proofs

A.1. Proof of Proposition 4.5. In this subsection we provide a proof of Proposition 4.5 by constructing a specific submodel of the Lyapunov model and then inductively showing that the determinant of the submatrix Σ iK,jK does not vanish for a trek T with endpoints i and j where the conditioning set K = V ( T ) \ { i, j } .

Let T be a trek with endpoints i and j such that the length of the left branch, whose sink is i , is ℓ and the length of the right branch, whose sink is j , is r . Let K = V ( T ) \{ i, j } be the set of all other nodes on T . Let M ⊆M ′ T T of the Lyapunov model on T obtained from the trek rule in Proposition 4.3 by setting every edge parameter m qp = 1 for all p → q ∈ E ( T ). In other words

$$\begin{array} {$$

̸

̸

Before we begin, note that every entry of every Σ ∈ M ′ T is divisible by 2 ζ . In the following sequence of lemmas, we will instead consider a slight alteration of the above parameterization given by

This follows directly from Proposition 4.3 since every trek monomial m T = ∏ e ∈ T m e = 1. Thus we can think of any matrix Σ ∈ M ′ T as a matrix whose entries are univariate polynomials in ζ . Our goal is now to show that for almost all Σ ∈ M ′ T , | Σ iK,jK | = 0. To do this, we will focus on computing the lowest degree term in | Σ iK,jK | as a polynomial in ζ . Formally, the lowest degree term of a univariate polynomial f = ∑ n k = k 0 f ζ k k with f k 0 = 0 is low( f ) := f k 0 ζ k 0 .

$$\mathcal { M } _ { \mathcal { T } } ^ { \prime \prime }$$

̸

This simply corresponds to factoring out the common 2 ζ from each entry of Σ which does not change whether or not a minor of Σ vanishes. Thus if Σ exhibits a conditional dependence then 2 ζ Σ will exhibit the same conditional dependence. To prove that our original model M T satisfies i ⊥ ⊥ j | K , it suffices to prove it for a matrix in M ′′ T .

The following easy facts are used throughout the proofs in this section:

- · The lowest degree term is multiplicative: low( fg ) = low( f ) low( g ).
- · Since T is a trek, there is a unique shortest trek between any two vertices and thus there is a well-defined notion of distance on the trek.
- · Let the shortest trek between i and j have length ( ℓ, r ). Since it is unique, we have low( σ ij ) = ( ℓ + r ℓ ) ζ ℓ + r for every Σ = Σ( ζ ) ∈ M ′′ T .
- · An important special case of this formula occurs when i and j are on the same branch of T . In that case ℓ = 0 or r = 0 and the binomial coefficient is 1.

We begin with the following lemma concerning the lowest degree terms of principal minors.

Lemma A.1. Let T be a trek and Σ ∈ M ′′ T . For any K ⊆ V we have low( Σ | K | ) = 1 .

Proof. Under the parametrization of M ′′ T , any diagonal entry σ ii has low( σ ii ) = 1 coming from the empty trek on i , while every off diagonal entry has deg(low( σ ij )) ≥ 1. The result then follows immediately from the Leibniz formula for the determinant. □

In the following, we compute minors of the form | Σ iK,jK | which are decisive for conditional independence by Corollary 3.4. If K = V \ { i, j } , then this minor is, up to sign, the adjoint entry adj(Σ) ij . To facilitate inductive proofs via Schur's formula (which computes | Σ iK,jK | using entries of adj(Σ K )), it is important to not only show that these minors are nonzero, but also to consistently compute the signs.

Convention A.2. Throughout this section we assume the following numbering of the vertices V = { 1 2 , , . . . , n } of the trek: the leaf on the left branch is labeled 1, its parent 2, and so on, until the top node is labeled ℓ +1; from there we descend down with ℓ +2 and so on until the leaf on the right branch is labeled n = ℓ + r +1. Any minor | Σ I,J | is computed by ordering the rows and columns of the matrix with respect to this ordering.

Lemma A.3. Let T be a trek and Σ ∈ M ′′ T . Let i, j ∈ V ( T ) be distinct and ( ℓ, r ) be the length of the shortest trek between them. For any K ⊆ V ( T ) \ { i, j } , either | Σ iK,jK | = 0 or deg(low( Σ | iK,jK | )) ≥ ℓ + r . If i and j are adjacent, then low( Σ | iK,jK | ) = ζ .

Proof. Consider the Leibniz formula

$$| \Sigma _ { i K, j K } | = \sum _ { \substack { \pi \colon i K \to j K \\ \text{bjective} } } \text{sgn} ( \pi ) \prod _ { k \in i K } \sigma _ { k \pi ( k ) }.$$

Note that Σ iK,jK is a matrix with a particular ordering on its rows and columns. This ordering fixes a bijection between iK and jK with respect to which sgn( π ) is well-defined.

It suffices to show that each summand is either zero or has lowest term of degree ≥ ℓ + . r Fix a bijection π : iK → jK . The image π i ( ) either equals j or is an element k 1 ∈ K . In the latter case, we may apply π again to get π k ( 1 ) = π 2 ( ) i which either equals j or k 2 ∈ K . Continuing in this way, collect all iterated images of i into a set A ⊆ iK including i and excluding the j which is eventually reached. Observe that

$$\deg \log \left ( \prod _ { k \in A } \sigma _ { k \pi ( k ) } \right ) = \sum _ { k \in A } \deg \log \left ( \sigma _ { k \pi ( k ) } \right ) \geq \log \left ( \sigma _ { i j } \right ) = \ell + r,$$

since the sequence low( σ iπ i ( ) ) , low( σ π i π ( ) 2 ( ) i ) , . . . , low( σ π m ( ) i j ) given by A encodes an undirected path from i to j which is at least as long as the shortest trek.

The lower bound in (A.1) can only be attained with equality if A is completely contained in the shortest trek T ′ connecting i and j and if the successive images π i , π ( ) 2 ( ) i , . . . are ordered such that π i ( ) is the closest node to i , π 2 ( ) the second closest and so on; while all i elements of K \ A are fixed points of π and do not contribute to the degree. In this case, the σ kk , k ∈ K \ A , have lowest term 1 and the shortest treks between adjacent members of the sequence i, π ( ) i , π 2 ( ) i , . . . , j give a disjoint decomposition of T ′ hence achieving degree exactly ℓ + r in total. If i and j are adjacent, then ℓ + r = 1 and the only bijection π which satisfies this sends i to j directly and fixes every other k ∈ K . Since i and j are adjacent, they appear in the same place in the ordering of iK and jK , respectively, hence π corresponds to the identity permutation of the rows and columns and thus we compute low( Σ | iK,jK | ) = low( σ ij ) = ζ . □

As the previous proof shows, the lower degree bound ℓ + r is achievable and the Leibniz formula always has summands of this degree. Whether or not these summands cancel to

. . .

Figure 8. The four cases in Lemma A.5.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00583v1-figure-6833f5f0-2948-45ae-9614-d904546e345b.jpg)

cause a bump in lowest degree depends on the locations of i and j on the trek and to each other, as well as the conditioning set K .

Convention A.4. In the following proofs, we will use the Schur complement expansion of the determinant:

$$\ s g n _ { K } ( i, j ) | \Sigma _ { i K, j K } | = \sigma _ { i j } | \Sigma _ { K } | - \Sigma _ { i, K } \text{ adj} ( \Sigma _ { K } ) \Sigma _ { K, j } = \sigma _ { i j } | \Sigma _ { K } | - \sum _ { p, q \in K } \sigma _ { i p } \hat { \sigma } _ { p q } \sigma _ { q j }, \quad ( A. 2 )$$

̸

where ˆ σ pq = adj(Σ K pq ) . To use Schur's formula to expand the determinant by the ( i, j )-entry of Σ iK,jK as indicated, we have to permute the i th row as well as the j th column to be the first row and column, respectively. The symbol sgn K ( i, j ) denotes the sign change in the determinant incurred by this. This sign depends only on the difference of the positions of i in iK and j in jK . We will assume below that K contains all interior nodes on the shortest trek T ′ between i and j . In this case and if i = , sgn j K ( i, j ) can be computed as

$$\text{sgn} _ { K } ( i, j ) = ( - 1 ) ^ { n ^ { \prime } } = ( - 1 ) ^ { | j - i | + 1 },$$

where n ′ is the number of nodes on T ′ , which equals | j - | i +1 under Convention A.2.

Lemma A.5. Let T be a trek and Σ ∈ M ′′ T . Let i, j ∈ V ( T ) be distinct, non-adjacent vertices on the same branch of T . Denote the length of the shortest path T ′ between i and j by ℓ &gt; 1 . If K ⊆ V ( T ) \ { i, j } contains all interior nodes on T ′ , then either | Σ iK,jK | = 0 or deg low( Σ | iK,jK | ) &gt; ℓ .

Proof. Without loss of generality, we may suppose by symmetry that both i and j are on the left branch of the trek and i &lt; j . Then the shortest path connecting them is also the unique shortest trek connecting them, which has length ( ℓ, 0).

We perform induction on ℓ ≥ 2. The case ℓ = 2 can be handled using the methods of Lemma A.3. The shortest trek T ′ between i and j is the form i → k → j where k ∈ K by assumption. Using the terminology of the proof of Lemma A.3, the only two bijections iK → jK which contribute to the ζ 2 term are:

. . .

↦

- (1) π 1 sending i → j and fixing every other element. Since k is between i and j in the ordering, the element i in iK occupies the same place as k in jK ; and k in iK the same place as j in jK . Hence π 1 corresponds to an adjacent transposition of the rows and columns and has sgn( π 1 ) = -1. This summand is thus -ζ 2 .

↦

- (2) π 2 sending i → k → j and fixing every other element. By similar reasoning as in the previous case, this corresponds to the identity permutation of the rows and columns and thus the contribution is ζ 2 .

↦

The two summands cancel and therefore either the entire determinant vanishes or the lowest term has degree &gt; 2.

For the induction step, we use the Schur complement expansion (A.2). By Lemma A.1 and the definition of M ′′ T , we know low( σ ij | Σ K | ) = ζ ℓ . The goal is to show that for each choice of p, q ∈ K , either σ ip ˆ σ pq σ qj = 0 or deg low( σ ip ˆ σ pq σ qj ) ≥ ℓ and in case of equality to compute the coefficient of the ζ ℓ term. We distinguish four cases (up to exchanging p and q or switching branches) depending on the locations of p and q with respect to i and j :

- (a) In case of Figure 8a, notice that the shortest trek from i to p strictly includes the entire path from i to j , so deg(low( σ ip ˆ σ pq σ qj )) ≥ deg(low( σ ip )) &gt; ℓ .
- (b) In case of Figure 8b, the reason is similar. Lemma A.3 shows deg(low( σ ip ˆ σ pq )) &gt; ℓ .
- (c) In case of Figure 8c, we know by Lemma A.3 that deg(low(ˆ σ pq )) &gt; ℓ .
- (d) The last case is depicted in Figure 8d. Then low( σ ip ) = ζ a and low( σ qj ) = ζ b and a + b + c = . So ℓ

$$\deg ( \text{low} ( \sigma _ { i p } \hat { \sigma } _ { p q } \sigma _ { q j } ) ) = \ell \iff \deg ( \text{low} ( \hat { \sigma } _ { p q } ) ) = c.$$

Notice that p and q are on the same branch of T . If they coincide ( c = 0), then low(ˆ σ pq ) = 1 by Lemma A.1; if they are adjacent ( c = 1), Lemma A.3 yields ˆ σ pq = -| Σ K \{ } q ,K \{ p } | = -ζ . Otherwise, the induction hypothesis ensures ˆ σ pq = 0 or deg low(ˆ σ pq ) &gt; c , which do not contribute to the lowest term. This implies that

$$\Phi \left ( \sum _ { p, q } \sigma _ { i p } \hat { \sigma } _ { p q } \sigma _ { q j } \right ) & = \sum _ { \substack { p, q \\ c = 0 } } \left ( \zeta ^ { a } \cdot 1 \cdot \zeta ^ { b } \right ) + \sum _ { \substack { p, q \\ c = 1 } } \left ( \zeta ^ { a } \cdot \left ( - \zeta \right ) \cdot \zeta ^ { b } \right ) \\ & = \ell \zeta ^ { \ell } - ( \ell - 1 ) \zeta ^ { \ell } = \zeta ^ { \ell }.$$

Thus the only case which contributes to the ζ ℓ term in | Σ iK,jK | is (d). Its contribution cancels with low( σ ij | Σ K | ) = ζ ℓ in (A.2), implying the result. □

̸

Proposition A.6. Let T be a trek and i, j ∈ V ( T ) be nodes on different branches of T (in particular neither of them is the top node). Denote the length of the shortest path T ′ between i and j by ( ℓ, r ) . If K ⊆ V ( T ) \ { i, j } contains all interior nodes on T ′ , then there exists Σ ∈ M ′′ T such that i ⊥ ⊥ j | K .

Proof. Without loss of generality, we may assume that i &lt; j in the ordering, so that i is on the left branch and j on the right. We claim that low( Σ | iK,jK | ) = ( -1) ℓ + +1 r ( ℓ + r -2 ℓ -1 ) ζ ℓ + r which shows that | Σ iK,jK | is non-zero and yields conditional dependence for all but a finite number of values of ζ .

Figure 9. The three cases in Proposition A.6.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00583v1-figure-fb363930-f9d4-4e3d-a9ac-a3ba9f4be86f.jpg)

We prove this statement by induction on ℓ, r ≥ 1, again using the Schur complement expansion (A.2),

$$( - 1 ) ^ { \ell + r + 1 } | \Sigma _ { i K, j K } | = \sigma _ { i j } | \Sigma _ { K } | - \Sigma _ { i, K } \text{ adj} ( \Sigma _ { K } ) \Sigma _ { K, j } = \sigma _ { i j } | \Sigma _ { K } | - \sum _ { p, q \in K } \sigma _ { i p } \hat { \sigma } _ { p q } \sigma _ { q j }$$

with ˆ σ pq = adj(Σ K pq ) . By Lemma A.1 and the definition of M ′′ T , we have low( σ ij | Σ K | ) = ( ℓ + r ℓ ) ζ ℓ + r . The other part of the expansion is a sum over σ ip ˆ σ pq σ qj which we deal with term by term, so fix p, q ∈ K . If p is strictly closer to j than q is (or, symmetrically, q is closer to i than p is), then deg low( σ ip σ qj ) &gt; deg low( σ ip σ pj ) ≥ deg low( σ ij ) = ℓ + . Additionally, r by Lemma A.3, p and q must be contained in V ( T ′ ), i.e. i &lt; p ≤ q &lt; j in our ordering of vertices by Convention A.2. Hence, to find the lowest-order terms, we only need to consider

$$p, q \in V ( \mathcal { T } ^ { \prime } ) \, \smallsetminus \, \{ i, j \} \subseteq K \text{ such that } i < p \leq q < j.$$

With these remarks, the basis of induction is easy to show. For ℓ = r = 1 the only possible option in (A.4) is that p = q = top node of T . Then low( σ ip ˆ σ pq σ qj ) = ζ · 1 · ζ = ζ 2 . With low( σ ij | Σ K | ) = 2 ζ 2 we get low( Σ | iK,jK | ) = ( -1) 3 · (2 ζ 2 -ζ 2 ) = -ζ 2 , as required.

In the general case, the reduction via (A.4) allows us to identify all relevant choices of p, q with decompositions of T ′ into three disjoint parts: the trek i to p , the trek p to q and the trek q to j . There are three cases:

- (a) Suppose that p and q are both on the left branch of T ′ and suppose that they are not both equal to the top node. Then the trek decomposes into four paths as in Figure 9a. We have low( σ ip σ qj ) = ( c + r c ) ζ a + + c d . By Lemma A.5, the lowest exponent of ˆ σ pq in this case is b if and only if b = 0 or b = 1. Otherwise it is larger than b and a + b + c + d &gt; ℓ + r disqualifying those cases from contributing to the lowest terms of the overall expression, so assume b = 0 or b = 1. Then low(ˆ σ pq ) = ( -1) b ζ b by Lemma A.3. In total, this case contributes ∑ ℓ -1 c =1 ( c + r c ) when b = 0 and -∑ ℓ -2 c =0 ( c + r c ) when b = 1 to the lowest-order terms. (The different summation ranges for c reflect that we must have a ≥ 1 and a + c = ℓ and also that if b = 0, in order for p = q not to be the top node, we must have c ≥ 1 whereas if b = 1, we only need c ≥ 0.)

Together, this gives ( ℓ + r -1 ℓ -1 ) -1. The case where p and q are both on the right branch is analogous and contributes ( ℓ + r -1 r -1 ) = ( ℓ + r -1 ℓ ) -1.

- (b) Suppose that p is on the left branch and q is on the right branch but neither of them is the top node. Then the trek decomposes into four paths as in Figure 9b and low( σ ip σ qj ) = ζ a + d . Here the induction hypothesis applies and computes | Σ pL,qL | = ( -1) q -p +1 ( b + c -2 b -1 ) ζ b + c , where L = K \ { p, q } . (Note that L contains all interior nodes on the shortest trek between p and q since K had the same property for i and j .) However, there is another sign to account for coming from the cofactor: ˆ σ pq = adj(Σ K pq ) = ( -1) q -p | Σ pL,qL | = -( b + c -2 b -1 ) ζ b + c . The sign only depends on q -p , again because L contains all vertices between them in the ordering. In total, this case contributes -∑ ℓ -1 b =1 ∑ r -1 c =1 ( b + c -2 b -1 ) = 1 -( ℓ + r -2 ℓ -1 ) to the lowest term.
- (c) There is a special case when p and q both coincide with the top node of the trek, as in Figure 9c. In this case, they are both simultaneously on the left and the right branch but this case was not counted in the previous bullet points. As in the base case of the induction, the contribution is low( σ ip ˆ σ pq σ qj ) = ζ ℓ + r .

In total,

$$& \text{In total}, \\ & \quad ( - 1 ) ^ { \ell + r + 1 } \text{low} ( | \Sigma _ { i K, j K } | ) \\ & \quad = \left ( \underbrace { \binom { \ell + r } { \ell } - \binom { \ell + r - 1 } { \ell - 1 } - \binom { \ell + r - 1 } { \ell } + 2 - \left ( 1 - \binom { \ell + r - 2 } { \ell - 1 } \right ) - 1 \right ) } _ { = 0 } + 2 - \left ( 1 - \binom { \ell + r - 2 } { \ell - 1 } \right ) - 1 \right ) \zeta ^ { \ell + r } \\ & \quad = \binom { \ell + r - 2 } { \ell - 1 } \zeta ^ { \ell + r }, \\ & \text{as claimed}.$$

as claimed.

□

The last remaining case needed to complete the proof of Proposition 4.5 is the case when T is a trek such that ℓ = 0 or r = 0, meaning that T is actually a path. In this case Lemma A.5 shows that the lowest degree term can have degree strictly higher than ℓ + . For r this reason, we instead analyze the highest degree term in ζ which appears in the expansion of | Σ iK,jK | . Let high( f ) denote the highest degree non-zero term of a univariate polynomial f .

̸

Lemma A.7. Let T be a directed path from i to j and K = V ( T ) \ { i, j } be the remaining vertices. Then there exists Σ ∈ M ′′ T such that i ⊥ ⊥ j | K .

̸

̸

Proof. To show that there exist Σ ∈ M ′′ T for which i ⊥ ⊥ j | K , it is enough to show that | Σ iK,jK | as a univariate polynomial in ζ has a non-zero coefficient. To this end, we compute the determinant of high(Σ iK,jK ), which is the matrix obtained by taking the highest degree term of each entry of Σ iK,jK , at ζ = 1. If this value is non-zero, it follows that | high(Σ iK,jK ) | = 0 is the highest degree term of | Σ iK,jK | which is thus non-zero as well.

We claim that | high(Σ iK,jK ) | = 1 at ζ = 1. Recall that M ′′ T is defined by the parametrization which maps σ ij to the sum of ζ ℓ + r ( ℓ + r ℓ ) over all possible treks between i and j . Hence, the entries of high(Σ iK,jK ) only take into account the longest treks between nodes (with the highest value of ℓ + ). For simplicity, we number the vertices from 1 to r n with i = 1 and j = n . Now, the longest and unique trek between 1 and k for any k ≤ n is precisely the

directed path from 1 to k . However, the longest trek between a and k for any 1 &lt; a &lt; k is one where the left and right branches are directed paths from 1 to a and 1 to k , respectively. This implies that the term with highest degree for σ 1 k has coefficient 1 (as ℓ = 0) and for that of σ ak has coefficient ( a + k -2 a -1 ) . This gives us that the ( a, b ) entry of high(Σ { 1 2 3 , , ,...,n - } { 1 , 2 3 , ,...,n } ) at ζ = 1 is precisely ( a + b -2 a -1 ) . However, using the identity in [8, Equation (1)], we can obtain an LU decomposition of this matrix where the diagonal terms of both L and U are 1. This implies that | high(Σ iK,jK ) | = 1 at ζ = 1, which proves our claim. □

## A.2. Proof Details for Lemma 4.7.

Proof Details for Lemma 4.7. In the proof of Lemma 4.7, we solve the following six equations resulting from the block-matrix partition for the corresponding blocks of Σ ( m ) for m →∞ :

$$( \mathbf A ) \ M _ { U } ^ { * } \Sigma _ { U } ^ { ( m ) } + \Sigma _ { U } ^ { ( m ) } M _ { U } ^ { * T } + 2 I _ { k - 1 } = 0$$

This equation is uniquely solved by Σ ( m ) U = Σ ∗ U due to (a).

$$( \mathbf B ) \ M _ { U } ^ { * } \Sigma _ { U, k } ^ { ( \bar { m } ) } + m \Sigma _ { U, k - 1 } ^ { ( m ) } - \bar { m } \Sigma _ { U, k } ^ { ( m ) } = 0 \, \stackrel { \cdot } { \lambda }$$

The matrix M ∗ U -mI k -1 is invertible due to the structure of M ∗ , so Lemma A.8 together with (A) yields

$$\lim _ { m \to \infty } \Sigma _ { U, k } ^ { ( m ) } = \lim _ { m \to \infty } \Sigma _ { U, k - 1 } ^ { ( m ) } = \Sigma _ { U, k - 1 } ^ { * }.$$

$$( \text{C} ) \ M _ { U } ^ { * } \Sigma _ { U, W } ^ { ( m ) } + \Sigma _ { U, W } ^ { ( m ) } M _ { \psi ( W ) } ^ { * } + \Sigma _ { U, k } ^ { ( m ) } ( m _ { k, k - 1 } ^ { * }, 0, \dots, 0 ) = 0$$

Taking the limit of the left-hand side for m →∞ yields

$$M _ { U } ^ { * } \Sigma _ { U, W } ^ { ( m ) } + \Sigma _ { U, W } ^ { ( m ) } M _ { \psi ( W ) } ^ { * } + \Sigma _ { U, k } ^ { ( m ) } ( m _ { k, k - 1 } ^ { * }, 0, \dots, 0 )$$

$$\xrightarrow { ( B ) } M _ { U } ^ { * } \left ( \lim _ { m \to \infty } \Sigma _ { U, W } ^ { ( m ) } \right ) + \left ( \lim _ { m \to \infty } \Sigma _ { U, W } ^ { ( m ) } \right ) M _ { \psi ( W ) } ^ { * } + \Sigma _ { U, k - 1 } ^ { * } ( m _ { k, k - 1 } ^ { * }, 0, \dots, 0 ),$$

$$M _ { U } ^ { * } \Sigma _ { U, W } ^ { ( m ) } + \Sigma _ { U, W } ^ { ( m ) } M _ { \psi ( W ) } ^ { * } + \Sigma _ { U, k } ^ { ( m ) } ( m _ { k, k - 1 } ^ { * }, 0, \dots, 0 ) \\ \overrightarrow { m \to \infty } \, M _ { U } ^ { * } \left ( \lim _ { m \to \infty } \Sigma _ { U, W } ^ { ( m ) } \right ) + \left ( \lim _ { m \to \infty } \Sigma _ { U, W } ^ { ( m ) } \right ) M _ { \psi ( W ) } ^ { * } + \Sigma _ { U, k - 1 } ^ { * } ( m _ { k, k - 1 } ^ { * }$$

where we employed (B) in the last summand. This last summand can further be rewritten as

$$\Sigma _ { U, k - 1 } ^ { * } ( m _ { k, k - 1 } ^ { * }, 0, \dots, 0 ) = \Sigma _ { U } ^ { * } M _ { \psi ( W ), U } ^ { * }$$

by attaching a ( k -2) × ( k -1) zero matrix to the existing row vector and thereby extending it to M ∗ ψ W ,U ( ) T . Setting the limit to zero yields equation (b) which is solved

$$( D ) \ 2 \,.$$

$$\begin{array} { c } \lambda & \psi _ { W } w _ { J }, \upsilon & \lambda \\ \text{by $\Sigma_{U,\psi(W)}^{*}$} \\ \text{$\text{$This gives us} \lim_{m\to\infty} \Sigma_{U,W}^{(m)} = \Sigma_{U,\psi(W)}^{*}$} \\ \text{(D)} & 2 - 2 m \Sigma _ { k, k } ^ { ( m ) } + m \left ( \Sigma _ { k, k - 1 } ^ { ( m ) } + \Sigma _ { k - 1, k } ^ { ( m ) } \right ) = 0 \end{array}$$

Solving for Σ ( m ) k,k , letting m go to infinity, and using the result from (B) yields

$$\lim _ { m \to \infty } \Sigma _ { k, k } ^ { ( m ) } = \lim _ { m \to \infty } \left ( \frac { 1 } { m } + \frac { \Sigma _ { k, k - 1 } ^ { ( m ) } + \Sigma _ { k - 1, k } ^ { ( m ) } } { 2 } \right ) = \frac { \Sigma _ { k - 1, k - 1 } ^ { * } + \Sigma _ { k - 1, k - 1 } ^ { * } } { 2 } = \Sigma _ { k - 1, k - 1 } ^ { * }.$$

$$( \text{E} ) \ m \Sigma _ { k - 1, W } ^ { ( m ) } - m \Sigma _ { k, W } ^ { ( m ) } + \Sigma _ { k, k } ^ { ( m ) } ( m _ { k, k - 1 } ^ { * }, 0, \dots, 0 ) + \Sigma _ { k, W } ^ { ( m ) } M _ { \psi ( W ) } ^ { * } = 0$$

The matrix M ∗ ψ W ( ) T -mI n -k is invertible due to the structure of M ∗ , so the transposed version of Lemma A.8 together with (C) and (D) yields

$$\lim _ { m \to \infty } \Sigma _ { k, W } ^ { ( m ) } = \lim _ { m \to \infty } \Sigma _ { k - 1, W } ^ { ( m ) } = \Sigma _ { k - 1, \psi ( W ) } ^ { * }.$$

$$\begin{array} { c c } ( \mathbf F ) & M _ { \psi ( W ) } ^ { * } \Sigma _ { W } ^ { ( m ) } + \Sigma _ { W } ^ { ( m ) } M _ { \psi ( W ) } ^ { * } + 2 I _ { n - k } \\ & + ( m _ { k, k - 1 } ^ { * }, 0, \dots, 0 ) ^ { T } \Sigma _ { k, W } ^ { ( m ) } + \Sigma _ { W, k } ^ { ( m ) } ( m _ { k, k - 1 } ^ { * }, 0, \dots, 0 ) = 0 \end{array}$$

We use that Σ ( m ) k,W → Σ ∗ k -1 ,ψ W ( ) as well as Σ ( m ) W,k → Σ ∗ ψ W ,k ( ) -1 for m →∞ due to (E) and take the limit on both sides. This yields for the left-hand side

$$\text{and take the main on your sites. } I m s y e r s o r o n e t e r e t e n a s i n s u e \\ & M _ { \psi ( W ) } ^ { * } \Sigma _ { W } ^ { ( m ) } + \Sigma _ { W } ^ { ( m ) } M _ { \psi ( W ) } ^ { * } + 2 I _ { n - k } \\ & + ( m _ { k, k - 1 } ^ { * }, 0, \dots, 0 ) ^ { T } \Sigma _ { k, W } ^ { ( m ) } + \Sigma _ { W, k } ^ { ( m ) } ( m _ { k, k - 1 } ^ { * }, 0, \dots, 0 ) \\ \frac { ( E ) } { m \to \infty } \, M _ { \psi ( W ) } ^ { * } \left ( \lim _ { m \to \infty } \Sigma _ { W } ^ { ( m ) } \right ) + \left ( \lim _ { m \to \infty } \Sigma _ { W } ^ { ( m ) } \right ) M _ { \psi ( W ) } ^ { * } + 2 I _ { n - k } \\ & + ( m _ { k, k - 1 } ^ { * }, 0, \dots, 0 ) ^ { T } \Sigma _ { k - 1, \psi ( W ) } ^ { * } + \Sigma _ { \psi ( W ), k - 1 } ^ { * } ( m _ { k, k - 1 } ^ { * }, 0, \dots, 0 ) \\ = \, M _ { \psi ( W ) } ^ { * } \left ( \lim _ { m \to \infty } \Sigma _ { W } ^ { ( m ) } \right ) + \left ( \lim _ { m \to \infty } \Sigma _ { W } ^ { ( m ) } \right ) M _ { \psi ( W ) } ^ { * } + 2 I _ { n - k } \\ & + M _ { \psi ( W ), U } ^ { * } \Sigma _ { U, \psi ( W ) } ^ { * } + \Sigma _ { \psi ( W ), U } ^ { * } M _ { \psi ( W ), U } ^ { * }, \\ \text{where we memoized the same technique of rewriting the last two commands as for } ( C )$$

where we employed the same technique of rewriting the last two summands as for (C). Setting the limit to zero yields equation ?? that is solved by Σ ∗ ψ W ( ) . We conclude

$$\lim _ { m \rightarrow \infty } \Sigma _ { W } ^ { ( m ) } = \Sigma _ { \psi ( W ) } ^ { * }.$$

□

Lemma A.8. Let A ( m ) and C ( m ) be convergent sequences of n × l matrices, and let D,E,F ∈ R n × n be constant matrices such that E -mI n is invertible for all m ∈ N &gt; 0 . If a further matrix sequence B ( m ) satisfies that

$$m \cdot A ^ { ( m ) } - m \cdot B ^ { ( m ) } + D \cdot A ^ { ( m ) } + E \cdot B ^ { ( m ) } + F \cdot C ^ { ( m ) } = 0$$

for all m ∈ N &gt; 0 , then

$$\lim _ { m \to \infty } B ^ { ( m ) } = \lim _ { m \to \infty } A ^ { ( m ) }.$$

Proof. By collecting the terms and dividing by m , we have

$$\left ( \frac { 1 } { m } D + I _ { n } \right ) A ^ { ( m ) } + \left ( \frac { 1 } { m } E - I _ { n } \right ) B ^ { ( m ) } + \frac { 1 } { m } F \cdot C ^ { ( m ) } = 0.$$

Since ( 1 m E -I n ) is invertible, we can write

$$B ^ { ( m ) } = - \left ( \frac { 1 } { m } E - I _ { n } \right ) ^ { - 1 } \left ( \left ( \frac { 1 } { m } D + I _ { n } \right ) A ^ { ( m ) } - \frac { 1 } { m } F \cdot C ^ { ( m ) } \right ).$$

As matrix inversion is continuous and A ( m ) and C ( m ) are convergent, letting m go to infinity on both sides yields the result. □

| Department of Mathematics, KTH Royal Institute of Technology, Sweden Email address : post@taboege.de                                                                        |
|-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| TUM School of Computation, Information and Technology, Technical University of Munich, and Munich Center for Machine Learning, Germany Email address : mathias.drton@tum.de |
| TUM School of Computation, Information and Technology, Technical University of Munich Email address : benjamin.hollering@tum.de                                             |
| TUM School of Computation, Information and Technology, Technical University of Munich Email address : sarah.lumpp@tum.de                                                    |
| TUM School of Computation, Information and Technology, Technical University of Munich Email address : pratik.misra@tum.de                                                   |
| TUM School of Computation, Information and Technology, Technical University of Munich Email address : daniela.schkoda@tum.de                                                |
Tobias Boege Mathias Drton Benjamin Hollering Sarah Lumpp Pratik Misra Daniela Schkoda
Conditional Independence in Stationary Diffusions
Stationary distributions of multivariate diffusion processes have recently
been proposed as probabilistic models of causal systems in statistics and
machine learning. Motivated by these developments, we study stationary
multivariate diffusion processes with a sparsely structured drift. Our main
result gives a characterization of the conditional independence relations that
hold in a stationary distribution. The result draws on a graphical
representation of the drift structure and pertains to conditional independence
relations that hold generally as a consequence of the drift's sparsity pattern.