本發(fā)明屬于神經(jīng)影像學(xué)和中藥學(xué)領(lǐng)域，涉及基于功能磁共振成像fMRI數(shù)據(jù)的大腦功能網(wǎng)絡(luò)研究領(lǐng)域，具體涉及一種檢測養(yǎng)血清腦丸改善認(rèn)知功能的腦網(wǎng)絡(luò)分析方法。

背景技術(shù)：

養(yǎng)血清腦丸是天津天士力制藥股份有限公司研制、生產(chǎn)的國家三類中藥新藥，以宋代中醫(yī)名著《太平惠民和劑局方》的“四物湯”為基礎(chǔ)加味并采用最新工藝研制而成的標(biāo)本兼治的制劑，具有“養(yǎng)血活血、平肝潛陽”的作用。其主要成分是當(dāng)歸、川芎、白芍、熟地黃、鉤藤、雞血藤、夏枯草、決明子、珍珠母、延胡索、細(xì)辛。藥理研究表明，養(yǎng)血清腦丸具有改善軟腦膜微循環(huán)、增加腦血流量、緩解血管痙攣、止痛等藥理作用，現(xiàn)已成為治療頭暈、頭痛、腦供血不足、高血壓、眩暈、失眠等疾病的中藥特效藥。更有一些研究表明，養(yǎng)血清腦丸可以改善由局部缺血或精神疾病等引起的認(rèn)知功能損害，同時(shí)也能夠改善亞健康人群的心理狀態(tài)。目前，雖然對養(yǎng)血清腦丸的藥理作用已有一定認(rèn)識，但養(yǎng)血清腦丸改善亞健康人群認(rèn)知功能的潛在神經(jīng)生理機(jī)制仍是未知的。

功能磁共振成像(fMRI)技術(shù)是一種基于磁共振成像的、非介入的、無創(chuàng)傷且非常有效的腦功能成像技術(shù)，其原理是通過磁共振造影量測神經(jīng)元活動引起的血液動力的改變。由于它具有可靠的理論基礎(chǔ)、良好的時(shí)空分辨率等諸多特性，近年來，從神經(jīng)影像中發(fā)現(xiàn)對中藥敏感的生物標(biāo)記物和功能連接特性已成為一個(gè)新的研究熱點(diǎn)。

在fMRI數(shù)據(jù)分析方法上，腦功能網(wǎng)絡(luò)分析已經(jīng)發(fā)展成為fMRI數(shù)據(jù)分析領(lǐng)域的一個(gè)重要的分支，它能夠在全腦水平上進(jìn)行功能連接分析并獲取腦功能網(wǎng)絡(luò)的拓?fù)鋵傩?。除腦功能網(wǎng)絡(luò)分析外，目前還存在有基于種子點(diǎn)的功能連接分析、局部一致性分析、獨(dú)立成分分析等功能連接方法，但是這些方法都不能實(shí)現(xiàn)對網(wǎng)絡(luò)的拓?fù)鋵傩赃M(jìn)行分析。相比較而言，腦網(wǎng)絡(luò)分析方法能夠繪制同一時(shí)間內(nèi)的所有大腦單元間的連接模式，并能夠用于探討該連接模式干預(yù)前后的變化。最常用的腦網(wǎng)絡(luò)分析方法是傳統(tǒng)腦網(wǎng)絡(luò)分析，它能有效地在全腦范圍內(nèi)評估節(jié)點(diǎn)間連接的組間差異。除此之外，基于數(shù)據(jù)挖掘技術(shù)對腦網(wǎng)絡(luò)進(jìn)行分析并從中發(fā)現(xiàn)規(guī)律，也已成為新的研究趨勢。頻繁子圖挖掘作為一種分析方法已逐漸開始被應(yīng)用于腦影像領(lǐng)域，該方法主要用于尋找一個(gè)圖集中的共有子圖，并根據(jù)組間的子圖頻度差對數(shù)據(jù)進(jìn)行特征分類。

技術(shù)實(shí)現(xiàn)要素：

本發(fā)明要解決的技術(shù)問題是：利用傳統(tǒng)腦網(wǎng)絡(luò)分析和頻繁子圖挖掘方法分析fMRI數(shù)據(jù)，進(jìn)而得到養(yǎng)血清腦丸改善認(rèn)知功能的神經(jīng)生理機(jī)制。

為了解決上述技術(shù)問題，本發(fā)明采用的技術(shù)方案為：

