二維離散數據的相關性檢驗

相關性檢驗

 

% 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