因为自己收的一批数据来自不同中心医院,而且每家医院的样本量也比较小,于是我采用了combat的数据融合方法。
处理过程是这样的:
1.先把各中心数据用DPABI进行预处理,再获得ReHo值、fALLF值。
2.生成所有被试的90%groupmasks
3.将每个被试ReHo、fALLFmap转换为矩阵,用reshape函数生成一列数值,再生成groupmasks范围内的矩阵,最后将得到特征值x被试数目的矩阵。去除为0的行。
4.将不同的被试所在的医院设为不同的batch,age、sex、头动、疾病类型设为mod,方法1经典贝叶斯、2非参数,进行combat处理
5.combat后获得到的矩阵补0,再转为nii文件。最后再放入DPABI进行统计分析。
如果有步骤不对的地方,欢迎大家交流指正。