SPM學(xué)習(xí)筆記

原創(chuàng)作者:xiaole
來(lái)源:https://www.douban.com/note/637152103/

在用spm進(jìn)行影像數(shù)據(jù)分析的時(shí)候敞临,如何解讀結(jié)果是困擾很多研究者的問(wèn)題金蜀,本文的主要目的就是試圖梳理spm結(jié)果表中的不同標(biāo)準(zhǔn)及其意義入桂。在實(shí)際使用的時(shí)候生音,當(dāng)然是滿足越多越嚴(yán)格的標(biāo)準(zhǔn)越好啦~

先放出spm示例結(jié)果鎮(zhèn)樓窿冯,再逐一進(jìn)行解析:


image.png

本文默認(rèn)未添加Mask枣抱,以全腦為Search Volume嫩舟;結(jié)果解讀從右向左:

0. mm mm mm

即為該結(jié)果在MNI坐標(biāo)下的x y z 值贷痪,查詢此結(jié)果對(duì)應(yīng)的腦區(qū)可通過(guò)耶魯大學(xué)的在線工具:

http://sprout022.sprout.yale.edu/mni2tal/mni2tal.html

1. peak-level

在peak-level欄目中城侧,有p-FWE-corr,q-FDR-corr咱筛,T搓幌,(Z_),p-uncorr 五欄迅箩,分別代表:
?? 該peak點(diǎn)經(jīng)過(guò)FWE校正的p值
?? 該peak點(diǎn)經(jīng)過(guò)FDR校正的p值
?? 該peak點(diǎn)的T值
?? 該peak點(diǎn)經(jīng)過(guò)Z變換后的T值
?? 該peak點(diǎn)未校正的p值

image.png

1.1 p-FWE-corr

??在結(jié)果圖中溉愁,左下角FWEp即為 voxel-wise FWE 的標(biāo)準(zhǔn),通過(guò) FWEp 標(biāo)準(zhǔn)的點(diǎn)代表其可經(jīng)過(guò)利用隨機(jī)場(chǎng)理論的多重比較校正饲趋。

1.2. p<0.05 voxel-wise FDR

??在結(jié)果圖中拐揭,左下角FDRp即為 voxel-wise FDR 的標(biāo)準(zhǔn)。這種閾值方法假設(shè)Voxel之間是獨(dú)立的篙贸,單一Voxel 滿足 voxel-wise FDR 的標(biāo)準(zhǔn)即可投队,而實(shí)際的fMRI數(shù)據(jù)經(jīng)過(guò)平滑之后,voxel之間是不獨(dú)立的爵川,因此對(duì)于一個(gè)激活的腦區(qū)來(lái)說(shuō)敷鸦,用這種voxel-wise FDR校驗(yàn)得到的錯(cuò)誤發(fā)現(xiàn)率就很高。

1.3 T

??無(wú)需解釋

1.4 (Z_)

??無(wú)需解釋

1.5 p-uncorr

??此處即為不校正的結(jié)果寝贡,一般p<0.001 uncorrected 算是能接受的最低閾值了扒披,低于此閾值的結(jié)果只能進(jìn)行探索性分析。

2. cluster-level

在cluster-level欄目中圃泡,有p-FWE-corr碟案,q-FDR-corr,k-E颇蜡,p-uncorr 四欄价说,分別代表:
??該Cluster經(jīng)過(guò)FWE校正的p值,
??該Cluster經(jīng)過(guò)FDR校正的p值风秤,
??該Cluster的大小鳖目,
??該Cluster未校正的p值,

2.1. p-FWE-corr

??在結(jié)果圖中缤弦,左下角FWEc即為 cluster-wise FWE 的標(biāo)準(zhǔn)领迈,通過(guò) FWEc 標(biāo)準(zhǔn)的Cluster代表其可經(jīng)過(guò)利用隨機(jī)場(chǎng)理論的多重比較校正。

2.2. q-FDR-corr

