Volume空间的功能连接基本都明白了。
ROI to brain
定义一个ROI,做它与全脑的功能连接,那么取这个ROI的时间序列平均值,然后把这个平均值和全脑 体素做FC,这个ROI的值是1,其他体素的值就是相关性系数的值。这样我们就得到了全脑的R值。
统计的时候,把多个被试的全脑R值放进去,两组进行对比,每个体素的R值就会有多个被试那么多,所以经过统计,就能得到每个体素的t值。如下图:
(未完待续)
ROI与ROI的对比就更简单了。
取ROI_1的平均时间序列,取ROI_2的平均时间序列,把他们进行对比,就能得到一个R值。多个被试就有多个R值。作图的时候可以显示平均R值。
干预前后,就是把干预前多个R值和干预后多个R值进行对比,就得到了1个t值和p值,那么这个数就可以代表干预有没有显著性,是显著上升还是显著下降。
这个明白了,都可以用R语言去做置换检验计算。我推荐用置换检验,尤其是对十多个被试的小样本来说。当然,其实是t检验,置换都做,哪个有结果用哪个。
ROI to ROI是不用多重校正的,因为只有1to1,如果是3个ROI,那就要两两比较的多重校正,这都符合统计常识。
DPABISurf
皮层的计算,其实原理和Volume差不太多。不同的一点是,多了一步,就是把Volume的体素投射到重建的皮层上,所以这个时候根据角度啥的,相应的把体素值给分解,然后贴到皮层上,简单理解就是:还是体素值,但是坐标系不同。
因此计算皮层的功能连接,比如定义一个ROI,给出坐标,那么这个ROI就是在皮层坐标上,然后取ROI的平均时间序列,去做FC,后面都是一样的操作,然后计算以后,把t值透射到皮层去。
皮层的结构像用fsaverage这个皮层空间坐标,皮层的功能像用fsaverage5这个皮层空间坐标。fsaverage5和MNI,这两个空间坐标需要记住。
今天遇到dparsf的报错,有点无语。
报错信息:XXX is not defined for values of class matlab.ui.figure。
报错了3次,很烦,后来排查。
解决方法:其实就是点run的时候,不能停留在dparsf的界面,需要点一下matlab的界面,回到matlab的command window。
原理:如果不回到command window去,那么check picture的时候,没办法自动跳过,就卡bug了,后续的reho,FC分析就会停止。只有回到window窗口,check picture才能自动跳过。
发这篇文,主要是为了记录这一个报错,附带着说一下皮层分析的原理,皮层肯定是更精细的,是趋势所在,所以明白理解了volume以后,还是需要早日进入surface空间去分析。