第九章秩和检验课件.ppt

《第九章秩和检验课件.ppt》由会员分享，可在线阅读，更多相关《第九章秩和检验课件.ppt（80页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、Statistics 医学统计学医学统计学 (Medical Statistics)(Medical Statistics)院医室楼院医室楼院医室楼院医室楼 4 4楼楼楼楼 Tel:Tel:0759-23885900759-2388590E-mail:E-mail:gdmyfyxgdmyfyx cn医学统计学教研室医学统计学教研室医学统计学教研室医学统计学教研室 DepartmentDepartment of ofStatistical MedicineStatistical MedicineStatistics主讲：孔丹莉主讲：孔丹莉秩和检验秩和检验秩和检验秩和检验2/45rank sum

2、test第十一章第十一章非参数检验非参数检验简单复习一下前面所学内容简单复习一下前面所学内容简单复习一下前面所学内容简单复习一下前面所学内容统计描述统计描述统计推断统计推断适用条件适用条件计量资料计量资料(定量资料定量资料)频数分布频数分布集中趋势集中趋势离散趋势离散趋势统计图表统计图表抽样误差、抽样误差、标准误标准误tuF检验检验秩和检验秩和检验正态分布正态分布方差齐性方差齐性计数资料计数资料(定性资料定性资料)(分类资料分类资料)(无序分类无序分类)频数分布频数分布相对数及相对数及其标准化其标准化统计图表统计图表u 2检验检验分布类型分布类型样本例数样本例数理论频数理论频数等级资料等级资料

3、(有序分类有序分类)秩和检验秩和检验 3/45 复习一下正态分布的判断方法复习一下正态分布的判断方法复习一下正态分布的判断方法复习一下正态分布的判断方法 Statistics4/45利用频数表或频数图进行判断利用频数表或频数图进行判断根据专业知识判断：疾病的潜伏期、住院天数和根据专业知识判断：疾病的潜伏期、住院天数和 临床生化指标大多为偏态临床生化指标大多为偏态正态性检验正态性检验正态分布的经验判断正态分布的经验判断 若若 ，则有理由怀疑资料呈偏态分布，则有理由怀疑资料呈偏态分布 若若 ，可认为资料呈偏态分布，可认为资料呈偏态分布参数统计（参数统计（parametricstatistics）参

4、数统计要求样本来自某种特定分布的参数统计要求样本来自某种特定分布的总体，而该分布中的某些参数未知，统总体，而该分布中的某些参数未知，统计分析的目的就是对这些未知参数进行计分析的目的就是对这些未知参数进行估计或推断，如估计或推断，如t检验、方差分析都要求检验、方差分析都要求原始数据来自于正态分布的总体，且方原始数据来自于正态分布的总体，且方差齐性。差齐性。非参数统计非参数统计（non-parametricstatistics）有许多资料不符合参数统计的要求，不有许多资料不符合参数统计的要求，不能用参数统计的方法进行检验，而需要能用参数统计的方法进行检验，而需要一种不依赖于总体分布类型，也不对总一

5、种不依赖于总体分布类型，也不对总体参数进行统计推断的假设检验，称为体参数进行统计推断的假设检验，称为非参数检验。非参数检验。非参数分析方法的优缺点非参数分析方法的优缺点优点：优点：不受总体分布条件的限制，适用范围广；某不受总体分布条件的限制，适用范围广；某些不便准确测定，只能以严重程度，好坏优劣，些不便准确测定，只能以严重程度，好坏优劣，次第先后等作记录的资料也可应用。次第先后等作记录的资料也可应用。对资料没有特殊要求，例如对资料没有特殊要求，例如 不受分布的影响（偏态、分布不明的资料）不受分布的影响（偏态、分布不明的资料）不受方差齐性的限制不受方差齐性的限制 不受变量类型的影响不受变量类型的

