ktest0.m

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

function [Dn,h]=ktest0(xi,Dnalpha)
%KTEST0 柯尔莫哥洛夫检验 
%参考文献：《应用数理统计》（叶慈南编）p143.
%柯尔莫哥洛夫检验
%xi为随机样本
%Dnalpha为柯尔莫哥洛夫检验临界值，可以通过查表得到

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

n=length(xi);
Dna=max((1:n)/n-F0(xi));
Dnb=max(F0(xi)-(0:n-1)/n);
Dn=max(Dna,Dnb);
if Dn<Dnalpha 
    h='接受H0，数据与理论分布F0(x)拟合得好';
else
    h='拒绝H0，数据与理论分布F0(x)拟合得不好';
end
function F=F0(x)
%F0(x)为检验的理论分布，它的具体形式需要根据具体情况相应地改变
if x>0 
    F=1-e-x/2;
else
    F=0;
end