医学统计学等级资料的秩和检验课件

1、等级资料的秩和检验,Department of Epidemiology (B组的实际秩和应接近理论秩和n2(N+1)/2). 或相差不大，差值很大的概率应很小 。 如果相差较大，超出了预定的界值，则可认为H0不成立。,A组 B组 和 实际秩和 25 53 78 理论秩和 n1(N+1)/2 n2(N+1)/2 N(N+1)/2 39 39 78 差值 14 14 0 抽样误差？ 如果H0成立，则理论秩和与实际秩和之差纯粹由抽样误差造成。,基本思想,两样本秩和检验 T 界值,n1=6,n2-n1=0 双侧 单侧 2850 0.10 0.05 26 52 0.05 0.025 24 54 0.0

2、2 0.01 23 55 0.01 0.005,间距 22 26 30 32,6(12+1)/2=39(理论值),检验结果,如果H0成立，则按0.05水准， A组秩和之界值为2652。 现A组的实际秩和为25，在界值之外，故拒绝H0，接受H1，认为两组的分布位置不同。 昆明治疗癫痫病医院 昆明治疗癫痫病医院 昆明治疗癫痫病医院 昆明癫痫病治疗,秩和检验的结论判断,A组的实际秩在界值之外， (小于或等于下界，大于或等于上界) 则拒绝H0，接受H1。 A组的实际秩在界值之内， (大于下界且小于上界) 则不拒绝H0。,例题,【例】 某实验室观察在缺氧条件下猫和兔的生存时间，结果见表，试检验在缺氧条件

3、下猫和兔的生存时间有无差别？,例题,检验假设 H0 ：猫和兔在缺氧条件下生存时间总体分布相同 ； H1 ：猫和兔在缺氧条件下生存时间总体分布不同 。 =0.05。 检验统计量T值 n1=8，n2=12，检验统计量T=127.5 确定P值和作出推断结论 查附表10得T界值是58110。则双侧P0.05，按 =0.05水准，拒绝H0，接受H1，差异有统计学意义，故可认为在缺氧条件下猫的生存时间较兔长。,例题,例8.2 用复方猪胆胶囊治疗老年性慢性支气管炎患者403例，疗效见表8.1第(1) (3)栏。问该药对此两型支气管炎疗效是否相同？,例题（page89）,例题（ page90 ）,例题（ pa

4、ge90 ）,建立检验假设 H0：两型老慢支疗效分布相同 ； H1：两型老慢支疗效分布不同； =0.05。 检验统计量 n1=182, n2=221，检验统计量T=40682.5 。,例题（ page90 ）,样本含量较大时，用 u 检验,相同秩次较多时需要校正,确定P值和作出推断结论 P0.01，按 =0.05水准，拒绝H0 ，接受H1，差异有统计学意义。可认为复方猪胆胶囊治疗老年性慢性支气管炎喘息型与单纯型的疗效有差别。,例题（ page90 ）,4 多组比较的秩和检验,Kruskal-Wallis法 先对所有数据编秩，求各组秩和； 计算 H 统计量； 查 H 界值表，或用近似 2 检验，

5、 计算 P 值； 界定 P 值，作出结论。,【例】 某医生在研究再生障碍性贫血时，测得不同程度再生障碍性贫血患者血清中可溶性CD8抗原水平（U/ml），问不同程度再生障碍性贫血患者血清中可溶性CD8抗原水平有无差别？,例题,例题 (假设不满足参数检验的要求),例题,多组等级比较的检验假设,建立检验假设 H0 ：各组总体的等级分布相同； H1 ：各组总体的等级分布不同或不全相同。 =0.05。,计算检验统计量H值,H近似服从 =k1的2分布。,多组等级比较的检验假设,确定P值和作出推断结论 本例k=3，查附表12, 得P0.001。按 = 0.05水准拒绝H0，接受H1，差异有统计学意义，可认为

6、不同程度再生障碍性贫血患者血清中可溶性CD8抗原水平有差别。,多组等级比较的检验假设,H 的校正(page91 公式8.4请更正),当有相同秩次时，H 需校正：,例8.3 某医院用三种复方小叶枇杷治疗老年性慢性支气管炎，数据见表8.2第(1)(4)栏，试比较其疗效有无差异。,例题（page91）,例题（ page91 ）,建立检验假设 H0：三药疗效总体分布相同 ； H1：三药疗效总体分布不同或不全相同 ； = 0.05。,例题（ page91 ）,计算检验统计量H值,例题（ page91 ）,由于此资料的相同秩次很多，须校正：,例题（page120）,确定P值和作出推断结论 故按=k1=31

7、=2查2界值表，得20.005,2 =10.60，Hc20.005,2，P0.005。按 =0.05水准拒绝H0，接受H1，差异有统计学意义。认为三药疗效有差别。,例题（page120）,多组间的两两比较,如果多组等级比较拒绝H0，认为组间存在差异，则需进行两两比较。 两两比较的方法很多。,5 配对符号秩和检验,Wilcoxon符号秩和检验 计算等级之差值，对差值进行编秩，并冠以差值的符号。 查 T 界值表，或用近似 u 检验，计算 P 值； 界定 P 值，作出结论。,符号秩和检验的假设,H0：差值的总体中位数为0； H1：差值的总体中位数不为0。 =0.05。 当n50时，查界值表 当n50

