在计算任意两个脑区时间序列(T*N,T为时间点数,N为脑区数)的偏相关时,当T<N时,结果出现NaN,请问出现这个结果的原因是啥?用的是matlab自带函数partialcorr
因为自由度不够,你的变量数量(脑区个数)比样本量(时间点)还多。
谢谢,那有什么解决方法么
如果脑区太多的话,我感觉可以先对要控制的脑区的时间序列做个PCA(因为不同脑区的时间序列可能高度相关,有很多冗余),这样就可以避免维度的问题(不过我没看到文献中有这种做法)。另一方面,我感觉如果是做功能连接的话,偏相关似乎并不常用,所以用偏相关是否真的更好,也许是要找更多文献支持的。
确实是的,在fMRI中计算FC的常用方式还是Pearson correlation,这里是为了做一个对比,谢谢您的解答!