目前有病人组前后两次的静息态MRI数据和健康对照组的静息态数据，请问怎样进行统计呢？

Submitted by wisepig119 on Sun, 04/13/2014 - 17:08

 我们对病人组在治疗前后均扫描了静息态数据，也采集了健康对照组的静息态数据，在对ALFF用REST统计上犯了难，目前的试了一下，我的步骤和问题如下：
1.病人前后和对照共三组数据做ANOVA，得到了一个脑区，再在这个脑区内用REST做两两比较（该脑区为mask），发现病人治疗前该脑区内大部分voxel比对照高（REST双样本T检验），而治疗后该脑区大部分voxel比治疗前更高（用了REST的配对T检验），但是治疗后的影像与对照比却没有得到阳性发现（REST双样本T检验，跟前面是同样的多重比较矫正水平），对于这一点我就觉得纳闷，为什么治疗前比对照高，治疗后更高，治疗后比对照反而没有影像学的阳性发现？
2、我也知道不是独立样本不适合做配对T检验和ANOVA，但是有什么更合适的办法吗？或者现在我的方法有没有可以被接受的可能？
请老师们指导一下了，谢谢

Submitted by 精进不休 on Wed, 04/16/2014 - 10:00

你好。
我也不清楚你治疗前后结果倒转的原因，我觉得最有可能是数据处理出现问题。
介于这个检验的方法，我刚过一些治疗效果的文章，大多都是选择治疗前的被试还有另外一批治疗后的被试。显然你这个是独立样本和配对样本的结合，我认为你的检验是可以的。

Submitted by ZangYF on Wed, 04/30/2014 - 11:53

1、我觉得这种情况下做ONE-way ANOVA有一些问题：病人前后是within-subject，而与正常人比较，则是between-subject。我不清楚有没有专门的这种ANOVA。由于你的“治疗”前后比较，并没有设置对照（安慰剂或者对照治疗），这种差异，只能说是前后的差异。这种差异的解释可以有很多可能性：比如，（1）第一次扫描时正常、第二次变得不正常；（2）第一次扫描不正常、第二次扫描时变得正常；等等。但不管怎样，无法说清楚是治疗变化、还是自然变化。
2、“不是独立样本不适合做配对T检验和ANOVA”，病人前后比较，最好做paired t-test；病人与正常人的比较，则用two-sample t-test。
3、具体到你目前的设计（如果你已经提供了全面的信息的话），我建议先比较病人第一次与正常人，得到异常脑区；第二步，在这些脑区内部，比较前后差异。相对来说，病人第一次扫描与正常人的比较，这种设计问题不大；病人前后的比较，无法说清楚是治疗变化、还是自然变化。

Submitted by wisepig119 on Wed, 05/07/2014 - 15:39

 明白了，谢谢臧老师的细致解答，我正是先比较了正常人与病人治疗前的差异，再比较病人治疗后的变化，不过如果局限于病人被发现的“脑区”，治疗后并没有统计学意义的变化，恰恰病人两次扫描变化的脑区在之前并没有被发现异常，但是病人这个脑区的改变恰恰跟症状的恢复有相关，我推测是由于这种药物治疗的作用——并非去逆转异常脑区，而是在正常脑区上做文章，接下来我也在搜关于这个药物在大脑作用机制的文献来支持。

Submitted by ZangYF on Fri, 05/09/2014 - 06:49

 你提到的这此脑区，即治疗之前未显示明显异常（与正常人比较无发现明显差异）、然后治疗之后发生了明显变化，并且这个脑区的变化（治疗前后的差值）与症状恢复（治疗前后的症状评分的差值）有关。这些脑区也的确是很值得关注的。