研究生-统计学讲义-第8讲-非参数检验与Ridit分析.pptVIP

第四节 完全随机设计多样本比较的秩和检验(H检验) 完全随机分组资料比较, H0为各总体分布相同. H0成立时, 容量分别为ni (1≤i≤k) 的k个样本来自同一总体. 定量资料编秩时, 同组相同数据取顺序秩次, 不同组相同数据取平均秩次. 分类资料编秩时, 同一等级取平均秩次. 记N=∑ni ,Ti为容量为ni样本的秩和． 当 ni 增大时, Ti构成的 H 统计量近似于自由度为k－1的卡方分布，即 ～x2，df＝k－1． 第三十页，编辑于星期六：七点 五十九分。 这种方法称完全随机分组秩和检验（Kruskal－Wallis法）, 在ni 较小时, 可查统计用表13, 用H值与H界值进行比较． 在ni 较大时，可用 H 统计量近似作 x2 检验，并在相同秩次太多时校正，即 例为比较不同批号药盒检验结果是否一致，将三个批号各四个药盒一一测定了某一标本。结果如表11-4，试问不同批号测定结果有无差别？ H0:不同批号测定值总体分布位置同；H1:不同批号的测定值总体分布位置不同或不全相同.α＝0.05。 第三十一页，编辑于星期六：七点 五十九分。 1号 1.80 1.92 2.00 2.01 2.11 n1=5 2号 2.11 2.21 2.25 2.57 n2=4 3号 2.75 3.10 3.27 3.27 n3=4 rank 1 2 3 4 5.5 T1=15.5 rank 5.5 7 8 9 T2=29.5 rank 10 11 12.5 12.5 T3=46 n1+n2+n3=N=13 H =10.449，P 值0.005；故拒绝H0，接受H1，可认为不同批号的测定值总体的分布位置不同或不全相同 第三十二页，编辑于星期六：七点 五十九分。 第三十三页，编辑于星期六：七点 五十九分。 输出结果：由Kruskal-Wallis H中的Test Statistica,b得χ2=10.449，近似P值=0.005；由Median Test中的Test Statisticb得χ2=8.976， P=0.011，故单侧检验均拒绝H0，接受H1，可认为不同批号的测定值总体的分布位置不同或不全相同 第三十四页，编辑于星期六：七点 五十九分。 例9.4 用祖传及一般针炙疗法治疗哮喘病人46例及28例判断祖传针灸疗法的疗效是否高于一般针炙疗法 疗效 一般 祖传 合计 范围 平均 一般秩和 祖传秩和 无效 5 3 8 1～8 4.5 22.5 13.5 好转 14 15 29 9～37 23 322 345 显效 5 16 21 38～58 48 240 768 痊愈 4 12 16 59～74 66.5 266 798 合计 n1＝28 n2＝46 N＝74 T1＝850.5 T2=1924.5 分类资料，单向有序且样本容量较小，不宜Ridit分析 第三十五页，编辑于星期六：七点 五十九分。 H0：两总体分布相同，H1：两总体分布不同 计算各等级合计数，秩次范围，平均秩次，秩和。如疗效为“无效”者合计8例，平均秩次为(1＋8)/2＝4.5，一般疗法组的秩和为4.5×5＝22.5 确定T＝850.5，t1＝8、t2＝29、t3＝21、t4＝16 ＝2.3305 由uC＞u0.05/2，双侧P＜0.05，以?＝0.05水准的双侧检验拒绝H0，两总体的分布不同 第三十六页，编辑于星期六：七点 五十九分。 由T1＜T2，可以认为祖传针灸疗法的疗效高于一般针炙疗法 第四节 完全随机分组秩和检验(K-W法) ⑴ 完全随机分组秩和检验，H0：各总体分布相同 ni增大时,Ti构成H统计量近似df=k－1卡方分布 定量资料同组相同数据取顺序秩,不同组相同数据取平均秩。分类资料同一等级取平均秩 第三十七页，编辑于星期六：七点 五十九分。 ni较小时查统计用表13，用H值与H界值进行比较 ni较大时用H统计量近似x2检验,相同秩多时校正 ⑵ 在各总体分布不全相同结论下多重比较，H0：第i、j个总体分布相同 df＝N－k 第三十八页，编辑于星期六：七点 五十九分。 在无相同数据时 在有相同数据时 例9.5 用三种药物杀灭钉螺，每批用200只活钉螺，计算死亡率，问三种药物杀灭钉螺的效果有无差异 第三十九页，编辑于星期六：七点 五十九分