NeurIPS 2026/Temporal causal discovery

SCDM: Scalable Causal Discovery in
Nonlinear Temporal Systems
with Meta-Learning

A shared-parameter framework for nonlinear temporal causal discovery.

University of Tennessee at Chattanooga

arXiv: 2610.29992 Project: Page Code: GitHub Data: 3 benchmarks License: MIT

Explore the core implementation and runnable examples on GitHub.

The idea at a glanceConceptual illustration
01Many temporal variables

Each target variable defines a prediction task.

02One shared initialization
Θ = {θ, W}predictor + structure

Sample targets, adapt on support, update on query.

03Directed graph scores

Read out continuous scores after training.

Share the predictor. Keep the target-specific structural signal.Explore the method ↓
01 / Motivation

A quadratic problem.
A shared learning strategy.

SCDM treats the prediction of each target variable as a related task, sharing a temporal predictor and a causal-strength matrix across the system.

When the number of variables grows, so does the space of possible directed relations. Target-wise fitting and repeated conditional testing can become expensive. SCDM processes a subset of target tasks per episode, allowing information from those tasks to update a shared initialization.

Task sampling reduces per-episode target processing. The dense structural matrix and complete graph readout still depend on system dimension.

Explore the candidate edge spaced(d − 1)
999,000possible directed pairs

Counts exclude self-edges. This calculation illustrates graph size, not measured runtime.

Read the abstract

02 / Method

Adapt locally.
Learn collectively.

Temporary task adaptation and persistent shared learning play different roles. Follow one SCDM episode below.

Original paper overview
Each sampled target starts from the shared parameters. Support adaptation creates a temporary copy; query losses update the persistent initialization. Graph scoring follows training.
03 / Controlled recovery

Recover directed structure.
Across controlled systems.

Known directed graphs let us evaluate edge-ranking performance. Results vary with the dataset and the metric.

Evaluation metric

Mean ± 95% CI · higher is better

All benchmark results

Bold: highest reported mean
Controlled graph recovery, mean plus or minus 95 percent confidence interval

Values transcribed from the manuscript’s controlled-recovery table. AUROC and AUPRC are threshold-independent and unitless; self-edges are excluded. Rankings use reported means and do not imply statistical significance.

04 / Visual evidence

From brain networks
to traffic systems.

Controlled simulations and real-world observations answer different questions. Explore the figures and their evaluation roles.

100ROI-level simulated fMRI
1,000Voxel-level observed variables
3,000Sensors in LargeST-GLA-3000

These are benchmark dimensions. The voxel simulator includes latent regional drivers; LargeST has no complete directed causal ground truth.

05 / Theory & scope

A recovery principle,
with explicit conditions.

For a fully observed delayed system, the graph can be recovered when estimation errors remain small relative to the population risk margin.

Sufficient condition
εstat + εmeta + εreadout < Δmin / 2
εstat
Finite-sample error under temporal dependence
εmeta
Residual from finite episodic meta-training
εreadout
Difference between implemented and ideal readout

Under the stated assumptions, a threshold inside the separating margin recovers the delayed graph with high probability.

Interpretation matters

The guarantee assumes causal sufficiency, a correct history order, temporal regularity, a positive population margin, and controlled estimation errors. It does not cover arbitrary hidden confounding. Predictive readout and geographic locality are diagnostics; neither alone verifies an individual causal effect.

06 / Read & build on this work

The paper and
what comes next.

DATA & REPRODUCIBILITY

Download the original simulations.

Complete benchmark arrays with known directed graphs, preserved sequence partitions, and SHA-256 checksums.

These are simulated datasets, not recordings from human participants. fMRI and EEG retain the original train / validation / test sequence partitions. The voxel simulator includes latent regional drivers. The 64-variable quick start is a separate lightweight example. A single-seed run can differ from the paper's multi-seed results.

Citation

Reference this work.

Citation uses the paper title and author list.

Download image ↓