??在結(jié)果圖中,左下角FDRc即為 cluster-wise FDR 的標(biāo)準(zhǔn)狸捅,意思即在當(dāng)前發(fā)現(xiàn)中衷蜓,Cluster Size大于此值才能考慮顯著。若結(jié)果為FDPc: Inf 尘喝,則代表沒(méi)有結(jié)果可以通過(guò) cluster-wise FDR磁浇。默認(rèn)算法是基于高斯隨機(jī)場(chǎng)理論,認(rèn)為voxel之間不是獨(dú)立的瞧省,超過(guò)T threshold的voxel 構(gòu)成了一個(gè)個(gè)cluster扯夭,然后對(duì)這些cluster進(jìn)行FDR校驗(yàn),這樣獨(dú)立單元(原來(lái)是voxel總個(gè)數(shù)鞍匾,而現(xiàn)在就是cluster的個(gè)數(shù))就大大減小了。這種校驗(yàn)的錯(cuò)誤發(fā)現(xiàn)率是小于0.05的骑科,不論T threshold 設(shè)為多大橡淑。

在SPM12中,默認(rèn) FDR算法為 topological FDR咆爽,我目前的理解是 topological 算法既可以進(jìn)行 voxel-wise FDR 亦可以進(jìn)行 cluster-wise FDR 的校正梁棠。若要改為傳統(tǒng)算法,則需要修改 spm_defaults 這個(gè)文件斗埂,其中 defaults.stats.topoFDR 這個(gè)參數(shù)默認(rèn)為1符糊,改成0即切換成了傳統(tǒng)算法(voxel-wise,但我不清楚此處的voxel-wise 和我們實(shí)際用到的2標(biāo)準(zhǔn)的區(qū)別)

2.3 k-E

無(wú)需解釋

2.4 p-uncorr

無(wú)需解釋

3. set-level

??在set-level欄目中呛凶,有p和c兩欄男娄,p代表當(dāng)前結(jié)果不是隨機(jī)出現(xiàn)的概率,c代表當(dāng)前結(jié)果中共有多少個(gè)cluster漾稀。

個(gè)人認(rèn)為set-level可參考性不大模闲,因?yàn)槿绻轻槍?duì)全腦進(jìn)行的分析,這個(gè)p值只是告訴你結(jié)果不是隨機(jī)出現(xiàn)的崭捍,但并不代表有意義尸折,比如:結(jié)果全在腦室部分。

  1. 表尾信息

大都是一些基本的參數(shù)殷蛇,不再一一敘述实夹,如有不明白的地方可以留言。

以上信息為個(gè)人理解粒梦,如有錯(cuò)誤還請(qǐng)指正亮航。在具體進(jìn)行數(shù)據(jù)分析時(shí),雖然選擇哪一個(gè)標(biāo)準(zhǔn)都可以谍倦,但我認(rèn)為peak-level 和 cluster-level的結(jié)果均很重要塞赂,需要參考。如果結(jié)果實(shí)在不甚顯著昼蛀,可以利用ROI的方法縮小Search Volume的體積宴猾,這在很多分析中不失為一個(gè)好辦法圆存。

參考鏈接:

  1. http://blog.sciencenet.cn/blog-479412-572049.html

  2. https://www.jiscmail.ac.uk/cgi-bin/webadmin?A2=spm;d940a597.1201

  3. http://www.fil.ion.ucl.ac.uk/spm/doc/manual/results.htm

  4. http://www.scholarpedia.org/article/Statistical_parametric_mapping_(SPM)

