1、等级资料的秩和检验易洪刚Department of Epidemiology & Biostatistics, School of Public Health Nanjing Medical UniversityRank Sum Test of Ranked DataBiostatistics2资料的分类 数值变量资料 分类资料 二分类 多分类无序多分类有序多分类(等级资料)Biostatistics3内容提要：o 1 医学研究中的等级资料o 2 秩次与秩和 o 3 两样本比较的秩和检验o 4 多样本比较的秩和检验o 5 配对设计的秩和检验o 6 配伍组比较的秩和检验o 7 秩和检验的正确应用
2、 Biostatistics41 医学研究中的等级资料o 疗 效：痊愈、显效、有效、无效、恶化o 化验结果：、+、+o 体格发育：下等、中下、中等、中上、上等o 心功能分级：I、II、IIIo 文化程度：小学、中学、大学、研究生o 营养水平：差、一般、好Biostatistics5等级资料的特点o 既非呈连续分布的定量资料，也非仅按性质归属于独立的若干类的定性资料；o 比“定量”粗，而比一般的“定性”细；o 等级间既非等距，亦不能度量。 Biostatistics62 秩次与秩和 o Wilcoxon在1945年首先提出了比较两个总体分布函数的秩和检验。秩和检验以及其它的秩检验法，都是建立在秩
3、及秩统计量基础上的非参数方法。o 秩次(rank)，秩统计量 是指全部观察值按某种顺序排列的位序；o 秩和(rank sum) 同组秩次之和。Biostatistics7例1 编秩(page86) A组： 、+、+、+、+B组： +、+、+、+、+、+A： + + + + B： + + + + + + 1 2 3 4 5 7 6 8 9 10 11 12 1 2 4.5 4.5 4.5 8.5 4.5 8.5 8.5 8.5 11 12 1 2 4.5 4.5 4.5 7 4.5 8 9 10 11 12 1 2 4.5 4.5 4.5 8.5 4.5 8.5 8.5 8.5 11.5 11.
4、5 尿白细胞：秩次相同(tie)取平均秩次！Biostatistics8 秩和o A组： 、+、+、+、 + 秩和： 1 2 4.5 4.5 4.5 8.5 TA25 o B组： +、+、+、+、+、+ 秩和： 4.5 8.5 8.5 8.5 11.5 11.5 TB53 TA+TB=N(N+1)/2=78Biostatistics9o 秩次：在一定程度上反映了等级的高低；o 秩和：在一定程度上反映了等级的分布 位置。o 对等级的分析，转化为对秩次的分析。o 秩和检验就是通过秩次的排列求出秩和，进行假设检验。2 秩次与秩和 Biostatistics10Sum of ranks = 37Sum
5、 of ranks = 41769213458101112假如两组等级分布相同In this case the sum of ranks for the two samples will be close to one another.Populations2 秩次与秩和 Biostatistics11Sum of ranks = 29Sum of ranks = 49769213458101112假如两组等级分布不相同2 秩次与秩和 Biostatistics12 3 两样本比较的秩和检验o 检验假设 H0 ：A、B两组等级分布相同； H1 ：A、B两组等级分布不同。 =0.05。 Dist
6、ribution of two populations when their locations are samePopulation LocationsThe location of popn 1 is to the left of the location of popn 2The location of popn 1 is to the right of the location of popn 2population 1population 2population 2population 1Biostatistics15H0: The two population locations
7、are the same.H1: The location of population 1 is different from the location of population 2 3 两样本比较的秩和检验Biostatistics16基本思想o 如果H0 成立，即两组分布位置相同， 则A组的实际秩和应接近理论秩和n1(N+1)/2; (B组的实际秩和应接近理论秩和n2(N+1)/2). 或相差不大，差值很大的概率应很小 。o 如果相差较大，超出了预定的界值，则可认为H0不成立。Biostatistics17 A组 B组 和 o 实际秩和 25 53 78 o 理论秩和 n1(N+1)/2
8、 n2(N+1)/2 N(N+1)/2 39 39 78o 差值 14 14 0 抽样误差？o如果H0成立，则理论秩和与实际秩和之差纯粹由抽样误差造成。 基本思想Biostatistics18两样本秩和检验 T 界值 n1=6,n2-n1=0 双侧 单侧 2850 0.10 0.05 26 52 0.05 0.025 24 54 0.02 0.01 23 55 0.01 0.005 间距 22 26 30 326(12+1)/2=39(理论值) Biostatistics19 检验结果o 如果H0成立，则按0.05水准， A组秩和之界值为2652。o 现A组的实际秩和为25，在界值之外，故拒绝
9、H0，接受H1，认为两组的分布位置不同。Biostatistics20 秩和检验的结论判断o A组的实际秩在界值之外， (小于或等于下界，大于或等于上界) 则拒绝H0，接受H1。o A组的实际秩在界值之内， (大于下界且小于上界) 则不拒绝H0。Biostatistics21例题 【例】 某实验室观察在缺氧条件下猫和兔的生存时间，结果见表，试检验在缺氧条件下猫和兔的生存时间有无差别？ Biostatistics22例题Biostatistics23例题Biostatistics24o检验假设 H0 ：猫和兔在缺氧条件下生存时间总体分布相同 ； H1 ：猫和兔在缺氧条件下生存时间总体分布不同 。
10、 =0.05。o检验统计量T值 n1=8，n2=12，检验统计量T=127.5 o确定P值和作出推断结论 查附表10得T界值是58110。则双侧P0.05，按 =0.05水准，拒绝H0，接受H1，差异有统计学意义，故可认为在缺氧条件下猫的生存时间较兔长。例题Biostatistics25例8.2 用复方猪胆胶囊治疗老年性慢性支气管炎患者403例，疗效见表8.1第(1) (3)栏。问该药对此两型支气管炎疗效是否相同？例题（page87）Biostatistics26例题（ page87 ）Biostatistics27例题（ page87 ）Biostatistics28例题（ page87 ）
11、Biostatistics29o建立检验假设 H0：两型老慢支疗效分布相同 ； H1：两型老慢支疗效分布不同； =0.05。o检验统计量 n1=182, n2=221，检验统计量T=40682.5 。 例题（ page88 ）Biostatistics30样本含量较大时，用 u 检验相同秩次较多时需要校正Biostatistics32o确定P值和作出推断结论 P0.01，按 =0.05水准，拒绝H0 ，接受H1，差异有统计学意义。可认为复方猪胆胶囊治疗老年性慢性支气管炎喘息型与单纯型的疗效有差别。例题（ page88 ）Biostatistics334 多组比较的秩和检验o Kruskal-W
12、allis法o 先对所有数据编秩，求各组秩和；o 计算 H 统计量；o 查 H 界值表，或用近似 2 检验， 计算 P 值；o 界定 P 值，作出结论。 Biostatistics34 【例】 某医生在研究再生障碍性贫血时，测得不同程度再生障碍性贫血患者血清中可溶性CD8抗原水平（U/ml），问不同程度再生障碍性贫血患者血清中可溶性CD8抗原水平有无差别？ 例题Biostatistics35例题 (假设不满足参数检验的要求)Biostatistics36例题Biostatistics37多组等级比较的检验假设o 建立检验假设 H0 ：各组总体的等级分布相同；H1 ：各组总体的等级分布不同或不全
13、相同。 =0.05。o 计算检验统计量H值 Biostatistics38H近似服从 =k1的2分布。 多组等级比较的检验假设Biostatistics39o 确定P值和作出推断结论 本例k=3，查附表12, 得P20.005,2，P0.05。按双侧0.05水准，不拒绝H0，差异无统计学意义，故尚不能认为两法测定血清中谷-丙转氨酶含量有差别。查表时，若T在T上、下界值范围外，则P。 例题Biostatistics616 配伍组比较的秩和检验(page93)o Friedman法o 分别对每一配伍组中的观察值进行编秩；o 计算Friedman M 统计量；o 查M界值表，或用近似 2 检验，计算
14、 P 值；o 界定 P 值，作出结论。 Biostatistics62葡萄酒的等级例： 白 兰 地评委WXYZA四二一三B四一二三C三一二四D四二三一E三一二四Ri (4)(2)(1)(3)(4)(1)(2)(3)(3)(1)(2)(4)(4)(2)(3)(1)(3)(1)(2)(4)1871015Biostatistics63葡萄酒等级资料的分析o H0：对四种葡萄酒评判结果的总体分布相同； H1：评判结果的总体分布不同或不全相同。 =0.05。 o 式中b为配伍组数，k为处理组数。本例b=5，k=4。Biostatistics64检验结果o 本例b=5，k=4，查M界值表得：M0.05=7
15、.80，M0.01=9.96。M0.05MM0.01，0.01P50mg”50mg”等等 ) )缺点：方法比较粗糙，对于符合参数检验条件者，采用非参数检验会损失部分信息，其检验效能较低；样本含量较大时，两者结论常相同。Biostatistics68秩和检验用于定量资料o 定量资料中：n 极度偏态资料，或个别数值偏离过大n 各组离散度相差悬殊n 资料中含有不确定值o 大于5年 o 0.001o 1:1024以上 n 兼有等级和定量性质的资料Biostatistics69心肌坏死面积例 用二种食物配方饲养大白鼠，观察心肌坏死面积。分析二组间的差异。A组：(n=29, mean 3.61) 0, 0
16、, 0, 0, 0, 0, 0, 0, 0, 0, 0.2, 0.3, 0.4, 0.4, 0.6, 1, 1.6, 2.2, 2.6, 3.3, 4.3, 5.1, 5.4, 5.5, 6.1, 6.2, 9.7, 13.8, 36B组： (n=28, mean1.06) 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.2, 0.2, 0.2 0.3, 0.4, 0.4, 0.9, 0.9, 1.3, 1.7, 2.8, 7.4, 13Biostatistics70特点：o0特别多，无论用什么变换均不可能改变分布的偏性；o0与其它数的区别是质的区
17、别。Biostatistics71分析一：二组心肌坏死率比较，按四格表作2检验： 无心肌坏死 有心肌坏死 合计 心肌坏死率 甲组 10 19 29 65.5% 乙组 15 13 28 46.4% 2= 2.1081, P = 0.147Biostatistics72分析二：二组平均心肌坏死面积的比较二组平均心肌坏死面积的 t 检验： t =1.7755，P = 0.0814Biostatistics73分析三：按等级资料处理：两组秩和检验： n 秩和 理论秩和 A组29 968 841B组28 685 812合计57 1653 1653u = 2.119， P =0.0341。Biostati
18、stics74成组设计两样本比较o 如资料满足 t 检验的条件,应该用 t 检验进行分析。此时，如果对这类资料用Wilcoxon秩和检验，实际上是将观察单位的具体数值舍弃不用，只保留了秩次的信息，使检验功效降低；尤其样本含量较小时，降低更加明显。 o 如资料不满足 t 检验的条件，而用了t 检验，同样降低了检验效能。Biostatistics75秩变换检验o 基本思想：n 将等级转换为秩次；n 对秩次作效应的参数检验(t检验，方差分析，相关等)。o 在样本含量较大时，秩变换检验与秩和检验结果相同或相近。Biostatistics76Biostatistics77例 不同时间空气中CO2的含量的
19、秩和检验课前课中课后含量含量含量(1)(3)(5)0.484.452.950.534.733.070.554.773.180.554.823.200.584.893.300.625.004.45Biostatistics78例 不同时间空气中CO2的含量的秩和检验课前课中课后含量秩次含量秩次含量秩次(1)(2)(3)(4)(5)(6)0.4814.4512.52.9570.5324.73143.0780.5534.77153.1890.5544.82163.20100.5854.89173.30110.6265.00184.4512.22192.557.5R1R2R3Biostatistics79H0 ：各组总体的等级分布相同；H1 ：各组总体的等级分布不同或不全相同。 =0.05。例 不同时间空气中CO2的含量的秩和检验Biostatistics80n例数较少时可以查附表H界值表得到界值H0.05=5.801，P0.05。n按0.05水准，拒绝H0，接受H1，可以认为三个不同时间空气中CO2的含量不同。例 不同时间空气中CO2的含量的秩和检验
