第四章 符号和检验法 函数 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.2 x=[20.3 23.5 22 19.1 21 24.7 16.1 18.5 21.9 24.2 23.4 25]; y=[18 21.7 22.5 17 21.2 24.8 17.2 14.9 20 21.1 22.7 23.7]; [p,h,stats]=signtest(x,y) p = 0.3877 h = 0 stats = sign: 4 length(find((x-y)>0)) ans = 8 2*(1-binocdf(7,12,0.5)) ans = 0.3877 p =0.3877与书上算的不一样,书上算错了。 符号检验接受原假设。 W+ 个数 取的值 0 1 不取任何数 1 1 1 2 1 2 3 2 3、(1,2) 4 2 4、(1,3) 5 3 5、(1+4)、(2+3)) 6 4 6、(1,5)、(2,4)、(1,2,3) 7 5 7、(1,6)、(2,5)、(3,4)、(1,2,4) 8 6 8、(1,7)、(2,6)、(3,5)、(1,2,5)、(1,3,4) 9 8 9、(1,8)、(2,7)、(3,6)、(4,5)、(1,2,6)、(1,3,5)、(2,3,4) 10 10 10、(1,9)、(2,8)、(3,7)、(4,6)、(1,2,7)、(1,3,6)、(2,3,5)、(4,5,1)、(1、2、3、4) 11 12 11、(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.02685546875000 h = 1 stats = signedrank: 11 P43,表 4.5中: 3+2+1+5=11,12*13/2-11 ans = 67 a = 1 1 1 2 2 3 4 5 6 8 10 12 p=2*sum(a)/2^12 p = 0.02685546875000 在显著性水平0.05下,拒绝原假设。 符号秩和检验应用的条件:假设总体服从对称分布,而符号检验不需要。 习题四 1. x1=[22....