请教一个相关分析的问题

Submitted by Chen0075 on Wed, 03/28/2012 - 17:29

老师，您好！
参照REST的视频讲解，我先分别做了两组被试者的组内分析（within-group analysis），然后做了“并集”计算得出一个MASK，用于2-sample t-test。得到了两组间有差异的脑区。
我的问题是：如果我继续做correlation analysis（与行为量表相关），应该是把分析的范围限定在“并集MASK”内，还是限定在存在差异的脑区？那种做法会比较好？
另外，请问SPM可以直接做correlation analysis吗？

Submitted by YAN Chao-Gan on Sat, 03/31/2012 - 02:25

1. 取决于你的问题。如果你是要看差异脑区在病人组是否同疾病程度相关，应该选用组上差异显著的cluster。
2. 不太清楚。

Submitted by Chen0075 on Wed, 04/04/2012 - 20:39

但是，在已经做了组间比较后，已经得知了哪些脑区存在异常，那么你再去观察其他没有“异常的脑区”，是不是没有太大意义？！换句话说，在这种情况下，做“独立的全脑范围相关分析”，其得到的结果应该如何解释呢，这里能否给出一个“好文章”的例子！
其实，该问题在ANOVA的分析中，显示的更“彻底”，因为我们做post hoc组间分析时，往往是把范围设定在ANOVA分析后的significant脑区内；前一段时间一个专业的statistical reviewer就提出说把correlation analysis 扩张到“非异常的脑区”，很明显会降低statistical power！
另外，顺便请教一下：我做了组间比较得到4个clusters（并用REST提出其中的FC值）, 同时我有3个行为学量表，这样我做相关分析是就要做12次，请问这里的多次比较校正应该怎么做呢，要用什么软件，什么方法，能否给一个比较详细的操作步骤！谢谢！

Submitted by ZHANG_RESTadmin on Wed, 04/04/2012 - 22:50

我个人认为，在单样本t检验mask里看相关，或者是更全面地在全脑看相关，没有人非说这是错的。

一般可以不用校正，校正我一般用严格的p/N (N为比较次数，p为校正后阈值)。
ANOVA分析中的posthoc分析可以采用Sidak校正。不知能否应用在你的研究中。

Submitted by Chen0075 on Thu, 04/05/2012 - 10:02

P/N的校正是不是就是Bonferroni校正？
我的ANOVA后来是参考了你们那篇2011年Human Brain Mapping上的文章（就是ALFF在MCI、AD中的应用）：把Post hoc的范围限制在ANOVA后存活的脑区内；但是这种做法似乎reviewer觉得并不是最好的，他认为：Was a random effects model fit and contrasts estimated for determining whether pairs of groups differed (preferred method), or were two-sample t-tests for all pairwise comparisons performed? If the latter, note that contrasts from the ANOVA model should be used。
你说的Sidak校正，是否就是reviewer建议的做法（之一），REST可以做吗？
谢谢！

Submitted by ZHANG_RESTadmin on Thu, 04/05/2012 - 14:15

审稿人说的那段话我没有看太理解，好像说的是ANOVA，或者pairwise ttest。

校正是在统计完了以后进行的。

Sidak的方法你可以google一下，目前REST没有。Sidak比Bonforroni校正略宽松。它实际上已经考虑到了你的假设：就是ANOVA做出来显著的地方，一定是更加可能出现pairwise difference的。

Submitted by ZangYF on Thu, 04/05/2012 - 15:16

我也认为，审稿人的这句话太含糊，是否有上下文供进一步猜测？或者，你直接回信，请editor问一下审稿人，这句话到底是什么意思，显然不清楚。关于ANOVA之后的post hoc analysis的校正问题，似乎脑影像学软件都没提供这个功能。我个人认为，多重比较校正，不必过分关心，讲清楚就可以了，对于相对较大的p值，下结论时不要太绝对就是了。

Submitted by Chen0075 on Thu, 04/05/2012 - 19:21

十分感谢各位老师的细心解答！
我想确认一下P/N的校正是不是就是Bonferroni校正，我在你们另一篇文章（Xia Wu；2011）中看到这个表达：
To correct for the possible inflated type-I errors, multiple comparisons were corrected using Bonferroni at P = 0.05.
是不是就是对该校正的描述（为了写成英文回答reviewer）。
谢谢！

Submitted by Chen0075 on Mon, 04/09/2012 - 19:43

张老师，你好，能回答一下么：
P/N的校正是不是就是Bonferroni校正？我在你们另一篇文章（Xia Wu；2011）中看到这个表达：
To correct for the possible inflated type-I errors, multiple comparisons were corrected using Bonferroni at P = 0.05.
是不是就是对该校正的描述？