首先感謝大神空里流霜耐心的講解,這篇筆記內容主要是整理他的諄諄教導,雖然他也看不到><
所有數據都要經過平滑。
Paired t-test雖然在2nd-level analysis中,但是由於本次實驗每個被試只有一張圖(自己內部已經沒有東西可以比較了),所以可以直接在2nd-level analysis中做檢驗。
P.S. 如何判定需不需要做1st-level analysis呢?簡單地說,就是每個被試只有一張圖的時候不需要再做1st-level analysis。
>Specfy 2nd-level >選目錄,design下選擇paired t-test即可,注意一下下面的masking中在explicit里選掩模,運行
>生成設計矩陣,保存在SPM.mat中(將來需要看的時候在menu的Review打開就好啦)。
>Estimate >選擇剛剛生成的SPM.mat,默認運行即可
這時候會發現在SPM.mat所在的文件夾下生成了一系列新的img文件。其中beta_0001.nii,beta_0002.nii...是設計矩陣中所有系數β;ResMS.nii就是多元回歸殘差。
>Result >Define new contrast(在蠟筆趙欣SPM8視頻中,如果在1st-level analysis下定義了Factorial design那么這里會自動生成一系列contrast) >(1 -1),t-contrast > mask:選image >inclusive(注意這里選擇mask只是作用在顯示上 ,如果沒有在前面Specify 2nd-level那里選擇masking的話依然會產生腦外激活點) >校正方式,選擇none就好了,FWE是一種比較嚴格的校正選了很可能什么都得不到...(據說修改SPM某個腳本的話可以顯示出FDR方式),選擇未校正則將來會在p處顯示unc >p = 0.001(比較寬松是0.005)>最后一項填10就可以(我也不知道什么意思...蠟筆在他的視頻里是這么說的...)
>>>>>>>>> 然后就出結果了!!!可以看到顏色越深的點表示激活程度越高。xjview可以幫助更好地觀察腦區。