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)是多体素-volumetric像素3X3 (Y变量)中随时间变化的大脑活动(X变量)。

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

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

有没有想过如何修复代码，这样我就可以得到所有体素之间的部分相关性？

Answer 1

该警告来自于Matlab在求协方差矩阵时出现的问题。

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