The recent works on causal discovery have followed a similar trend of learning partial ancestral graphs (PAGs) since observational data constrain the true causal directed acyclic graph (DAG) only up to a Markov equivalence class. This limits their application in the majority of downstream tasks, as uncertainty in causal relations remains unresolved. We propose a new refinement framework, CausalSAGE, for converting PAGs to DAGs while respecting the underlying causal relations. The framework expands discrete variables into state-level representations, constrains the search space using structural knowledge and soft priors, and applies a unified differentiable objective for joint optimization. The final DAG is obtained by aggregating the optimized structures and enforcing acyclicity when necessary. Our experimental evaluations show that the obtained DAGs preserve the underlying causal relations while also being efficient to obtain.
Figures & tables
Figure 1: Why PAGs are not enough. (a) Most downstream tasks require a fully specified causal DAG. (b) FCI handles latent confounders but outputs partially oriented graphs (PAGs). (c) Ambiguous parent sets (e.g., Pa(X) ) prevent computing interventional effects.
Figure 2: CausalSAGE . The framework expands discrete variables into state representations, constrains the search space using random/LLM priors, and finally applies a unified differentiable objective for joint optimization to output a DAG.
Figure 3: SHD reduction via refinement (Q1). Bars show Δ SHD = (Baseline − Refined), where the positive values indicate improvement.
Figure 4: Elimination of Orientation Ambiguity. FCI/RFCI leaves a large fraction of edges unresolved (46%–86%) across datasets. CausalSAGE reduces this ratio to 0% in all cases, effectively breaking Markov equivalence and producing fully directed DAGs.
Figure 5: Comparison with standard DAG learners (10k samples). SHD measures structural error, while F1 measures directional accuracy. Points closer to the top-left corner indicate stronger overall performance.
Figure 6: Effect of sample size on SHD. As the number of observational samples increases, structural error consistently decreases across datasets.
Dataset
Variables
Sample Size
Runtime (s)
Sachs
11
855
2.8
Child
20
5000
2.6
Insurance
27
10000
38
Alarm
37
10000
32
Hailfinder
56
10000
59
Win95pts
76
10000
71
Table 1: CausalSAGE is efficient. The runtime as a function of graph size shows the efficiency of our method.
Appendix figures & tables10 assets
Supplementary material from the paper’s appendix.
Appendix
Dataset
Nodes
Edges
Samples
Type
Sachs
11
17
855
Continuous
Child
20
25
5000
Discrete
Insurance
27
52
10000
Discrete
Alarm
37
46
10000
Discrete
Hailfinder
56
66
10000
Discrete
Win95pts
76
112
10000
Discrete
Appendix
Table 2: Benchmark datasets from the bnlearn repository.
Causal discovery has shifted from data-centric methods to hybrid strategies that integrate semantic knowledge from experts or large language models (LLMs). Such external information is vital for identifying causal structures beyond the Markov Equivalence Class (MEC), which data alone cannot resolve. However, expert availability is often limited, and LLMs frequently misidentify causal directions in specialized domains. To overcome such shortcomings, we propose a tag-based approach that leverages semantically meaningful labels while deriving causal directionality directly from data. Using variable-level tag assignments from available sources (e.g., LLMs), our tags for DAGs method learns from identifiable data structures to extract higher-level causal relations. These are then used to orient undirected edges, enabling causal discovery to move beyond the MEC without reliance on fallible external knowledge.
Florian Peter Busch, Moritz Willig, Florian Guldan +2
Computer Science Department, Technical University of Darmstadt · Hessian Center for AI (hessian.AI) · German Research Center for AI (DFKI) +2
Bayesian causal discovery seeks to determine the posterior distribution of causal theories, which are interpreted as directed acyclic graphs (DAGs) that explain the observed data. The resulting posterior allows systematic reasoning regarding epistemic uncertainty within these theories. Nonetheless, finding such graphs is difficult due to identifiability problems and limited observational data. Furthermore, precisely approximating posterior over graphs is challenging given vast range of potential DAGs. Recent Bayesian approaches have addressed some of these challenges, yet they remain limited as they fail to encode dependencies between edges, and lack principled ways to incorporate domain knowledge as inductive biases during the search process. To overcome these limitations, we propose SVI-DAG, a structured variational inference approach to Bayesian causal discovery using observational data and prior beliefs that uses normalizing flows to model dependencies between edges, supporting expressive and multimodal posterior learning over DAGs. To mitigate mode seeking behaviour in evidence lower bound optimization and promote mode coverage, we use stein variational gradient descent to update the node potentials using a kernel in acyclicity space. We evaluate SVI-DAG against 5 state-of-the-art Bayesian DAG learning methods and demonstrate superior performance in uncertainty quantification while remaining competitive in terms of structural accuracy.
Shrenik Zinage
Massachusetts Institute of Technology · This work was conducted while the author was a PhD student at Purdue University, independently of the author’s dissertation research.
Causal discovery from observational data is a fundamental yet challenging task in scientific research. While existing approaches are primarily based on conditional independence tests, structure scores, or restrictive functional assumptions, we propose Decoupled Causal Discovery (DCD), a novel decoupling-based perspective that does not rely on these methodologies. DCD directly identifies the Markov boundary (MB) by decoupling non-target variables via weighting functions, such that only variables within the MB preserve dependence with the target under the decoupled distribution. Building on this, DCD iteratively constructs the Completed Partially Directed Acyclic Graph (CPDAG) by exploiting structural asymmetries within the MBs. We establish the theoretical identifiability, soundness, and completeness of DCD. Empirical evaluations demonstrate that DCD achieves strong performance, particularly excelling in challenging noise regimes.
Zhengkang Guan, Fei Wu, Kun Kuang
College of Computer Science and Technology Zhejiang University