Causal-Audit: A Framework for Risk Assessment of Assumption Violations in Time-Series Causal Discovery
Authors: Marco Ruiz, Miguel Arana-Catania, David R. Ardila, Rodrigo Ventura
Organizations: Institute for Systems and Robotics, Instituto Superior Técnico, 1049-001 Lisbon, Portugal · Digital Scholarship at Oxford, University of Oxford, Oxford, OX1 3BG, United Kingdom · Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA 91109, USA
Time-series causal discovery methods rely on assumptions such as stationarity, regular sampling, and bounded temporal dependence. When these assumptions are violated, structure learning can produce confident but misleading causal graphs without warning. We introduce Causal-Audit, a framework that formalizes assumption validation as calibrated risk assessment. The framework computes effect-size diagnostics across five assumption families (stationarity, irregularity, persistence, nonlinearity, and confounding proxies), aggregates them into four calibrated risk scores with uncertainty intervals, and applies an abstention-aware decision policy that recommends methods (e.g., PCMCI+, VAR-based Granger causality) only when evidence supports reliable inference. The semi-automatic diagnostic stage can also be used independently for structured assumption auditing in individual studies. Evaluation on a synthetic atlas of 500 data-generating processes (DGPs) spanning 10 violation families demonstrates well-calibrated risk scores (AUROC > 0.95), a 62% false positive reduction among recommended datasets, and 78% abstention on severe-violation cases. On 21 external evaluations from TimeGraph (18 categories) and CausalTime (3 domains), recommend-or-abstain decisions are consistent with benchmark specifications in all cases. An open-source implementation of our framework is available.
Figures & tables
Family
Assumption
Consequence if Violated
Error
@BTrule []Stationarity
No trends or regime changes
Spurious links from trends
FP
Irregularity
No gaps or missing values
Selection bias, incorrect lags
Bias
Persistence
Low autocorrelation
Inflated significance
FP
Nonlinearity
Linear relationships
Missed nonlinear effects
FN
Confounding proxies
Stable model coefficients
Spurious links from instability
FP
Table 1: Assumption families and consequences of violations. Causal-Audit computes diagnostics for all five families and aggregates four (all except nonlinearity) into calibrated risk scores.
Figure 4: Detailed flowcharts for each pipeline stage. (a) Stage I: Diagnostic Auditing computes five diagnostic families from input X , producing the diagnostic vector d=[d1,…,d5] . (b) Stage II: Risk Estimation transforms diagnostics into calibrated risk scores via logistic aggregation, isotonic calibration, and bootstrap uncertainty quantification. (c) Stage III: Decision Policy evaluates thresholds to output Recommend m∗ or Abstain .
VIF Range
Linear (VAR/Granger)
PCMCI+ (ParCorr)
Transfer Entropy
@BTrule [] <10
Reliable
Reliable
Reliable
10 to 100
Caution
With regularization
Reliable
100 to 104
Unreliable
Pseudo-inverse needed
Caution
>104
Unreliable
Consider alternatives
Caution
Table 2: Expected method-specific VIF tolerance, derived from the statistical properties of each estimator class. The VIF <10 threshold follows the widely adopted rule of thumb [ 70 ] , though recent work suggests context-dependent interpretation [ 71 ] . Higher VIF ranges reflect the known sensitivity of OLS-based and partial-correlation-based estimators to collinear inputs; information-theoretic estimators avoid numerical instability from collinearity but suffer from degraded conditional distribution estimates in finite samples.
Risk Dimension
Model
Key Features
@BTrule [] Rnonstat
σ(−2.0+1.0⋅xbreak_mag+2.0⋅xdrift
Break magnitude and drift
+0.5⋅xadf+0.5⋅xkpss)
slope dominate
@BTrule [] Rirreg
σ(−1.5+2.5⋅xgap_cv+2.0⋅xmissing
Gap regularity is primary
+1.5⋅xseasonal_miss)
contributor
@BTrule [] Rpersist
σ(−1.5+(−3.0)⋅xTeff/T
Negative weight: lower
+2.0⋅xτint)
Teff/T→ higher risk
Table 3: Risk-specific logistic regression models. Weights are initialized from domain knowledge and refined via hierarchical Bayesian calibration.
Parameter
Value
@BTrule []Calibration / validation split
396(79%) / 100 (21%)
Stratification
By violation family
Isotonic binning
10 quantile bins
Stability assessment
5-fold and 10-fold stratified CV
Table 4: Calibration configuration.
Method
Rcomp
Rnonstat
Rirreg
Rconfound
Rpersist
@BTrule []Granger (VAR)
≤0.30
Hard ≤0.60
Hard ≤0.50
Hard ≤0.60
Soft ≤0.75
PCMCI+
≤0.70
Hard ≤0.80
Soft ≤0.60
Soft ≤0.80
Soft ≤0.85
LPCMCI †
≤0.80
Hard ≤0.80
Soft ≤0.70
Soft ≤0.90
Soft ≤0.85
Transfer Entropy †
≤0.90
Soft ≤0.85
Soft ≤0.80
Soft ≤0.95
Soft ≤0.90
Table 5: Method-assumption constraint mapping. Per-dimension thresholds reflect the relative robustness of each method to each violation type, informed by documented method properties (see text); composite thresholds are validated to ensure expected failure rates below 20% on the calibration corpus.
Table 6: Synthetic DGP Atlas: 10 violation families with severity gradations. Each family contains 50 datasets whose violation-controlling parameters span the severity levels listed in the rightmost column. Families F1–F9 use N∈{5,6,7,8} variables and T∈{500,750,1000} observations; F10 uses wider ranges ( N∈{3–12} , T∈{200–2000} ) to stress-test boundary conditions. Families F2–F5 correspond directly to the four calibrated risk dimensions. The remaining families target additional assumptions or combine multiple violations.
ISR-Lisbon, Instituto Superior Técnico, Av. Rovisco Pais 1, Lisbon, 1049-001, Portugal · Digital Scholarship at Oxford, University of Oxford, Oxford, OX1 3BG, United Kingdom · Jet Propulsion Lab., Caltech, Pasadena, CA 91109, California, USA