求助

Submitted by sunshine on

我按照你们的那个视频教程处理了一组静息态数据，用求reho的方法做的，里面有一步是将得到的smRehoMap——图像减1，我想向各位老师请教一下在什么情况下能够得到reho的值为1？使得减1后为0.谢谢！

祝好！

smRehoMap，通常是指将计算得到的ReHo图除以自己的全脑均值，在这种情况下，一般来说直接进入二次统计就可以了，比如说两组间的比较，或者进行one-sample t-test (通常为one-sided，也就是说，只关心大于1的那些voxels)。但有些软件进行one-sample t-test时，可能不提供与1进行比较，而只能与0进行比较，在这种情况下，"smrehoMap"就只能减1，然后进行one-sample t-test。至于，“在什么情况下能够得到reho的值为1？使得减1后为0”，在没除以全脑均值以前，ReHo（准确地说是KCC）理论最大值是1，除以全脑均值后，全脑均值为1，这样，全脑的voxels中，KCC值接近于1的voxels就会非常非常多。如果进行两组间的比较，建议不再对smRehoMap减1。

我是按照视频教程Multimedia Course: Data Processing of Resting-State fMRI (静息功能磁共振脑影像数据处理视频教程)做的，上面说要减1.我用spm5做的单样本和双样本t检验，需要减1吗？十分感谢！

YAN Chao-Gan

Mon, 02/23/2009 - 06:10

Permalink

SPM5做单样本T检验的时候，只有与0比的功能，因此需要减1才能进行统计。相当于是一种权宜的办法。

而AFNI做单样本T检验，是提供与1比的功能的，因此可以不需要减1。

正如视频里面所说，在做双样本T检验的时候，是不需要减1的。因为在单样本T检验中，我们关心的是显著大于全脑均值的区域（即mReHo大于1），而双样本T检验（或配对T检验）中，我们关心的问题变成了一组显著不同于另一组（或一种条件显著不同于另一种条件）的区域。

哦，是这样啊。谢谢！

我在做时间滤波时经常出现内存溢出的现象，请问该怎么解决？

matlan显示的错误是：

 Load mask "Default".
  Build band pass filtered mask. Wait...
  Band Pass Filter working. Wait.......................
  ReConstructing 3D+time Dataset. Wait...
Exception occured. (MATLAB:nomem)
 Error using ==> zeros
Out of memory. Type HELP MEMORY for your options.
 178#line,  rest_bandpass, in "D:\soft\spm\rest\REST2007V1.2_080926\REST2007V1.2_080926\rest_bandpass.m"
 1220#line,  BandPass, in "D:\soft\spm\rest\REST2007V1.2_080926\REST2007V1.2_080926\reho_gui.m"
 <a href="error:D:\soft\spm\rest\REST2007V1.2_080926\REST2007V1.2_080926\reho_gui.m,641,0">641#line</a>,  btnBandPass_Callback, in "D:\soft\spm\rest\REST2007V1.2_080926\REST2007V1.2_080926\reho_gui.m"
 75#line,  gui_mainfcn, in "D:\Program Files\MATLAB71\toolbox\matlab\uitools\gui_mainfcn.m"
 33#line,  reho_gui, in "D:\soft\spm\rest\REST2007V1.2_080926\REST2007V1.2_080926\reho_gui.m"
>>

YAN Chao-Gan

Mon, 02/23/2009 - 11:40

Permalink

通常内存溢出可能与你的内存容量有关，也会与你内存分配的复杂程序有关。有以下几个建议：

1、扩大内存。

2、在滤波时尽量不要运行其他程序，甚至可以先重启计算机以得到干净的内存。

3、不推荐使用MATLAB7.0以下的版本。

4、打开WINXP的3GB功能，可以分配到更大的内存，详情请参考：http://www.microsoft.com/whdc/system/platform/server/PAE/PAEmem.mspx。

YAN Chao-Gan

Mon, 02/23/2009 - 14:54

Permalink

你好！对于REST来说，三个人的内存需要与一个人的内存需要是一样的，因为它在每一时刻只处理一个人的数据。

但你的电脑以前可以正常处理，现在会内存溢出，原因比较复杂，估计应该是MATLAB或是操作系统，或是硬件有了些变化。这个问题我也说不好，有的时候，也许重装系统可能会解决一些莫名其妙的问题。

sunshine

Wed, 02/25/2009 - 11:46

Permalink

根据视频教程Multimedia Course: Data Processing of Resting-State fMRI (静息功能磁共振脑影像数据处理视频教程)，在计算完reho之后要进行求平均mreho,还要平滑得到smreho.我想向各位老师请教一下问什么还要进行平滑？谢谢！

 

求平均mreho，相当于大家都用全脑均值做了一下规一化，在不同个体相比时更有可比性一些。当然，对于这一点，也有观点认为不需要做规一化，直接相比。这需要一些进一步的研究，不过目前的文献都是做了规一化的。

处理fMRI数据，基本都是要进行平滑的，这有助于减少配准精度不够带来的误差。由于平滑会影响ReHo的计算 （ReHo是计算时间序列之间的肯德尔和谐系数），因此我们选择在计算完ReHo之后再进行平滑。

Forums
Discuss