我正在使用以下代码来获取部分相关矩阵(来自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。
有关如何修复代码的任何想法,以便我得到所有体素之间的部分相关性?
答案 0 :(得分:1)
警告来自Matlab,在反演协方差矩阵时出现问题。
一种解决方案可能是尝试pinv() http://www.mathworks.com/help/techdoc/ref/pinv.html