医学统计学秩转换的非参数检验课件.ppt

1、如果总体分布为已知的数学形式，对其总体参数作假设检验。如：t 检验和 F 检验。对总体分布不作严格假定，又称任意分 布检验（distribution-free test），它直接对总体分布作假设检验。参数检验非参数检验271 推断一个总体表达分布位置的中位数M（非参数）和已知M0、两个或多个总体的分 布是否有差别。先将数值变量从小到大，或等级从弱到强转换成秩后，再计算检验统计量。秩转换的非参数检验特点：秩转换的非参数检验假设检验的结果对总体分布的形状差别不敏感，只对总体分布的位置差别敏感。371应用范围：对于计量资料：1.不满足正态和方差齐性条件的小样本资料；2.分布不明的小样本资料；3.一端

2、或二端是不确定数值（如0.5、5.0等）的资料（必选）；对于等级资料：若选行列表资料的 检验，只能推断构成比差别，而选秩转换的非参数检验，可推断等级强度差别。2471注意：注意：如果已知其计量资料满足（或近似如果已知其计量资料满足（或近似满足）满足）检验或检验或 检验条件，当然选检验条件，当然选 检检验或验或 检验，因为这时若选秩转换的非检验，因为这时若选秩转换的非参数检验，会降低检验效能。参数检验，会降低检验效能。ttFF571本章内容：第一节 配对样本比较的Wilcoxon符号秩检验第二节 两个独立样本比较的秩和检验（成组设计）第三节 完全随机设计多个样本比较的Kruskal-Wallis

3、 H 检验 （成组设计）第四节 随机区组设计多个样本比较的Friedman M检验671第一节配对样本比较的Wilcoxon符号秩检验 7711配对样本差值的中位数和0比较 目的是推断配对差值的总体中位数是否和0有差别，即推断配对的两个相关样本所来自的两个总体中位数是否有差别。方法步骤见例8-1。871 例8-1 对12份血清分别用原方法（检测时间20分钟）和新方法（检测时间10分钟）测谷-丙转氨酶，结果见表8-1的（2）、（3）栏。问两法所得结果有无差别？971107111710H：差值的总体中位数0dM 1H：0dM 0.05 1.建立检验假设，确定检验水平2.求检验统计量T值 省略所有差

4、值为0的对子数，令余下的有效对子数为n n，见表8-1第（4）栏，本例 n=12-1=11；检验步骤 若多个差值为若多个差值为0，可通过提高测量工具的精，可通过提高测量工具的精度来解决（度来解决（减少差值为减少差值为0 的个数的个数）。）。1271 按差值的绝对值从小到大编秩，然后分别冠以正负号。遇差值绝对值相等则取平均秩，称为相同秩（ties）（样本较小时，如果相同秩较多，检验结果会存在偏性，因此应提高测量精度，尽量避免出现较多的相同秩）,表8-1第（4）栏差值的绝对值为2有2个，其秩依次应为1，2，皆取平均秩为1.5，见表8-1第（5）、（6）栏；任取正秩和或负秩和为T，本例取T=11.5

5、。13713.确定P值，作出推断结论 当n50时，查T界值表（附表9,p824），判断原则：内大外小。1471若当n50，超出附表9范围，可用正态近似法作u检验。3(1)/4()(1)(21)2448jjTn nuttn nn式中jt（j=1,2,）为第j个相同秩的个数 假定相同秩（即平均秩）中有 2 个 1.5，5 个 8，3 个 14，则 12t，25t，33t，33()(2jjtt 332)(55)(33)150 Z 检验1571 配对等级资料采用符号秩检验最好选用大样本。注意(相同秩次多)16712单个样本中位数和已知总体中位数比较 目的：推断样本所来自的总体中位数M和某个已知的总体中

