二维离散数据的相关性检验

相关性检验

 

% x 为赋值，y为均值，N为分组数量
% mn为二维直方图统计，m，n为一维直方图统计

function [mn m n]=hist_2(x, y, N)

cout=length(x);

m=hist(x,N);
n=hist(y,N);
mn=zeros(N);

% 区间间隔
d_x=(max(x)-min(x))/N;
d_y=(max(y)-min(y))/N;

% 区间位置数
X=ceil(x/d_x-min(x)/d_x);
Y=ceil(y/d_y-min(y)/d_y);

% 找出0元素，赋值为1
for i=1:cout
   if X(i)>0
       X(i)=X(i);
   else X(i)=1;
   end
   if Y(i)>0;
       Y(i)=Y(i);
   else Y(i)=1;
   end
     
end

% 二维区间计数
for i=1:cout
    mn(X(i),Y(i))=mn(X(i),Y(i))+1;
end
    
end