对于不需要回归协变量的组间统计对比,可参照TBSS维基百科公布的事例即可,即如下
假设数据包含7个control (命名为Con*)和11个患者(命名为PAT*)
代码:
cd …/stats
design_ttest2 design 7 11
randomise -i all_FA_skeletonised -o tbss -m mean_FA_skeleton_mask -d design.mat -t design.con -n 500 —T2
对于需要回归协变量的统计对比,对于二分类数据,如性别可用0和1来表示。对于连续变量数据,如年龄,受教育程度,则需要对数据进行归一化处理,即个体数据值减去所有被试均值。
设计矩阵:
- 先用无需回归协变量代码得到design.con和design.mat
- 用文本打开以上两个文件,进行相应修改。
- 如下对design.con做修改,这里我回归了3个协变量,所以NumWaves由2改为5,接下来两行修改为/NumPoints 2,我这里做的两组统计对比,所以是2;接下来的矩阵依据协变量数目补0即可。
- 如下对design.mat做修改,同样NumWaves修改为5,删除第三行,矩阵补充三列协变量即可。