北京大学医学部医学统计学基础第5讲 非参数检验.pptVIP

非参数检验 公卫学院流行病与统计学系 易伟宁 讲师yiweining@163.com 概念 t检验和方差分析，都要求个体变量值服从正态分布，或根据中心极限定理，当样本较大时样本均数服从正态分布，两组或多组比较时，还要求方差齐性。这样的检验方法称为参数检验。 而非参数检验是对样本来自的总体分布不作要求的一类假设检验方法。比较常用的非参检验方法是秩和检验。其优点是适用资料范围广，缺点是功效稍低。 常用的参数检验方法和条件 检验方法名称 主要应用条件 单样本t检验 已知总体均数 样本均数服从正态 成组双样本t检验 两样本均数服从正态 两总体方差相等 配对双样本t检验 差值均数服从正态 成组多样本比较的 方差分析 各样本均数服从正态 各组的总体方差相等 部分参数和非参数检验 资料类型 参数检验 非参数检验 成组双样本 成组t检验 Wilcoxon秩和检验 配对双样本 配对t检验 符号秩检验 成组多样本 方差分析 Kruskal-Wallis检验 Wilcoxon秩和检验 比较实验组(n1=10)和对照组(n2=12)带瘤小鼠的生存天数有无差别 考虑作秩和检验的原因：实验组有存活90天以上小鼠，存在末端不确定数据，根据经验，生存天数有可能不服从正态分布，而且未能通过方差齐性检验 Wilcoxon秩和检验思想：对总共22个观察值统一排顺序(秩)，计算秩和。当H0成立时，两组平均秩次应该相差不大 查表法的检验步骤 H0: 两组小鼠生存天数的总体分布相同 H1: 两组的总体分布不同 α=0.05 两组秩和分别为T1=170，T2=83，课本 取例数较小组的秩和，即T1为检验统计 量T，查T界值表。由于170在84~146的 范围之外，所以P0.05。按照α=0.05的 检验水准，拒绝H0，接受H1，认为… R统计软件的程序 X=1:22 T=c() for (i in 1:100000 ) { T[i]=sum(sample(x,12,replace = F )) } hist(T) 符号秩检验 A、B两种方法测定10名健康人的尿汞值，问所得结果有无差别 考虑作秩和检验原因：不清楚方法B减去方法A的差值均数是否服从正态分布 思想：分别计算正负秩和，任取一个作为检验统计量T。若H0成立，随着n增大，T的分布将逼近均数为n(n+1)/4，方差为n(n+1)(2n+1)/24的正态分布 检验步骤 H0: 差值的总体中位数为零 H1:差值的总体中位数不为零 α=0.05 正、负秩和分别为18.5和26.5，取18.5为 检验统计量T，查T界值表。由于18.5在 5~40的范围之内，所以P0.05。按照 α=0.05的检验水准，尚不能拒绝H0， 认为差值的总体中位数为零，即两法无 差别(例数小，SPSS正态近似效果不好) R统计软件的程序 x=1:9 T=c() for (i in 1:100000 ) { T[i]=sum(x * (rbinom(9,1,0.5)) ) } hist(T) sum(T=18.5)/100000 *2 Kruskal-Wallis检验 比较小白鼠接种三种不同菌型的伤寒杆菌后的存活天数有无差别 考虑作秩和检验的原因：根据经验，生存天数有可能不服从正态分布 检验步骤 H0: 三组小鼠生存天数的总体分布相同 H1: 三组的总体分布不全相同 α=0.05 检验统计量H为9.94，自由度为k-1=2，可证明，当例数较大时H近似服从自由度为2的卡方分布，故渐进P值为0.0070.05。按照α=0.05的检验水准，拒绝H0，接受H1，认为… 两两比较的问题 当Kruskal-Wallis检验得出P0.05的结论后，如果想进一步推断是哪两组的总体分布不同，可用Nemenyi法检验。 软件没有提供Nemenyi法作两两比较的途径，需要手工计算。或者也可以考虑将全部观察值统一排秩后，对秩次作单因素方差分析，并进行两两比较。 Nemenyi法检验 H0：第1组和第2组生存天数的总体分布相同 H1：两组生存天数的总体分布不同 α=0.05 经查卡方界值表，P0.05，接受H1，认为…. 同理可得另两组卡方值为8.15，0.016 相同秩 个数t 2.0 3 4.5 2 7.0 3 10.5 4 15.5 6 21.0 5 26.5 2 28.5 2 两组等级资料比较 24名正常人痰液标本的白细胞检查结果 - + ++ +++ 11 10 3 0 44名慢性支气管炎病人痰液标本的结果 - + ++ +++ 5 18 16 5 问二者是否有不同？ 多组等级资料比较 三种不同药物治疗单纯性慢性支气管炎，结果为 药物 治愈 显效 好转 无效 A 17 51 33 7 B