一、背景
通常我们会根据一个分区模板将全脑划分为N个脑区,然后提取每个脑区的平均时间序列,计算任意两个脑区时间序列的相关系数,这样就得到了一个N*N
的功能连接矩阵。这N个脑区可能属于M个网络(比如Yeo的七网络),根据N*N
的功能连接矩阵和网络标签,我们可以得到一个M*M
的网络功能连接矩阵,这个矩阵里的对角线元素表示一个网络内部的平均功能连接,非对角线元素表示任意两个网络之间的平均功能连接。这个Matlab小程序的目的就是根据功能连接矩阵和每个脑区的网络标签,计算网络功能连接矩阵。
二、程序的使用和输入输出
-
在附件的压缩文件中,包含两个名为Network_FC_GUI.m和Network_FC_Core.m的文件,前者用来生成一个图形界面,后者用来进行实际的计算过程。还包含一个名为ExampleData的文件夹,该文件夹下有一个名为FC的子文件夹,存放着每个被试的功能连接矩阵;还有一个名为network_index.txt的文本文件,文件内容为每个脑区所属的网络标签(用不同的数字来表示不同的网络)。功能连接矩阵数据来源于NBS工具包。
-
整个程序的分析过程是:读入每个被试的FC矩阵(矩阵中的每个元素表示任意两个脑区的皮尔逊相关系数)、对FC进行Fisher’s Z变换、根据网络标签计算每个网络内部的平均功能连接和任意两个网络之间的平均功能连接。
-
在Matlab中,将当前路径切换到GUI程序所在目录(或者把GUI程序目录添加到搜索路径中),在命令行窗口中输入Network_FC_GUI,会弹出如下图所示的界面,分别选择所有被试的FC文件、网络标签文件、保存分析结果的输出目录以及输出文件的前缀(比如,输入文件是subject01.txt,前缀是NetFC,输出文件是NetFC_subject01.txt),然后选择Run Calculation开始分析。
三、其他
这个程序我是在Linux系统下测试的,不确定在Windows下是否有问题。对Matlab版本应该没有要求。这个程序仅用于学习交流,对于可能存在的错误,我不能负责。如果遇到问题,欢迎回帖进行交流。
Network_FC.zip (37.8 KB)