Questions in AlphaSim in REST

Submitted by frog on Mon, 10/25/2010 - 22:38

 抱歉有两个不是问题的问题请教：

1. AFNI里AlphaSim还有关于user specified activate region情况，请问REST有考虑将来的版本升级呢？
    另外请教，相应参数Zsep (separation between the mean of the noise distribution and the mean of the signal distribution)如何设定呢？

2. REST里AlphaSim关于xthr (xthr = sd*zthr + mean)计算，sd和mean都是每迭代一次，要累积前面迭代的值，我不是很明白。即如下：
    count=count+nxyz;   为什么不每次迭代count都赋值为nxyz呢？

    suma=sum(fimca)+suma;
    sumsq=sum(fimca.*fimca)+sumsq;
    mean=suma/count;
    sd = sqrt((sumsq - (suma * suma)/count) / (count-1));

谢谢啊。由于以前不懂这个内容，看了两天好多资料，脑子有些糊涂，请不吝赐教啊。

Submitted by dongzy08 on Tue, 10/26/2010 - 11:13

第一个问题，对于AFNI，目前我还不太会用，Douglas在2000年的那篇文章也没给例子，弄明白用法，REST应该会更新这部分内容。

第二个问题，是按照AFNI的源程序的写法，参见AlphaSim.c，其目的是计算累计的均值和P值，并不是针对某一次。
 if (*count < 1.0e+09)
    {
      *count += nxyz;
      for (ixyz = 0;  ixyz < nxyz;  ixyz++)
           {
               *sum += fim[ixyz];
               *sumsq += fim[ixyz] * fim[ixyz];
            }
    }

Submitted by frog on Sat, 10/30/2010 - 23:19

 我是一点都没用过AFNI，今天我下载了，看了一下这个函数代码，确实也是用的累积值，我猜想，这样可能是希望对mean,sd估计更准确吧。
但是因为不同的迭代轮，会导致对mean,sd不同的样本估计值，我是没有想清楚这样是否会对模拟产生偏差。