Why do I get NaN from corrcoef12?
1 次查看(过去 30 天)
显示 更早的评论
Why does this give me NaN as the result?
load('Varablesuv.mat') %see attached
rho = corrcoef12(invnormcdf(u),invnormcdf(v));
Here is the code for corrcoef12:
function xy = corrcoef(x,y)
%CORRCOEF Correlation coefficients.
% CORRCOEF(X) is a matrix of correlation coefficients formed
% from array X whose each row is an observation, and each
% column is a variable.
% CORRCOEF(X,Y), where X and Y are column vectors is the same as
% CORRCOEF([X Y]).
%
% If C is the covariance matrix, C = COV(X), then CORRCOEF(X) is
% the matrix whose (i,j)'th element is
%
% C(i,j)/SQRT(C(i,i)*C(j,j)).
%
% See also COV, STD.
% J. Little 5-5-86
% Revised 6-9-88 LS 2-13-95 BJ
% Copyright (c) 1984-98 by The MathWorks, Inc.
% $Revision: 5.5 $ $Date: 1997/11/21 23:23:34 $
switch nargin
case 1
c = cov(x);
case 2
c = cov(x,y);
otherwise
error('Not enough input arguments.');
end
d = diag(c);
xy = c./sqrt(d*d');
xy = xy(1,2);
From debugging this I can see that the problem is because: c = cov(x,y); gives NaN's however I'm not clear why cov is giving NaN's?
2 个评论
Walter Roberson
2017-4-6
Are you getting invnormcdf from https://github.com/beckel/nilm-eval/blob/master/Matlab/lib/lightspeed/invnormcdf.m ?
采纳的回答
Walter Roberson
2017-4-7
Your data contains 1.0 values. The invnormcdf formula is x = erfinv(2*p-1)*sqrt(2); which for 1.0 would be erfinv(2*1-1) which is erfinv(1) which is inf. Therefore you would be taking cov() of vectors with infinite values, and that is leading to NaN.
0 个评论
更多回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Special Functions 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!