Resting-state fMRI topics

请教关于DPARSFA_v2.2_Parallel Works报错问题

 我是6核电脑,运行dparsfa时调用4个CPU, 33个功能图像,再完成FunImg和T1Img转换以及去除10个时间点后,报错如下;盼老师予以解答,万分感谢!Configuration:

  Crash Decoding  : Disabled
  Default Encoding: GBK
  MATLAB Root     : C:\Program Files\MATLAB\R2012a
  MATLAB Version  : 7.14.0.739 (R2012a)
  Operating System: Microsoft Windows 7

REST 软件相关分析

老师您好,根据视频中讲的,用REST Correlation Analysis,将功能连接图与一种评分做相关分析,得到R图,请问怎么从这个R图中得到结果呢?比如怎么看哪些脑区的连接强度与评分的相关具有显著性,以及得到这些脑区的信息。盼复,谢谢。

请教一个关于做ALFF低频振幅全脑相关的问题

老师好:
    我想将低频振幅的结果和一个问卷得分进行相关,看看是不是大脑自发活动和问卷得分有相关的脑区,方法是SPM中的多元回归,回归掉性别、年龄等因素,有下面几个问题:
      1、选择mALFF-1的数据是否妥当
      2、使用被试3D结构像分割后做成的被试最优灰质Mask作为数据分析的Explicit Mask是否妥当
      3、使用这个Mask从mALFF-1文件中提取每个被试的低频振幅值,作为协变量协变掉是否需要

请教一个关于Alphasim矫正的问题

       请问一下,在rest里面进行alphasim矫正时,如果是自己算,那么FWHM那个地方是按着预处理平滑盒大小来输的吗?比如预处理是8mm,那么这个地方是不是也是8呢?
       另外在选择mask的时候,能否将自己两组单样本的并集mask作为mask呢?(因为一般来说alphasim是基于先验假设进行的一种矫正方法,这个地方的mask应该选择自己假设中的mask,但这样选并集mask有没有问题呢?)
       最后,当我们在选择第一次出现小于0.05的时候,包不包括0.05呢?还是说必须是第一次出现小于0.05的值我们才要呢?
       非常谢谢~~~

time points DPARSFa 2.2

I am running a test subject on the newest version 2.2 dparsfa. THe program gives the error at the normalization step:
Moving Normalized Files:A4671_001 OK
Generating the pictures for checking normalization: A4671_001 OK
    'Error in Normalize, time point number doesn't match: A4671_001'

DPARSF2.2运行出错

刚开始运行就报错。安装的是spm8和rest1.8

??? Error using ==> parallel_function at 598
Error in ==> DPARSFA_run>(parfor body) at 222
Index exceeds matrix dimensions.

Error in ==> DPARSFA_run at 218
    parfor i=1:AutoDataProcessParameter.SubjectNum

Error in ==> DPARSFA>pushbuttonRun_Callback at 1569
    [Error]=DPARSFA_run(handles.Cfg);

Error in ==> gui_mainfcn at 96
        feval(varargin{:});

关于alphasim校正新的思考

最近又研究了一下alphasim校正,但是不清楚自己想的对不对,所以来这里和大家探讨下。

其实个人感觉alphasim并不是那么弱,只是很多研究者可能在错误的使用它。最近通过看一些文章发现:alphasim中的平滑核大小并不是我们进行预处理时smooth平滑核的大小。这个大小是需要第三方重新进行估计的。如: 3dFWHM。所以rest中给出的平滑核大小和需要的cluster size就失去了本身的意义。

比如:预处理平滑的时候用的8 8 8,最后会出来如下这个,那么我们alphasim校正时的应该是14.3 13.8 13.4 mm mm mm,我的理解对吗?但是我没有想通为什么这个平滑核和预处理时的平滑核相差会这么大呢?而且下面这三行该如何解读?resel代表了什么?

How to convert from voxel index to MNI coordinates and vice versa.

Please see the figure below (aal.nii opened by MRIcroN).
The MNI coordinates of the current position is -66 -92 -40, and its value is 0. (-66x-92x-40= 0 on the titile).
The voxel index is 25 34 32. (X 25 Y 34 Z 32).