? 本文版權(quán)歸 xiaole 所有,任何形式轉(zhuǎn)載請(qǐng)聯(lián)系作者仇哆。

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末沦辙,一起剝皮案震驚了整個(gè)濱河市,隨后出現(xiàn)的幾起案子讹剔,更是在濱河造成了極大的恐慌油讯,老刑警劉巖,帶你破解...
    沈念sama閱讀 216,544評(píng)論 6 501
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件延欠,死亡現(xiàn)場(chǎng)離奇詭異陌兑,居然都是意外死亡,警方通過(guò)查閱死者的電腦和手機(jī)由捎,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 92,430評(píng)論 3 392
  • 文/潘曉璐 我一進(jìn)店門兔综,熙熙樓的掌柜王于貴愁眉苦臉地迎上來(lái),“玉大人狞玛,你說(shuō)我怎么就攤上這事软驰。” “怎么了锭亏?”我有些...
    開封第一講書人閱讀 162,764評(píng)論 0 353
  • 文/不壞的土叔 我叫張陵,是天一觀的道長(zhǎng)硬鞍。 經(jīng)常有香客問(wèn)我慧瘤,道長(zhǎng),這世上最難降的妖魔是什么膳凝? 我笑而不...
    開封第一講書人閱讀 58,193評(píng)論 1 292
  • 正文 為了忘掉前任碑隆,我火速辦了婚禮,結(jié)果婚禮上蹬音,老公的妹妹穿的比我還像新娘上煤。我一直安慰自己,他們只是感情好著淆,可當(dāng)我...
    茶點(diǎn)故事閱讀 67,216評(píng)論 6 388
  • 文/花漫 我一把揭開白布劫狠。 她就那樣靜靜地躺著,像睡著了一般永部。 火紅的嫁衣襯著肌膚如雪独泞。 梳的紋絲不亂的頭發(fā)上,一...
    開封第一講書人閱讀 51,182評(píng)論 1 299
  • 那天苔埋,我揣著相機(jī)與錄音懦砂,去河邊找鬼。 笑死,一個(gè)胖子當(dāng)著我的面吹牛荞膘,可吹牛的內(nèi)容都是我干的罚随。 我是一名探鬼主播,決...
    沈念sama閱讀 40,063評(píng)論 3 418
  • 文/蒼蘭香墨 我猛地睜開眼羽资,長(zhǎng)吁一口氣:“原來(lái)是場(chǎng)噩夢(mèng)啊……” “哼淘菩!你這毒婦竟也來(lái)了?” 一聲冷哼從身側(cè)響起屠升,我...
    開封第一講書人閱讀 38,917評(píng)論 0 274
  • 序言:老撾萬(wàn)榮一對(duì)情侶失蹤潮改,失蹤者是張志新(化名)和其女友劉穎,沒(méi)想到半個(gè)月后腹暖,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體汇在,經(jīng)...
    沈念sama閱讀 45,329評(píng)論 1 310
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 37,543評(píng)論 2 332
  • 正文 我和宋清朗相戀三年脏答,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了趾疚。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點(diǎn)故事閱讀 39,722評(píng)論 1 348
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡以蕴,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出辛孵,到底是詐尸還是另有隱情丛肮,我是刑警寧澤,帶...
    沈念sama閱讀 35,425評(píng)論 5 343
  • 正文 年R本政府宣布魄缚,位于F島的核電站宝与,受9級(jí)特大地震影響,放射性物質(zhì)發(fā)生泄漏冶匹。R本人自食惡果不足惜习劫,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,019評(píng)論 3 326
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望嚼隘。 院中可真熱鬧诽里,春花似錦、人聲如沸飞蛹。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,671評(píng)論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽(yáng)卧檐。三九已至墓懂,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間霉囚,已是汗流浹背捕仔。 一陣腳步聲響...
    開封第一講書人閱讀 32,825評(píng)論 1 269
  • 我被黑心中介騙來(lái)泰國(guó)打工, 沒(méi)想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留,地道東北人榜跌。 一個(gè)月前我還...
    沈念sama閱讀 47,729評(píng)論 2 368
  • 正文 我出身青樓闪唆,卻偏偏與公主長(zhǎng)得像,于是被迫代替她去往敵國(guó)和親斜做。 傳聞我的和親對(duì)象是個(gè)殘疾皇子苞氮,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 44,614評(píng)論 2 353

推薦閱讀更多精彩內(nèi)容