2025-11-24T14:22:17.661777

Robust Causal Discovery in Real-World Time Series with Power-Laws

Tusoni, Masi, Coletta et al.

Exploring causal relationships in stochastic time series is a challenging yet crucial task with a vast range of applications, including finance, economics, neuroscience, and climate science. Many algorithms for Causal Discovery (CD) have been proposed, but they often exhibit a high sensitivity to noise, resulting in misleading causal inferences when applied to real data. In this paper, we observe that the frequency spectra of typical real-world time series follow a power-law distribution, notably due to an inherent self-organizing behavior. Leveraging this insight, we build a robust CD method based on the extraction of power -law spectral features that amplify genuine causal signals. Our method consistently outperforms state-of-the-art alternatives on both synthetic benchmarks and real-world datasets with known causal structures, demonstrating its robustness and practical relevance.

academic

Robuste Kausaldiskriminierung in realen Zeitreihen mit Potenzgesetzen

Grundlegende Informationen

Paper-ID: 2507.12257
Titel: Robust Causal Discovery in Real-World Time Series with Power-Laws
Autoren: Matteo Tusoni, Giuseppe Masi, Andrea Coletta, Aldo Glielmo, Viviana Arrigoni, Novella Bartolini
Klassifizierung: cs.LG physics.data-an stat.ML stat.OT
Veröffentlichungsdatum: 12. Oktober 2025 (arXiv v2)
Paper-Link: https://arxiv.org/abs/2507.12257

Zusammenfassung

Die Erforschung von Kausalbeziehungen in stochastischen Zeitreihen ist eine herausfordernde, aber entscheidende Aufgabe mit breiter Anwendung in Finanzen, Wirtschaft, Neurowissenschaften und Klimawissenschaften. Obwohl zahlreiche Algorithmen zur Kausaldiskriminierung (CD) vorgeschlagen wurden, sind diese häufig hochgradig lärmempfindlich und führen bei der Anwendung auf reale Daten zu irreführenden Kausalschlüssen. Diese Arbeit beobachtet, dass das Frequenzspektrum typischer realer Zeitreihen einer Potenzgesetzverteilung folgt, was hauptsächlich auf das inhärente selbstorganisierende Verhalten zurückzuführen ist. Basierend auf dieser Erkenntnis konstruieren wir eine robuste Kausaldiskriminierungsmethode, die auf der Merkmalsextraktion von Potenzgesetz-Spektren basiert und echte Kausalsignale verstärken kann. Unsere Methode übertrifft konsistent modernste Alternativen bei synthetischen Benchmarks und realen Datensätzen mit bekannter Kausalstruktur und demonstriert ihre Robustheit und praktische Relevanz.

Forschungshintergrund und Motivation

Problemdefinition

Diese Forschung befasst sich mit dem Problem der Kausaldiskriminierung in Zeitreihendaten, d. h. der Identifikation von Kausalbeziehungen zwischen Variablen aus Beobachtungsdaten. Traditionelle Kausaldiskriminierungsmethoden, insbesondere solche auf Basis der Granger-Kausalität, zeigen bei komplexen realen Daten folgende Einschränkungen:

Lärmempfindlichkeit: Traditionelle Methoden sind hochgradig empfindlich gegenüber nicht-gaußschem Lärm, Nichtstationarität und nichtlinearen Störungen
Annahmebeschränkungen: Abhängigkeit von strengen Annahmen wie Lärmstationarität und einzelner Merkmalsskalen
Erkennung falscher Beziehungen: Neigung, Lärmkorrelationen als Kausalbeziehungen zu interpretieren

Forschungsmotivation

Die Autoren beobachten, dass reale Systeme allgemein Potenzgesetz-Spektraleigenschaften aufweisen, die aus folgenden Gründen entstehen:

Selbstorganisierendes Verhalten mehrerer wechselwirkender Einheiten
Skalierungsinvarianz aufgrund fehlender externer Koordination
Fraktale Eigenschaften und Langzeitkorrelationen des Systems