一種檢測養(yǎng)血清腦丸改善認(rèn)知功能的腦網(wǎng)絡(luò)分析方法，包括以下步驟：

(1)實(shí)驗(yàn)設(shè)計(jì)：采用靜息態(tài)fMRI實(shí)驗(yàn)，以亞健康人群為研究對象；

(2)采集靜息狀態(tài)下的fMRI數(shù)據(jù)；

(3)針對采集到的fMRI數(shù)據(jù)進(jìn)行預(yù)處理；

(4)節(jié)點(diǎn)、邊的定義及閾值的選擇；

(5)選擇傳統(tǒng)腦網(wǎng)絡(luò)分析方法對fMRI數(shù)據(jù)進(jìn)行分析；

(6)選擇頻繁子圖挖掘方法對fMRI數(shù)據(jù)進(jìn)行分析；

(7)分析、比較并整合兩種方法的結(jié)果，得到養(yǎng)血清腦丸改善認(rèn)知功能的神經(jīng)機(jī)制。

本發(fā)明所述步驟(1)還包括以下步驟：

(1.1)篩選符合要求的亞健康人員組成一組；

(1.2)對入組的亞健康人員進(jìn)行第一次靜息態(tài)fMRI掃描和認(rèn)知測評，記錄第一次靜息態(tài)fMRI掃描和認(rèn)知測評結(jié)果；

(1.3)對入組亞健康人員進(jìn)行藥物干預(yù)；

(1.4)干預(yù)結(jié)束后對入組亞健康人員進(jìn)行第二次靜息態(tài)fMRI掃描和認(rèn)知測評，記錄第二次靜息態(tài)fMRI掃描和認(rèn)知測評結(jié)果。

本發(fā)明所述步驟(2)的具體方法為：使用德國西門子公司生產(chǎn)的3T全身磁共振成像系統(tǒng)，收集被試的全腦fMRI數(shù)據(jù)，采集回波平面成像序列，橫斷位從下到上隔層掃描，32層，層厚＝4毫米，層間隔＝0毫米，脈沖重復(fù)時(shí)間TR＝2500毫秒，脈沖回波時(shí)間TE＝30毫秒，采集矩陣＝64×64，翻轉(zhuǎn)角度FA＝90度，掃描視野FOV＝240mm，針對每個(gè)被試，采集212個(gè)時(shí)間點(diǎn)，掃描持續(xù)時(shí)間8分57秒。

本發(fā)明所述步驟(3)還包括以下步驟：

(3.1)所有的預(yù)處理操作都是使用統(tǒng)計(jì)參數(shù)圖SPM軟件包完成；

(3.2)為避免磁共振機(jī)器每個(gè)session啟動時(shí)的勻場效應(yīng)及受試者初入環(huán)境的不適應(yīng)可能對結(jié)果造成一定影響，刪除每個(gè)被試每個(gè)session前十幀功能圖像；

(3.3)剔除fMRI數(shù)據(jù)中頭部挪動幅度超過2.5mm或三維旋轉(zhuǎn)角度超過2.5°的被試；

(3.4)針對fMRI數(shù)據(jù)進(jìn)行層間時(shí)間校正；

(3.5)針對fMRI數(shù)據(jù)進(jìn)行頭動校正；

(3.6)針對fMRI數(shù)據(jù)進(jìn)行空間標(biāo)準(zhǔn)化處理；

(3.7)利用4mm半寬高的高斯核函數(shù)對fMRI數(shù)據(jù)進(jìn)行空間平滑處理；

(3.8)利用0.01-0.08Hz頻率區(qū)間對fMRI數(shù)據(jù)進(jìn)行時(shí)域?yàn)V波；

(3.9)針對fMRI數(shù)據(jù)進(jìn)行回歸頭動參數(shù)和腦脊液平均信號。

本發(fā)明所述步驟(4)還包括以下步驟：

(4.1)節(jié)點(diǎn)的定義：使用標(biāo)準(zhǔn)腦解剖圖譜AAL模板作為腦節(jié)點(diǎn)的定義方法，對全腦進(jìn)行區(qū)域級別的節(jié)點(diǎn)定義，共定義90個(gè)節(jié)點(diǎn)；

