《非参数统计》与MATLAB编程第4章符号秩和检验法.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.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.32 25.76 24.23 21.35 23.43 26.97 18.36 20.75 24.07 26.43 25.41 27.22]; x2=[21.25 23.97 24.77 19.26 23.12 26.00 19.40 17.18 22.23 23.35 24.98 25.90] 符号秩和 [p,h,stats]=signrank(x1,x2) p = 0.01220703125000 h = 1 stats = signedrank: 8 在0.05显著性下，拒绝原假设，有差异。 2． 符号检验法： [p,h,stats]=signtest(x1,x2) p = 0.03857421875000 h = 1 stats = sign: 2 [h,p,ci,tstat]=ttest(x1,x2) h = 1 p= 0.00888227075133 ci = 0.37992410311173 2.101