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

বাস্তব-বিশ্বের সময় শ্রেণী তে শক্তি-আইন সহ শক্তিশালী কার্যকারণ আবিষ্কার

মৌলিক তথ্য

পেপার আইডি: 2507.12257
শিরোনাম: বাস্তব-বিশ্বের সময় শ্রেণীতে শক্তি-আইন সহ শক্তিশালী কার্যকারণ আবিষ্কার
লেখক: ম্যাটিও তুসোনি, গিউসেপ মাসি, আন্দ্রিয়া কোলেটা, আল্ডো গ্লিয়েলমো, ভিভিয়ানা অ্যারিগোনি, নোভেল্লা বার্টোলিনি
শ্রেণীবিভাগ: cs.LG physics.data-an stat.ML stat.OT
প্রকাশনার সময়: ২০২৫ সালের অক্টোবর ১২ (arXiv v2)
পেপার লিঙ্ক: https://arxiv.org/abs/2507.12257

সারসংক্ষেপ

স্টোকাস্টিক সময় শ্রেণীতে কার্যকারণ সম্পর্ক অন্বেষণ একটি চ্যালেঞ্জিং কিন্তু অত্যন্ত গুরুত্বপূর্ণ কাজ, যা আর্থিক, অর্থনৈতিক, স্নায়ুবিজ্ঞান এবং জলবায়ু বিজ্ঞান সহ বিভিন্ন ক্ষেত্রে ব্যাপক প্রয়োগ রয়েছে। যদিও অনেক কার্যকারণ আবিষ্কার (CD) অ্যালগরিদম প্রস্তাব করা হয়েছে, তারা প্রায়শই শব্দের প্রতি অত্যন্ত সংবেদনশীল এবং বাস্তব ডেটায় প্রয়োগ করার সময় বিভ্রান্তিকর কার্যকারণ অনুমান তৈরি করে। এই পেপারটি পর্যবেক্ষণ করে যে সাধারণ বাস্তব-বিশ্বের সময় শ্রেণীর ফ্রিকোয়েন্সি বর্ণনা শক্তি-আইন বিতরণ অনুসরণ করে, যা প্রধানত তাদের অন্তর্নিহিত স্ব-সংগঠিত আচরণের কারণে। এই অন্তর্দৃষ্টির উপর ভিত্তি করে, আমরা শক্তি-আইন বর্ণনামূলক বৈশিষ্ট্য নিষ্কাশনের উপর ভিত্তি করে একটি শক্তিশালী কার্যকারণ আবিষ্কার পদ্ধতি তৈরি করেছি, যা প্রকৃত কার্যকারণ সংকেত প্রসারিত করতে পারে। আমাদের পদ্ধতি সংশ্লেষিত বেঞ্চমার্ক এবং পরিচিত কার্যকারণ কাঠামো সহ বাস্তব-বিশ্বের ডেটাসেটগুলিতে ধারাবাহিকভাবে অত্যাধুনিক বিকল্প পদ্ধতিগুলিকে ছাড়িয়ে যায়, এর শক্তিশালীতা এবং ব্যবহারিক প্রাসঙ্গিকতা প্রদর্শন করে।

গবেষণা পটভূমি এবং প্রেরণা

সমস্যা সংজ্ঞা

এই গবেষণা সময় শ্রেণী ডেটায় কার্যকারণ আবিষ্কার সমস্যা সমাধানের জন্য নিবেদিত, অর্থাৎ পর্যবেক্ষিত ডেটা থেকে ভেরিয়েবলগুলির মধ্যে কার্যকারণ সম্পর্ক চিহ্নিত করা। ঐতিহ্যবাহী কার্যকারণ আবিষ্কার পদ্ধতি, বিশেষত গ্রেঞ্জার কার্যকারণের উপর ভিত্তি করে, বাস্তব-বিশ্বের জটিল ডেটার মুখোমুখি হওয়ার সময় নিম্নলিখিত সীমাবদ্ধতা প্রদর্শন করে:

শব্দ সংবেদনশীলতা: ঐতিহ্যবাহী পদ্ধতি অ-গাউসিয়ান শব্দ, অ-স্থিরতা এবং অ-রৈখিক বিঘ্নের প্রতি অত্যন্ত সংবেদনশীল
অনুমান সীমাবদ্ধতা: শব্দ স্থিরতা এবং একক বৈশিষ্ট্য স্কেলের মতো কঠোর অনুমানের উপর নির্ভর করে
মিথ্যা সম্পর্ক সনাক্তকরণ: শব্দ সম্পর্ককে কার্যকারণ সম্পর্ক হিসাবে ভুল করার প্রবণতা

