indenpendencetest0.m

来自「统计学习的Matlab程序:包括各种假设检验,模型识别和线性回归等.」· M 代码 · 共 27 行

function [chi2,chi21,h]=indenpendencetest0(nij,alpha) 
%INDENPENDENCETEST0 独立性检验--列联表（contingency table)方法
%参考文献：《应用数理统计》（叶慈南编）p139. 
%独立性检验（列联表（contingency table)方法) 
%p141例3.5.4的数据：nij=nij=[120 30;479 371];alpha=0.01 
%p142例3.5.5的数据：nij=[32 40 59 ;133 121 106];alpha=0.05 
%调用格式：[chi2,chi21,h]=indenpendencetest0(nij,alpha) 

%By：Ji Lin
%Email: [e-mail]
%Blog: http://linji526.spaces.live.com
%Date: 2007，2008/08/18

if nargin<2 
    alpha=0.05; 
end 
[r,s]=size(nij); 
n=sum(sum(nij)); 
ni=sum(nij.'); 
nj=sum(nij); 
chi2=n*sum(sum(((nij-ni'*nj/n).^2)./(ni'*nj))); 
chi21=chi2inv(1-alpha,(r-1)*(s-1)); 
if chi2>=chi21 
    h='X与Y有关'; 
else 
    h='X与Y独立'; 
end