← All publications
BETA This research page is an early version. Explanatory text, figure selection, and presentation are still being revised.
NeurIPS 2023 — workshop on Optimal Transport and Machine Learning (OTML)

Causal Discovery via Monotone Triangular Transport Maps

Sina Akbari · Luca Ganassali · Negar Kiyavash
NeurIPS 2023 Workshop on Optimal Transport and Machine Learning · 2023
Keywords.
causal discovery optimal transport conditional independence transport maps

TL;DR.

Monotone triangular transport maps provide a common engine for both constraint-based and score-based causal discovery. We use them to build noise-agnostic CI tests, recover graphs up to Markov equivalence even with latent variables, and—under additive-noise or post-nonlinear structure—score candidate graphs for unique recovery.

Why this matters.

Transport maps let causal discovery move beyond Gaussian or narrowly parametric CI tests while keeping one coherent representation of the full data distribution. The same fitted object supports sparsity discovery and stronger model-based orientation.

Abstract

We study the problem of causal structure learning from data using transport maps. Specifically, we first provide a constraint-based method which builds upon lower-triangular monotone parametric transport maps to design conditional independence tests which are agnostic to the noise distribution. We provide an algorithm for causal discovery up to Markov Equivalence for general structural equations and noise distributions, which allows for settings with latent variables. Our approach also extends to score-based causal discovery by providing a novel means for defining scores. This allows us to uniquely recover the causal graph under additional identifiability and structural assumptions, such as additive noise or post-nonlinear models. We provide experimental results to compare the proposed approach with the state of the art on both synthetic and real-world datasets.

Main message. Transport maps are not only distribution estimators: they provide a unified toolkit for conditional-independence testing, graph recovery, and structural scoring.

The paper at a glance

An informal guide to the problem, the idea, and the main results.

optimal transport transport maps conditional independence causal discovery

1. The problem

Many causal-discovery tools are tightly tied to simple distributional assumptions, such as linear-Gaussian structure. We want discovery machinery that remains flexible while still being statistically and computationally useful.

2. The key idea

Monotone triangular transport maps give a rich representation of distributions. Once we can fit them well, we can use the representation to test conditional independence or score candidate graphs in a causally meaningful way.

3. Noise-agnostic constraint-based discovery

Method + theory

We use lower-triangular monotone transport maps to construct conditional-independence tests for general continuous distributions and integrate them into a PC-style causal-discovery procedure. The method recovers the graph up to Markov equivalence and can be used in the presence of latent variables.

4. Unique recovery under structural models

Theory

For additive-noise and post-nonlinear structural equation models, we derive transport-map characterizations that are necessary and sufficient for membership in these model classes. These characterizations yield scores that can distinguish candidate DAGs within a Markov equivalence class.

5. How it works

  1. Fit monotone triangular transport maps to the observed distribution.
  2. Use the fitted maps to test conditional independences or evaluate graph scores.
  3. Run the corresponding discovery routine.
  4. Under stronger assumptions, leverage the same machinery for sharper identification of the graph.

What the numerical experiments show

Main empirical message.

In non-Gaussian synthetic experiments, PC-OT outperforms Gaussian-PC and correlation-based Grow-Shrink once the sample size is sufficiently large, reducing both misoriented edges and overall graph loss. In the additive-noise experiment, the transport-based score clearly separates the true causal ordering from the other DAGs in the same Markov equivalence class.

Where this helps

Good fit

  • You want causal discovery with flexible distributional modeling.
  • You are interested in optimal-transport ideas beyond pure density fitting.
  • You want one framework that can support both testing and score-based discovery.

Keep in mind

  • The same transport-map fit supports both constraint-based and score-based discovery, so sparsity learning and stronger causal orientation live in one framework.
  • The constraint-based component is designed to be agnostic to the noise distribution, while stronger ANM/PNL structure can be used when available to obtain unique recovery.

Cite this paper

@inproceedings{akbari2023transport,
  title     = {Causal Discovery via Monotone Triangular Transport Maps},
  author    = {Sina Akbari and Luca Ganassali and Negar Kiyavash},
  booktitle = {NeurIPS 2023 Workshop on Optimal Transport and Machine Learning},
  year      = {2023},
  url       = {https://openreview.net/forum?id=PMGGrKTIii}
}