《非参数统计》与MATLAB编程第三章符号检验法

第三章 符号检验法 §3.1 符号检验 函数 signtest 格式 p = signtest(x) 原假设为 x的中位数为 0，显著性水平为 0.05的双侧检验。 p = signtest(x,m) 原假设为 x的中位数为 m，显著性水平为 0.05的双侧检验。 p = signtest(x,m，alpha) 原假设为 x的中位数为 m，显著性水平为alpha的双侧检验。 [p,h,stats] = signtest(x,m,alpha)，当样本小于 100时，stats只会显示 sign,取 x小于假设中位数 m 的个数与大于假设中位数 m 的个数的较小值，当样本容量大于或等于 100时，stats还将显示'zval'，即正态统计量值。 45.02)1(5.0nnkpnpnpkZ k为大于中位数的个数，当 k大于 n/2时，±取负号，当k小于 n/2时，±号取正号。 [p,h] = signtest(...,'alpha', alpha) 例：[p,h] = signtest(...,'alpha', 0.01) [p,h,stats] = signtest(...,'method', ‘exact’) 用精确的方法 [p,h] = signtest(...,'method', ‘approximate’) 用正态近似的方法 [p,h,stats]=signtest(x,y,'alpha',0.01,'method','exact') [p,h,stats]=signtest(y1,y2,0.01,'method',' approximate ') 例3.1某市劳动和社会保障部门的资料说明，1998年高级技师的年收入的中位数为21700，该市某个行业有一个由50名高级技师组成的样本，这些高级技师的年收入如下: 23072 24370 20327 24296 22256 19140 25669 22404 26744 26744 23406 20438 24890 24815 24556 18472 21514 22516 25112 23480 26522 24074 18064 22590 25261 21180 26188 21625 24333 23146 18324 3598 26040 20846 20438 19474 19214 23072 26744 23443 24630 26893 26485 18138 20179 26744 23554 25706 21588 17990 21700)(:21700)(:10XMHXMHee 建立一个50行和1列的向量x，即上面那个数据。 a=[]; %把上表数据粘到 a中 b=size(a); x=reshape(a,b(1)*b(2),1); format long [p,h,stats]=signtest(x,21700) p = 0.06490864707227 h = 0 stats = sign: 18 sign值的计算： length(x(find(x>21700))) ans = 32 length(x(find(x<21700))) ans = 18 取两者较小值，18。 P值的计算：2*(1-binocdf(31,50,0.5)) 得：0.06490864707229 P值也可按：2*(binocdf(18,50,0.5)) 得：0.06490864707227 P值大于缺失的显著性水平 0.05， 接受原假设。 注意：符号检验...