Stable Differentiable Causal Discovery
About
Inferring causal relationships as directed acyclic graphs (DAGs) is an important but challenging problem. Differentiable Causal Discovery (DCD) is a promising approach to this problem, framing the search as a continuous optimization. But existing DCD methods are numerically unstable, with poor performance beyond tens of variables. In this paper, we propose Stable Differentiable Causal Discovery (SDCD), a new method that improves previous DCD methods in two ways: (1) It employs an alternative constraint for acyclicity; this constraint is more stable, both theoretically and empirically, and fast to compute. (2) It uses a training procedure tailored for sparse causal graphs, which are common in real-world scenarios. We first derive SDCD and prove its stability and correctness. We then evaluate it with both observational and interventional data and on both small-scale and large-scale settings. We find that SDCD outperforms existing methods in both convergence speed and accuracy and can scale to thousands of variables. We provide code at https://github.com/azizilab/sdcd.
Related benchmarks
| Task | Dataset | Result | Rank | |
|---|---|---|---|---|
| Causal Discovery | Synthetic Data | Runtime1.65e+3 | 57 | |
| Causal Discovery | Synthetic (n=100, |E|=400, sample size=1000) | mAP65.7 | 36 | |
| Causal Discovery | Synthetic n=1000, |E|=2000, sample size=1000 | mAP59.6 | 32 | |
| Causal Discovery | Sachs real data d=11 | SHD13 | 19 | |
| Causal Discovery | Semantic Causal Environment observation-only | F1 Score47.9 | 15 | |
| Causal Discovery | Synthetic Data Observation-only (1000 samples) | Rank10.1 | 15 | |
| Causal Discovery | K562 Perturb-seq CausalBench | W Score18.6 | 13 | |
| Causal Discovery | SACHS p = 11, s = 20, n = 100 (real flow cytometry) | F1 Score41 | 13 | |
| Causal Discovery | RPE1 Perturb-seq CausalBench | W Score12.6 | 13 | |
| Causal Discovery | Synthetic Data Mixed-interventional (800 observational + 200 interventional samples) | Rank4.9 | 9 |