大腦狀態(tài)動力學(xué)研究分析

劉亮穎

摘要：基于圖論分析的腦網(wǎng)絡(luò)研究表明腦網(wǎng)絡(luò)具有小世界特性，該文介紹基于最小生成樹的K-Means算法將全腦按功能劃分腦區(qū)，并利用隱馬爾科夫模型揭示隱藏的大腦狀態(tài)以及狀態(tài)間轉(zhuǎn)換的規(guī)律性，為了衡量網(wǎng)絡(luò)內(nèi)聚力的程度，采用奇異值分解與滑動窗口結(jié)合的算法。揭示了在靜息態(tài)下，網(wǎng)絡(luò)間的轉(zhuǎn)換是非任意的，呈現(xiàn)出兩個不同的亞穩(wěn)態(tài)網(wǎng)絡(luò)系統(tǒng)。同時提供了一種衡量網(wǎng)絡(luò)內(nèi)部凝聚力的方法。

關(guān)鍵詞： 復(fù)雜腦網(wǎng)絡(luò); 最小生成樹;K-Means算法; 隱馬爾科夫模型; 網(wǎng)絡(luò)內(nèi)聚力

中圖分類號：TP311? ? ? ? 文獻(xiàn)標(biāo)識碼：A? ? ? ? 文章編號：1009-3044（2018）36-0179-01

復(fù)雜腦網(wǎng)絡(luò)研究證明腦網(wǎng)絡(luò)具有“經(jīng)濟(jì)的”小世界特性，即這種網(wǎng)絡(luò)具有大的聚類系數(shù)和小的特征路徑長度，因此小世界網(wǎng)絡(luò)可以用圖來表征，“節(jié)點(diǎn)”和“邊”是構(gòu)成圖的兩個要素，這里將每一個腦區(qū)（ROI）作為節(jié)點(diǎn)，腦區(qū)間的連接作為邊。腦區(qū)內(nèi)部節(jié)點(diǎn)高密度連接，腦區(qū)間存在稀疏的長連接。K-Means算法是一種經(jīng)典的劃分聚類算法，將性質(zhì)相似的數(shù)據(jù)聚集在一起，不同的簇之間差別較大，因此提供了劃分ROI的一種思想。但由于K的值難以預(yù)先確定，故本文提出用最小生成樹算法估計(jì)K值，改進(jìn)K-Means算法的不足。

腦功能連接的概念最早出現(xiàn)在 EEG研究中，它度量空間上分離的不同腦區(qū)間時間上的相關(guān)性和功能活動的統(tǒng)計(jì)依賴關(guān)系，是描述腦區(qū)之間協(xié)同工作模式的有效手段之一。從功能性的角度分析，研究發(fā)現(xiàn)網(wǎng)絡(luò)之間的轉(zhuǎn)換是非隨機(jī)的，某些網(wǎng)絡(luò)更可能在其他網(wǎng)絡(luò)之后發(fā)生。隱馬爾科夫模型（HMM）提供了一種概率（生成）模型，通過一個貝葉斯推理的單一過程，對時間周期和時間周期包含的方式進(jìn)行建模，從而揭示不同網(wǎng)絡(luò)間的轉(zhuǎn)換規(guī)律。

網(wǎng)絡(luò)間節(jié)點(diǎn)的緊密程度越高，代表進(jìn)行活動的效率越高，代價越小。為了衡量網(wǎng)絡(luò)內(nèi)聚力的程度，采用奇異值分解與滑動窗口結(jié)合的算法，產(chǎn)生酉矩陣，其對角線上的數(shù)據(jù)代表每個變量的方差。

本文主要提出劃分腦區(qū)的一種方法，并介紹現(xiàn)有的關(guān)于揭示網(wǎng)絡(luò)間轉(zhuǎn)化規(guī)律的研究工具，以及衡量網(wǎng)絡(luò)內(nèi)聚力的手段，全面分析當(dāng)前大腦狀態(tài)動力學(xué)的有關(guān)研究。

1 基于最小生成樹的K-Means算法

使用K-Means算法需要預(yù)先確定K的值，故本文利用最小生成樹算法估計(jì)K值，再利用K-Means算法將數(shù)據(jù)進(jìn)行分類，以達(dá)到劃分腦區(qū)或ROI的效果。

1.1 Kruskal算法

1） 將全部邊e1、e2、…、en-1按照權(quán)值由小到大排序。

2）按順序（邊權(quán)由小到大的順序）考慮每條邊，若這條邊ei和我們已經(jīng)選擇的邊不構(gòu)成圈，就保留這條邊，否則放棄這條邊。

3） 重復(fù)2），直至選擇n-1條邊。

1.2K-Means算法

1） 在n個數(shù)據(jù)點(diǎn)中隨機(jī)選取k個對象作為初始聚類中心，每個對象代表一個類別的中心。