(4.2)邊的定義：使用皮爾遜相關(guān)系數(shù)作為網(wǎng)絡(luò)中邊的定義，對于所得到的每個(gè)腦的時(shí)間序列，進(jìn)行多元線性回歸，以去除頭動影響，之后，利用所得到的殘差，兩兩之間進(jìn)行皮爾遜相關(guān)計(jì)算，得到90×90的關(guān)聯(lián)矩陣，最后，根據(jù)預(yù)先設(shè)定的閾值，將關(guān)聯(lián)矩陣轉(zhuǎn)化為二值的鄰接矩陣；

(4.3)閾值的選擇：使用稀疏度S進(jìn)行閾值設(shè)置，稀疏度為網(wǎng)絡(luò)中實(shí)際存在的邊數(shù)與網(wǎng)絡(luò)中可能存在的最大邊數(shù)的比值。利用小世界屬性這一腦網(wǎng)絡(luò)的基本特性進(jìn)行閾值空間的設(shè)定，以滿足在保證小世界屬性的同時(shí)最大程度的去掉偽連接，根據(jù)此規(guī)則，計(jì)算得到閾值空間為S∈(5％，40％)，隨后的網(wǎng)絡(luò)分析均在此空間中進(jìn)行，且間隔設(shè)置為0.01。

本發(fā)明所述步驟(5)還包括以下步驟：

(5.1)對于每個(gè)閾值，分析三個(gè)節(jié)點(diǎn)指標(biāo)：度、節(jié)點(diǎn)效率和中間中心度；

(5.2)計(jì)算在閾值范圍內(nèi)的曲線下面積，為每個(gè)節(jié)點(diǎn)指標(biāo)提供一個(gè)總標(biāo)量；

(5.3)使用非參數(shù)置換檢驗(yàn)分析比較用藥前后兩組數(shù)據(jù)的組間差異，使用多元線性回歸去除統(tǒng)計(jì)過程中性別和年齡的影響。

本發(fā)明所述步驟(6)使用gSpan算法進(jìn)行頻繁子圖挖掘。

本發(fā)明所述的步驟(1.3)的藥物干預(yù)為：采用由天津天士力公司生產(chǎn)的養(yǎng)血清腦丸進(jìn)行藥物干預(yù)，商品規(guī)格是2.5g×15袋，批準(zhǔn)文號是國藥準(zhǔn)字Z20063808，入組亞健康人員每日3次，每次1袋，口服，連續(xù)服用12周。

本發(fā)明一種檢測養(yǎng)血清腦丸改善認(rèn)知功能的腦網(wǎng)絡(luò)分析方法具有下述優(yōu)點(diǎn)：本發(fā)明將傳統(tǒng)腦網(wǎng)絡(luò)分析和頻繁子圖挖掘兩種腦網(wǎng)絡(luò)分析方法結(jié)合起來，從全腦網(wǎng)絡(luò)和頻繁子網(wǎng)絡(luò)兩個(gè)角度探討藥物干預(yù)前后被試腦功能網(wǎng)絡(luò)的改變，進(jìn)而探尋養(yǎng)血清腦丸改善亞健康人群認(rèn)知功能的神經(jīng)生理機(jī)制，具有原理簡單、流程規(guī)范、實(shí)現(xiàn)簡便、被試內(nèi)結(jié)果穩(wěn)定性高等優(yōu)點(diǎn)。

附圖說明

圖1為本發(fā)明實(shí)施方法的基本流程示意圖；

圖2為本發(fā)明實(shí)施例中，分別使用藥物干預(yù)前后采集的fMRI數(shù)據(jù)構(gòu)建的腦功能網(wǎng)絡(luò)圖；

圖3為本發(fā)明實(shí)施例中，使用傳統(tǒng)腦網(wǎng)絡(luò)分析方法發(fā)現(xiàn)的藥物干預(yù)前后腦功能網(wǎng)絡(luò)節(jié)點(diǎn)屬性的差異圖；

圖4為本發(fā)明實(shí)施例中，使用頻繁子圖挖掘方法發(fā)現(xiàn)的藥物干預(yù)前后頻繁子網(wǎng)絡(luò)的差異圖；

圖5為本發(fā)明實(shí)施例中，全腦網(wǎng)絡(luò)和子網(wǎng)絡(luò)分析的結(jié)果整合圖。

具體實(shí)施方式

如圖1－圖5所示，本實(shí)施例檢測養(yǎng)血清腦丸改善認(rèn)知功能的腦網(wǎng)絡(luò)分析方法，包括以下步驟：

(1)實(shí)驗(yàn)設(shè)計(jì)：采用靜息態(tài)fMRI實(shí)驗(yàn)，以亞健康人群為研究對象；

