Questions in AlphaSim in REST

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

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));

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

第一个问题,对于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];
            }
    }

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