第八章秩转换的非参数检验(正式).ppt

* 本例n1=10, T=141.5，查附表10确定T概率 界值： 假设成立T1理论秩和=115 76 79 84 89 115 141 146 151 154 0.05 0.025 0.01 0.005 本例 p＜0.05。在α=0.05水准， p＜0.05，拒绝H0，两组 变量值分布差别有统计意义，肺癌组平均秩次（141.5/10=14.15）高于矽肺0期组平均秩次（111.5/12=9.29）。 * 若n110或n2-n110，可用正态近似法作u检验，令N=n1+n2 ，按下式计算 值。。 式中tj（j=1,2,…）为第j个相同秩的个数。 * 例4 39名吸烟工人和40名不吸烟工人的碳氧血红蛋白HbCO(%)含量见表8-6。问吸烟工人的HbCO(%)含量是否高于不吸烟工人的HbCO(%)含量？ 本例两样本资料为等级资料，现用Wilcoxon秩和检验。 2．频数表资料和等级资料的两样本比较 * 表 吸烟工人和不吸烟工人的HbCO(%)含量比较 含 量 吸烟 工人 不吸烟 工人 合计 秩范围 平均秩 秩 和 吸烟工人 不吸烟工人 (1) (2) (3) (4) (5) (6) (7)=(2)(6) (8)=(3)(6) 很低 1 2 3 1～3 2 2 4 低 8 23 31 4～34 19 152 437 中 16 11 27 35～61 48 768 528 偏高 10 4 14 62～75 68.5 685 274 高 4 0 4 76～79 77.5 310 0 合 计 39(n1) 40(n2) 79 ─ ─ 1917(T1) 1243(T2) * H0：两组观察值的总体分布位置相同 H1：两组观察值的总体分布位置不同 α=0.05 方法同前：本例n1=39,n2=40,n1+n2=79, T=1917（例数较小组） 查u界值表得单侧P0.005，按α＝0.05水准拒绝H0，接受H1，可认为吸烟工人的HbCO含量高于不吸烟工人的HbCO含量。 * H0：两总体分布相同。 Wilcoxon秩和检验的基本原理 编秩 编秩 样本观察值既不会集中于高秩，也不会集中于低秩。 * H0：两总体分布位置相同。 只需要保证经编秩后，总体分布能够转换成相同的均匀分布。 编秩转换只对分布的位置敏感，而对分布的形态不甚敏感。 从而，可以将H0（两总体分布相同）这一较严格的前提，转换为H0（两总体分布位置相同）这一较宽松的前提。 * 第三节 完全随机设计多个样本比较的Kruskal-Wallis H 检验 * Kruskal和Wallis在1952年设计的一种类似Wilcoxon秩和检验的方法，以进行多个独立样本比较的非参数检验，简称K－W检验或H检验。 其基本原理与Wilcoxon秩和检验相同。该方法与总体具体属于什么分布无关，将多组样本混合起来按大小编秩，计算两组的平均秩和，比较各组分布的中心位置是否不同。 * 一、多个独立样本比较的 Kruskal-Wallis H 检验 * 多个独立样本资料的比较 1.编秩方法同两组资料（不配对）相同，即，混合编秩； 2.计算检验公式为H 3.确定P值，包括查表法(总例数N≤15时)和大样本法(N＞15时, 用H值与卡方界值做比较) 4.如有统计意义，做组间的两两比较(Nemenyi检验) * 例5 用三种药物杀灭钉螺，每批用200只活钉螺，用药后清点每批钉螺的死亡数、再计算死亡率（%），结果见表8-9。问三种药物杀灭钉螺的效果有无差别？ 本例两样本资料为百分率资料，不满足正态分布，现用Kruskal-Wallis H检验。 1．原始数据的两样本比较 * 甲药 乙药 丙药 死亡率 秩 死亡率 秩 死亡率 秩 32.5 10 16.0 4 6.5 1 35.5 11 20.5 6 9.0 2 40.5 13 22.5 7 12.5 3 46.0 14 29.0 9 18.0 5 49.0 15 36.0 12 24.0 8 Ri 63 ─ 38 ─ 19 ni