请教rest做协方差分析时去协变量的问题？

Submitted by LiangHong on Thu, 11/07/2013 - 11:43

您好！
     四组被试用rest做ANCOVA，性别年龄做协变量，请问去除协变量时，是每组被试各自去除完该组对应的性别年龄协变量之后进行方差分析，还是四组被试一起去除性别年龄协变量之后进行方差分析？
    谢谢

Submitted by ZangYF on Thu, 11/07/2013 - 14:46

这是一个非常好、非常普遍存在的问题。理论上说，好像默认是一起回归掉的。但我觉得各有各的道理。
在AFNI中的多元回归分析（我忘是哪一条命令了）中，由用户自己定义。如果放在一排，表示一起回归掉，如果用两排（空的时候补0），表示分别回归掉。或许在SPM和REST中也可以自己定义，但我实在说不好，请高手回答。

Submitted by ZHANG_RESTadmin on Fri, 11/08/2013 - 00:04

接着臧老师的来说，你提的问题是个很好的问题。我觉得各自去协变量和一起去协变量的结果肯定不同。
因为，一起去协变量，默认四组来自同一个总体，所以可以把协变量在四组上接起来做一个回归子，得到一个beta值。
分别去的话，默认四组不同，所以分别作为协变量，其他地方补0，得到4个beta值。

在REST里面，是一起去的协变量。
在SPM里，对于SPM oneway ANOVA的covariate输入来看，是一次输所有人的covariate，从设计矩阵来看，也是一列全部的形式。所以就one way ANOVA来说，两者一样。


Submitted by lixiangcx on Fri, 01/03/2014 - 20:05

 谢谢各位老师，那请问，在预处理后 regress nuisance variate这一步（6个头动变量、CSF、WM、global mean），与t-test、ANOVA统计中的去除协变量是不是相同的？也是将不同组作为一个整体回归？前面做过后面统计就不需要再做协变量回归了吧？

Submitted by ZangYF on Mon, 11/11/2013 - 15:40

当时设计REST的utilities的时候，的确讨论过这个问题，但由于一般都是默认一起回归的，所以，REST也是按一起回归掉的设计的。但我不确定，如果加两个协变量是否可以达到分别去协变量的目的。我的意思是：第一个协变量中，对应第一组的数据全部为0，对应第二组的数据为实际数值；第二个协变量则在对应第一组的数据的部分为实际数值，而在对应第二组的数据的部分全部设0. 不知这样是否可行？

如果可行，需要明确说明这个做法与一起回归掉的结果是不同的，因为这种做法不常见。

Submitted by ZHANG_RESTadmin on Mon, 11/11/2013 - 18:56

臧老师，这样在REST里行不通。REST里是要求一组一组来输协变量，先放一组的人的所有协变量，再放另一组人的。
在程序中，是将这些协变量按组连接起来做的，就是说形成了一批包含所有组所有人的协变量。
如果需要改，可以在程序中修改，提供给用户要求用哪种方式。

Submitted by LiangHong on Thu, 11/14/2013 - 10:24

谢谢各位老师的回复。
年龄，性别和受教育程度作为协变量，对四组被试的ALFF图像做ANCOVA，阈值0.01，20下得到一些cluser, 将这些cluster作为ROI，提取每个ROI的ALFF值，我想提取每个ROI的去除协变量后的fitted value，这该怎么做? 请各位老师指导。

Submitted by ZHANG_RESTadmin on Fri, 11/15/2013 - 14:26

 目前没有办法通过REST实现，不过你可以提取ROI的原始值以后，再通过matlab编程实现，
具体可以用regress这个函数：

[b,bint,r] = regress(y,X)

其中y是所有人的ROI的值，X是协变量（列的形式），r就是fitted值。