医学统计学精品教学汕头大学第十章非参数检验

1、2021-11-141第十章第十章 非参数秩和检验非参数秩和检验吴吴 库库 生生汕头大学医学院预防医学教研室2021-11-14秩和检验2参数统计与非参数统计1、参数统计（参数统计（parametric statistics） 样本所来自的总体分布具有某个已知样本所来自的总体分布具有某个已知的函数形式（如正态分布），而其中的函数形式（如正态分布），而其中有有的的参数参数是未知的，统计分析的目的是对是未知的，统计分析的目的是对这些未知的参数进行估计或检验。这些未知的参数进行估计或检验。2021-11-14秩和检验3但实际上有些医学研究资料并但实际上有些医学研究资料并不符合上述不符合上述条件条件，

2、即使用变量变换的方法也难以达，即使用变量变换的方法也难以达到应用参数方法的要求。到应用参数方法的要求。这就需要一种这就需要一种不依赖于总体分布不依赖于总体分布的具体形的具体形式的统计方法来解决此类问题。式的统计方法来解决此类问题。由于这类方法不受总体参数的影响，故称由于这类方法不受总体参数的影响，故称非参数检验法。非参数检验法。它检验的是它检验的是分布分布，而非，而非参数。参数。2021-11-14秩和检验42、非参数统计（、非参数统计（ non-parametric statistics ） 不考虑研究对象总体分布的具体形式，也不考虑研究对象总体分布的具体形式，也不对总体参数进行统计推断，而

3、是通过检验不对总体参数进行统计推断，而是通过检验样本所代表的样本所代表的总体分布形式是否一致总体分布形式是否一致来得出来得出统计结论。由于这类方法不受总体参数的限统计结论。由于这类方法不受总体参数的限制，故称非参数统计法。有时也称为任意分制，故称非参数统计法。有时也称为任意分布检验布检验(distribution-free statistics)(distribution-free statistics)。2021-11-14秩和检验53、非参数统计的优点非参数统计的优点 适用范围广，不论样本所代表的总体分适用范围广，不论样本所代表的总体分布形式如何，甚至是未知的，都可适用；布形式如何，甚至是

4、未知的，都可适用； 资料可用资料可用“等级等级”、“符号符号”表示，收表示，收集方便；集方便； 多数非参数检验多数非参数检验方法简便，易于理解和方法简便，易于理解和掌握。掌握。2021-11-14秩和检验64、非参数统计的缺点非参数统计的缺点 对符合用参数检验的资料，如用非参数对符合用参数检验的资料，如用非参数统计，会丢失信息，导致检验效率下降，统计，会丢失信息，导致检验效率下降，犯第犯第iiii类错误的可能性比参数检验大类错误的可能性比参数检验大, ,即即值值； 非参数检验的有些问题的计算仍显繁冗。非参数检验的有些问题的计算仍显繁冗。2021-11-14秩和检验7何时使用非参数检验何时使用非

5、参数检验?1.1.等级资料（半定量资料）等级资料（半定量资料） 当指标只能用严重程度当指标只能用严重程度( (轻、中、重轻、中、重) )，优劣，优劣等级（好、中、差），治疗转归（治愈、好转、等级（好、中、差），治疗转归（治愈、好转、无效）等形式表达时，原始资料并非定量的，这无效）等形式表达时，原始资料并非定量的，这些数据比些数据比“定量定量”粗，而比一般的粗，而比一般的“定性定性”细，细，是相继的若干阶梯。但毗邻的阶梯之间并非等距是相继的若干阶梯。但毗邻的阶梯之间并非等距离，亦不能度量。离，亦不能度量。2.2.偏态分布资料偏态分布资料2021-11-14秩和检验8何时使用非参数检验何时使用非参

6、数检验?3.3.个别数据偏离过大资料个别数据偏离过大资料 这里指的是随机的偏离而不属于这里指的是随机的偏离而不属于“过失误差过失误差”所致。还有一种情况就是数据的某一端无确定数所致。还有一种情况就是数据的某一端无确定数值，如值，如“50mg”等等，只给一等等，只给一个上限或下限，而没有具体数值。个上限或下限，而没有具体数值。4.4.各组离散程度相差悬殊各组离散程度相差悬殊( (方差不齐方差不齐) )的资料的资料 如经变量变换达到方差齐性，则可用参数检如经变量变换达到方差齐性，则可用参数检验；否则的话就要用非参数法处理。验；否则的话就要用非参数法处理。2021-11-14秩和检验9何时使用非参数