গবেষণা প্রেরণা

লেখকরা পর্যবেক্ষণ করেছেন যে বাস্তব-বিশ্বের সিস্টেমগুলি সর্বজনীনভাবে শক্তি-আইন ফ্রিকোয়েন্সি বর্ণনা বৈশিষ্ট্য প্রদর্শন করে, যা উদ্ভূত হয়:

পারস্পরিক ক্রিয়াশীল ইউনিটের স্ব-সংগঠিত আচরণ
বাহ্যিক সমন্বয়কারীর অনুপস্থিতি দ্বারা সৃষ্ট স্কেল অপরিবর্তনীয়তা
সিস্টেমের ফ্র্যাক্টাল বৈশিষ্ট্য এবং দীর্ঘ-পরিসীমা সময় সম্পর্ক

এই পর্যবেক্ষণের উপর ভিত্তি করে, এই পেপারটি আরও শক্তিশালী কার্যকারণ আবিষ্কারের জন্য শক্তি-আইন বর্ণনামূলক বৈশিষ্ট্য ব্যবহার করার প্রস্তাব করে।

মূল অবদান

PLaCy ফ্রেমওয়ার্ক প্রস্তাব: শক্তি-আইন বর্ণনামূলক বৈশিষ্ট্যের উপর ভিত্তি করে একটি উপন্যাস কার্যকারণ আবিষ্কার পদ্ধতি
তাত্ত্বিক গ্যারান্টি: ফ্রিকোয়েন্সি ডোমেইন রূপান্তরের অধীনে কার্যকারণ গ্রাফ কাঠামোর অপরিবর্তনীয়তা প্রমাণ করা (উপপাদ্য 1)
পরীক্ষামূলক যাচাইকরণ: সংশ্লেষিত এবং বাস্তব ডেটাসেটে ব্যাপক মূল্যায়ন, উচ্চতর শক্তিশালীতা প্রদর্শন করে
পদ্ধতি সর্বজনীনতা: অন্যান্য কার্যকারণ আবিষ্কার অ্যালগরিদমের উপর বর্ণনামূলক পূর্ব-প্রক্রিয়াকরণের উন্নতিমূলক প্রভাব প্রদর্শন করে

পদ্ধতি বিস্তারিত

কাজের সংজ্ঞা

বহু-ভেরিয়েট সময় শ্রেণী $x \in \mathbb{R}^{L \times d}$ দেওয়া, লক্ষ্য হল নির্দেশিত গ্রাফ $G = (V, E)$ অনুমান করা, যেখানে:

$V = \{1, 2, ..., d\}$ সিস্টেম ভেরিয়েবল প্রতিনিধিত্ব করে
$E \subseteq V \times V$ কার্যকারণ প্রান্ত সেট প্রতিনিধিত্ব করে
নির্দেশিত প্রান্ত $(i,j)$ বিদ্যমান যখন এবং শুধুমাত্র যখন $x_i$ হল $x_j$ এর কার্যকারণ কারণ

মডেল আর্কিটেকচার

1. স্লাইডিং উইন্ডো বিভাজন

প্রতিটি সময় শ্রেণীকে দৈর্ঘ্য $l$ এর ওভারল্যাপিং উইন্ডোতে বিভক্ত করুন, ধাপ আকার $s$ : $w_i^k = (x_i(k \cdot s), ..., x_i(k \cdot s + l - 1))$

2. বর্ণনামূলক বৈশিষ্ট্য নিষ্কাশন

প্রতিটি উইন্ডোতে বিচ্ছিন্ন ফুরিয়ার রূপান্তর প্রয়োগ করুন: $\phi(k) = \sum_{t=0}^{L-1} x(t) e^{-i2\pi k t/L}$

বর্ণনামূলক প্রশস্ততা গণনা করুন: $A(f_k) = |\phi(k)|$

3. শক্তি-আইন ফিটিং

লগ-লগ স্থানে রৈখিক মডেল ফিট করুন: $\log A(f) = a - \lambda \log f$

যেখানে $a$ হল ইন্টারসেপ্ট পরামিতি, $\lambda > 0$ হল বর্ণনামূলক সূচক।

4. কার্যকারণ বিশ্লেষণ