Basierend auf dieser Beobachtung schlägt die Arbeit vor, Potenzgesetz-Spektralmerkmale für robustere Kausaldiskriminierung zu nutzen.

Kernbeiträge

Vorschlag des PLaCy-Rahmens: Eine neuartige Kausaldiskriminierungsmethode basierend auf Potenzgesetz-Spektralmerkmalen
Theoretische Garantien: Beweis der Invarianz der Kausaldiagrammstruktur unter Frequenzbereichstransformation (Theorem 1)
Experimentelle Validierung: Umfassende Bewertung auf synthetischen und realen Datensätzen mit überlegener Robustheit
Methodische Universalität: Demonstration der Verbesserungseffekte der Spektralvorverarbeitung auf andere Kausaldiskriminierungsalgorithmen

Methodische Details

Aufgabendefinition

Gegeben eine multivariate Zeitreihe $x \in \mathbb{R}^{L \times d}$ , besteht das Ziel darin, einen gerichteten Graphen $G = (V, E)$ abzuleiten, wobei:

$V = \{1, 2, ..., d\}$ Systemvariablen darstellt
$E \subseteq V \times V$ die Menge der Kausalkanten darstellt
Eine gerichtete Kante $(i,j)$ existiert genau dann, wenn $x_i$ eine Kausalursache von $x_j$ ist

Modellarchitektur

1. Schiebefenster-Segmentierung

Segmentierung jeder Zeitreihe in überlappende Fenster der Länge $l$ mit Schrittweite $s$ : $w_i^k = (x_i(k \cdot s), ..., x_i(k \cdot s + l - 1))$

2. Spektrale Merkmalsextraktion

Anwendung der diskreten Fourier-Transformation auf jedes Fenster: $\phi(k) = \sum_{t=0}^{L-1} x(t) e^{-i2\pi k t/L}$

Berechnung der Spektralamplitude: $A(f_k) = |\phi(k)|$

3. Potenzgesetz-Anpassung

Anpassung eines linearen Modells im Log-Log-Raum: $\log A(f) = a - \lambda \log f$

wobei $a$ der Abschnittparameter und $\lambda > 0$ der Spektralexponent ist.

4. Kausalanalyse

Anwendung des multivariaten Granger-Kausalitätstests auf die extrahierten Spektralparameter-Zeitreihen $(a_i, \lambda_i)$ , um die Vorhersagekraft von $(\lambda_i, a_i)$ auf $\lambda_j$ zu bewerten.

Algorithmus-Ablauf (PLaCy)

Eingabe: Zeitreihe x = (x₁, ..., xₐ), Fenstergröße l, Schrittweite s
Ausgabe: Kausaldiagramm G

1. Segmentiere jedes xᵢ in ⌊(L-l)/s⌋+1 Schiebefenster wᵢᵏ
2. for each i ∈ {1, ..., d} do
3.   for each k ∈ {0, ..., ⌊(L-l)/s⌋} do
4.     Wende DFT auf wᵢᵏ an, um φᵢᵏ zu erhalten
5.     Erhalte (aᵢᵏ, λᵢᵏ) durch Anpassung von Gleichung (2)
6.   Verkette (aᵢᵏ, λᵢᵏ) zu Zeitreihe (aᵢ, λᵢ)
7. for each i,j ∈ {1, ..., d}, i ≠ j do
8.   Gᵢ,ⱼ ← Granger-Kausalitätstest mit (aᵢ,λᵢ) als Ursache, λⱼ als Ergebnis
9. return G

Technische Innovationen

Frequenzbereich-Kausaldiskriminierung: Erste systematische Nutzung von Potenzgesetz-Spektralmerkmalen für Kausalschluss
Adaptive Fensterauswahl: Automatische Auswahl der optimalen Fensterlänge durch p-Wert-Kriterium
Lärmrobustheit: Spektrale Anpassung als natürlicher Entstörungsschritt, verbesserte Robustheit gegenüber nicht-gaußschen Schwankungen
Theoretische Grundlage: Theoretischer Beweis der Invarianz des Kausaldiagramms unter Spektraltransformation