6、位数M0是否有差别。用样本各变量值和M0的差值，即推断差值的总体中位数和0是否有差别。方法步骤见例8-2。1771 例8-2 已知某地正常人尿氟含量的中位数为45.30 。今在该地某厂随机抽取12名工人，测得尿氟含量见表8-2第（1）栏。问该厂工人的尿氟含量是否高于当地正常人的尿氟含量？mol/L187119712071检验步骤 0H：尿氟含量的总体中位数45.30M 1H：45.30M 0.05 据表8-2第（3）、（4）栏，取T=1.5。任取正秩和或负秩和为T2171第二节两个独立样本比较的Wilcoxon秩和检验目的：用于推断两个独立样本所来自的目的：用于推断两个独立样本所来自的两个总体

7、分布是否有差别。两个总体分布是否有差别。22710:两个总体分布位置相同两个总体分布位置相同23710 0:两个总体分布位置相同两个总体分布位置相同1 1:两个总体分布位置不相同两个总体分布位置不相同 （或左、右一些）（或左、右一些）24711原始数据的两样本比较 例8-3 对10例肺癌病人和12例矽肺0期工人用X光片测量肺门横径右侧距RD值（cm），结果见表8-5。问肺癌病人的RD值是否高于矽肺0期工人的RD值？2571肺癌病人 矽肺 0 期工人 RD 值 秩 RD 值 秩 2.78 1 3.23 2.5 3.23 2.5 3.50 4 4.20 7 4.04 5 4.87 14 4.15

8、6 5.12 17 4.28 8 6.21 18 4.34 9 7.18 19 4.47 10 8.05 20 4.64 11 8.56 21 4.75 12 9.60 22 4.82 13 4.95 15 5.10 16 110n 1141.5T 212n 2111.5T 表8-5 肺癌病人和矽肺0期工人的RD值（cm）比较 本例两样本资料经方差齐性检验，推断得两总体方差不等（0.01P）26710H：肺癌病人和矽肺 0 期工人的RD 值总体分布位置相同 1H：肺癌病人的 RD 值高于矽肺 0 期工人的 RD 值 0.05 检验 步骤求检验统计量T值：2771确定P值，作出推断结论：2871

9、13123(1)/2()(1)112jjTn Nuttn n NNN式中(1,2,)jtj L为第j个相同秩的个数 29712频数表资料和等级资料的两样本比较 例8-4 39名吸烟工人和40名不吸烟工人的碳氧血红蛋白HbCO(%)含量见表8-6。问吸烟工人的HbCO(%)含量是否高于不吸烟工人的HbCO(%)含量？计量资料为频数表资料，是按数量区间分组；等级资料是按等级分组。现以等级资料为例，方法步骤见例8-4。3071 表8-6 吸烟工人和不吸烟工人的HbCO(%)含量比较 31710H：吸烟工人和不吸烟工人的HbCO含量总体分布位置相同 1H：吸烟工人的HbCO含量高于不吸烟工人的HbCO

10、含量 0.05 先确定各等级的合计人数、秩范围和平均秩，见表8-6的（4）栏、（5）栏和（6）栏，再计算两样本各等级的秩和，见（7）栏和（8）栏；本例T=1917（此组样本量较小）；327131917 39(79 1)/23.702339 40(79 1)52230(1)127979u139n ，240n，394079N 计算u值333333()(33)(3131)(2727)(1414)(44)52230jjtt3371第三节完全随机设计多个样本比较的Kruskal-Wallis H 检验3471一、多个独立样本比较的 Kruskal-Wallis H 检验357136711原始数据的多个样

11、本比较 例8-5 用三种药物杀灭钉螺，每批用200只活钉螺，用药后清点每批钉螺的死亡数、再计算死亡率（%），结果见表8-9。问三种药物杀灭钉螺的效果有无差别？3771甲药 乙药 丙药 死亡率 秩 死亡率 秩 死亡率 秩 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 iR 63 38 19 in 5 5 5 表8-9 三种药物杀灭钉螺的死亡率（%）比较 3871求检验统计量H值：212()3(1)(1)iiRHNN Nn397

