用深度柯尔莫哥洛夫算子发现非线性系统的因果关系
Deep Koopman operator framework for causal discovery in nonlinear dynamical systems
- 基于深度学习自动寻找最优可观测变量
- 在希尔伯特空间中通过动态距离量化因果强度
- 适用于气候等复杂非线性系统,可解释性强
我们提出一种基于深度柯尔莫哥洛夫算子理论的新型因果发现算法Kausal。传统统计方法如格兰杰因果无法有效处理非线性动力系统中的复杂反馈、时标混合和非平稳性问题,难以揭示真实因果机制。而柯尔莫哥洛夫算子方法能将非线性动力系统近似映射到可观测变量的线性空间中。Kausal利用深度学习自动推断最优可观测变量,并在再生核希尔伯特空间中评估因果关系,定义为效应变量的边缘动态与因果-效应联合动态之间的距离。数值实验表明,相较于依赖预设可观测变量的方法,Kausal在识别和刻画因果信号方面具有显著优势。最后,我们将算法应用于厄尔尼诺-南方涛动(El Niño-Southern Oscillation)观测数据,验证了其在真实世界现象中的适用性。代码已开源。
原文摘要 · Abstract (English)
We use a deep Koopman operator-theoretic formalism to develop a novel causal discovery algorithm, Kausal. Causal discovery aims to identify cause-effect mechanisms for better scientific understanding, explainable decision-making, and more accurate modeling. Standard statistical frameworks, such as Granger causality, lack the ability to quantify causal relationships in nonlinear dynamics due to the presence of complex feedback mechanisms, timescale mixing, and nonstationarity. This presents a challenge in studying many real-world systems, such as the Earth's climate. Meanwhile, Koopman operator methods have emerged as a promising tool for approximating nonlinear dynamics in a linear space of observables. In Kausal, we propose to leverage this powerful idea for causal analysis where optimal observables are inferred using deep learning. Causal estimates are then evaluated in a reproducing kernel Hilbert space, and defined as the distance between the marginal dynamics of the effect and the joint dynamics of the cause-effect observables. Our numerical experiments demonstrate Kausal's superior ability in discovering and characterizing causal signals compared to existing approaches of prescribed observables. Lastly, we extend our analysis to observations of El Niño-Southern Oscillation highlighting our algorithm's applicability to real-world phenomena. Our code is available at https://github.com/juannat7/kausal.
Thank you to arXiv for use of its open access interoperability. PaperDance 不是 arXiv 官方产品;中文卡片由大模型生成,请以原文为准。