0

我正在使用以下代码来获取偏相关矩阵（来自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 - time course）。我收到以下警告-

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

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

4

1 回答 1

1

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

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

于 2011-10-23T21:33:36.203 回答