---
title: Estimate Collapsibility in Causal Graphs
url: https://www.emergentmind.com/topics/estimate-collapsibility
type: topic
---

# Estimate Collapsibility in Causal Graphs

Estimate collapsibility is a notion of estimator invariance under marginalization in causal Bayesian networks and completed partially directed acyclic graphs (CPDAGs). For a causal Bayesian network \(\mathcal{B}=(G,\mathcal{P}(G))\), it is defined by the equality \(\widehat{\mathcal{P}(G_A)}=\widehat{\mathcal{P}(G)_A}\), meaning that the maximum-likelihood estimate obtained from the marginal model induced by a subgraph equals the marginal of the maximum-likelihood estimate from the full model. In the CPDAG setting, the concept is used to reduce a graph around a treatment–outcome pair while preserving the causal effect estimator before and after marginalization; the central structural object is the strong d-convex hull, which characterizes the minimal collapsible set [2606.08941].

## 1. Formal setting and scope

The 2026 formulation places estimate collapsibility within three collapsibility notions for causal Bayesian networks: **CI-collapsibility**, where the independence model is preserved after marginalization; **model collapsibility**, where the distribution family is preserved after marginalization; and **estimate collapsibility**, where the maximum-likelihood estimators are preserved after marginalization [2606.08941]. This separation is central because the target is not merely preservation of conditional independences or model class, but preservation of the estimator actually used for causal effect computation.

The ambient graphical setting is that of DAGs and CPDAGs. A CPDAG \(\mathcal{G}\) represents a Markov equivalence class of DAGs with the same independence model,
\[
I(G_1)=I(G_2).
\]
For a causal Bayesian network, the factorization is
\[
P(x_V;\theta)=\prod_{v\in V} P(x_v\mid pa_G(x_v);\theta_v).
\]
The paper studies how to identify, for a target pair \(\{X,Y\}\), a smallest induced subgraph on which estimation can be performed without changing the resulting estimator [2606.08941].

This formulation makes estimate collapsibility a graph-reduction principle for causal estimation. A plausible implication is that it serves as a local-computation analogue of classical collapsibility in graphical modeling, but specialized to estimator preservation rather than to parameter or independence preservation.

## 2. Causal estimate collapsibility

For general estimation in a causal Bayesian network \(\mathcal{B}=(G,\mathcal{P}(G))\), estimate collapsibility onto a subset \(A\) is defined by
\[
\widehat{\mathcal{P}(G_A)} = \widehat{\mathcal{P}(G)_A}.
\]
This states that fitting the reduced graph \(G_A\) and then estimating is equivalent to fitting the full graph and marginalizing afterward [2606.08941].

For causal effect estimation in a DAG, the paper defines a stronger causal version. Let \(X\) and \(Y \notin pa_G(X)\) be non-adjacent, and let \(R\subseteq V\) contain \(X,Y\). Then \(G\) is **causal estimate collapsible onto \(R\)** if, for every non-empty valid back-door adjustment set \(Z\) in \(G_R\),
\[
\hat{P}(Y \mid do(X=x)) = \sum_z \hat{P}_{G_R}(Y \mid X=x, Z=z)\,\hat{P}_{G_R}(Z=z).
\]
The paper also gives the analogous empty-adjustment case: if \(Z=\emptyset\) is valid in \(G_R\), then causal estimate collapsibility holds when \(R\) is the strong d-convex hull of \(X\) and \(Y\) in \(G_{\bar X}\) [2606.08941].

The causal estimand is the average causal effect
\[
ACE = E(Y\mid do(X=1)) - E(Y\mid do(X=0)),
\]
and the back-door formula used in the reduced graph is
\[
P(Y\mid do(X=x))=\sum_z P(Y\mid X=x,Z=z)P(Z=z).
\]
The framework therefore preserves not only local likelihood-based estimation but also the causal intervention estimate built from valid adjustment in the reduced graph [2606.08941].

A key point is that estimate collapsibility is defined at the level of the estimator. This distinguishes it from criteria that preserve only structural or probabilistic properties of the graph.

## 3. Strong d-convexity and minimal collapsible sets

The main structural result is that minimal collapsible sets are exactly strong d-convex hulls. A subset \(R \subseteq V\) is **strongly d-convex** if:

1. **No inducing path for \(R\)** exists in \(G\).
2. \(Ch_G(M)\cap An_G(R)\) is **linearly ordered with respect to \(R\)**, where \(M=V\setminus R\),
   \[
   Ch_G(M)=ch_G(M)\cup M,\qquad An_G(R)=an_G(R)\cup R.
   \]

The paper defines the linearly ordered condition as follows: for any node in that set, any two distinct parents are adjacent unless both parents lie in \(R\) [2606.08941].

