《非参数统计》与MATLAB编程第5章两样本问题.docVIP

第五章 两样本问题 §5.1 a=[hygepdf(0,22,11,12) hygepdf(1,22,11,12) hygepdf(2,22,11,12) hygepdf(3,22,11,12)] a = 0 0.00001701085292 0.00093559691083 0.01403395366244 sum(a) ans = 0.01498656142619 hygecdf(3,22,11,12) ans = 0.01498656142619 §5.2.5 P72 x1=[20.6 19.9 18.6 18.9 18.8 20.2 21 20.5 19.8 19.8 19.2 20.5]; x2=[21.3 17.6 17.4 18.5 19.7 21.1 17.3 18.8 17.8 16.9 18 20.1]; format long [p,h,stats]=ranksum(x1,x2) p = 0.04317169310436 h = 1 stats = zval: 2.02204509752701 ranksum: 1.855000000000000e+002 zval的算法： n=12;m=12;N=24;me=n*(N+1)/2,de=n*m*(N+1)/12-n*m*3*(2^3-2)/(12*N*(N-1)) me = 150 de = 2.996086956521739e+002 z=(185.5-me-0.5)/sqrt(de) z = 2.02204509752701 P值的计算： p=2*(1-normcdf(z)) p = 0.04317169310436 精算算法： [p,h,stats]=ranksum(x1,x2,alpha,0.05,method,exact) p = 0.04004280818118 h = 1 stats = ranksum: 1.855000000000000e+002 习题五 1. 女职工男职工2850030650397003370031000350503325036300228003560031800372503235026900382003395030450313503080037750382002895032250367003410032900380503610030150313003480026550335503135032750392002735035700388004100025200359002990040400320503520037400355002655030450　　 x =[ 28500 31000 22800 32350 30450 38200 34100 30150 33550 27350 25200 32050 26550 30650 35050 35600 26900 31350 28950 32900 31300 31350 35700 35900 35200 30450] y =[ 39700 33250 31800 38200 30800 32250 38050 34800 32750 38800 29900 37400 33700 36300 37250 33950 37750 36700 36100 26550 39200 41000 40400 35500] ①me=median([x;y]) me = 33400 [length(find(xme)) length(find(xme));length(find(yme)) le