Matlab的逆问题 - fMRI数据 - 部分相关算法

Question

我用下面的代码获取部分相关矩阵（从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 - 时间过程）。 我收到以下警告 -

警告：矩阵接近单数或严重缩放。 结果可能不准确。 RCOND = 4.998365e-022。

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

Answer 1

该警告来自Matlab，它有一个反演协方差矩阵的问题。

一种解决方案是尝试PINV（） http://www.mathworks.com/help/techdoc/ref/pinv.html