6、影响 不受样本量的影响不受样本量的影响非参数分析方法的优缺点非参数分析方法的优缺点缺点：缺点：适用于参数检验的资料若用非参数检验会造适用于参数检验的资料若用非参数检验会造成信息损失，导致检验效率较低。即当成信息损失，导致检验效率较低。即当H0H0不真，不真，非参数检验可能不如参数检验能较灵敏地拒绝非参数检验可能不如参数检验能较灵敏地拒绝H0 H0，犯第二类错误的概率比参数检验大犯第二类错误的概率比参数检验大。若资料适用于参数方法（如正态分布，方差若资料适用于参数方法（如正态分布，方差齐性等，级就尽可能不用非参数方法）齐性等，级就尽可能不用非参数方法）第一类错误与第二类错误第一类错误与第二类错误

7、 客观实际客观实际拒绝拒绝H0不拒绝不拒绝H0H0成立成立 第一类错误第一类错误,（误诊率，假阳误诊率，假阳性率性率）推断正确推断正确H0不成立推断正确推断正确第二类错误第二类错误（漏诊率，假（漏诊率，假阴性率）阴性率）优点点：对资料的信息利用充分对资料的信息利用充分统计分析的效率高统计分析的效率高 缺点：缺点：对资料的要求高对资料的要求高适用范围有限适用范围有限Statistics7/45参数分析方法的优缺点参数分析方法的优缺点 常用的非参数检验方法常用的非参数检验方法秩和检验秩和检验Ridit分析分析卡方检验卡方检验K-S法法秩相关分析（等级相关分析）秩相关分析（等级相关分析）cpd法法(

8、积差交叉法）积差交叉法）游程检验游程检验 Statistics9/45秩秩和和检检验验秩和检验（秩和检验（ranksumtest）是最常用的的非）是最常用的的非参数检验方法，也称参数检验方法，也称秩转换秩转换（ranktransform-ation），该方法在非参数检验中占有重要地位。），该方法在非参数检验中占有重要地位。原理：原理：是首先将原始数据从小到大，或是首先将原始数据从小到大，或等级从弱到等级从弱到强转换成秩后，再对基于强转换成秩后，再对基于秩次秩次的统计量（如秩的统计量（如秩和），进行检验，做出统计推断。故又称基于秩和），进行检验，做出统计推断。故又称基于秩次的非参数检验。次的非参

9、数检验。秩号：秩号：将各原始数据从小到大排列，分别给每个将各原始数据从小到大排列，分别给每个数据一个顺序号，也就是秩号数据一个顺序号，也就是秩号(rank)。如：967.513 秩号：3124 秩和：秩和：用秩次代替原始数据，在计算各组秩次之用秩次代替原始数据，在计算各组秩次之和；和；秩和检验：秩和检验：基于秩和提供的信息，对不同总基于秩和提供的信息，对不同总体的平均水平进行假设检验。体的平均水平进行假设检验。秩和检验常用方法秩和检验常用方法 l 配对设计资料的符号秩和检验配对设计资料的符号秩和检验l两独立样本差别的秩和检验两独立样本差别的秩和检验l 完全随机设计多组差别的秩和检验完全随机设计

10、多组差别的秩和检验l 多个样本间两两比较的秩和检验多个样本间两两比较的秩和检验第一节第一节 WilcoxonWilcoxon 符号秩和检验符号秩和检验一、配对设计的两样本比较一、配对设计的两样本比较一、配对设计的两样本比较一、配对设计的两样本比较 Statistics10/45 什么叫配对设计？什么叫配对设计？是是将将实验对象按一定条件配成对子，如将年龄、实验对象按一定条件配成对子，如将年龄、体重一致的动物配成对子，再随机分配每对中的两体重一致的动物配成对子，再随机分配每对中的两个对象接受不同的处理方式；或同一样品用两种方个对象接受不同的处理方式；或同一样品用两种方法（或仪器等）检验的结果；或

11、同一受试对象两个法（或仪器等）检验的结果；或同一受试对象两个部位的数据。配对设计可增强组间的均衡性，提高部位的数据。配对设计可增强组间的均衡性，提高实验效率。实验效率。WilcoxonWilcoxon符号秩和检验符号秩和检验（WilcoxonWilcoxon 配对法，或配对法，或WilcoxonWilcoxon signed rank signed rank testtest）是推断其差值是否来自中位数为零的总体是推断其差值是否来自中位数为零的总体的方法，可用于计量配对设计资料差值的比的方法，可用于计量配对设计资料差值的比较和单一样本与总体中位数的比较较和单一样本与总体中位数的比较 Stati

12、stics11/45请看请看课本课本P152P152例题例题11-111-1为观察血浆置换法治疗出凝血功能异常为观察血浆置换法治疗出凝血功能异常的临床疗效的临床疗效,某医师治疗了某医师治疗了11例出凝血功例出凝血功能异常患者，置换前后各患者的凝血酶原能异常患者，置换前后各患者的凝血酶原时间见表时间见表11-1。该医师采用两样本均数比。该医师采用两样本均数比较的检验，结果，差异有统计学意义，由较的检验，结果，差异有统计学意义，由此认为血浆置换治疗前后凝血酶原时间有此认为血浆置换治疗前后凝血酶原时间有差别。差别。Statistics11/45问题问题l（1）该资料为何种类型资料？）该资料为何种类型

13、资料？l（2）该研究属于何种设计方案？）该研究属于何种设计方案？l（3）该医师所选用的统计分析方法是否正）该医师所选用的统计分析方法是否正确？为什么？确？为什么？l（4）该资料应采用何种统计方法进行分）该资料应采用何种统计方法进行分析？析？其步骤如何？其步骤如何？Statistics11/45Statistics11/45（1）该资料为计量资料。）该资料为计量资料。（2）该研究属自身配对设计方案。）该研究属自身配对设计方案。（3）该该医医师师所所选选用的用的统计统计分析方法不正确。分析方法不正确。原因在于：原因在于：自身配自身配对设计资对设计资料一般料一般选选用配用配对对检验，而该医师采用完全

14、随检验，而该医师采用完全随机设计两样本均数比较的机设计两样本均数比较的t检验，统计检验，统计方法与设计类型不符；方法与设计类型不符；该医师未考该医师未考虑资料是否满足虑资料是否满足t检验的应用条件。检验的应用条件。建立检验假设，确定检验水准建立检验假设，确定检验水准 H H0 0：差值的总体中位数为零，即：差值的总体中位数为零，即M Md d=0=0 H H1 1：M Md d0 0 =0.05=0.05 计算检验统计量计算检验统计量v 求各对测量值的差值求各对测量值的差值v 编秩：编秩原则编秩：编秩原则 12/45按照差值的绝对值从小到大编秩按照差值的绝对值从小到大编秩v差值为差值为 0 0

15、 者不参加编秩者不参加编秩v绝对值相等，符号相同时顺次编绝对值相等，符号相同时顺次编秩秩v绝对值相等，符号相反时取平均秩次绝对值相等，符号相反时取平均秩次 v 求秩和、确定检验统计量求秩和、确定检验统计量T值值分别求出正、负秩和，正秩和以分别求出正、负秩和，正秩和以T+T+表示，负表示，负秩和以秩和以T-T-表示。任取正秩和或负秩和为统计表示。任取正秩和或负秩和为统计量量T值。其中值。其中 T+T-=n(n+1)/2T+T-=n(n+1)/2，n n 为有效对子数为有效对子数 确定确定P P 值、作出统计推断结论值、作出统计推断结论v查表法：查表法：用于有效对子数用于有效对子数 n50 (P3

16、45)n50 (P345)若若T值在上、下界值范围内，则值在上、下界值范围内，则P0.05若若T值在上、下界值上或范围外值在上、下界值上或范围外，则，则P0.05 14/45StatisticsStatistics15/45v正态近似法：正态近似法：用于有效对子数用于有效对子数 n n5050 1)1)当差值绝对值相同的个数当差值绝对值相同的个数n25%n25%时时 差值绝对值相同差值绝对值相同即指相同秩。指差值的绝对值即指相同秩。指差值的绝对值相等者取平均秩。也就是秩次相等的个数，可用相等者取平均秩。也就是秩次相等的个数，可用tjtj表示。如表示。如t1t1指第一个出现相同秩次的个数。本例中