8、时，用u近似,扁平足疗效例,用配对设计观察两种方法治疗扁平足效果记录如下，问那种方法好。 病例号 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 A法 好 好 好 好 差 中 好 好 中 差 好 差 好 中 好 中 B法 差 好 差 中 中 差 中 差 中 差 好 差 中 差 中 差,差值 2 0 2 1 -1 1 1 2 0 0 0 0 1 1 1 1,符号秩和检验,病例号 A法 B法 差值 秩次 (1)(2)(3)(4)(5) 131210 2330 331210 43214.5 512-1-4.5 62114.5 73214.5 831210 9220 10

9、110 11330 12110 133214.5 142114.5 153214.5 162114.5,T+ =61.5 T- = 4.5,符号秩和检验的基本思想,总秩和为TN(N+1)/2 如H0成立，则正负各半，T+ 与 T 均接近 N(N+1)/4。 如果相差太大，超出了事先规定的界值， 则H0不成立。,符号秩和检验 T 界值表,N=11 双侧 单侧 1353 0.10 0.05 10 56 0.05 0.025 7 59 0.02 0.01 5 61 0.01 0.005,间距 40 46 52 56,11(11+1)/4=33(理论值),本例T=61.5，T=4.5，已超出附表10中

10、双侧P0.01相对应的界值561，故得P 0.01。按 =0.05水准拒绝H0，接受H1，认为两种方法疗效的总体分布不同。结合资料可认为A法优于B法。,扁平足疗效例,u 的校正,当重复的秩次较多时，u 需要校正：,【例】 某医院检验科试用新旧两种方法检测谷-草转氨酶，新方法的检测时间由20分钟缩短为10分钟。用两种方法测一份血清，结果见表第（2）、（3）栏，问两法所得结果有无差别？,例题,例题,建立检验假设 H0：两种方法所测值的总体分布相同（差值总体中位数Md = 0）； H1：两种方法所测值的总体分布不同（差值总体中位数Md 0）； = 0.05。 求检验统计量T值 T+ =18.5，T-

11、 =36.5，其和为55， n(n+1)/2=10(10+1)/2=55。,例题,符号秩和检验 T 界值表,N=10 双侧 单侧 1045 0.10 0.05 8 47 0.05 0.025 5 50 0.02 0.01 3 52 0.01 0.005,间距 40 46 52 56,10(10+1)/4=27.5(理论值),确定P值和作出推断结论 本例n=10，T=18.5，查附表10，得 =0.05时的T界值为847，T在双侧界值范围内，故P 0.05。按双侧0.05水准，不拒绝H0，差异无统计学意义，故尚不能认为两法测定血清中谷-丙转氨酶含量有差别。 查表时， 若T在T上、下界值范围外，则

12、P。,例题,6 配伍组比较的秩和检验,Friedman法 分别对每一配伍组中的观察值进行编秩； 计算Friedman M 统计量； 查M界值表，或用近似 2 检验，计算 P 值； 界定 P 值，作出结论。,葡萄酒的等级例：,白 兰 地 评委WXYZ A四二一三 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) 1871015,葡萄酒等级资料的分析,H0：对四种葡萄酒评判结果的总体分布相同； H1：评判结果的总体分布不同或不全相同。 =0.05。 式中b为配伍组数

13、，k为处理组数。本例b=5，k=4。,检验结果,本例b=5，k=4，查M界值表得：M0.05=7.80，M0.01=9.96。M0.05MM0.01，0.01P0.05。按 =0.05水准拒绝H0，接受H1，差异有统计学意义，即可认为四种酒的等级有差别。,7 秩和检验的正确应用,主要对等级资料进行分析； 秩和检验可用于任意分布(distribution free)的资料； 非参数统计的主要优点是不受总体分布的限制，适用范围广。 非参数统计的主要缺点是符合参数检验的资料（如两样本均数比较的t检验），如用非参数检验，因没有充分利用资料提供的信息，检验效率降低。,已知总体分布类型，对未知参数进行统计

14、推断,依赖于特定分布类型，比较的是参数,参数检验 （parametric test）,非参数检验 （nonparametric test）,对总体的分布类型不作严格要求,不受分布类型的影响，比较的是总体分布位置,优点：方法简便、易学易用，易于推广使用、应用范围广；可用于参数检验难以处理的资料(如等级资料，或含数值“50mg”等 ),缺点：方法比较粗糙，对于符合参数检验条件者，采用非参数检验会损失部分信息，其检验效能较低；样本含量较大时，两者结论常相同。,秩和检验用于定量资料,计量资料中： 极度偏态资料，或个别数值偏离过大 各组离散度相差悬殊 资料中含有不确定值 大于5年 0.001 1:102

15、4以上 兼有等级和定量性质的资料,心肌坏死面积例,用二种食物配方饲养大白鼠，观察心肌坏死面积。分析二组间的差异。 A组：(n=29, mean 3.61) 0, 0, 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, 36 B组： (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,

16、0.9, 1.3, 1.7, 2.8, 7.4, 13,特点：,0特别多，无论用什么变换均不可能改变分布的偏性； 0与其它数的区别是质的区别。,分析一：二组心肌坏死率比较，按四格表作2检验：,无心肌坏死 有心肌坏死 合计 心肌坏死率 甲组 10 19 29 65.5% 乙组 15 13 28 46.4% C2= 1.404, P = 0.24,分析二：二组平均心肌坏死面积的比较,二组平均心肌坏死面积的 t 检验： t =1.7755，P = 0.0814,分析三：按等级资料处理：,两组秩和检验： n 秩和 理论秩和 A组29 968 841 B组28 685 812 合计57 1653 1653 u = 2.11