• 
    

    
    

      99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看

      ?

      機(jī)器學(xué)習(xí)方法在fMRI數(shù)據(jù)中的應(yīng)用研究

      2023-07-14 14:50:15黃慶坤賀政剛楊鵬
      計(jì)算機(jī)應(yīng)用文摘 2023年12期
      關(guān)鍵詞:復(fù)雜網(wǎng)絡(luò)統(tǒng)計(jì)分析機(jī)器學(xué)習(xí)

      黃慶坤 賀政剛 楊鵬

      摘要:文章通過(guò)機(jī)器學(xué)習(xí)方法對(duì)功能性核磁共振成像(fMRI)數(shù)據(jù)進(jìn)行分類研究.旨在探討精神分裂患者的差異性腦網(wǎng)絡(luò)結(jié)構(gòu)。在特征選擇上,采用了復(fù)雜網(wǎng)絡(luò)分析來(lái)度量腦網(wǎng)絡(luò)的拓?fù)鋵傩?,并且以年齡為協(xié)變量通過(guò)雙樣本£檢驗(yàn)構(gòu)建顯著差異性特征集;在分類算法上,采用粒子群優(yōu)化支持向量機(jī)分類算法,在10折交叉驗(yàn)證下進(jìn)行建模分析。結(jié)果表明精神分裂患者(SCH)組和正常被試(NC)組在全局屬性和節(jié)點(diǎn)屬性上存在顯著差異,在分類模型下能夠得到最佳準(zhǔn)確率93%。

      關(guān)鍵詞:機(jī)器學(xué)習(xí);功能性核磁共振成像;統(tǒng)計(jì)分析;復(fù)雜網(wǎng)絡(luò)

      中圖法分類號(hào):TPI83 文獻(xiàn)標(biāo)識(shí)碼:A

      1 引言

      人在處理事件或完成某項(xiàng)任務(wù)時(shí),大腦需要依賴各個(gè)功能腦區(qū)的相互協(xié)作,盡管在休息時(shí),大腦仍需要消耗大量的能量來(lái)進(jìn)行信息處理。從各個(gè)功能腦區(qū)相互作用所形成的復(fù)雜腦網(wǎng)絡(luò)這一角度研究大腦,能進(jìn)一步闡述腦內(nèi)信息加工機(jī)制,對(duì)相關(guān)腦疾病診斷具有較大的理論價(jià)值。近年來(lái),FMRI 技術(shù)已廣泛應(yīng)用于腦功能網(wǎng)絡(luò)的研究,Rubinov 等[1] 通過(guò)復(fù)雜網(wǎng)絡(luò)分析來(lái)量化大腦網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu);Wang 等[2] 通過(guò)靜息態(tài)fMRI 數(shù)據(jù)研究腦功能網(wǎng)絡(luò)拓?fù)渲笜?biāo)并對(duì)TRT 參數(shù)進(jìn)行評(píng)估;Onias 等[3] 概述了復(fù)雜網(wǎng)絡(luò)的統(tǒng)計(jì)參數(shù)并對(duì)癲癇患者的腦網(wǎng)絡(luò)結(jié)構(gòu)進(jìn)行度量;王靜等[4] 對(duì)抑郁癥患者大腦的復(fù)雜網(wǎng)絡(luò)研究進(jìn)行綜述等。在先前的研究中[5] ,研究了大腦感興趣區(qū)域(ROI)的功能連接(FC),將其作為特征向量,通過(guò)PCA 算法構(gòu)建特征集并進(jìn)行分類實(shí)驗(yàn),其本質(zhì)上屬于種子點(diǎn)分析法,僅通過(guò)建模分類來(lái)說(shuō)明在精神分裂患者中存在差異性腦區(qū)或腦網(wǎng)絡(luò)機(jī)制,但未深入研究大腦網(wǎng)絡(luò)結(jié)構(gòu)的拓?fù)鋵傩砸约案鱾€(gè)節(jié)點(diǎn)之間存在的差異性。對(duì)此,本文首先通過(guò)復(fù)雜網(wǎng)絡(luò)統(tǒng)計(jì)量描述來(lái)度量SCH 組與NC 組的腦網(wǎng)絡(luò)結(jié)構(gòu);其次將年齡作為協(xié)變量,通過(guò)統(tǒng)計(jì)檢驗(yàn)方法討論精神分裂患者的差異性腦網(wǎng)絡(luò)結(jié)構(gòu)并構(gòu)建特征集;最后使用粒子群優(yōu)化SVM 分類模型進(jìn)行建模論證,通過(guò)機(jī)器學(xué)習(xí)方法進(jìn)一步解釋精神分裂患者的潛在患病機(jī)理。

      2 相關(guān)理論

      2.1 復(fù)雜網(wǎng)絡(luò)的統(tǒng)計(jì)描述

      在大腦中,將以ROI 為節(jié)點(diǎn)所構(gòu)成的復(fù)雜網(wǎng)絡(luò)中根據(jù)全局指標(biāo)(Global)和節(jié)點(diǎn)指標(biāo)(Node)作為網(wǎng)絡(luò)度量方法[6] 。全局指標(biāo)根據(jù)整個(gè)網(wǎng)絡(luò)結(jié)構(gòu)對(duì)相關(guān)統(tǒng)計(jì)進(jìn)行描述, 具體包括: 聚類系數(shù)( Clusteringcoefficient)、最短路徑長(zhǎng)度(Shortest path length)、小世界屬性(Small world)、全局效率(Global efficiency)、同配性(Assortativity)、同步性(Synchronization)、層次(Hierarchy)系數(shù)[7] 。節(jié)點(diǎn)指標(biāo)根據(jù)節(jié)點(diǎn)信息進(jìn)行計(jì)算,具體包括:節(jié)點(diǎn)聚類系數(shù)、節(jié)點(diǎn)最短路徑長(zhǎng)度、節(jié)點(diǎn)效率( Nodal efficiency )、度中心性( DegreeCentrality)、介數(shù)中心性(Betweenness Centrality)。

      2.2 粒子群算法優(yōu)化分類

      3 實(shí)驗(yàn)

      3.1 實(shí)驗(yàn)平臺(tái)

      實(shí)驗(yàn)基于Matlab 平臺(tái)R2021a 版本和Pycharm 平臺(tái)2021.2 版本。Matlab 中主要使用的工具箱為gretnav2. 0 版本[8] 。Pycharm 中主要使用Anaconda3,nilearn0.8.1 版本,sklearn0.24.2 版本以及數(shù)據(jù)處理庫(kù)和繪圖庫(kù)等。

      3.2 實(shí)驗(yàn)數(shù)據(jù)

      實(shí)驗(yàn)數(shù)據(jù)來(lái)源于生物醫(yī)學(xué)研究卓越中心提供的精神分裂癥數(shù)據(jù)集,實(shí)驗(yàn)選擇74 名NC 和72 名SCH的rs?fMRI(fmri_XXXXXXX.nii.gz)和被試采集記錄構(gòu)建樣本數(shù)據(jù)。通過(guò)校正、標(biāo)準(zhǔn)化、平滑、去線性趨勢(shì)、回歸協(xié)變量、濾波等程序?qū)Γ妫停遥?數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)預(yù)處理。

      3.3 實(shí)驗(yàn)流程

      實(shí)驗(yàn)通過(guò)Python 編程語(yǔ)言進(jìn)行代碼簡(jiǎn)寫(xiě),如圖1所示。(1) 劃分樣本數(shù)據(jù)。通過(guò)cobre 數(shù)據(jù)集中phenotypic_data.tsv 文件將fMRI 數(shù)據(jù)分為NC 組(74名正常被試NC.nii)和SCH 組(72 名精神分裂患者SCH.nii)2 類,并構(gòu)造類別標(biāo)簽集(NC.label 和SCH.label)。(2) 全局和節(jié)點(diǎn)網(wǎng)絡(luò)分析。實(shí)驗(yàn)使用解剖學(xué)自動(dòng)標(biāo)記圖譜ALL?90,通過(guò)gretna 計(jì)算相關(guān)全局與節(jié)點(diǎn)指標(biāo),具體見(jiàn)2.1 節(jié),其中使用Network Sparsity 閾值法(閾值被定義為實(shí)際連邊數(shù)與最大可能性連邊數(shù)的比率)進(jìn)行篩選,并且與隨機(jī)網(wǎng)絡(luò)進(jìn)行對(duì)比,以確定非隨機(jī)拓?fù)浣Y(jié)構(gòu),稀疏閾值設(shè)置為[0.05, 0.1, 0.15,0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5]。(3)統(tǒng)計(jì)檢驗(yàn)構(gòu)建特征集。將年齡作為協(xié)變量,通過(guò)雙樣本t 檢驗(yàn)對(duì)比NC 組和SCH 組的差異性指標(biāo)(p<0.05,FDR 矯正)并構(gòu)建特征集(NC.train 和SCH.train)。(4)建立SVM 模型。通過(guò)PSO 優(yōu)化算法在10 折交叉驗(yàn)證下對(duì)特征集和標(biāo)簽集進(jìn)行二分類實(shí)驗(yàn)。

      4 實(shí)驗(yàn)結(jié)果

      4.1 全局指標(biāo)

      通過(guò)曲下面積(AUC)分析發(fā)現(xiàn):聚類系數(shù)(Cp)有顯著差異p = 0.038,t = 2.09,NC 組大于SCH 組;最短路徑長(zhǎng)度(Lp)有顯著差異p =0.003,t =3.01,NC 組大于SCH 組;小世界指數(shù)(σ 有)顯著差異p =0.01,t =-2.58,NC 組小于SCH 組;全局效率(Eg)有顯著差異p =0.037,t = -2.1,NC 組小于SCH 組;層級(jí)系數(shù)(ab)有顯著差異p =0.017,t =-2.42,NC 組小于SCH 組,如圖2 所示;局部效率(Eloc)、同配性(ar)、同步性(as)沒(méi)有顯著差異(P>0.05)。

      4.2 節(jié)點(diǎn)指標(biāo)

      節(jié)點(diǎn)聚類系數(shù)(NCp)有顯著差異的節(jié)點(diǎn)為57(左側(cè)中央后回)。提取57 節(jié)點(diǎn)的值進(jìn)行雙樣本t 檢驗(yàn):節(jié)點(diǎn)57 的t =3.75,NC 組大于SCH 組,如圖3 所示。

      節(jié)點(diǎn)效率(NE)有顯著差異的節(jié)點(diǎn)為4(右側(cè)背外側(cè)額上回)、74(右側(cè)豆?fàn)顨ず耍?、77(左?cè)丘腦)、78(右側(cè)丘腦)。提取4 個(gè)節(jié)點(diǎn)的值進(jìn)行雙樣本t 檢驗(yàn):節(jié)點(diǎn)4 的t =-3.42,NC 組小于SCH 組;節(jié)點(diǎn)74 的t =-3.46,NC 組小于SCH 組;節(jié)點(diǎn)77 的t =-4.23,NC 組小于SCH 組;節(jié)點(diǎn)78 的t =-4.12,NC 組小于SCH 組,如圖4 所示。

      節(jié)點(diǎn)度中心性(NDc)有顯著差異的腦區(qū)為69(左側(cè)中央旁小葉)、74(右側(cè)豆?fàn)顨ず耍ⅲ罚罚ㄗ髠?cè)丘腦)、78(右側(cè)丘腦)。提?。?個(gè)腦區(qū)的值進(jìn)行雙樣本t 檢驗(yàn):節(jié)點(diǎn)69 的t =3.28,NC 組大于SCH 組;節(jié)點(diǎn)74 的t=-3.46,NC 組小于SCH 組;節(jié)點(diǎn)77 的t = -4.18,NC組小于SCH 組;節(jié)點(diǎn)78 的t = -3.58,NC 組小于SCH組,如圖5 所示。

      節(jié)點(diǎn)局部效率(NLe)沒(méi)有顯著差異;節(jié)點(diǎn)介數(shù)中心性(NBc)沒(méi)有顯著差異;節(jié)點(diǎn)最短路徑長(zhǎng)度(NLp)沒(méi)有顯著差異。

      4.3 分類結(jié)果

      通過(guò)顯著差異的全局指標(biāo)和節(jié)點(diǎn)指標(biāo)構(gòu)建特征集,如圖6 所示。實(shí)驗(yàn)對(duì)比高斯核函數(shù)和線性核函數(shù),在k =10 的交叉驗(yàn)證下(KFlod),rbf SVM 分類準(zhǔn)確率平均為76% (C = 130. 554 9,gam = 0. 001),linearSVM 分類準(zhǔn)確率平均為75.1% (C = 130.554 9)。其中,當(dāng)K =1 時(shí),linear SVM 分類準(zhǔn)確率最高為93.3%,當(dāng)K =10 時(shí),linear SMV 和rbf SVM 對(duì)于特征集的訓(xùn)練效果普遍較差,這是由于KFold 交叉驗(yàn)證劃分訓(xùn)練樣本不均所導(dǎo)致的。實(shí)驗(yàn)同時(shí)對(duì)比全局指標(biāo)和節(jié)點(diǎn)指標(biāo),其中節(jié)點(diǎn)指標(biāo)平均準(zhǔn)確率均高于全局指標(biāo)平均準(zhǔn)確率,同時(shí)低于特征集平均準(zhǔn)確率,如圖7 所示。

      5 討論

      本文在特征選擇上屬于種子點(diǎn)分析法,區(qū)別于傳統(tǒng)FC 分析和PCA 降維等方法在特征上選擇具有顯著差異性全局屬性和節(jié)點(diǎn)屬性。通過(guò)雙樣本t 檢驗(yàn)發(fā)現(xiàn),患者腦網(wǎng)絡(luò)中聚類系數(shù)和最短路徑長(zhǎng)度低于正常被試,在節(jié)點(diǎn)中央后回(PoCG.L)的節(jié)點(diǎn)聚類系數(shù)存在顯著差異(P<0.001),吳偉斌等[9] 解釋了該結(jié)果的生理學(xué)原理。全局效率患者高于正常被試,在節(jié)點(diǎn)背外側(cè)額上回(SFGdor.R)、豆?fàn)顨ず耍ǎ校眨裕?R)、丘腦(THA)的節(jié)點(diǎn)效率存在顯著差異(P<0.001),說(shuō)明精神分裂患者腦網(wǎng)絡(luò)的集團(tuán)化程度和信息傳輸速率存在差異,特別是各種感覺(jué)的傳導(dǎo)通路均在THA 內(nèi)更換神經(jīng)元后投射到大腦皮層,解釋了精神分裂患者存在幻想癥、思維混亂、極度激動(dòng)等行為,由于小世界屬性參數(shù)受到隨機(jī)網(wǎng)絡(luò)的聚類系數(shù)與最短路徑影響,因此出現(xiàn)小世界屬性SCH 組大于NC 組的現(xiàn)象。在實(shí)驗(yàn)中為了減少個(gè)體性差異對(duì)實(shí)驗(yàn)的影響,在進(jìn)行雙樣本t 檢驗(yàn)時(shí),均以年齡為協(xié)變量對(duì)比兩組被試的全局指標(biāo)和節(jié)點(diǎn)指標(biāo),考慮到性別所帶來(lái)的影響,實(shí)驗(yàn)對(duì)比僅發(fā)現(xiàn)男性患者(58 人)和女性患者(14 人)在節(jié)點(diǎn)度中心性的上距狀裂周圍皮層(CAL.R)腦區(qū)存在顯著差異(p =0.0002,t =3.92,FDR 矯正),具有較強(qiáng)的特征解釋性。在分類方法上屬于機(jī)器學(xué)習(xí)中的監(jiān)督學(xué)習(xí),實(shí)驗(yàn)通過(guò)SVM 算法解決fMRI 數(shù)據(jù)小樣本問(wèn)題,通過(guò)核函數(shù)解決了fMRI 數(shù)據(jù)高維度特性問(wèn)題,相較于之前的實(shí)驗(yàn)結(jié)果,最優(yōu)分類準(zhǔn)確率93%也進(jìn)一步證實(shí)從復(fù)雜網(wǎng)絡(luò)研究角度構(gòu)建特征來(lái)進(jìn)一步區(qū)別患者與正常被試的可能性,為計(jì)算機(jī)輔助診斷提供了參考。實(shí)驗(yàn)中仍有未充分考慮的其他干擾因素導(dǎo)致樣本泛化能力有待提高,如KFlod 交叉驗(yàn)證中正負(fù)樣本均等劃分以及大腦結(jié)構(gòu)成像中掃描的灰質(zhì)、白質(zhì)和腦脊液的形態(tài)結(jié)構(gòu)以及被試的性別、慣用手等因素,因此對(duì)于部分劃分的樣本訓(xùn)練效果較差。下一步工作可結(jié)合多模態(tài)MRI 分析,提取不同特征值和嘗試不同算法[10] ,以盡量減少由于個(gè)體差異性以及多變的生理學(xué)因素所帶來(lái)的干擾。

      參考文獻(xiàn):

      [1] RUBINOV M,SPORNS O.Complex network measures of brainconnectivity: uses and interpretations.[J].Neuroimage,2010,52(3):1059?1069.

      [2] WANG J H,ZUO X N,SURIL G,et al. Graph theoreticalanalysis of functional brain networks:test-retest evaluation onshort- and long - term resting - state functional MRI data[J].PloS one,2011,6(7):21976.

      [3] ONIAS H, VIOL A, PALHANO?FONTES F, et al. Braincomplex network analysis by means of resting state fMRI andgraph analysis: Will it be helpful in clinical epilepsy?[J].EPILEPSY AND BEHAVIOR,2014,38(1):71?80.

      [4] 王靜,孔令茵,雷炳業(yè),等.抑郁癥的腦復(fù)雜網(wǎng)絡(luò)研究進(jìn)展[J].中國(guó)醫(yī)學(xué)物理學(xué)雜志,2020,37(6):6.

      [5] 黃慶坤,楊鵬.基于神經(jīng)網(wǎng)絡(luò)的精神分裂fMRI 數(shù)據(jù)分類[J].電子技術(shù)與軟件工程,2020(3):2.

      [6] SILVA T C,ZHAO L.Machine Learning in Complex Networks[M].Springer International Publishing,2016.

      [7] ERZSéBET Ra, LáSZLó B A. Hierarchical organization incomplex networks [ J]. Physical Review E, 2003, 67 ( 2):026112.

      [8] WANG J, WANG X, XIA M, et al. GRETNA: a graphtheoretical network analysis toolbox for imaging connectomics[J].Frontiers in Human Neuroscience,2015,9(386):386.

      [9] 吳偉斌,傅禮洪,上官文博,等.青少年精神分裂癥BOLD?fMRI 特征與認(rèn)知功能的關(guān)系[J].中國(guó)CT 和MRI 雜志,2022,20(6):4?7.

      [10] AHMAD F,AHMAD I,GUERRERO?S?NCHEZ Y.Classificationof schizophrenia?associated brain regions in resting?state fMRI[J].The European Physical Journal Plus,2023,138(1):58.

      作者簡(jiǎn)介:黃慶坤(1994—),碩士,助教,研究方向:計(jì)算機(jī)技術(shù)。

      賀政剛(1994—),碩士,助教,研究方向:應(yīng)用數(shù)學(xué)。

      楊鵬(1995—),本科,助教,研究方向:臨床醫(yī)學(xué)。

      猜你喜歡
      復(fù)雜網(wǎng)絡(luò)統(tǒng)計(jì)分析機(jī)器學(xué)習(xí)
      基于圖熵聚類的重疊社區(qū)發(fā)現(xiàn)算法
      基于機(jī)器學(xué)習(xí)的圖像特征提取技術(shù)在圖像版權(quán)保護(hù)中的應(yīng)用
      基于網(wǎng)絡(luò)搜索數(shù)據(jù)的平遙旅游客流量預(yù)測(cè)分析
      如何發(fā)揮新時(shí)期統(tǒng)計(jì)工作的作用之我見(jiàn)
      前綴字母為特征在維吾爾語(yǔ)文本情感分類中的研究
      基于復(fù)雜網(wǎng)絡(luò)理論的通用機(jī)場(chǎng)保障網(wǎng)絡(luò)研究
      企業(yè)固定資產(chǎn)分析存在問(wèn)題及對(duì)策分析
      以統(tǒng)計(jì)分析為基礎(chǔ)的房地產(chǎn)稅收優(yōu)化分析
      基于支持向量機(jī)的金融數(shù)據(jù)分析研究
      城市群復(fù)合交通網(wǎng)絡(luò)復(fù)雜性實(shí)證研究
      科技視界(2016年20期)2016-09-29 11:19:34
      休宁县| 昂仁县| 华坪县| 三原县| 曲阳县| 屯门区| 青阳县| 大埔县| 武定县| 梅州市| 镇沅| 西和县| 周至县| 穆棱市| 云龙县| 平顺县| 汝阳县| 中西区| 新竹市| 武平县| 阳山县| 道真| 泌阳县| 金沙县| 岗巴县| 大田县| 沅江市| 工布江达县| 鲁甸县| 陆川县| 乐至县| 淮阳县| 富顺县| 浠水县| 东安县| 邓州市| 龙门县| 顺平县| 东平县| 无锡市| 长海县|