(2)采集靜息狀態(tài)下的fMRI數(shù)據(jù)；

(3)針對采集到的fMRI數(shù)據(jù)進(jìn)行預(yù)處理；

(4)節(jié)點(diǎn)、邊的定義及閾值的選擇；

(5)選擇傳統(tǒng)腦網(wǎng)絡(luò)分析方法對fMRI數(shù)據(jù)進(jìn)行分析；

(6)選擇頻繁子圖挖掘方法對fMRI數(shù)據(jù)進(jìn)行分析；

(7)分析、比較并整合兩種方法的結(jié)果，得到養(yǎng)血清腦丸改善認(rèn)知功能的神經(jīng)機(jī)制。

本實(shí)施例步驟(1)的詳細(xì)步驟包括：

(1.1)篩選符合要求的亞健康人員組成一組；

(1.2)對入組的亞健康人員進(jìn)行第一次靜息態(tài)fMRI掃描和認(rèn)知測評，記錄第一次靜息態(tài)fMRI掃描和認(rèn)知測評結(jié)果；

(1.3)對入組亞健康人員進(jìn)行藥物干預(yù)；

(1.4)干預(yù)結(jié)束后對入組亞健康人員進(jìn)行第二次靜息態(tài)fMRI掃描和認(rèn)知測評，記錄第二次靜息態(tài)fMRI掃描和認(rèn)知測評結(jié)果。

本實(shí)施例中，步驟(1.1)具體是篩選符合以下五項(xiàng)要求的30名亞健康人員：一、經(jīng)中國亞健康學(xué)術(shù)研討會制定的亞健康狀態(tài)的檢測標(biāo)準(zhǔn)檢測，處于亞健康狀態(tài)，即在一年內(nèi)持續(xù)超過1個(gè)月表現(xiàn)為以下6類綜合癥狀群中任意1項(xiàng)以上者：(1)身體癥狀：疲勞乏力，頸、肩、上肢或腿麻木僵硬，頭暈，耳鳴，頭痛，咽部有異物感；(2)心理癥狀：記憶力減退，心煩意亂，孤單、孤獨(dú)，焦慮、煩躁，多夢易醒，休息欠佳，注意力不集中；(3)活力：活力不足或減退，對自己周圍事物和環(huán)境不感興趣，情緒不良；(4)適應(yīng)社會能力：工作能力下降或自覺吃力，同事關(guān)系不和諧；(5)免疫力：容易患病，如感冒或其他疾病；(6)到醫(yī)院看病：感覺自己有身體某方面的疾病，但卻無法被確診。二、男女均可，18-55歲，右利手，初中及以上文化程度，能配合完成相關(guān)研究項(xiàng)目。三、無重大軀體疾病及顱腦外傷史，且兩系三代無精神病史或神經(jīng)病史。四、符合磁共振掃描條件：體內(nèi)無金屬材料，如假牙，心臟起搏器等；無幽閉恐懼癥。五、對本研究知情同意。

本實(shí)施例中，步驟(1.3)具體是采用由天津天士力公司生產(chǎn)的養(yǎng)血清腦丸(Cerebralcare Granule)進(jìn)行藥物干預(yù)，商品規(guī)格是2.5g×15袋，批準(zhǔn)文號是國藥準(zhǔn)字Z20063808。入組志愿者每日3次，每次1袋，口服，連續(xù)服用12周。

本實(shí)施例中，步驟(1.2)和(1.4)具體是采用韋氏成人智力測驗(yàn)和韋氏記憶測驗(yàn)進(jìn)行認(rèn)知測評。

本實(shí)施例中，步驟(2)具體是使用德國西門子公司生產(chǎn)的3T全身磁共振成像系統(tǒng)，收集被試的全腦fMRI數(shù)據(jù)，在掃描過程中，要求被試閉上雙眼，保持清醒狀態(tài)，盡量不進(jìn)行任何思考。使用專門的泡沫墊固定被試頭部以避免頭動過大影響數(shù)據(jù)質(zhì)量，給每位被試佩戴專門的防噪音耳機(jī)以降低機(jī)器噪聲對被試產(chǎn)生的影響。掃描結(jié)束后，詢問每位被試是否在掃描過程中睡著。采集回波平面成像序列，橫斷位從下到上隔層掃描，32層，層厚＝4毫米，層間隔＝0毫米，脈沖重復(fù)時(shí)間TR＝2500毫秒，脈沖回波時(shí)間TE＝30毫秒，采集矩陣＝64×64，翻轉(zhuǎn)角度FA＝90度，掃描視野FOV＝240mm。針對每個(gè)被試，采集212個(gè)時(shí)間點(diǎn)，掃描持續(xù)時(shí)間8分57秒。

