关于rest的correlation analysis

Submitted by liufeng on Sun, 10/16/2011 - 11:44

管理员好！
最近我们尝试用rest的correlation analysis做量表与reho全脑的相关分析，想看下哪个脑区和量表相关比较显著。然后因为spm中的multiple regression也可以做全脑与量表算相关，又用spm算了下。 我们算的时候将年龄、性别等作为了协变量。
1、可是当我们做完发现，rest的结果和spm的结果基本上不一致，区域差别很明显。就是说基本上没有重叠的脑区，而且存在某个重叠的脑区在rest里是正相关、spm中是负相关。我们研究了好久，还是没有找出原因。只知道rest是用partialcorr做的去协变量。spm是用regress去的协变量。其实partialcorr里面也是用到了回归。我们就是想问下，为什么两个软件做的都是同一件事情可是得到的结果为什么差异这么明显呢？那么哪个结果做出来的比较真实客观呢？
2、我手头上有个问题就是两组人（病人vs正常人）的alff做双样本t检验，把年龄和教育年限作为了协变量进行了回归（就是在spm 2次统计中将他们添加为covariates）。得到差异区域后，我想提取出差异区域的值与量表（HAMD）做corr相关分析。我想问的第二个问题是，当我提取出两组人的这个区域的平均alff值的向量后，也应该先做去协变量然后再和量表做相关吧，去协变量我使用regress函数然后用corr函数得出相关值，还是直接用partialcorr得出相关系数呢？

谢谢解答！

Submitted by YAN Chao-Gan on Mon, 10/17/2011 - 06:34

multiple regression 会给出beta值。
而REST或partialcorr是要给出r值。r值具有统计意义，它与beta值具有一定的关系，统计书上会给出相应的关系。
至于使用哪个值，可能需要考虑一下。个人一般使用R值做为统计。

Submitted by liufeng on Mon, 10/17/2011 - 14:13

每一个r值都会对应了一个p值是吧，因为做corr的话，r和p是一起出来的。那么rest做r图校正的原理是直接用的一起出来的p值。还是先把r值转化到了t值，然后将t值做了校正呢？我感觉是不是可以直接用一起出来的p值做校正，直接做p<0.05+cluster size这样的校正呢？我指的一起出来的p值是这个[RHO,PVAL] = partialcorr(...)；RHO是r值；PVAL是相对应的p值。 Rest里是如何对r图做统计的呢？是我说的两种方法的哪一种呢？
谢谢！

Submitted by YAN Chao-Gan on Wed, 10/26/2011 - 03:10

每个r值都会有一个对应的p值。
REST Slice Viewer里面df按钮是设置自由度用的，你可以看到如何计算相关系数的自由度。网上google一下能找到r和p之间的转换公式。
REST->Statistical Anlyasis->Correlation Analysis会自动写入自由度信息。

Submitted by YAN Chao-Gan on Sun, 10/30/2011 - 09:29

是的，不需要转成t。