《非参数统计》与MATLAB编程--第四章--符号秩和检验法.doc

______________________________________________________________________________________________________________精品资料第四章符号和检验法函数signrank格式p=signrank(x)原假设为x的中位数为0，显著性水平为0.05的双侧检验。p=signrank(x,m)原假设为x的中位数为m，显著性水平为0.05的双侧检验。p=signrank(x,m，alpha)原假设为x的中位数为m，显著性水平为alpha的双侧检验。[p,h]=signrank(...,'alpha',alpha)例：[p,h]=signrank(...,'alpha',0.01)[p,h,stats]=signrank(...,'method',‘exact’)用精确的方法[p,h]=signrank(...,'method',‘approximate’)用正态近似的方法[p,h,stats]=signrank(x,y,'alpha',0.01,'method','exact')[p,h,stats]=signrank(y1,y2,0.01,'method','approximate')所P值除以2，得到相应单侧检验的P值。§4.2x=[20.323.52219.12124.716.118.521.924.223.425];y=[1821.722.51721.224.817.214.92021.122.723.7];[p,h,stats]=signtest(x,y)p=0.3877h=0stats=sign:4length(find((x-y)>0))ans=82*(1-binocdf(7,12,0.5))ans=0.3877p=0.3877与书上算的不一样，书上算错了。符号检验接受原假设。W+个数取的值01不取任何数111212323、（1，2）424、（1，3）535、（1+4）、（2+3））646、（1，5）、（2，4）、（1，2，3）757、（1，6）、（2，5）、（3，4）、（1，2，4）868、（1，7）、（2，6）、（3，5）、（1，2，5）、（1，3，4）989、（1，8）、（2，7）、（3，6）、（4，5）、（1，2，6）、（1，3，5）、（2，3，4）101010、（1，9）、（2，8）、（3，7）、（4，6）、（1，2，7）、（1，3，6）、（2，3，5）、（4，5，1）、（1、2、3、4）111211、（1，10）、（2，9）、（3，8）、（4，7）、（5，6）、（1，2，8）、（1，3，7）、（1，4，6）、（2，3，6）、（2，4，5）、（1，2，3，5）符号秩和检验：[p,h,stats]=signrank(x,y)p=0.02685546875000h=1stats=signedrank:11P43，表4.5中：3+2+1+5=11，12*13/2-11ans=67a=11122345681012p=2*sum(a)/2^12p=0.02685546875000在显著性水平0.05下，拒绝原假设。符号秩和检验应用的条件：假设总体服从对称分布，而符号检验不需要。习题四1．x1=[22.3225.7624.2321.3523.4326.9718.3620.7524.0726.4325.4127.22];x2=[21.2523.9724.7719.2623.1226.0019.4017.1822.2323.3524.9825.90]符号秩和[p,h,stats]=signrank(x1,x2)p=0.01220703125000h=1stats=signedrank:8在0.05显著性下，拒绝原假设，有差异。2．符号检验法：[p,h,stats]=signtest(x1,x2)p=0.03857421875000h=1stats=sign:2[h,p,ci,tstat]=ttest(x1,x2)h=1p=0.00888227075133ci=0.379924103111732.10174256355494tstat=tstat:3.17229233575613df:11sd:1.354973722209693．x1=[390390450380400390350400370430]x2=[270280350300300340290320280320]x3=x2+100[phstats]=signrank(x1,x3)p=0.26562500000000h=0stats=signedrank:9或者：[p,h,stats]=signrank(x1-x2,100)p=0.26562500000000h=0stats=sig