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