Causal discovery aims to understand the relationships between individual random variables. In many applications, such as brain imaging and climate modeling, it is more meaningful to consider interactions among groups of variables. Existing methods assume that knowledge of such groups or clusters is explicitly available when modeling interactions. However, in practice, these clusters as well as the causal relationships among them, are latent. In this paper, we present a novel approach based on variational inference to simultaneously infer both the latent clusters and causal structures. We learn an approximate posterior over clusters and graph-structure by considering variational distributions based on categorical and Bernoulli models respectively. We derive variational lower bounds and estimation techniques to learn variational and model parameters. The effectiveness of our proposed methods for cluster and causal discovery are demonstrated on both synthetic and real data sets.
Figures & tables
Figure 1: Example of a cluster-DAG
Figure 2: Graphical model of the generative process
Figure 3: Example of a) C-DAG b) Expanded Graph
Table 1: Results on linear synthetic datasets. ExG-SHD, E-SHD, and HD: lower is better; E-RI and AUC: higher is better. NA: Not Applicable.
Table 2: Results on nonlinear synthetic datasets. ExG-SHD and E-SHD: lower is better; E-RI and AUC: higher is better.
Data- sets
Methods
ExG- SHD
E- RI
AUC
E- SHD
HD
Protein Dataset
Baseline
3.0
0.55
1.0
0.3
NA
Lin Fact-C-E
27.4
0.57
1.0
1.0
1.0
Lin AR-C-Fact-E
0.0
0.56
1.0
0.0
0.7
Neur Fact-C-E
0.04
0.74
1.0
0.0
NA
Neur AR-C-Fact-E
0.0
0.56
1.0
0.0
NA
Neur Fact-C-Cond-E
0.0
0.80
1.0
0.0
NA
Table 3: Results on real-world datasets. NA means Not Applicable. ExG-SHD, E-SHD, HD lower is better. E-RI and AUC higher is better.
Appendix figures & tables8 assets
Supplementary material from the paper’s appendix.
Appendix
Figure 4: Data Generating CDAG for 4 variable 3 cluster synthetic dataset a) fork b) chain c) v-structure
Figure 5: Data Generating CDAG for 7 variable 3 cluster synthetic dataset a) fork b) chain c) v-structure
Dataset
C 1
C 2
4 variable 3 cluster
fork
Tanh
Tanh
chain
Sin
ReLU
v-structure
-
Cubic Polynomial
7 variable 3 cluster
fork
Sin
ReLU
chain
ReLU
Sin
v-structure
-
Tanh
Appendix
Table 4: Element-wise non-linearities applied after linear combination of parents
Figure 6: Sachs Dataset
Figure 7: Ground truth C-DAG for climate dataset
Figure 8: Autoregressive generation of cluster assignments using the linear model
Dataset
Hidden-State Dimension
4 variable 3 cluster
fork
32
chain
4
v-structure
4
7 variable 3 cluster
fork
8
chain
128
v-structure
48
Appendix
Table 5: Linear Autoregressive Model hidden-state dimension choices
Dataset
E-RI±std
E-SHD±std
4 variable 3 cluster
fork
0.7±0.1
1.08±0.51
chain
0.75±0.12
0.94±0.64
v-structure
0.72±0.08
0.95±0.21
7 variable 3 cluster
fork
0.68±0.07
1.05±0.22
chain
0.60±0.06
0.97±0.5
v-structure
0.65±0.06
1.0±0.0
Appendix
Table 6: Variation across random restarts for autoregressive linear cluster model for selected hidden dimension size
Latent variables pose a fundamental challenge to causal discovery and inference. Conventional local methods focus on direct neighbors but fail to provide macro level insights. Cluster level methods enable macro causal reasoning but either assume clusters are known a priori or require causal sufficiency. Moreover, directly applying single variable causal discovery methods to cluster level problems violates causal sufficiency and leads to incorrect results. To overcome these limitations, this paper proposes L2C (Local to Cluster Causal Abstraction), a unified framework that bridges local structure learning and cluster level causal discovery. Unlike prior work that requires a complete manual assignment of micro variables to clusters, L2C discovers the partition automatically from local causal patterns. Our solution leverages a cluster reduction theorem to reduce any cluster to at most three nodes without loss of causal information, applies local causal discovery to identify direct causes, effects, and V structures in the presence of latent variables, and performs macro level causal inference via cluster level calculus on the learned cluster graph. L2C does not assume causal sufficiency, as latent variables are handled through local discovery. Theoretical analysis shows that L2C ensures soundness, atomic completeness, and computational efficiency. Extensive experiments on synthetic and real world data demonstrate that L2C accurately recovers ground truth clusters and achieves superior macro causal effect identification compared to existing baselines.
Recent algorithmic advances have made directed acyclic graph (DAG) structure learning scalable for causal discovery. Yet, the currently available techniques assume a completely homogeneous population, precluding their application to clustered data where cluster-specific variations (e.g., patient-specific effects) are common. We address this issue by introducing a new approach that estimates a global structure while accounting for local cluster-level effects. The key idea is to extend the fixed- and random-effects framework of classical mixed models to the structure learning setting. Towards this end, we present a differentiable graph coupling mechanism that guarantees the union of the fixed- and random-effects graphs remains acyclic. Computationally, we provide a provably convergent first-order method and leverage efficient batched updates across clusters. Statistically, we establish identifiability of the model and show that our approach recovers the true structure asymptotically. In experiments on real and synthetic data, our proposal detects dependencies missed by alternative estimators, underscoring its value for structure learning in clustered settings.
Ryan Thompson, Matt P. Wand, Veerabhadran Baladandayuthapani
School of Mathematical and Physical Sciences University of Technology Sydney Ultimo, NSW 2007, Australia · Department of Biostatistics University of Michigan Ann Arbor, MI 48105, USA
Causal discovery aims to uncover the underlying causal relationships given data generated from a system. The goal, however, is not merely to predict causal edges given data, but also to be able to interpret and explain either observed or hypothesized phenomena, such as a particularly large causal effect. We consider this task of conditional causal discovery and cast it as a Bayesian inference problem, in which we target the posterior over causal graphs and parameters conditional on an event such as a causal-effect constraint. Unfortunately, this poses a computational challenge: existing approaches to Bayesian causal discovery struggle when the event has small posterior mass. To address this, we adapt rare-event estimation techniques to perform inference the joint graph-parameter space. Our method gradually drives a particle population toward the constrained region while maintaining samples that approximate the conditional posterior. Empirical evaluation on synthetic graphs validates the accuracy of our approach at small and large scales, and we show in a case study on the Sachs protein dataset how our method can be used to aid scientific exploration by providing pathway-level summaries.
Cixuan Zhang, Guy Van den Broeck, Benjie Wang
Computer Science Dept., Yale University, New Haven, Connecticut, USA · Computer Science Dept., University of California, Los Angeles, Los Angeles, California, USA