SPM一阶水平分析

在进行一阶水平分析时要保证你的数据是经过预处理的。

点击Specify 1st-level

SPM一阶水平分析_第1张图片

Directory是选择输出文件夹 ,units for design有seconds和scans两个选项,这里选择scans,因为scans比较简单。

interscan interval: 就是TR,Microtime resolution:是切片数量(slice number),MIcroctime oneset:参考切片(reference slice)

Data Design:双击添加数据,

conditions:有几组对比条件就添加几个,会作为GLM模型的回归因子(\beta

SPM一阶水平分析_第2张图片

SPM一阶水平分析_第3张图片

 Scans进行分析的选择fmri文件,双击conditions创建一个condition

name:起一个名字

Onesets:每个刺激块开始的时间

Durations :每个刺激块持续的时间

Mutiple conditions:可以添加经过头动矫正时获得的rp开头的txt文件。

SPM一阶水平分析_第4张图片

 以我的数据来举例说明:

一个run有5个block,每个block持续30s,TR为2s,所以每个block有15个scans,每个刺激快之后是30s的休息,第一个刺激快的开始时间切片为16,总共有165个全脑,故Onesets填16:30:165

Duration,为15

SPM一阶水平分析_第5张图片

完成之后,在menu窗口:点击Estimate,进行一阶水平分析 

 SPM一阶水平分析_第6张图片

选择上一步生成的SPM.mat文件。 点击RUN

SPM一阶水平分析_第7张图片

完成之后,在menu窗口点击Results,和上一步中一样选择,SPM.mat文件,出现如下窗口,

然后点击Define new contrast,定义新的对比分析

SPM一阶水平分析_第8张图片

 SPM一阶水平分析_第9张图片

 name:输入自定义对比分析的名字,比如我输入的是task-rest,即任务态和静息态的对比

type:选择默认的t检验

contrast :输入需要对比的\beta因子系数,因为我只定义了一个condition,所以只在第一个\beta因子输入1,其他头动因子为零,这样就是任务态和静息态的对比。当你有多个conditions时,可以使用不同的系数,来比较不同条件的对比。比如你有两个条件,positive,negetive,和6个头动因子,总共有8个回归因子,系数可以定义为1  -1  0 0 0 0 0 0 ,即两个条件间的对比

然后点击OK

 是否使用mask,这里选择none

SPM一阶水平分析_第10张图片

 FWE是比较严格条件,这里选择FWE,p值为0.05, voxels大小的阈值为30,避免噪声影响。

SPM一阶水平分析_第11张图片

这里的体素大小太了,以为设置的条件太严格,把FWE去掉,选择p值为0.001

可以发现体素的集合变大了

SPM一阶水平分析_第12张图片

 SPM一阶水平分析_第13张图片

 在这个页面,可以看到t值为5.36

选择overlays-->sections-->你的T1结构像,可以叠加到T1结构像上去看

SPM一阶水平分析_第14张图片

定义两个conditions的方法:

第一完成第一个seesion,之后点击replicateSession(1),得到第二个Session,需要更改第二个session中的Scans,conditions中的name,multiple regressor 注意对应相应的文件

SPM一阶水平分析_第15张图片

点击运行之后会有得到下图:左上一个为1 Back,右下为2 Back

 SPM一阶水平分析_第16张图片

 之后的操作与上面一致。

定义contrasts,会有些不同, 不需要写完所有的0,SPM会自动补齐0

这里定义的是2Back与1Back的对比

SPM一阶水平分析_第17张图片

 还可以定义1&2Back  ,这是任务态与静息态的对比

SPM一阶水平分析_第18张图片

 首先来看第一个对比:选择FWE,p: 0.05,voxels: 30

会发现只有两个很小的区域

SPM一阶水平分析_第19张图片

 再看第二个contrast:切换到第二个contrast

SPM一阶水平分析_第20张图片

 结果图:

SPM一阶水平分析_第21张图片

你可能感兴趣的:(经验分享,其他,matlab)