SPM paired t-test步骤

首先感谢大神空里流霜耐心的讲解,这篇笔记内容主要是整理他的谆谆教导,虽然他也看不到>< 

 

所有数据都要经过平滑。

 

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可以帮助更好地观察脑区。

 

原文地址:https://www.cnblogs.com/minks/p/4715329.html