2） 根據(jù)歐式距離，將數(shù)據(jù)對象分配給距離最小的聚類中心。

3） 重新計(jì)算聚類中心，計(jì)算每個類別中所有對象的均值作為該類別的新聚類中心，計(jì)算所有樣本到其所在類別聚類中心的距離平方和，即J（C）值。

4） 重復(fù)步驟2）3），直到達(dá)到迭代次數(shù)或J（C）值不再變化。

1.3 基于最小生成樹的K-Means算法

1） 對于n個數(shù)據(jù)點(diǎn)，利用kruskal算法生成一顆最小樹。

2） 選定最大距離即閾值，將最小生成樹分割成m棵子樹，即m個簇。

3） 將m個簇按數(shù)據(jù)量從大到小排列，選取前k個簇看作初始的簇。

4） 利用K-Means算法迭代計(jì)算。

2 基于隱馬爾科夫模型的大腦狀態(tài)轉(zhuǎn)換

為了研究大腦的大范圍的網(wǎng)絡(luò)動力學(xué)的時間組織，使用了一種設(shè)計(jì)用來發(fā)現(xiàn)隨時間重復(fù)的網(wǎng)絡(luò)的方法（稱為大腦狀態(tài)）。我們將網(wǎng)絡(luò)定義為表示概率分布，不僅具有不同的激活模式，還具有不同的功能連接模式。

1） 利用隱馬爾科夫模型（HMM）推斷出腦的隱藏狀態(tài)，HMM狀態(tài)表示在不同時間點(diǎn)重復(fù)出現(xiàn)的具有不同活動和功能連接的獨(dú)特大腦網(wǎng)絡(luò)（其中功能連接被定義為腦區(qū)域間隨時間的皮爾遜相關(guān)）。

2） 計(jì)算狀態(tài)轉(zhuǎn)移概率矩陣，它指定了從任何狀態(tài)轉(zhuǎn)換到另一個狀態(tài)的概率。

3） 計(jì)算訪問每個狀態(tài)的時間比例FO，得到FO相關(guān)性矩陣。

3 基于SVD算法的網(wǎng)絡(luò)內(nèi)聚力

1） 在Matlab中對每個網(wǎng)絡(luò)在10個TRs （23.4 s）的滑動時間窗口內(nèi)的時間序列數(shù)據(jù)執(zhí)行SVD，得到每個參與者、網(wǎng)絡(luò)和時間窗口的酉矩陣。

2） 酉矩陣由奇異值組成，矩陣對角線表示每個分量的方差，第一個奇異值S1，1，代表網(wǎng)絡(luò)的凝聚力，對角線上奇異值的總和[i=110S]i，i量化了網(wǎng)絡(luò)中的總方差，為每個網(wǎng)絡(luò)和時間窗口提供解釋的百分比方差（PVE）。

3） S1，1 / [i=110S]，得每個網(wǎng)絡(luò)的內(nèi)聚力。

4 結(jié)論

本文提出基于最小生成樹的K-Means算法劃分腦區(qū)或ROI。利用隱馬爾科夫模型揭示了大腦不同網(wǎng)絡(luò)之間的轉(zhuǎn)換不是完全隨機(jī)發(fā)生的，這種非隨機(jī)的大腦網(wǎng)絡(luò)序列本身是有層次結(jié)構(gòu)的，呈現(xiàn)出兩個不同的亞穩(wěn)態(tài)網(wǎng)絡(luò)系統(tǒng)?；赟VD算法計(jì)算出網(wǎng)絡(luò)內(nèi)聚力的大小，全面分析當(dāng)前大腦狀態(tài)動力學(xué)的有關(guān)研究。

參考文獻(xiàn)：

[1] 孫俊峰.復(fù)雜腦網(wǎng)絡(luò)研究進(jìn)展——結(jié)構(gòu)、功能、計(jì)算和應(yīng)用[J].復(fù)雜系統(tǒng)與復(fù)雜性科學(xué)，2010， 7（4）：74-90.

[2] 田麗霞.基于圖論的復(fù)雜腦網(wǎng)絡(luò)分析[J].北京生物醫(yī)學(xué)工程，2010（1）.

[3] 蔣田仔.腦網(wǎng)絡(luò)：從腦結(jié)構(gòu)到腦功能[J].生命科學(xué)，2009（2）.

[4] Diego Vidaurrea. Brain network dynamics are hierarchically organized in time [J].PANS，2017（5）.

[5] 王千.K-means 聚類算法研究綜述[J].電子設(shè)計(jì)工程，2012（7）：21-24.

[6] Christina B. Young. Dynamic Shifts in Large-Scale Brain Network Balance As a Function of Arousal [J].The Journal of Neuroscience，2017（3）.

[7] 肖雋瑋. 基于fMRI數(shù)據(jù)的腦功能網(wǎng)絡(luò)聚類研究[D]. 南京理工大學(xué)，2017.

[通聯(lián)編輯：唐一東]