REST-GCA

Submitted by zangzx on
Based on Resting State fMRI Data Analysis Toolkit (REST), we implemented GCA on MATLAB as a graphical user interface (GUI) toolkit. This toolkit, namely REST-GCA (as a component of REST), could output both the F on residual and the signed-path coefficient. REST-GCA also intergrates a programme that could transform the distribution of F on the residual to normal distribution and then permit parametric statistical inference at group level.

Please see the attachments for detailed discription and usage of REST-GCA.

Citation: Zang ZX, Yan CG, Dong ZY, Huang J, Zang YF (In Press) Granger Causality Analysis implementation on MATLAB: a graphic user interface toolkit for fMRI data processing. Journal of Neuroscience Methods. In Press.

Thank you for using REST!
by Zang Zhenxiang, 2010-05-03


The readme file and usage-guide ppt file of RESTv1.8-GCA have been uploaded. To view the changes of REST-GCA, please download the two recently uploaded files! Thank you for using REST!

老师我想请教一下关于GCA的问题，
我采用的是一个run的任务态的数据，选择了一个6mm的ROI来做voxel wise的GCA，然后我将所有被试的结果做one-sample t test，出来的结果是整个脑袋都亮了，我以为是因为smooth的原因，但是我又用没smooth的数据做，结果还是一样，请问我哪里出错了么？

 你好，
请问你是不是按照默认的参数运行的GCA？如果是，那么输出结果是回归系数图。这个图可以直接做t test。如果不是，那么你用的是基于残差的结果，是不能直接做t test的。

请问你设置的p值是多少，有没有进行校正？

非常感谢老师您的回答，我已经再重新跑了一次预处理，这一次的结果貌似合理很多，谢谢老师的帮助！
另外关于GCA我还有几个问题：
第一：是不是只能所有的被试统一的ROI，要是不同的ROI就只能每个被试分开做
第二：GCA的结果图在设定cluster size的时候，会报错，是因为我做的是基于voxel wise的GCA么，所以只能设定cluster为0？
第三：能不能直接将rest产生的基于个体的GCA的map用来做spm的单样本t检验？要是可以的话，在基于voxel-wise的分析得到的都是voxel的结果，没有办法选出cluster要怎么去解释自己的结果，因为得到的结果不能完全的集中到一定的区域，只能得到voxel的GCA的表格是么？
非常感谢老师的帮助！

你好
每个被试用不同的ROI要分开做。但我不太理解这样做的目的。
能否把报错的信息贴上来，我看一看是什么错误？
你的意思是能不能用spm做T检验？一组人可以做，但一个人的结果做不了。

就是我在用slice viwer查看GCA的组分析的结果的时候设置cluster size的时候会出现以下的报错信息，（比如我将cluster size设置为100） 我在想是不是因为分析是基于voxel wise的，所以才会出现这样的情况，以下是报错信息

??? Error using ==> iptcheckconn at 49
Function BWLABELN expected its second input argument, CONN,
to be a valid connectivity specifier.
A scalar connectivity specifier must be 1, 4, 6, 8, 18, or 26.

Error in ==> bwlabeln>parse_inputs at 97
    iptcheckconn(conn,mfilename,'CONN',2);

Error in ==> bwlabeln at 74
[A,conn] = parse_inputs(varargin{:});

Error in ==> rest_sliceviewer>ThrdOverlayCluster at 3058
 [theObjMask, theObjNum]=bwlabeln(Result,AConfig.Overlay.ClusterConnectivityCriterion); %DONG Zhang-Ye and YAN Chao-Gan 090711, make the
    Cluster Connectivity Criterion flexible. %[theObjMask, theObjN
Error in ==> rest_sliceviewer>ThresholdOverlayVolume at 3086
 AConfig.Overlay.VolumeThrd =ThrdOverlayCluster(AConfig, AConfig.Overlay.VolumeThrd);

Error in ==> rest_sliceviewer at 669
    REST_SliceViewer_Cfg.Config(theCardinal) =ThresholdOverlayVolume(REST_SliceViewer_Cfg.Config(theCardinal));
 
??? Error while evaluating uicontrol Callback

谢谢老师！

另外，关于针对个体所做的分析，我的想法是，因为我想用任务态的数据来做分析，但是每个被试可能激活的脑区会存在个体差异，所以想先通过每个个体的激活的区域定义每个个体的ROI，再用来做GCA，这样是不是就把因为激活存在的个体差异给考虑了。谢谢老师！

谢谢老师的答复，我这里选的是不同被试在同一个区域的不同的peak值，就是根据他们各自的激活的peak值来选择ROI，但是所有被试激活的区域大致相同。
另外老师我想补充一个问题：
就是我现在是在分析一个任务态的数据，两个run，两个不同的任务，但是在滤波的时候遇到问题，要是用不滤波的数据进行GCA的分析，voxel-wise的结果显示是全部都亮了。
但是要是采用静息态的滤波频段应该是不合适的，是么。我有尝试过用脑电的滤波频段1~30HZ滤波，但是那样出来的结果是啥都没有。我的任务是字体判断任务和语义判断任务，到底什么样的滤波较为合适，还是直接选择不要滤波？谢谢老师指教！

如果每个被试用个体的peak来划定ROI,那么ROI之间会存在很大的变异，这个没法进行组统计。
通常情况下，任务会选择高通滤波，选择的频点是0.008HZ，也就是每128s一周期。静息状态下一般选择带通滤波，具体选哪段，不好说，没有统一的标准。你用脑电的频段对fMRI数据滤波，这个没道理。通常情况下，fMRI数据全脑扫描高频率也不会超过1HZ，你取1HZ以上的频段，是没有信号的，自然什么结果都没有。

connectivity criterion 这个参数要填进去。卡cluster size需要这个参数。
你这里没有填这个参数才报错。

 老师好！
请教您一个问题，我做了voxel-wise的GCA，每个被试都有x2x_AR_1.nii, y2y_AR_1.nii, x2y.nii, y2x.nii四个文件，我现在想进行统计分析，请问应该选什么文件？我选了所有被试的x2y.nii文件做one-sample t-test，结果报错了。
麻烦老师了，谢谢！

老师您好！

请问可以给我发一份Read me和ppt吗（xiexinhui@gmail.com）？

本链接中的下载文件失效了

“The requested page "/sites/default/files/Usage_of_REST-GCA_1_0.ppt" could not be found.”

“The requested page "/sites/default/files/REST-GCA_readme_v1.8.doc" could not be found.”

谢谢！！

Forums
Discuss