নিষ্কাশিত বর্ণনামূলক পরামিতি সময় শ্রেণী $(a_i, \lambda_i)$ এ বহু-ভেরিয়েট গ্রেঞ্জার কার্যকারণ পরীক্ষা প্রয়োগ করুন, $(\lambda_i, a_i)$ এর $\lambda_j$ এর পূর্বাভাসমূলক ক্ষমতা মূল্যায়ন করুন।

অ্যালগরিদম প্রবাহ (PLaCy)

ইনপুট: সময় শ্রেণী x = (x₁, ..., xₐ), উইন্ডো আকার l, ধাপ আকার s
আউটপুট: কার্যকারণ গ্রাফ G

1. প্রতিটি xᵢ কে ⌊(L-l)/s⌋+1 স্লাইডিং উইন্ডোতে বিভক্ত করুন wᵢᵏ
2. for each i ∈ {1, ..., d} do
3.   for each k ∈ {0, ..., ⌊(L-l)/s⌋} do
4.     wᵢᵏ এ DFT প্রয়োগ করুন φᵢᵏ পেতে
5.     সমীকরণ(2) এর ফিটিং এর মাধ্যমে (aᵢᵏ, λᵢᵏ) পান
6.   (aᵢᵏ, λᵢᵏ) সংযুক্ত করুন সময় শ্রেণী (aᵢ, λᵢ) পেতে
7. for each i,j ∈ {1, ..., d}, i ≠ j do
8.   Gᵢ,ⱼ ← গ্রেঞ্জার কার্যকারণ পরীক্ষা, (aᵢ,λᵢ) কারণ হিসাবে, λⱼ ফলাফল হিসাবে
9. return G

প্রযুক্তিগত উদ্ভাবন পয়েন্ট

ফ্রিকোয়েন্সি ডোমেইন কার্যকারণ আবিষ্কার: প্রথমবারের মতো শক্তি-আইন বর্ণনামূলক বৈশিষ্ট্য ব্যবহার করে কার্যকারণ অনুমানের জন্য সিস্টেমেটিক পদ্ধতি
অভিযোজিত উইন্ডো নির্বাচন: p-মান মানদণ্ডের মাধ্যমে স্বয়ংক্রিয়ভাবে সর্বোত্তম উইন্ডো দৈর্ঘ্য নির্বাচন করুন
শব্দ শক্তিশালীতা: বর্ণনামূলক ফিটিং প্রাকৃতিক ডি-নয়েজিং ধাপ হিসাবে কাজ করে, অ-গাউসিয়ান ওঠানামার প্রতি শক্তিশালীতা উন্নত করে
তাত্ত্বিক ভিত্তি: বর্ণনামূলক রূপান্তরের অধীনে কার্যকারণ গ্রাফ অপরিবর্তনীয়তার তাত্ত্বিক প্রমাণ প্রদান করে

পরীক্ষামূলক সেটআপ

ডেটাসেট

সংশ্লেষিত ডেটাসেট

সাধারণীকৃত Ornstein-Uhlenbeck প্রক্রিয়ার উপর ভিত্তি করে চারটি পরিস্থিতি তৈরি করা: $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$ ): গুণক শব্দ ছাড়াই ভারসাম্যপূর্ণ অবস্থা
OU( $\sigma_g^m > 0$ ): গুণক শব্দ সহ ভারসাম্যপূর্ণ অবস্থা
ÔU( $\sigma_g^m = 0$ ): গুণক শব্দ ছাড়াই অ-ভারসাম্যপূর্ণ অবস্থা
ÔU( $\sigma_g^m > 0$ ): গুণক শব্দ সহ অ-ভারসাম্যপূর্ণ অবস্থা

বাস্তব ডেটাসেট

নদী ডেটাসেট: জার্মানির দক্ষিণে তিনটি জলবিজ্ঞান স্টেশনের নদীর জল স্তর এবং বৃষ্টিপাত ডেটা
বায়ু গুণমান ডেটাসেট: চীনের একাধিক শহরের PM2.5 দূষণ পর্যবেক্ষণ ডেটা

মূল্যায়ন মেট্রিক্স

F1 স্কোর: কার্যকারণ সম্পর্ক সনাক্তকরণের সামগ্রিক কর্মক্ষমতা পরিমাপ করে
সত্য নেতিবাচক হার (TNR): মিথ্যা সংযোগ বাদ দেওয়ার অ্যালগরিদমের ক্ষমতা মূল্যায়ন করে

তুলনা পদ্ধতি

