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 변수).

문제는 독립 변수 (x - time course)보다 종속 변수 (y - voxels)가 더 많다는 것입니다. 다음과 같은 경고 메시지가 표시됩니다.

경고 : 매트릭스가 단수이거나 축척이 잘못되었습니다. 결과가 정확하지 않을 수 있습니다. RCOND = 4.998365e-022.

모든 보셀 간의 부분 상관 관계를 얻을 수 있도록 코드를 수정하는 방법에 대한 의견이 있으십니까?

Answer 1

경고는 Matlab에서 공분산 행렬을 반전하는 데 문제가 있음을 나타냅니다.

하나의 해결책은 pinv() http://www.mathworks.com/help/techdoc/ref/pinv.html