17、指第一个出现相同秩次的个数。本例中有两个秩次为有两个秩次为2.5(2.5(指绝对值指绝对值)，则，则t1=2;t1=2;有两个秩有两个秩次为次为6.5,6.5,则则t2=2t2=2。当。当tjn25%tjn25%时，公式如下时，公式如下 2)2)当差值绝对值相同的个数当差值绝对值相同的个数n25%n25%时时上式中，上式中，tjtj 为为第第 j j 个差值绝对值相同的个数个差值绝对值相同的个数Statistics16/45=(2=(23 3-2)+(2-2)+(23 3-2)=12-2)=12如果还有相同的，再相应加上，依此类推如果还有相同的，再相应加上，依此类推（2 2）单一样本与总体中位

18、数比较）单一样本与总体中位数比较 （见课本（见课本P154P154例题例题11-211-2）已已知知某某地地正正常常人人尿尿汞汞含含量量的的中中位位数数为为2.50g/L，某某医医师师从从该该地地某某厂厂从从事事土土法法炼炼金金（汞汞齐齐法法）的的汞汞作作业业工工人人中中随随机机抽抽取取10名名工工人人，测测得得尿尿汞汞含含量量(g/L)为为11.01，2.13，2.56，2.79，12.95，3.12，3.56，4.37，5.13，18.90。采采用用单单样样本本检检验验对对资资料料进进行行分分析析，得得，按按水水准准，拒拒绝绝H0，接接受受H1，差差异异有有统统计计学学意意义义，可可以以认

