Converting Covariance Matrix to Correlation Matrix

2014 10 13
3 个回答
更新时间：2014 10 24
6 次查看（30 天）

0 个投票

I know of MATLAB function to convert a covariance matrix to a correlation matrix but I'm a teacher and want to demonstrate matrix manipulation. I can do the conversion using for loops as below but can't figure out how to do this more elegantly (using matrix operations).
d=diag(cov).^.5;
for i=1:length(cov)
for j=1:length(cov)
corr(i,j)=cov(i,j)./(d(i)*d(j));
end
end
Thanks in advance for any help.

请先登录，再回答此问题。

回答（3 个）

Peter Perkins
Peter Perkins 2014-10-14

0 个投票

Jim, if you mean, "matrix operations" as in "something out of a math book", then you can right and left multiply by diag(d), where d is 1./sqrt(diag(C)). If you mean "MATLAB matrix operations", then you can use ./ with d*d'.
You might find it interesting to look at the guts of corrcoef, where the code takes a bit more care.
Hope this helps.

2 个评论
显示 无 隐藏 无

Jim Moser
Jim Moser 2014-10-14
Peter, yes I saw that suggested elsewhere and it puzzled me. If the covariance matrix is 3x3, then left multiplying it by the transpose of the diagonal gives a 1x3 and multiplying the diagonal by that gives a 1x1. The code I listed gives a 3x3 correlation matrix.
And yes by "matrix operations" I was hoping for something that would match a standard text treatment.
Peter Perkins
Peter Perkins 2014-10-24
Sorry, I think I meant diag(1./sqrt(diag(C))). If C is a 3x3 cov matrix, then diag(C) is a vector, and diag(1./sqrt(diag(C))) is a 3x3 diagonal matrix.

请先登录，再进行评论。

James Tursa
James Tursa 2014-10-24
编辑：James Tursa 2014-10-24

0 个投票

A method using outer product and element-wise divide:
d = sqrt(diag(cov));
corr = cov./(d*d');
Matt J
Matt J 2014-10-24

0 个投票

corr=bsxfun(@rdivide,cov, d(:));
corr=bsxfun(@rdivide, corr,d(:).');

类别

在 帮助中心 和 File Exchange 中查找有关 Correlation and Convolution 的更多信息

标签

提问:

Jim Moser
Jim Moser
2014-10-13

回答：

Matt J
Matt J
2014-10-24

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by