Matlab逆问题 - fmri数据 - 偏相关算法

时间:2011-09-06 11:16:16

标签: matlab partial correlation inverse

我正在使用以下代码来获取部分相关矩阵(来自http://www.fmrib.ox.ac.uk/analysis/netsim/的原始代码)

ic=-inv(cov(ts1)); % raw negative inverse covariance matrix
r=(ic ./ repmat(sqrt(diag(ic)),1,Nnodes)) ./ repmat(sqrt(diag(ic))',Nnodes,1); % use diagonal to get normalised coefficients
r=r+eye(Nnodes); % remove diagonal 

我的原始矩阵(ts1)是多个体素中的大脑活动(X变量) - 体积像素3X3(Y变量)。

问题是,我有更多因变量(y -voxels)而不是自变量(x-时间过程)。 我得到以下警告 -

警告:Matrix接近单一或严重缩放。          结果可能不准确。 RCOND = 4.998365e-022。

有关如何修复代码的任何想法,以便我得到所有体素之间的部分相关性?

1 个答案:

答案 0 :(得分:1)

警告来自Matlab,在反演协方差矩阵时出现问题。

一种解决方案可能是尝试pinv() http://www.mathworks.com/help/techdoc/ref/pinv.html