19、认为为该该厂厂工工人人尿尿汞含量与正常人尿汞含量不同。汞含量与正常人尿汞含量不同。Statistics17/45Statistics17/45建立检验假设建立检验假设H0：工人的尿汞含量总体中位数等于工人的尿汞含量总体中位数等于2.15，M=2.50H1：M2.50单侧单侧=0.05计算检验统计量计算检验统计量v求差值：求差值：这里的差值是各观察值与已知这里的差值是各观察值与已知总体中位数（总体中位数（2.50）之差）之差v编秩：编秩原则同上编秩：编秩原则同上v求秩和、确定检验统计量求秩和、确定检验统计量T值：值：取取T+和和T-的任意一个作为检验统计量的任意一个作为检验统计量v确定确定P值、

20、作推断结论值、作推断结论 本本例例n=10,T=4.5，查查T界界值值表表（配配对对比比较较的的符符号号秩秩和和检检验验用用，附附表表10），得得P 20.05,3,故故P0.05，按按水水准准拒拒绝绝H0，接接受受H1，差差异异有有统统计计学学意意义义，可可认认为为斯斯康康杜杜尼尼和和利利多多卡卡因因麻麻醉醉效效果果不不相相同同，斯斯康康杜杜尼尼的的麻麻醉醉效果好于利多卡因。效果好于利多卡因。Statistics28/451建立检验假设，确定检验水准建立检验假设，确定检验水准H0：两种药物麻醉效果的总体分布相同两种药物麻醉效果的总体分布相同H1：两种两种药物麻醉效果的药物麻醉效果的总体分布不

21、同总体分布不同=0.052计算检验统计量：详见课本计算检验统计量：详见课本P158例题例题（代入公式如下）（代入公式如下）uc=4.79555/0.8881=5.40Statistics29/45最最后后与与1.96比比较较作作出出结结论论，P0.05,按按=0.05水水准准拒拒绝绝H0，接接受受H1，差差异异有有统统计计学学意意义义，可可以以认认为为斯斯康康杜杜尼尼和和利利多多卡卡因因应应用用于于深深牙牙周周袋袋刮刮治治术术的的局局部部麻麻醉醉效效果果的的总总体体分分布布不不同同，斯斯康康杜杜尼尼的的效效果果好好于于利多卡因。利多卡因。l 两样本比较的Wilcoxon秩和检验用于完全随机（成

22、组）设计两样本资料的比较，是针对呈偏态分布或方差不齐不满足t检验条件的计量资料、等级资料和其它不能精确测量的资料的比较。31/45第三节第三节 成组设计多个样本比较的秩和检验成组设计多个样本比较的秩和检验一、一、KruskalKruskal-Wallis-Wallis 法，又称法，又称H H 检验法检验法 该法是由该法是由KruskalKruskal和和WallisWallis在在WilcoxonWilcoxon秩和检秩和检验的基础上扩展的方法，又称为验的基础上扩展的方法，又称为K-WK-W检验或检验或H H检验。检验。本法利用多个样本的秩和来推断各样本分别代表本法利用多个样本的秩和来推断各样

