课件：卫生统计学第八版李晓松第十章基于秩的非参数检验.ppt

卫生统计学,贺 佳 第二军医大学 陶育纯 吉林大学,第十章 基于秩的非参数检验,目录,第一节：配对样本的比较,第二节：两组独立样本的比较,第三节：多组独立样本的比较,01,02,02,03,重点难点,符号秩和检验、两独立样本秩和检验和Kruskal-Wallis检验的基本思想 符号秩和检验、两独立样本秩和检验和Kruskal-Wallis检验的编秩规则 符号秩和检验、两独立样本秩和检验和Kruskal-Wallis检验中检验统计 量的计算方法,第一节 配对样本的比较,（一）单样本数据的符号秩和检验,第一节 配对样本的比较,例1 动物实验发现DON可导致家兔膝关节软骨和滑膜损伤，为研究大骨节病是否与粮食中DON含量有关，采集大骨节病高发地区面粉20份，测量面粉中DON含量，结果(g/g)如下：0, 0, 0, 0, 0, 12.4, 34.1, 69.0, 98.4, 129.5, 156.1, 163.5, 170.9, 177.6, 172.4, 180.3, 189.2, 192.2, 196.8, 205.3，中位数为142.8 g/g。根据前期研究发现，非大骨节病区面粉中DON含量平均水平（中位数）为18.9 g/g。是否可以认为大骨节病区与非大骨节病区面粉中DON含量不同？,正态分布拟合优度检验得2=369.011，P0.001，可认为该数据不服从正态分布。,基本思想 1.假设样本所对应的总体中位数与给定的总体中位数相同，H0：M1 = M0。 2.计算样本中所有数值与给定中位数的差值，根据所有差值绝对值进行编秩，得到正差值的秩和R+和负差值的秩和R-。 3.若H0成立，理论上，R+与R-的总体均数应相等，等于： ，总体标准差也应相等，等于： 。 4.若R+与R-相差悬殊，均远离M0，则有理由拒绝H0 。具体通过R+ 或R- 的抽样分布计算P值获得推断结论。,（一）单样本数据的符号秩和检验,第一节 配对样本的比较,检验步骤 (1) 建立检验假设，确定检验水准 H0：样本所对应总体的中位数等于已知总体中位数，M1 = M0 H1： 样本所对应总体的中位数不等于已知总体中位数，M1 M0 =0.05 (2) 求差值、编秩、求秩和 (3) 计算检验统计量 由表可知本例的检验统计量为R+=184或R-=26。,（一）单样本数据的符号秩和检验,第一节 配对样本的比较,检验步骤 (4) 确定P值，作出推断 样本量较大(n 50)则检验统计量近似服从正态分布，标准化变换获得Z值： P=0.004，拒绝H0，接受H1 ，样本与总体中位数的差异有统计学意义，可认为大骨节病病区面粉中DON含量高于非大骨节病地区。,（一）单样本数据的符号秩和检验,第一节 配对样本的比较,（一）单样本数据的符号秩和检验,第一节 配对样本的比较,例1数据的编秩结果,（二）配对样本数据的符号秩和检验,第一节 配对样本的比较,例2 采集10名正常成年男性志愿者的血清，分别用放射免疫法和酶联免疫法测量甲胎蛋白的含量(g/L)，结果见表。两种方法测量结果有无差异？,两种方法测量正常成年男性血清中甲胎蛋白含量(g/L)的结果,正态分布拟合优度检验得2=14.000，P=0.003，可认为该数据不服从正态分布。,基本思想 配对数据符号秩和检验的基本思想与单样本符号秩和检验是一致的。不同之处在于，配对数据中每个配对数值的差值可以看作是一个单独的样本，给定的总体中位数为0，即推断差值的单样本是否来自给定中位数为0的总体。其余部分则与单样本符号秩和检验并无差别。,（二）配对样本数据的符号秩和检验,第一节 配对样本的比较,检验步骤 (1) 建立检验假设，确定检验水准 H0：差值的总体中位数等于0， Md = 0 H1：差值的总体中位数不等于0， Md 0 =0.05 (2) 求差值、编秩、求秩和 首先计算每对数据的差值，并对差值进行编秩。分别计算正、负差值的秩和，得出 R+与R- ，如表所示。,（二）配对样本数据的符号秩和检验,第一节 配对样本的比较,检验步骤 (3) 计算检验统计量 本例的检验统计量为R+=21.5或R-=23.5。 (4) 确定P值，作出推断 本例P=0.880。在=0.05水准下不拒绝H0，差值的总体中位数与0的差异无统计学意义，尚不能认为放射免疫法与酶联免疫法测量正常成年男性血清甲胎蛋白的结果有差异。,（二）配对样本数据的符号秩和检验,第一节 配对样本的比较,（二）配对样本数据的符号秩和检验,第一节 配对样本的比较,两种方法测量正常成年男性血清中甲胎蛋白含量(g/L)的结果,第二节 两组独立样本的比较,（一）两组定量数据的比较,第二节 两组独立样本的比较,例3 在某小学随机采集12岁男童和女童各10名的头发样品，检测发样中钙(Ca)含量(g/g)，数据见下表。男童与女童头发中Ca含量有无差异？,12岁男、女童发样中Ca含量(g/g)比较,正态分布拟合优度检验，男童组2=2891.297，P0.001，女童组2=9290.446，P0.001。可以认为，男童组与女童组的数据均不服从正态分布。,基本思想 分别抽取样本量为n1和n2的两个样本，总例数N= n1+n2 。将全部数据统一编秩，取任意样本（如样本量为n1的样本）的秩和作为Wilcoxon秩和检验统计量W，假设两个总体分布相同(H0),则W的均数和标准差分别等于： 当W远离W，有理由拒绝 H0具体通过W的抽样分布计算P值获得推断结论。,（一）两组定量数据的比较,第二节 两组独立样本的比较,检验步骤 (1) 建立检验假设，确定检验水准 H0：男童与女童头发中Ca含量的总体分布相同 H1：男童与女童头发中Ca含量的总体分布不同 =0.05 (2) 编秩、求秩和 先将男童组与女童组发样中Ca含量的数值由小到大统一编秩，将两组秩分别相加得每组秩和。见表。,（一）两组定量数据的比较,第二节 两组独立样本的比较,(3) 计算检验统计量 本例W=77，Z=-2.117。 (4) 确定P值，作出推断 本例P=0.034，按=0.05 水准拒绝H0 ，接受H1 ，可以认为男童与女童的头发中Ca含量差异有统计学意义。男童组平均秩为77/10=7.7，女童组平均秩为133/10=13.3，可认为女童的头发中Ca含量高于男童。,（一）两组定量数据的比较,第二节 两组独立样本的比较,（一）两组定量数据的比较,第二节 两组独立样本的比较,12岁男童与女童发样中Ca含量(g/g)的比较,当两个样本例数较大时（n210或n2-n110 ），秩和检验统计量W将近似正态分布，可通过对W采取标准化变换获得值： 本例Z=2.078，P=0.032，结论同前。,（一）两组定量数据的比较,第二节 两组独立样本的比较,当数据包含相同秩时，秩和检验统计量W的精确分布会改变，同时，W的标准差必须进行如下调整： j 为出现相同秩的总次数， ti 为第i 次相同秩的个数,（一）两组定量数据的比较,第二节 两组独立样本的比较,例4 2003年10月在两个不同主食地区的中老年人手骨关节炎患病率随机抽样调查中，测得甲地区（主食大米）的175人和乙地区（主食面粉）的194人的手骨关节炎分值，数据见下表。现比较甲地区人群与乙地区人群的手骨关节炎分值是否具有统计学差异？,（一）两组定量数据的比较,第二节 两组独立样本的比较,甲地区与乙地区随机抽样人群的手骨关节炎分值比较,检验步骤 (1) 建立检验假设，确定检验水准 H0 ：甲地区人群与乙地区人群的手骨关节炎分值的总体分布相同 H1 ：甲地区人群与乙地区人群的手骨关节炎分值的总体分布不同 =0.05 (2) 编秩、求秩和（略）,（一）两组定量数据的比较,第二节 两组独立样本的比较,检验步骤 (3) 计算检验统计量 W = R1 =18642，Z = -13.447 (4) 确定P值，作出推断 P0.001，按 =0.05水准拒绝H0，接受H1，可以认为甲地区人群与乙地区人群的手骨关节炎分值的差异有统计学意义，可认为乙地区人群比甲地区人群的手骨关节炎严重。,（一）两组定量数据的比较,第二节 两组独立样本的比较,例5 2007年7月至11月对西藏拉萨市和山南地区大骨节病进行流行病学调查，获得了两个地区受检人群大骨节病不同临床分度的数据，见下表。现比较两个地区的大骨节病临床分度是否具有统计学差异？,（二）两组等级变量的比较,第二节 两组独立样本的比较,2007年拉萨市和山南地区大骨节病临床分度的比较,检验步骤 (1) 建立检验假设，确定检验水准 H0 ：两个地区的大骨节病临床分度的总体分布相同 H1 ：两个地区的大骨节病临床分度的总体分布不同 =0.05 (2) 编秩、求秩和（略） (3) 计算检验统计量 W = R1 =10608.5，Z = 0.145 （4）确定P值，作出推断 P0.05，按 =0.05水准不拒绝H0，尚不能认为拉萨市和山南地区的大骨节病临床分度的差异有统计学意义。,（二）两组等级变量的比较,第二节 两组独立样本的比较,第三节 多组独立样本的比较,（一）多组定量数据的比较,第三节 多组独立样本的比较,例6 为了解不同剂量的DON对新西兰家兔膝关节软骨和滑膜的损伤情况，将15只新西兰家兔按体重随机分为对照组、低剂量组和高剂量组，分别注射无菌生理盐水、0.05g/g和0.10g/g剂量DON毒素进行实验处理，实验期满后测定关节冲洗液中肿瘤坏死因子(TNF-)的水平(g/L)，获得数据见下表10-6。现比较3组家兔关节冲洗液TNF-测定结果是否具有统计学差异？,3组家兔关节冲洗液TNF-（g/L）测定结果,基本思想 Kruskal-Wallis检验的基本思想就是用所有观测值的秩代替原始观测值进行单因素方差分析。检验统计量为： 其实质是用秩计算组间变异，当组间变异较大时，有理由认为组间存在差异。具体通过H分布计算P值获得推断结论。 若有相同秩出现，则需要校正：,（一）多组定量数据的比较,第三节 多组独立样本的比较,检验步骤 (1) 建立检验假设，确定检验水准 H0：3组家兔关节冲洗液TNF-测定结果的总体分布相同 H1：3组家兔关节冲洗液TNF-测定结果的总体分布不全相同 =0.05 (2) 编秩、求秩和 先将三组TNF-的数值由小到大统一编秩，将各组秩分别相加得每组秩和Ri，见下表。,（一）多组定量数据的比较,第三节 多组独立样本的比较,检验步骤 (3) 计算检验统计量 本例H=11.18。 (4) 确定P值，作出推断 P0.005，按=0.05 水准拒绝H0，接受H1，可认为三组家兔关节冲洗液TNF-测定结果的差异有统计学意义。,（一）多组定量数据的比较,第三节 多组独立样本的比较,（一）多组定量数据的比较,第三节 多组独立样本的比较,3组家兔关节冲洗液TNF-（g/L）测定结果,（二）多组等级变量的比较,第三节 多组独立样本的比较,例7 2007年7月至11月西藏拉萨市和山南地区大骨节病流行病学调查获得了两个地区受检人群不同年龄段（岁）的大骨节病不同临床分度的数据，见下表。现比较不同年龄段的大骨节病临床分度是否具有统计学差异？,不同年龄段（岁）的大骨节病临床分度的比较,检验步骤 (1) 建立检验假设，确定检验水准 H0 ：不同年龄段的大骨节病临床分度的总体分布相同 H1 ：不同年龄段的大骨节病临床分度的总体分布不全相同 =0.05 (2) 编秩、求秩和（略） (3) 计算检验统计量 Hc = 1.330 (4) 确定P值，作出推断 本例P0.05，按 =0.05水准不拒绝H0，尚不能认为不同年龄段的大骨节病临床分度的差异有统计学意义。,（二）多组等级变量的比较,第三节 多组独立样本的比较,小结,1.数据符合参数检验条件，最好仍用参数检验，若采用非参数检验，容易导致检验效能下降。 2.假设样本所对应总体中位数与已知总体中位数相同(H0)，计算样本中所有数值与已知标准值的差值，根据所有差值的绝对值进行编秩，并分别计算正、负差值的秩和。根据正差值秩和或负差值秩和与总体平均秩的距离判断是否拒绝 。 3.配对数据符号秩和检验的基本思想与单样本符号秩和检验的相似，配对数据中每个配对数值的差值可以看作是一个