ঐতিহ্যবাহী পদ্ধতি: গ্রেঞ্জার কার্যকারণ, PCMCI, PCMCIΩ
অপ্টিমাইজেশন পদ্ধতি: DYNOTEARS, RCV-VarLiNGAM
গভীর শিক্ষা: Rhino
অ-রৈখিক পদ্ধতি: CCM-Filtering
ফ্রিকোয়েন্সি ডোমেইন পদ্ধতি: BCGeweke, DTF, GewekeNP

বাস্তবায়ন বিবরণ

স্লাইডিং উইন্ডো দৈর্ঘ্য: $l = 50$ (p-মান মানদণ্ড দ্বারা নির্বাচিত)
ধাপ আকার: $s = 1$
বিলম্ব পদ: 10
পরিসংখ্যানগত তাৎপর্য থ্রেশহোল্ড: $p = 0.05$

পরীক্ষামূলক ফলাফল

প্রধান ফলাফল

সংশ্লেষিত ডেটাসেটে কর্মক্ষমতা (N=5, $\sigma_g^a = 1.0$ ):

ডেটাসেট	PLaCy F1	সেরা বেসলাইন F1	PLaCy TNR	সেরা বেসলাইন 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

বাস্তব ডেটাসেট ফলাফল:

ডেটাসেট	PLaCy F1	PLaCy TNR	সেরা বেসলাইন F1	সেরা বেসলাইন TNR
নদী	0.51±0.10	0.75±0.13	0.47±0.07	0.74±0.05
বায়ু গুণমান	0.45±0.04	0.66±0.07	0.44±0.01	0.95±0.02

মূল অনুসন্ধান

গুণক শব্দ শক্তিশালীতা: PLaCy গুণক শব্দ বিদ্যমান পরিস্থিতিতে বিশেষভাবে চমৎকার কর্মক্ষমতা প্রদর্শন করে
অ-ভারসাম্যপূর্ণ অবস্থা অভিযোজনযোগ্যতা: অ-ভারসাম্যপূর্ণ প্রাথমিক শর্তের অধীনেও ভাল কর্মক্ষমতা বজায় রাখে
ফ্রিকোয়েন্সি ডোমেইন পদ্ধতির সুবিধা: সময় ডোমেইন পদ্ধতির তুলনায়, ফ্রিকোয়েন্সি ডোমেইন বিশ্লেষণ উন্নত শব্দ প্রতিরোধ প্রদর্শন করে
সর্বজনীন উন্নতি: PCMCI এবং অন্যান্য পদ্ধতিতে বর্ণনামূলক পূর্ব-প্রক্রিয়াকরণ প্রয়োগ করা কর্মক্ষমতা উল্লেখযোগ্যভাবে উন্নত করতে পারে

অ্যাবলেশন পরীক্ষা

উইন্ডো দৈর্ঘ্য এবং ধাপ আকার বিশ্লেষণ দেখায় যে:

ধাপ আকার 1 হলে কর্মক্ষমতা সর্বোত্তম, স্বল্পমেয়াদী কার্যকারণ নির্ভরতা ক্যাপচার করতে পারে
উইন্ডো দৈর্ঘ্য p-মান মানদণ্ড দ্বারা অভিযোজিত নির্বাচন সেরা ফলাফল দেয়
খুব ছোট বা খুব দীর্ঘ উইন্ডো উভয়ই কর্মক্ষমতা হ্রাস করে

উপসংহার এবং আলোচনা

প্রধান উপসংহার

PLaCy শক্তি-আইন বর্ণনামূলক বৈশিষ্ট্য ব্যবহার করে আরও শক্তিশালী কার্যকারণ আবিষ্কার অর্জন করেছে
পদ্ধতি সংশ্লেষিত এবং বাস্তব ডেটায় উভয়ই উচ্চতর কর্মক্ষমতা প্রদর্শন করে
ফ্রিকোয়েন্সি ডোমেইন বিশ্লেষণ সময় শ্রেণী কার্যকারণ আবিষ্কারের জন্য নতুন দৃষ্টিভঙ্গি প্রদান করে

সীমাবদ্ধতা

বর্ণনা পরিবর্তন ধীর সিস্টেম: যেসব সিস্টেমের বর্ণনামূলক পরামিতি খুব ধীরে পরিবর্তিত হয় তাদের জন্য সীমিত কার্যকারিতা
সংক্ষিপ্ত সময় শ্রেণী: স্থিতিশীল বর্ণনামূলক অনুমানের জন্য যথেষ্ট দীর্ঘ শ্রেণীর প্রয়োজন
গণনামূলক জটিলতা: সাধারণ পদ্ধতির তুলনায় অতিরিক্ত বর্ণনামূলক বিশ্লেষণ ওভারহেড

