计算两组或者两组以上共同激活脑区

如何计算两组或者两组以上共同激活脑区(mfALFF指标为例)。

想法1:把这些分组的数据放到同一个文件夹,先和数值1做图像运算,再用spm做单样本t检验。
想法2:把两组+对照组这三组数据,用单因素方差分析导入spm,输入对比矩阵:组1 vs 对照组,组2 vs 对照组,同时选中,然后选择Conjunction分析。
想法3:分别计算组1和对照组,组2和对照组激活的区域,做图像计算,(组1>0)+(组2>0)*2,得到数值1代表组1激活,数值2代表组2激活,数值3代表共同激活。
用前面三种想法都计算过,结果不一一致,请教应当用哪一种计算方法?或者有其他正确的方法。

(1)我觉得方法2和方法3结果应该是非常接近的,如果卡的阈值一样的话。
(2)方法1我不太懂和数值1做图像运算是什么意思?
(3)“激活”这个词我觉得不太恰当,因为一般任务态fMRI才会涉及到激活。你这里的“激活”应该是指的组间差异?

很奇怪,方法2和方法3的结果并不一致。方法2分析的时候,设计对比矩阵选择T检验,或者F检验,都尝试过。可能是我对conjunction analysis的理解有误。不知可否帮忙解释conjunction analysis?
方法1是基于单样本t检验和某一个常数比较,-1是表示ALFF值显著高于全脑平均水平?不过这里面没有引入对照组,所以结果和其他方法不一致也可以理解。
倒是用两个结果图做点乘,或者方法3(组1>0)+(组2>0)*2,这两个运算得到的结果是一致的。

conjunction analysis就是两个显著结果的交集,比如在你这里就是组1 > 对照组的显著区域组2 > 对照组的显著区域的交集。你可以参考Handbook of functional MRI data analysis的7.4节或者第125页:“Valid conjunction inference is obtained by thresholding each statistic image separately and then taking the voxel-level intersection of above-threshold voxels.”