7、检验何时使用非参数检验?5.5.分布形态不明的资料分布形态不明的资料 小样本，但又不趋向于正态分布；不小样本，但又不趋向于正态分布；不稳定，例数少，分布不明确的资料。稳定，例数少，分布不明确的资料。 符合参数检验条件 首选参数检验 不符合条件 非参数检验2021-11-14秩和检验10秩和检验秩和检验秩次秩次（rankrank）：就是将观察值按顺序由小）：就是将观察值按顺序由小到大排列到大排列( (排名次排名次) )，并用序号（如，并用序号（如1, 2, 1, 2, 3 3 ）代替变量值本身，秩次即通常意）代替变量值本身，秩次即通常意义上的序号。义上的序号。秩和：秩和：即秩次的和。即秩次的和。

8、秩转换秩转换：将数值变量转换成秩次，再计算：将数值变量转换成秩次，再计算检验统计量的过程。检验统计量的过程。秩和检验秩和检验：就是通过秩次的排列求秩和进：就是通过秩次的排列求秩和进行假设检验的方法。行假设检验的方法。2021-11-14秩和检验11contents第一节第一节 配对资料的符号秩和检验（配对资料的符号秩和检验（wilcoxonwilcoxon配对法）配对法）第二节第二节 两独立样本比较的秩和检验（两独立样本比较的秩和检验（wilcoxonwilcoxon两样本两样本法）法）第三节第三节 完全随机设计多个样本比较的秩和检验完全随机设计多个样本比较的秩和检验（kruskal-krus

9、kal-allisallis检验）检验）2021-11-14秩和检验12第一节第一节 两配对样本差值的符号秩和检验两配对样本差值的符号秩和检验(wilcoxon signed rank test) 符号秩和检验由符号秩和检验由wilcoxonwilcoxon于于19451945年提年提出，作为配对出，作为配对t t检验的替代方法。检验配检验的替代方法。检验配对资料的差值是否来自对资料的差值是否来自中位数中位数为为0 0的总体。的总体。 在数据满足配对在数据满足配对t t检验的要求时，符检验的要求时，符号秩和检验的功效是配对号秩和检验的功效是配对t t检验效能的检验效能的95%95%左右。左右。

10、2021-11-14秩和检验13符号秩和检验的符号秩和检验的基本思想基本思想：在在h h0 0成立的条件下（两配对处理效应相成立的条件下（两配对处理效应相同同 ），两配对样本的差数的正负是随），两配对样本的差数的正负是随机的，则正差数的秩和与负差数的秩和机的，则正差数的秩和与负差数的秩和应该相差不会太大，应该相差不会太大，即总体中位数为即总体中位数为0 0；反之，若两秩和相差太悬殊，则认为反之，若两秩和相差太悬殊，则认为h h0 0成立的可能性受到怀疑，从而拒绝成立的可能性受到怀疑，从而拒绝h h0 0，接受接受h h1 1。2021-11-14秩和检验14例例1 1：分别用氰化高铁法与光电比

11、色法测：分别用氰化高铁法与光电比色法测得得1010个血样的血红蛋白含量，问两法测个血样的血红蛋白含量，问两法测得结果有无差别？得结果有无差别？表1 两种方法测定血清血红蛋白含量血样号 氰化高铁法 光电比色法 1 15.0 12.5 2 12.8 12.9 3 13.5 13.0 4 11.3 12.8 5 12.3 15.4 6 12.1 12.9 7 13.7 13.0 8 12.9 14.1 9 12.9 12.1配对差值经检验不符合正态分布，不满足t检验条件2021-11-14秩和检验16符号秩和检验的基本步骤：1.1.建立检验假设，确定检验水准建立检验假设，确定检验水准 h h0 0

12、：md=0md=0（两处理效应相同）（两处理效应相同） h h1 1 ：mdmd0 0（两处理效应不相同）（两处理效应不相同） =0.05=0.052021-11-14秩和检验172. 求两组数据的差数求两组数据的差数d d，对差数的绝对，对差数的绝对值编秩号值编秩号 差数为差数为0 0的数据忽略不计，对子数减的数据忽略不计，对子数减1 1； 余下的余下的n-in-i个差数个差数按绝对值按绝对值由小到大排由小到大排秩号，但排好后的秩号秩号，但排好后的秩号要保持原差数的正负要保持原差数的正负号号； 差数绝对值相等时，则取其平均秩次。差数绝对值相等时，则取其平均秩次。表1 两种方法测定血清血红蛋白