For \(R\subseteq V\), an induced subgraph \(G_H\) is the **strong d-convex hull** of \(R\) if it is the **smallest** strongly d-convex subgraph containing \(R\). The central theorem states:
\[
H \text{ is the minimal collapsible set containing } R \iff G_H \text{ is the strong d-convex hull of } R.
\]
This provides the exact graph-theoretic characterization of estimate collapsibility [2606.08941].

Several auxiliary results organize this characterization. The paper gives the equivalence
\[
\hat{P}(G_R)=\hat{P}(G)_R \iff \exists G' \in [G]\text{ such that }V\setminus R \text{ is terminal in }G' \iff R \text{ is strongly d-convex},
\]
thereby linking estimator preservation, Markov-equivalent orientations, and strong d-convexity. It also states that for a DAG \(G\) and \(R\subseteq V\), the following are equivalent:

1. \(R\) is d-convex.
2. \(\mathcal{MF}_G(R)=\emptyset\), where \(\mathcal{MF}_G(R)\) is the set of marginalization forbidden pairs.
3. There are no inducing structures of \(R\) in \(G_R\).

In addition, strong d-convexity is hereditary inside a strongly d-convex induced subgraph, intersections of strongly d-convex supersets remain strongly d-convex, and both the d-convex hull \(ch(R)\) and strong d-convex hull \(sch(R)\) exist and are unique [2606.08941].

These results make the reduced estimation domain canonical rather than heuristic: once \(R\) is fixed, its minimal collapsible completion is uniquely determined.

## 4. Algorithms for identifying collapsible sets

The paper gives three algorithms for computing the relevant hulls in DAGs, later generalized to CPDAGs [2606.08941].

| Algorithm | Output | Complexity |
|---|---|---|
| CVM(\(G,R\)) | Vertices on shortest inducing paths between non-adjacent vertices in \(R\) | \(O(k^2|V|^2)\) |
| ICHA(\(G,R\)) | A d-convex hull \(H\) | \(O(k^2|V|^3)\) |
| ISCHA(\(G,R\)) | Strong d-convex hull \(H\) | \(O(k^2|V|^4)\) |

CVM(\(G,R\)) starts with \(M=V\setminus R\), computes
\[
R_1 = mb_G(M)\cap R,
\]
and, for each non-adjacent pair \(r_1,r_2\) in \(R_1\), forms the moralized ancestral graph
\[
G^* = (G_{An_G(\{r_1,r_2\})})^m,
\]
restricts it to \(\{r_1,r_2\}\cup M\), and collects the vertices on shortest paths. The resulting set records vertices on minimal inducing structures [2606.08941].

ICHA(\(G,R\)) iterates CVM: it initializes \(H\leftarrow R\), repeatedly updates
\[
\mathrm{AIP} \leftarrow \mathrm{CVM}(G,H),\qquad H \leftarrow H\cup \mathrm{AIP},
\]
and stops when no new vertices are added. This yields the d-convex hull. ISCHA(\(G,R\)) then strengthens ICHA by enforcing the linear-ordering condition. After updating \(H\) by ICHA, it computes
\[
PA \gets \{pa_G(w)\mid w\in Ch_G(V\setminus H)\cap An_G(H),\ w \text{ is not linearly ordered}\},
\]
adds those parents to \(H\), and repeats until \(PA=\emptyset\) [2606.08941].

The algorithmic significance is twofold. First, the minimal collapsible set is computable directly from graph structure. Second, the procedure is constructive enough to be embedded into causal effect estimation workflows rather than used as a purely structural diagnostic.

## 5. Reduction of CPDAG-based causal estimation

The strong d-convex hull is integrated with IDA, the framework “Intervention calculus when the DAG is Absent,” to form a reduced-graph estimation procedure for CPDAGs [2606.08941]. The modified method proceeds as follows:

1. Obtain an arbitrary DAG \(G\in[\mathcal{G}]\) using Meek’s rules.
2. Compute
   \[
   R_{\text{init}} = \mathrm{ISCHA}(G,\{X,Y\}).
   \]
3. For each DAG \(G\) consistent with \(\mathcal{G}\), set \(R=R_{\text{init}}\) and \(Z=pa_{G_R}(X)\).
4. If \(Z=\emptyset\), build \(G_{\bar X}\), recompute \(R\), and update \(Z\).
5. If \(Z=\emptyset\), estimate using \(P(Y\mid X)\); otherwise estimate by back-door adjustment:
   \[
   \sum_z P(Y\mid X,Z=z)P(Z=z).
   \]
6. Return the multiset of possible causal effects.