12、122212633819()3(15 1)9.7415(15 1)5H本例确定P值，作出推断结论：当样本个数3g 和每个样本例数5in 时，查H界值表（附表 11，P536）本例15N，1235nnn，查附 表11 得0.01P，按0.05水准拒绝0H，接 受1H，可认为三种药物杀灭钉螺的效果不同。4071当 各 样 本 数 据 存 在 相 同 秩 时，按 公 式（8-4）算 得 的H值 偏 小，按下 式 求 校 正CH值。/CHH C，331()/()jjCttNN 4171 例8-6 比较小白鼠接种三种不同菌型伤寒杆菌9D、11C和DSC1后存活日数，结果见表8-10。问小白鼠接种三种不同

13、菌型伤寒杆菌的存活日数有无差别？42719D 11C D SC1 存 活 日 数 秩 存 活 日 数 秩 存 活 日 数 秩 2 2 5 10.5 3 4.5 2 2 5 10.5 5 10.5 2 2 6 15.5 6 15.5 3 4.5 6 15.5 6 15.5 4 7 6 15.5 6 15.5 4 7 7 21 7 21 4 7 8 24 7 21 5 10.5 10 26.5 9 25 7 21 12 30 10 26.5 7 21 11 28.5 11 28.5 iR 84 169 212 in 10 9 11 iR 8.40 18.78 19.27 表8-10 小白鼠接种三种

14、不同菌型伤寒杆菌的存活日数比较 4371本例为时间资料，不是正态分布，现用Kruskal-Wallis H检验。0H：接种三种总体分布位置相同 1H：接种三种总体分布位置不全相同 0.05 4471312。查附表 8（2界值表）得0.0050.01P，按0.05水准拒绝0H，接受1H，可认为小白鼠接种不同菌型伤寒杆菌的存活日数有差别。1091130N。按公式（8-4）和公式（8-5）2221284169212()3(30 1)9.7730(30 1)10911H 333333333(3 3)(2 2)(3 3)(4 4)(6 6)(5 5)(2 2)(2 2)10.9830 30C 9.77/

15、0.98 9.97CH 4571 例8-7 四种疾病患者痰液内嗜酸性白细胞的检查结果见表8-11。问四种疾病患者痰液内的嗜酸性白细胞有无差别？2频数表资料和等级资料的多个样本比较 4671白细胞 支气管 扩张 肺水肿 肺癌 病毒性呼 吸道感染 合计 秩范围 平均秩 -0 3 5 3 11 111 6+2 5 7 5 19 1230 21+9 5 3 3 20 3150 40.5+6 2 2 0 10 5160 55.5 iR 739.5 436.5 409.5 244.5 in 17 15 17 11 60 iR 43.50 29.10 24.09 22.23 表8-11 四种疾病患者痰液内的

16、嗜酸性白细胞比较 47710H：四种疾病总体分布位置相同 1H：四种疾病总体分布位置不全相同 0.05 如表 8-11 第栏的秩和1R是用第栏各等级的频数与第（8）栏平均秩相乘再求和，即10(6)2(21)9(40.5)6(55.5)739.5R，仿此得表811iR下部行。检验步骤4871222212739.5436.5409.5244.5()3(60 1)14.2860(60 1)17151711H 33333(1111)(1919)(2020)(1010)10.926060C 14.28/0.9215.52CH413。查附表 8（2界值表，P533）得0.005P，按0.05水准拒绝0H，

17、接受1H，可认为四种疾病患者痰液内的嗜酸性白细胞有差别。4971两独立样本比较，若1n和2n较大，如为频数表资料或等级资料时，本章第二节介绍了用Wilcoxon 秩和检验的公式（8-2），也可以用本节介绍的 Kruskal-Wallis H检验的公式（8-4）或公式（8-5）。两者的关系是：H（或CH）2u。5071二、多个独立样本两两比较的Nemenyi法检验51715271 例8-8 对例8-6资料（表8-10）作三个样本间的两两比较。0H：任意两存活日数总体分布位置相同 1H：任意两存活日数总体分布位置不同 0.05 5371本例例 8-6 中已算得校正系数0.98C，故据表 8-10

