《医学]研究生医学统计学非参数检验》由会员分享，可在线阅读，更多相关《医学]研究生医学统计学非参数检验（70页珍藏版）》请在金锄头文库上搜索。
1、2019/1/14,1/62,郑大公卫统计教研室 平智广,第十章 基于秩次的非参数检验,2019/1/14,2/62,参数统计(parametric statistics) 是以样本来自某已知分布总体（如正态分布、t分布、F分布等）为假设基础，对总体参数（如总体均数、总体方差等）进行估计或检验的方法。,2019/1/14,3/62,在实践中常遇到以下一些资料，如需比较患者和正常人的血铁蛋白、血铅值、不同药物的溶解时间、实验鼠发癌后的生存日数、护理效果评分等，这类资料： 资料的总体分布类型未知； 或某些变量可能无法精确测量： 疼痛程度 无、轻度、中度、重度 疾病治疗结果 控制、显效、有效、无效
2、。,2019/1/14,4/62,非参数统计(nonparametric statistics) 是一种不依赖于总体分布类型，不考虑被研究对象为何种分布以及分布是否已知的统计推断方法，此类方法通常不涉及研究对象的参数。这时统计推断是比较分布而不是参数。,2019/1/14,5/62,非参数统计的基本思想,一组数据的最基本信息是次序，将数值按大小次序排队，每个数值在整个数据中所占的位置和次序，称为秩(rank)。在一定的假设下，这些秩及其统计量的分布是可以求出来的，且与原来的总体分布无关，可进行所需要的统计推断。,2019/1/14,6/62,非参数统计方法的优缺点,优点： 适用范围广：对变量的
3、类型和分布无特殊要求。 对数据要求不严：对某些指标不便准确测定，只能以严重程度，优劣等级，先后次序等作记录的资料也可应用。 对样本量无严格要求。 缺点： 对于符合参数检验的资料如果用非参数检验，由于没有充分利用资料提供的信息，故检验效能低于参数检验，若要使检验效能相同，往往需要更大的样本含量。,2019/1/14,7/62,应用范围：,对于计量资料： 不满足正态和方差齐性条件的小样本资料； 分布不明的小样本资料； 一端或两端是不确定数值（如或某一数值）的资料。 对于等级资料： 若选行列表资料的c2检验，只能推断构成比差别，而选秩转换的非参数检验，可推断等级强度差别。,2019/1/14,8/6
4、2,Wilcoxon于1945年提出的符号秩和检验(Wilcoxon singned-rank test)，亦称符号秩检验，可用以推断总体中位数是否等于某已知数值，也可推断配对样本差值的总体中位数是否为0。,第一节 单样本和配对设计资料的符号秩和检验,2019/1/14,9/62,一、单样本资料的符号秩和检验,Wilcoxon符号秩和检验常用于不满足t检验条件的单样本定量资料的比较。其目的是推断样本中位数与已知总体中位数是否相等。 例1：已知某地正常人尿铅含量的中位数为2.50mmol/L。今在该地随机抽取16名工人，测得尿铅含量，问该厂工人的尿铅含量是否高于当地正常人？,2019/1/14,
5、10/62,表1 16名工人与当地正常人到尿铅含量（mol/L）测定结果,2019/1/14,11/62,尿铅含量的直方图,2019/1/14,12/62,1.建立检验假设，确定检验水准 H0：差值的总体中位数等于0，即该厂工人的尿铅含量与正常人相同 H1：差值的总体中位数大于0，即该厂工人的尿铅含量高于正常人 a=0.05,2019/1/14,13/62,2.计算检验统计量T值 （1）求差值 （2）编秩 差值为0 / 差值的绝对值相同 （3）分别求正、负秩和 （4）确定检验统计量：任意取T+或T-,2019/1/14,14/62,表2 16名工人与当地正常人到尿铅含量（mol/L）测定结果,
6、2019/1/14,15/62,编秩： 将差值按绝对值大小从小到大编秩，并按差值的正负给秩次加上正负号。若差值的绝对值相等，则取其平均秩次，如差值的绝对值为1.88的有3个，它们的位次是11、12、13，取平均秩次为（11+12+13）/3=12。编秩时如遇差值为0，则舍去不计。,2019/1/14,16/62,3. 确定P值，并做出推断结论 （1）查表法 5n50，判断原则：内大外小 本例，n=16，T=28或T=108 ，查用T界值表，得0.01P0.025，按a=0.05检验水准，拒绝H0 ，接受H1，可认为该厂工人尿铅含量高于当地正常人。,2019/1/14,17/62,（2）正态近似
7、法（n50时）超出附表9范围，可用正态近似法作u检验。,2019/1/14,18/62,若多次出现相持现象（如超过25%），求得的Z值偏小，应计算校正的统计量值Zc。,式中tj为第j (j=1，2)个相同秩次的个数。如有2个差值的绝对值为2.5，则t1=2；有3个差值均为5，则t2=3。于是， (23-2)+(33-3),2019/1/14,19/62,二、配对设计资料的符号秩和检验,目的是推断配对样本差值的总体中位数是否与0有差别，即两个总体中位数是否有差别。 例2：对11份工业污水测定氟离子浓度(mg/L)，每份水样同时采用电极法及分光光度法测定，结果见表。问就总体而言，这两种方法的测定结
8、果有无差别？,表3 两法测定11份工业污水中氟离子浓度结果,2019/1/14,21/62,本例配对样本差值经正态性检验，推断得总体不服从正态分布，现用Wilcoxon符号秩检验。 建立检验假设，确定检验水准 H0：差值的总体中位数等于0 H1：差值的总体中位数不等于0 a=0.05 2. 求检验统计量T值 (1) 求差值 (2) 编秩 (3) 求秩和 (4) 确定统计量,2019/1/14,22/62,3. 确定P值，作出推断 查表法 5n 50，查T界值表（附表9） 任取正秩和或负秩和为T，本例为T=11.5或T=43.5 。 P0.10，按照a=0.05水准不拒绝H0 ，尚不能认为两法测
9、定结果有差别。,2019/1/14,23/62,第二节 两组独立样本比较的秩和检验,理论上零假设H0为两总体分布相同，即两样本来自同一总体；其对立假设H1为两总体分布不同。秩和检验对两总体分布形状的差别不敏感，对位置相同、形状不同但类似的两总体分布，推断不出两总体分布形状是否有差别，故H1不能为两总体分布不同，而是两总体分布位置不同，简化为两总体的中位数相等。,2019/1/14,24/62,一、定量变量两组独立样本的秩和检验,例3：用两种药物杀灭钉螺，每批用200300只钉螺，用药后清点钉螺的死亡数，并计算死亡率(%)，问两种药物杀灭钉螺的效果有无差别？,2019/1/14,25/62,表4
10、 两种药物杀灭钉螺死亡率（%）的比较,2019/1/14,26/62,建立检验假设，确定检验水准 H0：两种药物杀灭钉螺死亡率的总体中位数相等 H1：两种药物杀灭钉螺死亡率的总体中位数不相等 a=0.05,2019/1/14,27/62,2.计算检验统计量T值 (1) 编秩 把两样本数据混合从小到大编秩，同组/不同组； (2) 求各组秩和 以样本例数小者为n1，其秩和为T1， (3) 确定检验统计量T值： 若n1n2，则T=T1； 若两样本例数相等，可任取一样本的秩和 本例n1=n2，可取T=T1=71.5或T=T2=33.5 。,2019/1/14,28/62,表4 两种药物杀灭钉螺死亡率（
11、%）的比较,2019/1/14,29/62,3.确定P值，作出推断 (1)查表法 查T界值表（成组设计用） n1=7，n2-n1=0， T=71.5或T=33.5 0.01P0.02 按照a=0.05水准，拒绝H0 ，可以认为两种药物杀灭钉螺的效果有差别。,2019/1/14,30/62,(2)正态近似法 n110或 n2-n110,2019/1/14,31/62,二、有序分类变量两组独立样本的秩和检验,例4：某医科大学营养教研室为了解居民体内核黄素营养状况，于某年夏冬两个季节收集成年居民口服5mg核黄素后4小时的负荷尿，测定体内核黄素含量，试比较该地居民夏冬两个季节体内核黄素含量有无差别？
12、表5 某地居民夏冬两个季节体内核黄素营养状况比较,2019/1/14,32/62,建立检验假设，确定检验水准 H0：夏冬两个季节居民体内核黄素含量的总体中位数相等 H1：夏冬两个季节居民体内核黄素含量的总体中位数不相等 a=0.05 2. 计算检验统计量T值 (1) 编秩：先确定各等级的合计、秩范围和平均秩 (2) 求各组秩和 (3) 确定检验统计量T值 T1=16.510+48.514+74.516=2036 T2=16.522+48.518+74.54=1534 T=T1=2036,2019/1/14,33/62,3.确定P值，作出推断 本例n1=40，超出T界值表可查范围，用正态近似法。
13、 两独立样本的非参检验： Mann-Whitney U检验(SPSS)和Wilcoxon秩和检验是等价的。,2019/1/14,34/62,第三节 多组独立样本比较的秩和检验,一、定量变量多组独立样本的秩和检验 二、有序变量多组独立样本的秩和检验 三、多个独立样本间的多重比较 Kruskal-Wallis H检验,一、定量变量多组独立样本的秩和检验,例5：某医院用3种不同方法治疗15例胰腺癌患者，每种方法各治疗5例。治疗后生存月数见表，问这3种方法对胰腺癌患者的疗效有无差别？ 表6 3种方法治疗胰腺癌患者的生存月数比较,2019/1/14,36/62,建立检验假设，确定检验水准 H0：3种方法
14、治疗后患者生存月数的中位数相等 H1：3种方法治疗后患者生存月数的中位数不全相等 a=0.05 2.计算检验统计量H值 (1) 编秩 三组统一编秩，相同数据同一组，可顺次编秩；在不同组，求平均秩次。 (2) 求各组秩和,2019/1/14,37/62,(3) 确定检验统计量H值： 若相持较多（如超过25%），应计算校正值,2019/1/14,38/62,3确定P值，做出推断 (1)当组数k =3，ni5，可查附表11中的界值，H的临界值等于5.78，P0.05，拒绝H0，接受H1，故可认为3种方法治疗后胰腺癌患者的生存月数有差别。 (2)当不满足条件(1)时，近似地服从自由度为k-1的c2分布
15、，可查附表8的c2界值表得到P值。 Kruskal-Wallis H检验的基本思想类似于单因素方差分析。相当于将秩次看做原始数据进行比较。,2019/1/14,39/62,二、有序变量多组独立样本的秩和检验,例6：某医院用3种方法治疗慢性喉炎，结果见表7。问这3种方法的疗效是否有差别？ 表7 3种方法治疗慢性喉炎的疗效比较,2019/1/14,40/62,1. 建立检验假设，确定检验水准 H0：3种方法疗效的总体分布位置相同 H1：3种方法疗效的总体分布位置不全相同 a=0.05 2. 计算检验统计量H值 (1) 编秩：先计算各等级的合计，再确定秩次范围及平均秩次。 (2) 求各组秩和 如本例
16、.,2019/1/14,41/62,表7 3种方法治疗慢性喉炎的疗效比较,2019/1/14,42/62,(3) 计算检验统计量H 由于相持较多，故需校正。,2019/1/14,43/62,3确定P值并做出推断结论 k=3，各组例数均大于5，近似服从n=k-1=2的2分布。查2界值表，得P0. 05。按a=0.05检验水准，拒绝H0，接受H1。可以认为3种方法治疗慢性喉炎的效果有差别。,2019/1/14,44/62,三、多个独立样本间的多重比较,2019/1/14,45/62,此外还有Nemenyi法检验 计算的是D或 c2 值,2019/1/14,46/62,补充： 一、随机区组设计的秩和检验,(一) M检验 Friedman M检验，用于推断随机区组设计的多个相关样本所来自的多个总体分布是否有差别。检验假设和备择假设和多个独立样本比较的Kruskal-Wallis检验相同。,2019/1/14,47/62,例7：8名受试对象在相同实验条件下分别接受4种不同频率声音的刺激，他们