ভবিষ্যত দিকনির্দেশনা

অ-VAR কার্যকারণ আবিষ্কার পদ্ধতিতে সম্প্রসারণ
বর্ণনামূলক ঘনত্বের পরিসংখ্যানগত পরামিতি গভীরভাবে অধ্যয়ন করা
সম্ভাব্য মিশ্রিত কারণের প্রভাব পরিচালনা করা
আরও দক্ষ অনলাইন কার্যকারণ আবিষ্কার অ্যালগরিদম বিকাশ করা

গভীর মূল্যায়ন

শক্তি

শক্তিশালী উদ্ভাবনী: প্রথমবারের মতো শক্তি-আইন বর্ণনামূলক বৈশিষ্ট্য কার্যকারণ আবিষ্কারে সিস্টেমেটিকভাবে ব্যবহার করা
দৃঢ় তত্ত্ব: কঠোর তাত্ত্বিক বিশ্লেষণ এবং প্রমাণ প্রদান করে
ব্যাপক পরীক্ষা: একাধিক সংশ্লেষিত পরিস্থিতি এবং বাস্তব প্রয়োগ অন্তর্ভুক্ত করে
উচ্চ ব্যবহারিক মূল্য: শব্দ পরিবেশে উল্লেখযোগ্য সুবিধা প্রদর্শন করে

অপূর্ণতা

প্রয়োগের পরিসীমা: প্রধানত শক্তি-আইন বর্ণনা বৈশিষ্ট্য সহ সিস্টেমের জন্য প্রযোজ্য
পরামিতি নির্বাচন: উইন্ডো দৈর্ঘ্য এবং অন্যান্য পরামিতির নির্বাচনের জন্য অভিজ্ঞতা প্রয়োজন
গণনামূলক দক্ষতা: সাধারণ পদ্ধতির তুলনায় গণনামূলক ওভারহেড বেশি

প্রভাব

একাডেমিক অবদান: সময় শ্রেণী কার্যকারণ আবিষ্কারের জন্য নতুন গবেষণা দিকনির্দেশনা প্রদান করে
ব্যবহারিক মূল্য: আর্থিক, জলবায়ু এবং অন্যান্য শক্তি-আইন বৈশিষ্ট্য সহ ক্ষেত্রে ব্যাপক প্রয়োগ সম্ভাবনা
পুনরুৎপাদনযোগ্যতা: সম্পূর্ণ অ্যালগরিদম বর্ণনা এবং ওপেন সোর্স কোড প্রদান করে

প্রযোজ্য পরিস্থিতি

আর্থিক বাজার ডেটা বিশ্লেষণ
জলবায়ু সিস্টেম মডেলিং
স্নায়ুবিজ্ঞান গবেষণা
সামাজিক নেটওয়ার্ক বিশ্লেষণ
যেকোনো স্ব-সংগঠিত বৈশিষ্ট্য সহ জটিল সিস্টেম

সংদর্ভ

পেপারটি কার্যকারণ আবিষ্কার, সময় শ্রেণী বিশ্লেষণ এবং জটিল সিস্টেম সহ একাধিক ক্ষেত্রের গুরুত্বপূর্ণ কাজ কভার করে 51টি সম্পর্কিত সংদর্ভ উদ্ধৃত করে, গবেষণার জন্য একটি দৃঢ় তাত্ত্বিক ভিত্তি প্রদান করে।

সামগ্রিক মূল্যায়ন: এটি সময় শ্রেণী কার্যকারণ আবিষ্কার ক্ষেত্রে একটি উদ্ভাবনী পদ্ধতি প্রস্তাব করে একটি উচ্চ-মানের গবেষণা পেপার। বাস্তব-বিশ্বের সিস্টেমের শক্তি-আইন বর্ণনা বৈশিষ্ট্য দক্ষতার সাথে ব্যবহার করে, এটি কার্যকারণ আবিষ্কারের শক্তিশালীতা সফলভাবে উন্নত করেছে। তাত্ত্বিক বিশ্লেষণ কঠোর, পরীক্ষামূলক ডিজাইন যুক্তিসঙ্গত এবং ফলাফল প্রভাবশালী। এই কাজ জটিল সিস্টেমে কার্যকারণ অনুমানের জন্য নতুন সরঞ্জাম এবং দৃষ্টিভঙ্গি প্রদান করে।