fMRI model specifictaion

GLM based

包括以下步骤:①明确GLM设计矩阵;②用经典或贝叶斯方法估计GLM参数;③利用对比向量检查结果,生成统计参数图(SPMs)以及后验概率图(PPMs)。

Design matrix:每行对应一个scan,每列对应一个影响变量或解释变量(回归量或解释变量)。关于GLM(一般线性模型)网络上有大量相关文章,

对于多被试分析,SPM通常分为两个阶段(level):

① 1st-level analysis:个体分析(within-subject)。通常有多少被试就要进行多少1st-level分析。使用 Specify fi rst level 和 Estimate 选项。这些选项的结果可以被用于做 case study。

② 2nd-level analysis:RFX analysis(Random effect analysis)。从被试推断被试所属主体的情况。

1st level analysis:

①Conditions:

name(1*n cell):condition name;

onsets(1*n cell):每个cell中为一保存condition开始时间的时间向量,每列对应一个condition,每行对应其出现时间(补充,心心水滴csddzh:任务刺激启动的扫描数(如1:14:70 ,代表任务从第 1个TR 开始,每14个TR 为一个周期,共70个TR));

durations:事件时长。

一般这些信息都保存在condition.mat文件中,使用的时候直接load进来,使用Multiple Conditions.

②Regressors

设计矩阵中用来和血液动力响应做卷积的一列。之后的leave-one-out实验要用。

暂时还没搞清楚是不是要用Multiple Regressors。

③Factorial design

设定该选项则SPM 可以直接生成主要效应和交互作用之间必要的test:F-test(被试个体(1st level)与被试间(2nd level)的contrast image)和T-test。

根据蜡笔的视频:被试根据三种商品颜色选择买或者不买,那么就有两个factor:color,3;choice,2。

SPM 默认第一个因素变化最慢,第二个因素变化第二慢...以此类推。

Factorial design一般在有复杂情况的时候才用,简单的时候可以不用设置。

其他因素没怎么看,因为在蜡笔的视频里都是默认。

=============================================

fMRI model estimation

模型参数可以用经典(ReML - Restricted Maximum Likelihood)或贝叶斯算法估计。

蜡笔视频中这一步直接使用了上一步生成(?)的SPM.mat之后用默认运行。

生成一SPM.mat,包含设计和估计的模型参数。

看结果:主面板Result。

=============================================

Factorial design specification(2nd level analysis)

配置设计矩阵,描述统计分析所需要的GLM,数据及必要的参数(保存在GLM中)。用于为2nd level analysis做准备。

Specify 2nd level analysis。

05-04 00:27