23、本分别代表的总体分布有无差别。的总体分布有无差别。（1）原始数据的多个样本比较）原始数据的多个样本比较见见课本课本P159P159例题例题11-511-5Statisticsl例例11-511-5为探讨血管紧张素转化酶抑制剂开博通抗动脉粥为探讨血管紧张素转化酶抑制剂开博通抗动脉粥样硬化的机制，某医师将样硬化的机制，某医师将15只兔随机分为三组，即正常饮只兔随机分为三组，即正常饮食组、高脂饮食组及高脂饮食加开博通组，喂养食组、高脂饮食组及高脂饮食加开博通组，喂养14周后测周后测量兔血清总胆固醇量兔血清总胆固醇(TC)水平，结果如下：水平，结果如下：l正常饮食组：正常饮食组：0.860.781.3

24、80.881.70l高脂饮食组：高脂饮食组：40.3541.7944.0040.7227.05l高脂饮食加开博通组：高脂饮食加开博通组：41.2033.0440.6032.0540.4932/45Statistics该医师应用方差分析对资料进行分析，F=104.41，P0.05，差异有统计学意义，采用t检验进一步作两两比较，得高脂饮食组及高脂饮食加开博通组均与正常组差异有统计学意义；但高脂饮食组与高脂饮食加开博通组两组间差异无统计学意义，故认为高脂饮食组及高脂饮食加开博通组的兔血清总胆固醇显著高于正常组；开博通对血脂无影响。方法步骤：方法步骤：建立检验假设建立检验假设H0：三个总体的分布位置相

25、同三个总体的分布位置相同H1：三个总体的分布位置不全相同三个总体的分布位置不全相同=0.05lKruskal-Wallis法法计算检验统计量计算检验统计量v编秩：将各组数据统一由小到大排序并编秩，编秩：将各组数据统一由小到大排序并编秩，如遇有相等数值则取平均秩次。如遇有相等数值则取平均秩次。v求秩和：分别将各组秩次相加，得各组秩和求秩和：分别将各组秩次相加，得各组秩和v计算检验统计量计算检验统计量H值值当相同秩次的个数当相同秩次的个数n25%时，用上式时，用上式Ti为各组的秩为各组的秩和，和，ni为各组对应的例数。为各组对应的例数。33/45Statistics 2 2）当）当相同秩次的个数相

26、同秩次的个数n25%时时 tj为第为第j个相同秩次的个数个相同秩次的个数 34/45Statistics确定确定P值、作出统计推断结论值、作出统计推断结论v查表法：组数查表法：组数k=3，每组例数每组例数ni 5时，时，可查可查H界值表界值表P347近似法近似法若组数若组数k=3，每组例数每组例数ni5；或者或者k3或者或者ni5时，则时，则H分布近似分布近似服从自由度为服从自由度为k-1的的分布分布,可查可查界值表确定界值表确定P值值（2）频数表资料（或等级资料）的多个样本比较）频数表资料（或等级资料）的多个样本比较步骤同上步骤同上 （详见课本详见课本P160例题例题11-6，此略），此略）

27、35/45Statistics多个样本比较的秩和检验多个样本比较的秩和检验（Kruskal-Wallis法，即法，即H检验）检验）例例11-4某医院外科用三种手术方法治疗肝癌某医院外科用三种手术方法治疗肝癌患者患者15例，每组例，每组5例，进入各组的患者系用随例，进入各组的患者系用随机方法分配，每例术后生存月数如表机方法分配，每例术后生存月数如表11-6的第的第、栏。试问三种不同手术方法治疗肝栏。试问三种不同手术方法治疗肝癌的效果有无差别。癌的效果有无差别。11-6三种手术方法治疗肝癌患者的术后表生存月数三种手术方法治疗肝癌患者的术后表生存月数 1.建立假设：建立假设：H0：三个总体分布的中心

