关于rest的correlation analysis

Submitted by liufeng on
管理员好! 最近我们尝试用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得出相关系数呢? 谢谢解答!

YAN Chao-Gan

Sun, 10/16/2011 - 22:34

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

liufeng

Mon, 10/17/2011 - 06:13

In reply to by YAN Chao-Gan

每一个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图做统计的呢?是我说的两种方法的哪一种呢? 谢谢!
Forums