本實(shí)施例步驟(3)的詳細(xì)步驟包括：

(3.1)所有的預(yù)處理操作都是使用統(tǒng)計(jì)參數(shù)圖SPM軟件包完成，該軟件是由倫敦大學(xué)Frison等人在MATLAB平臺下開發(fā)的開源軟件；

(3.2)為避免磁共振機(jī)器每個(gè)session啟動時(shí)的勻場效應(yīng)及受試者初入環(huán)境的不適應(yīng)可能對結(jié)果造成一定影響，刪除每個(gè)被試每個(gè)session前十幀功能圖像；

(3.3)剔除fMRI數(shù)據(jù)中頭部挪動幅度超過2.5mm或三維旋轉(zhuǎn)角度超過2.5°的被試；

(3.4)針對fMRI數(shù)據(jù)進(jìn)行層間時(shí)間校正：在Slice Timing選項(xiàng)中進(jìn)行相應(yīng)的掃描層數(shù)、脈沖重復(fù)時(shí)間、掃描層順序、參考層等參數(shù)設(shè)置；

(3.5)針對fMRI數(shù)據(jù)進(jìn)行頭動校正：在Realign選項(xiàng)中選擇層間時(shí)間校正后的結(jié)果文件進(jìn)行頭動校正；

(3.6)針對fMRI數(shù)據(jù)進(jìn)行空間標(biāo)準(zhǔn)化處理：在Nomalise選項(xiàng)中使用EPI模板對校正后的結(jié)果文件進(jìn)行標(biāo)準(zhǔn)化處理；

(3.7)利用4mm半寬高的高斯核函數(shù)對fMRI數(shù)據(jù)進(jìn)行空間平滑處理：在Smooth選項(xiàng)中選擇標(biāo)準(zhǔn)化后的結(jié)果文件進(jìn)行平滑處理；

(3.8)利用0.01-0.08Hz頻率區(qū)間對fMRI數(shù)據(jù)進(jìn)行時(shí)域?yàn)V波；

(3.9)針對fMRI數(shù)據(jù)進(jìn)行回歸頭動參數(shù)和腦脊液平均信號。

本實(shí)施例步驟(4)的詳細(xì)步驟包括：

(4.1)節(jié)點(diǎn)的定義：使用標(biāo)準(zhǔn)腦解剖圖譜AAL模板作為腦節(jié)點(diǎn)的定義方法，對全腦進(jìn)行區(qū)域級別的節(jié)點(diǎn)定義，共定義90個(gè)節(jié)點(diǎn)；

(4.2)邊的定義：使用皮爾遜相關(guān)系數(shù)作為網(wǎng)絡(luò)中邊的定義，對于所得到的每個(gè)腦的時(shí)間序列，進(jìn)行多元線性回歸，以去除頭動影響，之后，利用所得到的殘差，兩兩之間進(jìn)行皮爾遜相關(guān)計(jì)算，得到90×90的關(guān)聯(lián)矩陣。最后，根據(jù)預(yù)先設(shè)定的閾值，將關(guān)聯(lián)矩陣轉(zhuǎn)化為二值的鄰接矩陣；

(4.3)閾值的選擇：使用稀疏度S進(jìn)行閾值設(shè)置，稀疏度為網(wǎng)絡(luò)中實(shí)際存在的邊數(shù)與網(wǎng)絡(luò)中可能存在的最大邊數(shù)的比值，利用小世界屬性這一腦網(wǎng)絡(luò)的基本特性進(jìn)行閾值空間的設(shè)定，以滿足在保證小世界屬性的同時(shí)最大程度的去掉偽連接。根據(jù)此規(guī)則，計(jì)算得到閾值空間為S∈(5％，40％)。隨后的網(wǎng)絡(luò)分析均在此空間中進(jìn)行，且間隔設(shè)置為0.01。

本實(shí)施例步驟(5)的詳細(xì)步驟包括：

(5.1)對于每個(gè)閾值，分析三個(gè)節(jié)點(diǎn)指標(biāo)：度、節(jié)點(diǎn)效率和中間中心度；

(5.2)計(jì)算在閾值范圍內(nèi)的曲線下面積，為每個(gè)節(jié)點(diǎn)指標(biāo)提供一個(gè)總標(biāo)量；