28、位置相同。三个总体分布的中心位置相同。H1：三个总体分布的中心位置不全相同。三个总体分布的中心位置不全相同。0.052.编秩：编秩：将各组数据将各组数据混合由小到大编秩次，混合由小到大编秩次，如遇相同数值时，若相同数值在不同组内，如遇相同数值时，若相同数值在不同组内，则取平均秩次。则取平均秩次。表表11-6三种手术方法治疗肝癌患者的术后生存月数三种手术方法治疗肝癌患者的术后生存月数 3.求各组秩和：求各组秩和：4.计算检验统计量计算检验统计量H值：值：分子为分子为H值，分母值，分母C为校正数为校正数，校正后，校正后，HcH，P值减小。值减小。4.确定确定P值，判断结果：值，判断结果：查附表查附

29、表12，H界值表，界值表，H0.055.78，按，按=0.05水水准，得准，得P0.05，拒绝拒绝H0，故三种手术方法术后生，故三种手术方法术后生存月数不全相等。存月数不全相等。（二）频数表法（二）频数表法:例例比较不同孕周产妇的乳量。比较不同孕周产妇的乳量。表表11-8分娩时孕周与乳量的关系分娩时孕周与乳量的关系1建立假设：建立假设：H0：三个总体分布的中心位置相同。三个总体分布的中心位置相同。H1：三个总体分布中心位置不全相同。三个总体分布中心位置不全相同。0.052编秩（第编秩（第5、6、7列）列）3求秩和（第求秩和（第8、9、10列）列）4计算检验统计量：计算检验统计量：5.确定确定P

30、值，判断结果：值，判断结果：当样本组数或样本例数超过当样本组数或样本例数超过H界值界值表的范围，则认为表的范围，则认为H值近似服从值近似服从 2分布，分布，查查 2界值表。界值表。本例：本例：查查 2界值表，界值表，=k 1=3 1=2，20.005，2=10.60，得，得P0.005，可认为分娩可认为分娩时的孕周对乳量有影响。时的孕周对乳量有影响。第四节第四节 随机区组设计的秩和检验随机区组设计的秩和检验36/45Statistics一、随机区组设计的秩和检验一、随机区组设计的秩和检验一、随机区组设计的秩和检验一、随机区组设计的秩和检验 什么叫随机区组设计？什么叫随机区组设计？也叫配伍组设计

31、，是配对设计的扩大，是也叫配伍组设计，是配对设计的扩大，是将将条件条件相近的实验对象（三个及以上）配成一组，每组中相近的实验对象（三个及以上）配成一组，每组中的不同对象分别给以不同的处理措施。此种设计增的不同对象分别给以不同的处理措施。此种设计增强了各组间的均衡性，可进一步控制混杂性偏倚。强了各组间的均衡性，可进一步控制混杂性偏倚。本本法法是是由由M-FriedmanM-Friedman在在符符号号检检验验的的基基础础上上提提出出来的，常称为来的，常称为FriedmanFriedman检验，又称检验，又称 M M 检验。检验。36/45Statistics 【例例11-711-7】某某医医师师

32、为为探探讨讨多多发发伤伤后后血血清清内内皮皮细细胞胞特特异异性性分分子子-1(ESM-1)-1(ESM-1)的的变变化化，将将1212窝窝小小白白鼠鼠（每每窝窝3 3只只）制制作作成成多多发发伤伤动动物物模模型型，再再把把每每窝窝3 3只只多多发发伤伤小小白白鼠鼠随随机机分分配配到到3 3个个实实验验组组，采采静静脉脉血血4ml4ml，分分离离血血清清，成成批批检检测测ESM-1(ng/ml)ESM-1(ng/ml)，结结果果见见表表11-1011-10。该该医医师师采采用用两两样样本本均均数数比比较较的的t t检检验验，结结果果三三次次两两两两比比较较的的结结果果均均为为，差差异异均均无无统

