麻烦各位老师能够帮忙解答一下问题，感谢！我的问题是
1.多中心不同机器，不同参数采集的静息态数据，比如有TR时间2s的，multiband0.8s的，0.75s的，这些数据用DPARSF预处理并计算了ALFF，后用SPM进行两组间比较，想回归掉多中心这个因素，协变量如何添加？能否设置为1,2,3····这样子？
2.多中心数据如果想用GIFT做ICA分析，该如何进行？能把多中心所有病人和健康对照一起做ICA分析吗？还是每个中心分开分析，然后找共同的网络统计分析？
- 如果想把TR作为协变量，可以直接把TR放入模型（即把实际的TR值作为一列放入设计矩阵中）。但是不同中心采集的数据往往不止TR不同（即使扫描参数完全相同，还有机器本身带来的站点效应），如果是用1/2/3表示不同站点，那么1/2/3表示不同类别，不能当作一个连续变量放入模型，要用分类变量的方式放入模型。根据我的文献了解，我个人更建议先用ComBat去除ALFF或者其他指标中的站点效应，再用去除站点效应以后的数据去做常规的统计分析。
- 我个人倾向于先用所有数据一起做ICA，在对每个被试的网络使用ComBat去除站点效应，最后用去除站点效应后的数据做常规统计分析。每个站点分开分析的问题是，寻找共同网络并非很容易的事情，容易被质疑。
1 个赞
感谢 Alex老师的回答！