(5.3)使用非參數(shù)置換檢驗(yàn)分析比較用藥前后兩組數(shù)據(jù)的組間差異，使用多元線性回歸去除統(tǒng)計(jì)過程中性別和年齡的影響。

本實(shí)施例步驟(5.1)的詳細(xì)步驟包括：

(5.1.1)根據(jù)公式(1)計(jì)算節(jié)點(diǎn)的度：

k_i＝∑_j∈Ga_ij (1)

式(1)中，a_ij表示二值網(wǎng)絡(luò)A的元素(i,j)，度是對一個(gè)節(jié)點(diǎn)與網(wǎng)絡(luò)中其他節(jié)點(diǎn)連通性的基本度量，被定義為一個(gè)節(jié)點(diǎn)的邊數(shù)，通常被看作獲取網(wǎng)絡(luò)上流動內(nèi)容的直接程度。

(5.1.2)根據(jù)公式(2)計(jì)算節(jié)點(diǎn)的節(jié)點(diǎn)效率：

式(2)中，d_ij表示節(jié)點(diǎn)i和節(jié)點(diǎn)j間的最短路徑長度，最短路徑長度表示在二值矩陣中兩節(jié)點(diǎn)在所有可能的連接路徑中的最短連接數(shù)量，節(jié)點(diǎn)效率表示給定節(jié)點(diǎn)i和網(wǎng)絡(luò)中其他節(jié)點(diǎn)間的信息傳遞能力。

(5.1.3)根據(jù)公式(3)計(jì)算節(jié)點(diǎn)的中間中心度：

式(3)中，σ_mn是從節(jié)點(diǎn)m到節(jié)點(diǎn)n的最短路徑的數(shù)量，σ_mn(i)是從節(jié)點(diǎn)m到節(jié)點(diǎn)n的最短路徑經(jīng)過節(jié)點(diǎn)i的數(shù)量。中間中心度是節(jié)點(diǎn)在所處的圖中的中心性的度量，表現(xiàn)了一個(gè)節(jié)點(diǎn)在其他節(jié)點(diǎn)進(jìn)行信息傳遞時(shí)的重要性。

本實(shí)施例步驟(6)的詳細(xì)步驟包括：

(6.1)定義G＝(V,E,L,l)為有標(biāo)簽無向圖，V為節(jié)點(diǎn)集合，為邊的集合。e＝{u,v}表示節(jié)點(diǎn)u和v之間的一條邊。L是標(biāo)簽集，l是一個(gè)將標(biāo)簽映射到V中節(jié)點(diǎn)和E中邊的映射函數(shù)。

(6.2)定義子圖：對于兩個(gè)有標(biāo)簽無向圖G_S＝(V_S，E_S，L_S，l_S)和G＝(V,E,L,l)，如果則Gs是G的一個(gè)子圖。

(6.3)定義圖同構(gòu)：一個(gè)圖G₁＝(V₁，E₁，L_V1，L_E1，l₁)同構(gòu)于另一個(gè)圖G₂＝(V₂，E₂，L_V2，L_E2，l₂)，當(dāng)且僅當(dāng)存在一個(gè)雙射f:V₁→V₂使得：

(i)

(ii)

(iii)

雙射f指G1和G2之間的一個(gè)同構(gòu)。

(6.4)定義子圖支持度：給定一個(gè)圖集G，一個(gè)子圖g_s的頻繁度定義為：

(6.5)定義頻繁子圖：給定一個(gè)圖集G和一個(gè)支持度參數(shù)s，一個(gè)子圖g_s是一個(gè)頻繁子圖當(dāng)且僅當(dāng)g_s至少存在圖集G中S·|G|個(gè)圖中。

(6.6)使用gSpan算法進(jìn)行頻繁子圖挖掘，它的偽代碼如下：

以上所述僅是本發(fā)明的優(yōu)選實(shí)施方式，本發(fā)明的保護(hù)范圍并不僅局限于上述實(shí)施例，凡屬于本發(fā)明思路下的技術(shù)方案均屬于本發(fā)明的保護(hù)范圍。應(yīng)當(dāng)指出，對于本技術(shù)領(lǐng)域的普通技術(shù)人員來說，在不脫離本發(fā)明原理前提下的若干改進(jìn)和潤飾，這些改進(jìn)和潤飾也應(yīng)視為本發(fā)明的保護(hù)范圍。