非参数检验和matlab程序.pptVIP

* * 【例8】 设由实验获得Ⅰ，Ⅱ两组样本值，列表如下： Ⅰ 2.36 3.14 7.52 3.48 2.76 5.43 6.54 7.41 Ⅱ 4.38 4.25 6.54 3.28 7.21 6.54 试问两总体是否同分布？（ ） 解 采用秩和检验法检验假设 其中F1(x)，F2(y)分别为总体Ⅰ，Ⅱ的分布函数 将两组样本观测值混在一起，按从小到大排序， 并计算相应的秩，列表如下： 一、秩和检验法 * * 编号 1 2 3 4 5 6 7 8 9 10 11 12 13 14 Ⅰ 2.36 2.76 3.14 3.48 5.43 6.54 7.41 7.52 Ⅱ 3.28 4.25 4.38 6.54 6.54 7.21 秩 1 2 3 4 5 6 7 8 10 10 10 12 13 14 其中观测值6.54出现3次，序号分别为9，10，11 统计量 因而其秩应为 一、秩和检验法 * * 查秩和检验表，得 由 T1=32 T2=58 由于 ，故接受H0 。 即认为两组样本对应的总体同分布。 秩和检验表只列到m,n≤10的情形 当其大于10时，统计量 于是可用U检验法求拒绝域。 一、秩和检验法 * * matlab实现 一、秩和检验法 函数 ranksum 格式 p = ranksum(x,y,alpha) % x,y为两个总体的样本, 可以不等长,alpha为显著性水平 [p,h] = ranksum(x,y,alpha) % h为检验结果, h=0表示X与Y的总体差别不显著； h=1表示X与Y的总体差别显著 [p,h,stats] = ranksum(x,y,alpha) %stats中包括: ranksum为秩和统计量的值 说明 :P为两个总体样本X和Y为一致的显著性概率, 若P接近于0,则不一致较明显. * * 一、秩和检验法 例8的Matlab实现： 输入： x=[2.36 3.14 7.52 3.48 2.76 5.43 6.54 7.41]; y=[4.38 4.25 6.54 3.28 7.21 6.54]; [p,h,stats] = ranksum(x,y,0.05) 输出： p = 0.66 h = 0 stats = ranksum: 49.00 * * 从总体X，Y中分别取容量均为N的样本： X1，X2，…，XN和Y1，Y2，…，YN 检验假设 H0:F1(x)=F2(y) H1:F1(x)≠F2(y) 将数据配对排好，列成表。 当xiyi时，取“+”号；当xiyi时，取“-”号；当xi=yi时，取“0”，并用n+和n-分别表示“+”号与“-”号的个数。 注意：样本容量必须相同 二、符号检验法 X x1 x2 … xN Y y1 y2 … yN 符号 * * 因此，当H0为真时，n+与n-以很大的概率取n/2附近的 整数值，如果min{n+,n-}比n/2小得多，就应该拒绝H0，故 选统计量 对于n和给定α，查符号检验表 可得相应的 当 时，则拒绝H0 认为两总体分布有显著差异。 二、符号检验法 定理：当H0为真时， 其中n=n++n-. * * [例9] 研究车间播放音乐对工人生产效率的影响。该车间有 10名工厂，播放音乐前与播放音乐后各30天平均日产量（件） 如下表所示，由此能否说明音乐有助于提高生产率？ （ ） 播放音乐前后平均日产量 （件） 不放音乐x 90 80 92 84 88 87 82 85 70 79 播放音乐y 99 85 97 83 81 94 72 85 82 89 符 号 — — — + + — + 0 — — 二、符号检验法 * * 解 要检验播放音乐对工人生产效率有无影响，就是检验假设 由上表可知 n+=3 n-=6 n=9 S=3 查临界值表，得 S0.05（9）=1 由于 S S0.05（9） 故接受H0 即认为播放音乐对生产率没有显著影响 二、符号检验法 * * 二、符号检验法 matlab实现 函数：signtest 格式：[p,h,stats]=signtest(x,y,alpha) P = SIGNTE