Experimentelle Einrichtung

Datensätze

Synthetische Datensätze

Basierend auf verallgemeinertem Ornstein-Uhlenbeck-Prozess mit vier Szenarien: $x(t+\Delta t) = x(t) + \frac{\Delta t}{\tau_c}(\mu - x(t)) + (\sigma_b \epsilon_b(t) + \sigma_g^a \epsilon_g^a(t) + \sigma_g^m \epsilon_g^m(t) \cdot x(t))\sqrt{\Delta t}$

OU( $\sigma_g^m = 0$ ): Gleichgewichtszustand ohne multiplikativen Lärm
OU( $\sigma_g^m > 0$ ): Gleichgewichtszustand mit multiplikativem Lärm
ÔU( $\sigma_g^m = 0$ ): Nichtgleichgewichtszustand ohne multiplikativen Lärm
ÔU( $\sigma_g^m > 0$ ): Nichtgleichgewichtszustand mit multiplikativem Lärm

Reale Datensätze

Rivers-Datensatz: Wasserstands- und Niederschlagsdaten von drei hydrologischen Stationen in Süddeutschland
AirQuality-Datensatz: PM2.5-Verschmutzungsüberwachungsdaten mehrerer chinesischer Städte

Bewertungsmetriken

F1-Score: Misst die Gesamtleistung der Kausalbeziehungserkennung
True Negative Rate (TNR): Bewertet die Fähigkeit des Algorithmus, falsche Assoziationen auszuschließen

Vergleichsmethoden

Traditionelle Methoden: Granger Causality, PCMCI, PCMCIΩ
Optimierungsmethoden: DYNOTEARS, RCV-VarLiNGAM
Tiefes Lernen: Rhino
Nichtlineare Methoden: CCM-Filtering
Frequenzbereich-Methoden: BCGeweke, DTF, GewekeNP

Implementierungsdetails

Schiebefensterlänge: $l = 50$ (durch p-Wert-Kriterium ausgewählt)
Schrittweite: $s = 1$
Verzögerungsterme: 10
Statistische Signifikanzschwelle: $p = 0.05$

Experimentelle Ergebnisse

Hauptergebnisse

Leistung auf synthetischen Datensätzen (N=5, $\sigma_g^a = 1.0$ ):

Datensatz	PLaCy F1	Beste Baseline F1	PLaCy TNR	Beste Baseline TNR
OU( $\sigma_g^m = 0$ )	0.77±0.17	0.61±0.18	0.94±0.05	0.99±0.02
OU( $\sigma_g^m > 0$ )	0.80±0.17	0.79±0.11	0.94±0.06	0.98±0.03
ÔU( $\sigma_g^m = 0$ )	0.70±0.17	0.58±0.18	0.88±0.09	0.99±0.02
ÔU( $\sigma_g^m > 0$ )	0.80±0.17	0.71±0.13	0.93±0.07	0.98±0.03

Ergebnisse auf realen Datensätzen:

Datensatz	PLaCy F1	PLaCy TNR	Beste Baseline F1	Beste Baseline TNR
Rivers	0.51±0.10	0.75±0.13	0.47±0.07	0.74±0.05
AirQuality	0.45±0.04	0.66±0.07	0.44±0.01	0.95±0.02

Schlüsselfunde

Robustheit gegenüber multiplikativem Lärm: PLaCy zeigt besonders hervorragende Leistung in Szenarien mit multiplikativem Lärm
Anpassungsfähigkeit an Nichtgleichgewichtszustände: Behält gute Leistung unter Nichtgleichgewichtsbedingungen bei
Vorteile der Frequenzbereich-Methode: Frequenzbereich-Analyse zeigt bessere Lärmresistenz im Vergleich zu Zeitbereich-Methoden
Universelle Verbesserung: Die Anwendung der Spektralvorverarbeitung auf Methoden wie PCMCI kann die Leistung erheblich verbessern

