协变量

Submitted by annie1990923 on Sun, 10/18/2015 - 21:57

 老师：
      您好！我想请教一下关于统计过程中去除协变量的顺序和方法的问题，有三种方案分别是：

  a、做ANVOA的时候去除所有协变量因素，然后以ANVOA的结果作为mask做两两事后检验，做两两事后检验的时候就不再添加协变量因素。

  b、做ANVOA的时候去除所有协变量因素，然后以ANVOA的结果作为mask做两两事后检验，做两两事后检验的时候仍然添加协变量因素。

  c、做ANVOA的时候不去除所有协变量因素，然后以ANVOA的结果作为mask做两两事后检验，做两两事后检验的时候添加协变量因素。
     
    哪种方案是正确的呢？谢谢。

Submitted by ZangYF on Wed, 10/21/2015 - 05:13

这个问题真的是有难度。我个人认为，b是最好的。即使如此，我并不清楚，post-hoc的t-test时加协变量与第一步ANOVA时加协变量达到的效果是否相同，但不管如何，b已经是最好的了。

Submitted by annie1990923 on Wed, 10/21/2015 - 14:32

 谢谢老师，还想请教一些问题：
           1. 做功能连接，平滑选4，6有什么讲究吗？对结果影响大吗？
           2. An analysis of covariance (ANCOVA) was performed on the ALFF maps to identify brain areas with significant differences among the three groups (voxel-levelp<0.01, cluster size>918 mm3/34 voxels, corresponding to a correctedp<0.05 as determined by AlphaSim correction). These areas were then extracted as a mask. A two-sample post hoct-test of the ALFF maps within this mask was performed between each pair of the three groups (DPD vs NDPD, DPD vs HC, NDPD vs HC). The clusters showing significant differences (voxel-level p<0.01, cluster size>162 mm3/6 voxels, corresponding to a correctedp<0.05 as determined by AlphaSim correction) in ALFF between DPD and NDPD group were extracted.  我想请问1. 这段话里可以体现出做了事后两两比较的多重比较校正吗，上面只描述了AlphaSim 校正？2. 我们也按照上面方法的描述，跑了数据， 但出现了问题，事后两组比较voxel-level p<0.01, cluster size>111 mm3/3 voxels, corresponding to a corrected p<0.05 as determined by AlphaSim correction，根据mask跑出来对应corrected p<0.05 的 体素值很小，只有3， 最后报出来的有差异的脑区体素值都很小，只有几，十几的样子，感觉结果都不能用了，这是什么原因呢?跟mask有关系吗？

Submitted by ZangYF on Thu, 10/22/2015 - 07:23

1. 好像没有任何特别的规定。我以前好像看到过，说是1.2-1.5倍体素大小什么的。但是，有时为了现符合random field theory，平滑可能需要比较大。

2.1. 是的，在我以前回答其它问题时提到过。该描述只提到了multiple voxels的多重比较校正，并没有考虑到3次比较这个多重比较的校正。就我所知，目前的软件没有考虑这个问题。
2.2. 既然是在ANOVA的基础上进行的多重比较校正，得到的voxel数目肯定要比之前少。解决这个问题的办法是进行post hoc比较时，不校正。事实上，在SPSS中，post hoc的多重比较校正有许多选项，其中，bonferroni是最严格的校正，也可以选择不校正。只要解释清楚，我认为是可以的，这主要是因为目前并没有公认的解决方案，我个人认为也不是很重要，但要清楚这个问题，并且在解释结果时，对于不是非常显著的区域，解释时不要过度。

Submitted by annie1990923 on Thu, 10/22/2015 - 08:31

 谢谢老师，您说的进行post hoc 比较时，不校正是什么意思啊，是指只做multiple voxels的多重比较校正，而不做3次比较的校正吗？
               我们就是这样跑数据的，得到的Voxel 数目仍很小，那把post hoc比较时的voxel-level p 从0.01 调成 0.05，校正后还是0.05，这样会不会使voxel 数目变大？可以                这样做吗？
               如果在做anova的时候，我把voxel-level p 从0.01 调成 0.05，校正后还是0.05，可以使voxel 数目变大吗？可行吗？
               voxel-level p 值定0.01 还是0.05 有特殊的规定吗？还是可以根据结果，自己选择更合适的？
               谢谢。