《非参数统计》与MATLAB编程第四章符号秩和检验法VIP专享

第四章 符号和检验法 函数 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....