Ablationsstudien

Analysen der Fensterlänge und Schrittweite zeigen:

Optimale Leistung bei Schrittweite 1, erfasst kurzfristige Kausalabhängigkeiten
Adaptive Auswahl der Fensterlänge durch p-Wert-Kriterium ist am effektivsten
Zu kurze oder zu lange Fenster reduzieren die Leistung

Schlussfolgerungen und Diskussion

Hauptschlussfolgerungen

PLaCy erreicht robustere Kausaldiskriminierung durch Nutzung von Potenzgesetz-Spektralmerkmalen
Die Methode zeigt überlegene Leistung auf synthetischen und realen Daten
Frequenzbereich-Analyse bietet neue Perspektiven für Zeitreihen-Kausaldiskriminierung

Einschränkungen

Systeme mit langsamen Spektralveränderungen: Begrenzte Effektivität bei Systemen mit sehr langsamen Spektralparameteränderungen
Kurze Zeitreihen: Erfordert ausreichend lange Sequenzen für stabile Spektralschätzung
Rechenkomplexität: Zusätzlicher Rechenaufwand durch Spektralanalyse im Vergleich zu einfachen Methoden

Zukünftige Richtungen

Erweiterung auf nicht-VAR-Kausaldiskriminierungsmethoden
Tiefere Untersuchung statistischer Parameter der Spektraldichte
Behandlung latenter Störfaktoren
Entwicklung effizienterer Online-Kausaldiskriminierungsalgorithmen

Tiefgehende Bewertung

Stärken

Hohe Innovativität: Erste systematische Anwendung von Potenzgesetz-Spektralmerkmalen auf Kausaldiskriminierung
Solide Theorie: Strenge theoretische Analyse und Beweise
Umfassende Experimente: Abdeckung mehrerer synthetischer Szenarien und realer Anwendungen
Hoher praktischer Wert: Signifikante Vorteile in lärmreichen Umgebungen

Schwächen

Begrenzte Anwendbarkeit: Hauptsächlich für Systeme mit Potenzgesetz-Spektraleigenschaften geeignet
Parameterauswahl: Auswahl von Parametern wie Fensterlänge erfordert Erfahrung
Recheneffizienz: Höherer Rechenaufwand im Vergleich zu einfachen Methoden

Einfluss

Akademischer Beitrag: Bietet neue Forschungsrichtung für Zeitreihen-Kausaldiskriminierung
Praktischer Wert: Breite Anwendungsperspektiven in Finanzen, Klima und anderen Bereichen mit Potenzgesetz-Eigenschaften
Reproduzierbarkeit: Vollständige Algorithmusbeschreibung und Open-Source-Code

Anwendungsszenarien

Finanzmarktdatenanalyse
Klimasystemmodellierung
Neurowissenschaftliche Forschung
Soziale Netzwerkanalyse
Beliebige komplexe Systeme mit selbstorganisierenden Eigenschaften

Literaturverzeichnis

Das Papier zitiert 51 relevante Referenzen, die wichtige Arbeiten in den Bereichen Kausaldiskriminierung, Zeitreihenanalyse und komplexe Systeme abdecken und eine solide theoretische Grundlage für die Forschung bieten.

Gesamtbewertung: Dies ist ein hochqualitatives Forschungspapier, das eine innovative Methode im Bereich der Zeitreihen-Kausaldiskriminierung vorschlägt. Durch geschickte Nutzung der Potenzgesetz-Spektraleigenschaften realer Systeme wird die Robustheit der Kausaldiskriminierung erfolgreich verbessert. Die theoretische Analyse ist streng, das Experimentdesign rational und die Ergebnisse überzeugend. Diese Arbeit bietet neue Werkzeuge und Perspektiven für Kausalschluss in komplexen Systemen.