抱歉有两个不是问题的问题请教：
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呢？
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));
谢谢啊。由于以前不懂这个内容，看了两天好多资料，脑子有些糊涂，请不吝赐教啊。
谢谢啊。由于以前不懂这个内容，看了两天好多资料，脑子有些糊涂，请不吝赐教啊。
Forums
Re
第一个问题，对于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];
}
}
Re:Re
但是因为不同的迭代轮，会导致对mean,sd不同的样本估计值，我是没有想清楚这样是否会对模拟产生偏差。