18、下部in行和iR行数据，得 221,2(8.4018.78)6.7230(301)11()0.9812109 同样可算得：21,38.15，22,30.016 54715571第四节随机区组设计多个样本比较的Friedman M检验5671一、多个相关样本比较的Friedman M检验5771 例8-9 8名受试对象在相同实验条件下分别接受4种不同频率声音的刺激，他们的反应率（%）资料见表8-12。问4种频率声音刺激的反应率是否有差别？5871 频率 A 频率 B 频率 C 频率 D 受试号 反应率 秩 反应率 秩 反应率 秩 反应率 秩 1 8.4 1 9.6 2 9.8 3 11.7 4

19、2 11.6 1 12.7 4 11.8 2 12.0 3 3 9.4 2 9.1 1 10.4 4 9.8 3 4 9.8 2 8.7 1 9.9 3 12.0 4 5 8.3 2 8.0 1 8.6 3.5 8.6 3.5 6 8.6 1 9.8 3 9.6 2 10.6 4 7 8.9 1 9.0 2 10.6 3 11.4 4 8 7.8 1 8.2 2 8.5 3 10.8 4 iR 11 16 23.5 29.5 表8-12 8名受试对象对4种不同频率声音刺激的反应率（%）5971本 例8n，4g，32Nng，为百 分 率资 料，不符 合正 态分 布，故 用Friedman M检验

20、。0H：4 种频率声音刺激的反应率总体分布位置相同 1H：4 种频率声音刺激的反应率总体分布位置不全相同 0.05 6071求检验统计量M值：将每个区组的数据由小到大分别编秩，遇数据相等者取平均秩；计算各样本的秩和 ，平均秩和为 ；iR(1)/2Rn g按下式求M值。2222222222()(1)/4(111623.529.5)84(4 1)/4199.5iiMRRRn g g 61716271212 199 514 968 4(4 1)4 13 查2界值表得0.005P。对例8-9，n=8，g=4，已算得M=199.5，则 若n15或g15时，超出附表 12 的范围，可用2近似法，按下式计算

21、2值。212(1)Mng gC，33()1()jjttCn gg 1g6371二、多个相关样本两两比较的检验 6471 例8-10 对例8-9资料（表8-12）作四个样本间的两两比较。0H：任意两反应率总体分布位置相同 1H：任意两反应率总体分布位置不同 0.05 设有g个相关样本，当区组个数n较多时，按下式求第i个样本和第j个样本比较的q值。6571 ijRRqnMS误差23(1)(21)11()612(1)(1)ijjng ggRttnMSng误差q的自由度(1)(1)ng，样本间跨度 a 指把g个样本秩和从小到大排列后ijRR和之间 涵盖的秩和个数（包括ijRR和自身在内）。6671 频

22、率 A 频率 B 频率 C 频率 D 受试号 反应率 秩 反应率 秩 反应率 秩 反应率 秩 1 8.4 1 9.6 2 9.8 3 11.7 4 2 11.6 1 12.7 4 11.8 2 12.0 3 3 9.4 2 9.1 1 10.4 4 9.8 3 4 9.8 2 8.7 1 9.9 3 12.0 4 5 8.3 2 8.0 1 8.6 3.5 8.6 3.5 6 8.6 1 9.8 3 9.6 2 10.6 4 7 8.9 1 9.0 2 10.6 3 11.4 4 8 7.8 1 8.2 2 8.5 3 10.8 4 iR 11 16 23.5 29.5 表8-12 8名受试对

23、象对4种不同频率声音刺激的反应率（%）6771本例根据表8-12有：n=8,g=422222111623.529.51799.5iR，33()226jjtt。8 4(4 1)(2 4 1)/6 1799.5/8 6/120.69(8 1)(4 1)MS 误差1,429.5 117.878 0.69q6871表8-15 表8-12相关样本的两两比较 样本间跨度,1776971(81)(41)21。据相关样本两两比较的q值及、a查附表 4（q界值表），所得P值见表 8-15。可认为频率 A、B和频率 C、D 声音刺激的反应率有差别；尚不能认为频率 A 和频率 B、频率 C 和频率 D 声音刺激的反应率有差别。7071练习题练习题P178 一、最佳选择题一、最佳选择题 全做全做 三、计算分析题三、计算分析题 全做全做71717271