fitpoisschi2test0.m

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

function [chi2,chi21,h]=fitpoisschi2test0(ni,alpha)
%FITPOISSCHI2TEST0 含有未知参数lambda的泊松分布拟合的chi2检验
%参考文献：《应用数理统计》（叶慈南编）p135.
%含有未知参数lambda的泊松分布拟合的chi2检验
%调用格式：[chi2,chi21,h]=fitpoisschi2test0(ni,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
k=length(ni);
n=sum(ni);
lenth=0:k-1;
lambda=lenth*ni'/n;
pj=poisspdf(0:k-2,lambda);
pj(k)=1-sum(pj);
% npj=n*pj;
chi2=sum((ni-n.*pj).^2./(n.*pj));
chi21=chi2inv(1-alpha,k-1);
if chi2>=chi21
    h='拒绝H0';
else
    h='接受H0';
end