2025-11-24T16:10:25.080119

Using Information Geometry to Characterize Higher-Order Interactions in EEG

Albers, Marriott, Tatsuno

In neuroscience, methods from information geometry (IG) have been successfully applied in the modelling of binary vectors from spike train data, using the orthogonal decomposition of the Kullback-Leibler divergence and mutual information to isolate different orders of interaction between neurons. While spike train data is well-approximated with a binary model, here we apply these IG methods to data from electroencephalography (EEG), a continuous signal requiring appropriate discretization strategies. We developed and compared three different binarization methods and used them to identify third-order interactions in an experiment involving imagined motor movements. The statistical significance of these interactions was assessed using phase-randomized surrogate data that eliminated higher-order dependencies while preserving the spectral characteristics of the original signals. We validated our approach by implementing known second- and third-order dependencies in a forward model and quantified information attenuation at different steps of the analysis. This revealed that the greatest loss in information occurred when going from the idealized binary case to enforcing these dependencies using oscillatory signals. When applied to the real EEG dataset, our analysis detected statistically significant third-order interactions during the task condition despite the relatively sparse data (45 trials per condition). This work demonstrates that IG methods can successfully extract genuine higher-order dependencies from continuous neural recordings when paired with appropriate binarization schemes.

academic

情報幾何を用いたEEGにおける高次相互作用の特性化

基本情報

論文ID: 2510.14188
タイトル: Using Information Geometry to Characterize Higher-Order Interactions in EEG
著者: Eric Albers, Paul Marriott, Masami Tatsuno
分類: q-bio.NC (神経と認知), q-bio.QM (定量的方法)
発表日時: 2025年10月16日 (arXivプレプリント)
論文リンク: https://arxiv.org/abs/2510.14188

要旨

本研究は、情報幾何(Information Geometry, IG)手法を従来の二値脊椎波列データから連続的な脳電図(EEG)信号分析へと拡張するものである。Kullback-Leibler散度と相互情報の直交分解を通じて、ニューロン間の異なる次数の相互作用を識別する。本研究は、運動想像実験における三次相互作用を識別するための3つの二値化手法を開発し、位相ランダム化代替データを用いて統計的有意性を評価した。前向きモデル検証法により手法の有効性を検証し、分析の各段階における情報減衰を定量化した。結果として、データが比較的疎である(条件あたり45試行)にもかかわらず、本手法はタスク条件下で統計的に有意な三次相互作用を検出できることが示された。

研究背景と動機

問題の定義

従来の神経科学研究は主に脳領域間の対関係(二次相互作用)に焦点を当ててきたが、複雑系としての脳は対関係を超えた高次相互作用を有する可能性がある。既存の機能的結合ネットワークは対相関に基づいて構築されており、脳の情報処理の複雑性を完全に捉えられない可能性がある。

重要性

理論的意義: 認知機能の実行に三次以上の相互作用が必要であるかどうかを理解すること
方法論的意義: 離散的な脊椎波データから連続的なEEG信号への情報幾何手法の拡張
応用価値: 脳機械インターフェースと神経疾患診断のための新しい分析ツールの提供

既存手法の限界

情報幾何手法: 主に二値脊椎波データに適用され、連続信号に対する効果的な離散化戦略が不足している
従来的なEEG分析: 主に対相関に基づいており、高次依存関係を無視している
統計推論: 疎なデータ条件下では、標準漸近ツール(χ²分布など)が適用できない可能性がある

研究動機

脊椎波分析に成功した情報幾何手法をEEGデータに拡張し、連続神経記録における真の高次依存関係を捉えるための適切な二値化戦略を開発すること。

核心的貢献

方法論的革新: 連続EEG信号を情報幾何分析に適した二値表現に変換するための3つの二値化手法(Sign、Diff、Power)を開発
検証フレームワーク: 位相ランダム化代替データに基づく統計的有意性検定手法を確立
前向きモデリング: 既知の二次および三次依存関係を持つ前向きモデルを実装し、分析過程における情報減衰を定量化
実証的発見: 運動想像EEGデータにおいて統計的に有意な三次相互作用を検出
理論的洞察: 理想化された二値ケースから振動信号への実装への移行時に最大の情報損失が発生することを明らかにした

方法の詳細

タスク定義

入力: マルチチャネルEEG連続信号出力: チャネル三つ組間の一次、二次、三次相互情報成分制約: 疎なデータ(条件あたり45試行)と連続信号の離散化の課題に対処

情報幾何理論の基礎

3つの二値変数X₁、X₂、X₃について、結合確率分布は8つの確率のベクトルとして表現できる:

p = (p₀₀₀, p₀₀₁, p₀₁₀, p₀₁₁, p₁₀₀, p₁₀₁, p₁₁₀, p₁₁₁)

期待パラメータηの座標系:

η₁, η₂, η₃: 周辺活性化率
η₁₂, η₁₃, η₂₃: 対活性化率
η₁₂₃: 三元活性化率

自然パラメータθの座標系は対数比により定義される。例えば:

θ₁₂₃ = log(p₀₀₁p₀₁₀p₁₀₀p₁₁₁)/(p₁₁₀p₁₀₁p₀₁₁p₀₀₀)

KL散度の直交分解

混合座標系を使用すると、KL散度は直交分解できる:

D[p : q] = D[p : p̄] + D[p̄ : p̃] + D[p̃ : q]

ここで:

Dp : p̄: 三元相互情報
Dp̄ : p̃: 対相互情報
Dp̃ : q: 活性化率調整情報

二値化手法

1. Sign手法

binary_signal = 1 if EEG_signal > 0 else 0

粗い位相情報を捉え、振幅を無視する。

2. Diff手法

diff_signal = diff(EEG_signal)
binary_signal = 1 if diff_signal > 0 else 0

位相転移パターンを捉える。

3. Power手法

power = EEG_signal²
envelope = moving_average(power, 30_samples)
z_scores = (envelope - mean) / std
binary_signal = 1 if z_scores > 1 else 0

高振幅期間を捉え、位相に依存しない。

統計的有意性検定

検定統計量を使用:

λ = 2N·D[p : p̄] ~ χ²(1)

データが疎であるため、χ²近似は不適切であり、IAAFT(反復振幅調整フーリエ変換)代替データに基づく非パラメトリック検定を採用した。

実験設定

データセット

OpenNeuro運動想像データセット (Triana-Guzman et al., 2022):

参加者: 32名の健康被験者(女性16名)
電極: 17個の電極、国際10-20システムに従って配置
サンプリングレート: 250 Hz
試行設計:
- 6ブロック(座位3ブロック、立位3ブロック)
- 各ブロック30試行(運動想像15試行、アイドル状態15試行)
- 条件あたり合計45試行

試行構造:

注視(4秒): 画面の十字を注視
観察(3秒): 実行予定のタスクを表示
想像(4秒): 心理的タスクを実行(運動想像またはアイドル状態)
休息(4秒): 自由活動

データ前処理

フィルタリング: 0.5 Hz高域フィルタ、58-62 Hz陷波フィルタ
アーティファクト除去: ASR(アーティファクト部分空間再構成)手法を使用
周波数帯域フィルタリング: Delta(0.5-4Hz)、Theta(4-8Hz)、Alpha(8-12Hz)、Beta(12-30Hz)、Gamma(30-60Hz)に分割
エポック抽出: 想像タスク開始前7秒から開始後4秒までの11秒エポック