33、统计计学学意意义义，故故认认为为多多发发伤伤后后三三组组的的血血清清内内皮皮细细胞胞特特异异性分子性分子-1(ESM-1)-1(ESM-1)无差别。无差别。（1）M检验检验（Friedman法）查表法法）查表法（用于区组数（用于区组数b15，处理组数处理组数k15）方法步骤：方法步骤：建立检验假设建立检验假设计算检验统计量计算检验统计量v编秩、求各处理组的秩和编秩、求各处理组的秩和 Ri：每个区组的数每个区组的数据由小到大分别编秩，相同数据取平均秩次据由小到大分别编秩，相同数据取平均秩次v求平均秩和求平均秩和b区组数区组数k处理组数处理组数v求检验统计量求检验统计量M值值 确定确定P值、作出统

34、计推断结论值、作出统计推断结论(P162例题例题11-7）37/45（2）F检验检验（用于区组数（用于区组数b15或处理组数或处理组数k15）方法步骤：方法步骤：建立检验假设建立检验假设计算检验统计量计算检验统计量v编秩：编秩原则编秩：编秩原则v求各处理组的秩和求各处理组的秩和 Riv计算所有数据秩次的平方和计算所有数据秩次的平方和A值值：若各区组内无相同秩次，可用若各区组内无相同秩次，可用：38/45Statistics39/45v计算计算B值值：v计算检验统计量计算检验统计量F 值：值：确定确定P值、作出统计推断结论值、作出统计推断结论查查F 界值表：界值表：1=k-1 2=(b-1)(k

35、-1)Statistics第五节第五节 多个样本两两比较的秩和检验多个样本两两比较的秩和检验成组设计多个样本两两比较的秩和检验成组设计多个样本两两比较的秩和检验 Nemenyi法法例题例题11-8（以例题（以例题11-5为例）为例）41/45l例例11-511-5为探讨血管紧张素转化酶抑制剂开博通抗动脉粥为探讨血管紧张素转化酶抑制剂开博通抗动脉粥样硬化的机制，某医师将样硬化的机制，某医师将15只兔随机分为三组，即正常饮只兔随机分为三组，即正常饮食组、高脂饮食组及高脂饮食加开博通组，喂养食组、高脂饮食组及高脂饮食加开博通组，喂养14周后测周后测量兔血清总胆固醇量兔血清总胆固醇(TC)水平，结果如

36、下：水平，结果如下：l正常饮食组：正常饮食组：0.860.781.380.881.70l高脂饮食组：高脂饮食组：40.3541.7944.0040.7227.05l高脂饮食加开博通组：高脂饮食加开博通组：41.2033.0440.6032.0540.49该该资资料料不不服服从从正正态态分分布布，在在进进行行多多组组间间比比较较时时采采用用KruskalKruskal-Wallis-Wallis 检检验验，推推断断结结论论为为拒拒绝绝，接接受受，这这与与前前面面所所介介绍绍的的多多组组比比较较方方差差分分析析类类似似，只只能能得得出出三三组组资资料料总总体体分分布布不不同同或或不不全全相相同同的

37、的结结论论，不不能能判判断断三三组组资资料料中中哪哪些些组组间间差差异异具具有有统统计计学学意意义义。若若要要对对每每两两个个总总体体分分布布做做出出有有无无不不同同的的推推断，需要做组间的两两比较。断，需要做组间的两两比较。41/45l多组秩和检验中做进一步的两两比较，需调整检验多组秩和检验中做进一步的两两比较，需调整检验水准水准。l检验水准的调整：对检验水准的调整：对k个样本组进行多次两两比较，个样本组进行多次两两比较，会增加第会增加第I类错误的概率。为保证第类错误的概率。为保证第I类错误的概率类错误的概率总共不超过，每次比较的第总共不超过，每次比较的第I类错误概率应采用式类错误概率应采用

38、式(11-12)进行调整。进行调整。l=/总的比较次数总的比较次数l（1）=2/k(k-1)(多组两两比较）多组两两比较）（2）=/k-1(实验组与同一对照组比）实验组与同一对照组比）41/45方法步骤：（方法步骤：（P165)建立检验假设建立检验假设确定检验水准确定检验水准计算检验统计量计算检验统计量两两比较两两比较:设为设为k k个样本组，和分别为比较的第个样本组，和分别为比较的第i i组和组和 第第j j组样本的秩和，其平均秩和分别为和。组样本的秩和，其平均秩和分别为和。(1)精确法：精确法：样样 本本 含含 量量 较较 小小 时时，应应 采采 用用 两两 样样 本本 秩秩 和和 检检

39、验验 的的 方方法法，求求得得统统计计量量的的数数值值后后，借借助助SAS和和SPSS软软件件的的“exact”功功能能得得到到相相应应的的P值值（包包括括该该数数值值对对应应的的概概率和更极端情形的概率）率和更极端情形的概率）。41/45(1)精确法：精确法：样样 本本 含含 量量 较较 小小 时时，应应 采采 用用 两两 样样 本本 秩秩 和和 检检 验验 的的 方方法法，求求得得统统计计量量的的数数值值后后，借借助助SAS和和SPSS软软件件的的“exact”功功能能得得到到相相应应的的P值值（包包括括该该数数值值对对应应的的概概率和更极端情形的概率）率和更极端情形的概率）。(2)正态近

40、似法：正态近似法：当当各各样样本本组组例例数数较较大大时时，按按下下式式(11-15)求求得得第第i组组和和第第j组样本比较的值。组样本比较的值。41/45平均秩次平均秩次，公式的分母为公式的分母为的标准误。的标准误。经多个样本比较的秩和检验，各组间有经多个样本比较的秩和检验，各组间有差别，试再进一步作两两比较。差别，试再进一步作两两比较。l表表11-12三组兔血清总胆固醇水平三组兔血清总胆固醇水平(mmol/L)的两两比较的两两比较l利用标准正态分布或利用标准正态分布或SAS和和SPSS统计软件求得统计量统计软件求得统计量数值所对应的数值所对应的P值。将值。将P值与调整以后的检验水准值与调整

41、以后的检验水准比比较，若较，若P0.051组和组和3组比较：组比较：P0.052组和组和3组比较：组比较：P0.054.确定确定P值，判断结果：值，判断结果：按按0.05水准，除了水准，除了2组与组与3组比较，组比较，拒绝拒绝H0，其余每二组间均无差别，可认其余每二组间均无差别，可认为乙法优于丙法。为乙法优于丙法。二、随机区组设计多个样本间的多重比较二、随机区组设计多个样本间的多重比较与与完完全全随随机机设设计计多多个个样样本本组组间间的的多多重重比比较较类类似似，对对于于随随机机区区组组设设计计资资料料，当当用用Friedman Friedman M M检检验验拒拒绝绝后后，同同样样需需要要

42、对对各各处处理理组组间间进进行行多多重重比比较较。采采用用调调整整检检验验水水准准，两两两两做做相相关关样样本本的的秩秩和和检检验验，即即配配对对wilcoxonwilcoxon符号秩和检验。符号秩和检验。42/45以例以例11-7说明随机区组设计资料多重比说明随机区组设计资料多重比较的基本步骤。（较的基本步骤。（P162)某某医医师师为为探探讨讨多多发发伤伤后后血血清清内内皮皮细细胞胞特特异异性性分分子子-1(ESM-1)-1(ESM-1)的的变变化化，将将1212窝窝小小白白鼠鼠（每每窝窝3 3只只）制制作作成成多多发发伤伤动动物物模模型型，再再把把每每窝窝3 3只只多多发发伤伤小小白白鼠

43、鼠随随机机分分配配到到3 3个个实实验验组组，采采静静脉脉血血4ml4ml，分分离离血血清清，成成批批检检测测ESM-1(ng/ml)ESM-1(ng/ml)，结果见表，结果见表11-1011-10。42/45小小 结结1.掌握非参数检验的概念及秩和检验的优掌握非参数检验的概念及秩和检验的优缺点。缺点。2.掌握配对设计资料、两组和多组完全随机掌握配对设计资料、两组和多组完全随机设计资料秩和检验的编次方法。设计资料秩和检验的编次方法。3.掌握如何确定配对设计资料、两组和多组掌握如何确定配对设计资料、两组和多组完全随机设计资料秩和检验的检验统计量。完全随机设计资料秩和检验的检验统计量。Thanks for listening!If you have questions,please ask us or postGdmcyfyx 44/45 Bye Bye!Statistics45/45