JOVANA
Explore Library Glossary Getting Started Three Levels Fields How it works Mission
Join the mission
All guides

有向與源空間連結:研究前沿

格蘭傑因果、頻譜有向量測、以波束成形前往源空間,以及對「為何網路宣稱如此難以信任」的誠實檢視。

誰在驅動誰?

對稱量測告訴你兩個腦區有耦合;更困難、也更有趣的問題是方向性——是 X 驅動 Y,還是相反?回答它意味著從功能性走向有效/有向連結,其主力框架是建立在自迴歸模型上的格蘭傑因果

格蘭傑因果與 VAR 模型

其核心是預測:若 X 的過去能在 Y 自身過去之外,進一步改善對 Y 的預測,則 X 格蘭傑導致 Y。你擬合一個向量自迴歸(VAR)模型,比較納入與不納入另一通道歷史時的殘差變異數。

\mathbf{x}(t) = \sum_{k=1}^{p} \mathbf{A}_k\,\mathbf{x}(t-k) + \mathbf{e}(t)

VAR(p) 模型:當前活動是過去 p 個延遲的線性函數加上創新項 e(t)。其係數 A_k 承載了所有方向性資訊。

向量自回歸(VAR)模型用所有通道近期過去的加權組合來預測每個通道的當前值。這些權重捕捉了訊號如何隨時間互相驅動——這是有向連結的基礎。

\mathbf{x}(t)
當前時刻的多通道訊號。
\mathbf{A}_k
延遲 k 的係數矩陣——誰預測誰。
p
模型階數——用了多少過去步。
\mathbf{e}(t)
創新項——無法被預測的部分。
F_{X\to Y} = \ln\frac{\operatorname{var}\!\left(e_Y \mid Y^{-}\right)}{\operatorname{var}\!\left(e_Y \mid Y^{-},\,X^{-}\right)}

格蘭傑因果 X→Y:不納入與納入 X 過去時,Y 預測誤差變異數的對數比。為零代表無有向影響。

若知道 X 的過去能在 Y 自身過去之外進一步幫助預測 Y,就說 X 格蘭傑造成 Y。此式比較納入與不納入 X 歷史時 Y 的預測誤差變異數;若加入 X 使誤差變小,對數比就為正——即有向影響的證據。

F_{X\to Y}
從 X 到 Y 的格蘭傑因果。
\operatorname{var}(e_Y\mid Y^{-})
只用 Y 自身過去時,Y 的預測誤差變異數。
\operatorname{var}(e_Y\mid Y^{-},X^{-})
也用上 X 過去時,Y 的誤差變異數。
\ln
對數比——為零代表無有向影響。

頻譜有向量測:PDC 與 DTF

神經驅動具頻率特異性,因此我們要逐頻帶的方向性。把 VAR 係數轉換到頻域,就得到偏有向相干(PDC)與有向轉移函數(DTF)。PDC 分離出直接j\to i 影響;DTF 則捕捉總體(含間接)的影響。

\mathrm{PDC}_{j\to i}(f) = \frac{\lvert \bar{A}_{ij}(f)\rvert}{\sqrt{\sum_{k}\lvert \bar{A}_{kj}(f)\rvert^{2}}}, \qquad \bar{A}(f)=I-\sum_{k=1}^{p}\mathbf{A}_k\,e^{-\,i\,2\pi f k}

PDC 以源通道 j 的流出量對頻域 VAR 係數作正規化,給出有向、逐頻率的量測。

偏有向同調(PDC)把 VAR 模型帶到頻域,並在每個頻率問:通道 j 的流出量有多少專門流向通道 i。以 j 的總流出量作歸一化,就得到一個有向、逐頻率的連結度量。

\mathrm{PDC}_{j\to i}(f)
在頻率 fji 的有向影響。
\bar{A}(f)
轉入頻域的 VAR 係數。
\bar{A}_{ij}(f)
該頻域矩陣中 j\to i 的元素。
\sqrt{\sum_k|\bar{A}_{kj}(f)|^{2}}
正規化因子——通道 j 的總流出量。

源空間連結與波束成形器

對場擴散更根本的修正,是根本不在感測器上計算連結。用一個空間濾波器——波束成形器——重建源時間序列,再估計源與源之間的連結。LCMV 波束成形器在時域中,於目標位置維持單位增益的約束下最小化輸出功率來達成此事;其頻域孿生體 DICS 波束成形器則鎖定特定頻帶以分析振盪耦合。

\mathbf{W} = \bigl(\mathbf{L}^{\top}\Sigma^{-1}\mathbf{L}\bigr)^{-1}\mathbf{L}^{\top}\Sigma^{-1}

LCMV 權重:在導程場 L 所定義的源上維持單位增益,並(透過資料共變異數 Σ)達成最小輸出功率。

為每個源建立一個虛擬感測器的波束成形器權重:在對目標源導程場保持單位增益的同時,最小化總輸出功率。在這些源空間訊號上做連結分析,能避開許多感測器層級的洩漏。

\mathbf{W}
波束成形器權重矩陣。
\mathbf{L}
導程場——源到感測器的映射。
\Sigma
感測器資料共變異數。
(\mathbf{L}^{\top}\Sigma^{-1}\mathbf{L})^{-1}
在源上強制單位增益的正規化因子。

這有幫助,卻無法根治問題:逆問題是不適定的,因此波束成形器重建出的源仍會彼此滲漏——殘餘的源洩漏。因此頂尖管線會再加上洩漏校正,例如在計算幅值耦合前,對源訊號作成對正交化。

開放問題與誠實的但書

前沿正走向何方?走向與動力系統模型融合、而非靜態圖的連結分析;走向多模態的錨定(以同步記錄或刺激來驗證所宣稱的連結);以及走向 BCI 應用——用網路生物標記驅動閉迴路神經調控,並運用比原始功率更能跨場次泛化的連結特徵。

兩半,一條管線

請留意,這條軌道的兩個主題是不可分割的。第 1 至 3 篇的空間濾波器——拉普拉斯、CSP、ICA、波束成形器——正是讓第 4 至 5 篇的連結分析可信的工具,因為它們在你追問源如何交談之前,先把源解混。空間濾波與功能性連結並非兩個主題,而是同一條管線的正反兩面:從原始通道走到站得住腳的腦網路。