请教一个相关分析的问题

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

YAN Chao-Gan

Fri, 03/30/2012 - 18:25

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

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

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

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

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可以做吗?
谢谢!

ZHANG_RESTadmin

Thu, 04/05/2012 - 06:15

In reply to by Chen0075

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

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

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

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

十分感谢各位老师的细心解答!
我想确认一下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)。
谢谢!

张老师,你好,能回答一下么:
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.
是不是就是对该校正的描述?
Forums