[發明專利]一種用于多被試fMRI數據分析的分組張量方法

有效

• 申請號: 201410126455.8
• 申請日: 2014-03-30
• 公開（公告）號: CN103870710B
• 公開（公告）日: 2017-01-18
• 發明（設計）人: 林秋華;鄺利丹;龔曉峰;叢豐裕
• 申請（專利權）人: 大連理工大學
• 主分類號: G06F19/00
• 分類號: G06F19/00
• 代理公司: 大連理工大學專利中心21200
• 代理人: 李寶元,梅洪玉

1.一種用于多被試fMRI數據分析的分組張量方法，其特征是，將M個被試fMRI數據構建的一個大張量x，以被試為單位分成K個子組每個子組包含N_k個被試，2≤N_k＜M；K個子組的被試間互相關系數平均值之和在M種分組方案中為最大值，為子組k內各被試的任務相關時間過程成分及腦空間激活區成分之間的互相關系數平均值，k＝1,,K；對各子組張量分別進行CP分解：得到各子組被試間共享的時間過程成分aj(k)和腦空間激活區成分sj(k)、被試間的強度差異c_j(k)及殘差_(k)；

互相關系數計算的對象是M個被試fMRI數據的任務相關時間過程成分)及腦空間激活區成分由ICA方法分離得到；求取和的協方差陣R_a和R_s，得到不同被試p、q，p,q∈{1,2,M}時間過程成分和腦空間激活區成分之間的互相關系數和計算綜合矩陣R＝0.5|R_a|+0.5|R_s|，進行兩相關被試{p,q}尋找：以初始被試p開始，將R中第p行元素置零，在第p列元素中尋找最大值所對應的行數，記為q，則找到與被試p最為相關的新被試q；重復{p,q}尋找（N_k-1）次，第一次尋找令初始被試p=未分組被試中編號最小的，之后尋找令p=q，則構建完成含有N_k個被試的一個子組；重復子組構建K次，構建完成K個子組，則一種分組方案構建完成；

對于M個被試的情況，改變第一子組內第一次{p,q}尋找的初始被試分別為1,2,...,M，一共構建M種分組方案，然后計算每種分組方案K個子組協方差矩陣R₁,R₂,,R_K上三角元素的平均值之和α₁,α₂,,α_m，和值最大的分組方案為最終分組結果。