原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
....
PAGE/NUMPAGES
非参数检验与matlab实现
Kolmogorov-Smirnovtest：
检验两个样本是否有相同分布
Kstest
Teststatistics：
[h,p,ksstat,cv]=kstest(x,CDF,alpha,type)
x：被测试的数据样本，以列向量输入（continuousdistributiondefinedbycumulativedistributionfunction）
CDF：被检验的样本cumulativedistributionfunction，缺省值为N(0,1)
Alpha：显著性水平，缺省时为0.05
Type：字符输入。unequal（缺省值）检验两者分布是否相同
larger检验x的CDF大于给定的CDF
smaller检验x的CDF小于给定的CDF
hh=0不拒绝原假设，即两个分布相同
h=1拒绝原假设，即两个分布不同
p：拒绝原假设的最小显著性水平
ksstat：假设为真时，满足student分布
cv：criticalvalue/cutoffvalue，determiningifksstatissignificant.
Kstest2：
[h,p,ks2stat]=kstest2(x1,x2,alpha,type)
详见ketest
Lillieforstest：
检验两个样本是否有相同分布
Teststatistics：
2-sidedgoodness-of-fittest
lillietest
[h,p,kstat,critval]=lillietest(x,alpha,distr,mctol)
各参数参见kstest，特别的，mctol为使用蒙特卡洛方法计算p值
Jarque-Beratest
检验样本是否来自均值和方差未知的正态分布
two-sidedgoodness-of-fittest
假设：x为正态分布
teststatistic：，wherenisthesamplesize,sisthesampleskewness,andkisthesamplekurtosis.
jbtest
[h,p,jbstat,critval]=jbtest(x,alpha,mctol)
各个参数意义详见lillietest。
Wilcoxon-Mann-WhitneyRanksTest
检验两个样本是否来自于同一分布。
以下functionf=Wilcoxon_Rank_Test(x,y,alpha)直接是秩和检验方法，另有functionf=Pre_Wilcoxon(x,y,alpha)是改进的秩和检验方法。区别在于样本的处理方式。
functionf=Wilcoxon_Rank_Test(x,y,alpha);
%x,yarevectors
%x,ycanbeatdifferentlength
Z=[x;y];
%z1forrankingfromsmalltolarge
%z2forplace
[z1,z2]=sort(Z,ascend);
fori=1:size(x,1)
forj=i:size(z1,1)
ifz1(j)==x(i)
X(i,1)=j;
end
end
end
%Xistherankofxinvector
m=size(x,1);
n=size(y,1);
Ex=m*(m+n+1)/2;
Varx=m*n*(m+n+1)/12;
c=sqrt(Varx)*norminv(1-alpha/2);
beta=sum(X,1);
ifabs(beta-Ex)c
disp(RejectNullHypothesis);
f=1;
elsedisp(FailtoRejectNullHypothesis);
f=0;
end
functionf=Pre_Wilcoxon(x,y,alpha);
%tofindwhetherdivisionisneeded
memo=zeros(size(x,1),size(y,1));
fori=1:size(x,1)
forj=1:size(y,1)
ifx(i,1)==y(j,1)
memo(i,j)=1;
end
end
end
%memo
文档评论（0）