13、含量血样号 氰化高 光电比 差值 负秩 正秩 铁法 色法 (1) (2) (3) (4)=(2)-(3) (5) (6) 1 15.0 12.5 2.5 9 2 12.8 12.9 -0.1 -1.5 3 13.5 13.0 0.5 3 4 11.3 12.8 -1.5 -8 5 12.3 15.4 -3.1 -10 6 12.1 12.9 -0.8 -5.5 7 13.7 13.0 0.7 4 8 12.9 14.1 -1.2 -7 9 12.9 12.1 0.8 5.5对差对差数的数的绝对绝对值编值编秩秩2021-11-14秩和检验193. 分别求正、负差数的秩和，以绝对值较分别求正、负差

14、数的秩和，以绝对值较小者为小者为t t 验算：验算：t t+ +t+t- -=n(n+1)/2=n(n+1)/24. 4. 确定确定p p值和作出推断结论值和作出推断结论小样本时（小样本时（n n50, 50, n n为为d d00的对子的对子数数）, ,通过查通过查p226p226附表附表8(t8(t界值表界值表) )得到得到p p值，值，原则原则是：若统计量是：若统计量t t处于某处于某t t界值范界值范围之内，其围之内，其p p值大于或等于相应概率；反值大于或等于相应概率；反之，则小于相应概率。当之，则小于相应概率。当t t值恰好等于附值恰好等于附表中的界值时，其确切概率常小于相应表中的

15、界值时，其确切概率常小于相应的概率水平。的概率水平。2021-11-14秩和检验20本例中：本例中：t=23 对子数对子数n=10，查，查p226附表附表8得得2310 45 10 45 故故 p0.10下结论：在下结论：在=0.05=0.05水准上不拒绝水准上不拒绝h h0 0 ，故，故不能认为两法测得血红蛋白含量有差别。不能认为两法测得血红蛋白含量有差别。2021-11-14秩和检验21 大样本时（大样本时（ 对子数对子数n50）通过公式进）通过公式进行行u u转化，采用正态近似检验：转化，采用正态近似检验： （式（式10-110-1） 式中式中0.50.5是连续性校正数是连续性校正数,

16、,因为因为t t值是值是不连续的不连续的, ,而而u u分布是连续的。分布是连续的。 24/ ) 12)(1(5 . 0|4/ ) 1(|nnnnntu2021-11-14秩和检验22当存在同一秩较多时，按（式当存在同一秩较多时，按（式10-110-1）计算）计算的的u u值要比实际的值要比实际的u u值小，应采用（公式值小，应采用（公式10-210-2）对）对u u加以校正加以校正: :其中，其中，t tj j是相同秩号的个数。是相同秩号的个数。48)(24) 12)(1(5 . 0|4/ ) 1(|3jjttnnnnntu2021-11-14秩和检验23符号秩和检验的spss实现 spss

17、(statistical package for social science,社会科学统计软件包)一直是国际上最流行且最具权威性的统计分析软件之一。spss最显著的特点是菜单和对话框操作方式，无需编程，且图形功能强大，易于掌握和操作。2021-11-14秩和检验242021-11-14秩和检验25符号秩和检验的spss实现2021-11-14秩和检验26符号秩和检验的spss实现2021-11-14秩和检验272021-11-14秩和检验28符号秩和检验的spss实现2021-11-14秩和检验29符号秩和检验的spss实现2021-11-14秩和检验30ranksranks5a4.6023

18、.005b6.4032.000c10negative rankspositive rankstiestotal光电比色法 - 氰化高铁法nmean ranksum of ranks光电比色法 氰化高铁法b. 光电比色法 = 氰化高铁法c. test statisticstest statisticsb b-.459a.646zasymp. sig. (2-tailed)光电比色法 -氰化高铁法based on negative ranks.a. wilcoxon signed ranks testb. 观测值x a 7 14 22 36 40 48 63 98 41.0029.81b 3 5 6 10 17 18 20 39 14.7511.73 phhh, ,则则p;pp；ii.ii.大样本情况：若大样本情况：若k k 3 3或或n ni i 5 5时，理论时，理论上，上，h h近似服从自由度为近似服从自由度为k1k1的的 分布，可查附表分布，可查附表7 7（ 界值表）得界值表）得出出p p值，最后按值，最后按p p值作出推断结论。值作出推断结论。222021-11-14秩和检验73 本例处理数本例处理数k=3k=3，各组样本含量均为，各组样本含量均为1010，超出了附表超出了