This procedure reduces the size of the graph on which IDA is run, the number of equivalent DAGs that must be enumerated, and the number of variables entering adjustment sets. The paper states that the local estimator is **exactly identical** to the global IDA estimator when the strong d-convex hull is used [2606.08941].

A further invariance theorem extends the construction from DAGs to CPDAGs: if \(G\) is any DAG in the Markov equivalence class of a CPDAG \(\mathcal{G}\), then for any \(R\supseteq\{X,Y\}\),
\[
R \text{ is causal estimate collapsible in } \mathcal{G} \iff R \text{ is causal estimate collapsible in } G.
\]
This result is essential because it shows that collapsibility is not an artifact of one representative DAG but a property stable across the entire equivalence class [2606.08941].

In effect, estimate collapsibility turns a global equivalence-class problem into a local induced-subgraph problem without changing the causal answer returned by IDA.

## 6. Relation to other collapsibility frameworks

Estimate collapsibility belongs to a broader family of collapsibility notions, but it is distinguished by its estimator-preservation target. In hierarchical log-linear graphical models, for example, collapsibility onto a target set \(A\) is defined by
\[
\hat p(i_A)=\hat p_A(i_A),
\]
and a graphical model is collapsible onto \(A\) if and only if \(A\) contains every minimal separator between every pair of non-adjacent vertices in \(A\); equivalently, \(A\) contains at least one minimal \(xy\)-separator for every non-adjacent pair \(x,y\in A\) [2510.09024]. This separator characterization is closely related in spirit to the strong d-convex hull characterization, but the objects preserved are marginal model estimates rather than causal effect estimators in CPDAGs.

In the counterfactual causal framework, collapsibility is defined for effect measures rather than graph-based estimators. The causal risk difference is collapsible over arbitrary baseline covariates with weights \(\Pr(V=v)\), the causal risk ratio \(RR(-)\) is collapsible with weights \(\Pr(V=v\mid Y^{a=0}=1)\), and no general weights exist for the odds ratio [1610.00033]. A geometric reformulation using Rothman diagrams states that a measure of association is collapsible if and only if all its contour lines are straight; in that framework, risk difference and risk ratio are collapsible, while odds ratio and cumulative hazard ratio are noncollapsible except at the null line [2506.23927]. These are properties of effect measures, not of reduced graphical estimators.

For categorical contingency tables, marginal log-linear models define collapsibility by equality of parameters across margins,
\[
\lambda_L^M(x_L)=\lambda_L^N(x_L),
\]
with strict collapsibility additionally requiring that all interactions involving removed variables vanish [1711.00680]. A further generalization is **average collapsibility**, under which the conditional average of an association measure equals the marginal association measure; this is weaker than ordinary collapsibility and is developed for the expectation dependence function, the mixed derivative of interaction, and the log-expectation dependence measure [1110.1740].

These neighboring frameworks clarify the specificity of estimate collapsibility. It is not primarily a statement about effect-measure algebra, contingency-table parameter equality, or separator structure alone. It is a statement that marginalization over graph vertices leaves a designated likelihood-based causal estimator unchanged.

## 7. Empirical behavior and significance

The empirical evaluation attached to the strong d-convex hull framework emphasizes exactness together with substantial reduction in computational burden [2606.08941].

| Domain | Setting | Reported outcome |
|---|---|---|
| Probabilistic reasoning | Hailfinder, Win95pts, Pathfinder, Munin1, Andes | 56 nodes reduced to about 14; 76 to 8.77; 109 to 4.68; 186 to 22.35; 223 to 68.83 |
| Probabilistic reasoning | Same benchmarks | KL divergence on the \(10^{-18}\) scale; node reduction up to **96%** |
| Synthetic causal estimation | Random DAGs with \(20\) to \(100\) nodes and densities \(0.05\) to \(0.15\) | Recall near 1; Precision is 1; speedup up to about **18×** |
| Real-world causal estimation | Sachs, Insurance, Alarm, Hepar2, Pathfinder, Munin1 | Precision stays at **1.0**; Recall typically \(>0.86\); speedup **1.13×** to **30.37×** |

The paper interprets the slight recall drop in real-world causal experiments as finite-sample noise rather than a theoretical failure [2606.08941]. That interpretation is consistent with the exact structural theorems: the reduced estimator is guaranteed to match the full estimator when the strong d-convex hull criterion is satisfied.

The broader significance of estimate collapsibility is that it provides a precise notion of causally sufficient locality in CPDAGs. Instead of treating graph reduction as an approximation, it identifies the smallest induced subgraph that preserves maximum-likelihood estimation and back-door-based causal estimation exactly. This makes the concept relevant both for theoretical characterizations of estimator invariance and for practical acceleration of equivalence-class causal effect estimation.

Source: https://www.emergentmind.com/topics/estimate-collapsibility