《医学统计学秩和检验幻灯片.ppt》由会员分享,可在线阅读,更多相关《医学统计学秩和检验幻灯片.ppt(80页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、医学统计学秩和检验第1页,共80页,编辑于2022年,星期五医学统计学医学统计学 第七章第七章 秩和检验秩和检验 第2页,共80页,编辑于2022年,星期五本章学习要点:本章学习要点:非参数检验的概念非参数检验的概念参数检验与非参数检验的区别参数检验与非参数检验的区别几种常用的秩和检验方法几种常用的秩和检验方法第3页,共80页,编辑于2022年,星期五参数统计和非参数统计的概念参数统计和非参数统计的概念参数统计概念参数统计概念 参参数数统统计计是是指指在在样样本本资资料料满满足足一一定定条条件件(正正态态分分布布,方方差差齐齐)的的基基础础上上,对对总总体体的的参参数数(总总体体均均数或率)进
2、行假设检验的方法,称为参数统计。数或率)进行假设检验的方法,称为参数统计。参数统计均有一定的条件要求,若不满足条件则参数统计均有一定的条件要求,若不满足条件则不能使用,可考虑进行变量变换或采用不要求条件的不能使用,可考虑进行变量变换或采用不要求条件的其它统计方法,如秩和检验。其它统计方法,如秩和检验。第4页,共80页,编辑于2022年,星期五非参数统计概念、优点和缺点非参数统计概念、优点和缺点 1 1、概念、概念非非参参数数统统计计是是指指不不考考虑虑资资料料的的分分布布类类型型,不不对对总总体体参参数数进进行行检检验验,而而对对资资料料的的分分布布是是否否相相同同进进行行检检验验,这这种种统
3、统计计方方法法称称非非参参数数检检验验。秩秩和和检检验验、符符号号检检验验、游游程程检检验验、RiditRidit分分析析、X X2 2检检验等均属于非参数检验方法。验等均属于非参数检验方法。第5页,共80页,编辑于2022年,星期五2、优点、优点(1)不受总体分布的限定,适用范围广;可用于)不受总体分布的限定,适用范围广;可用于各种统计资料,主要用于偏态分布资料、分布不明资各种统计资料,主要用于偏态分布资料、分布不明资料;料;(2)对数据要求不严格;可用于不能准确定量的资料,)对数据要求不严格;可用于不能准确定量的资料,主要用于等级资料,开口资料;主要用于等级资料,开口资料;(3)有些方法在
4、样本例数不多时,尚简便易行。)有些方法在样本例数不多时,尚简便易行。第6页,共80页,编辑于2022年,星期五3、缺点、缺点不不能能充充分分利利用用资资料料所所提提供供的的信信息息(仅仅考考虑虑位位次次大大小小),故故检检验验效效率率较较参参数数检检验验低低,犯犯第第二二类类错错误误的的概概率率较较参参数数检检验验大大,同同一一资资料料要要达达到到相相同同的的检检验验效效能能(1-),则则非非参参数数检检验验比比参参数数检检验验所所需需的的样样本本例例数多。数多。因因此此,在在进进行行统统计计分分析析时时,应应首首先先考考虑虑是是否否满满足足参参数检验,不满足参数检验时才考虑使用非参数检验。数
5、检验,不满足参数检验时才考虑使用非参数检验。第7页,共80页,编辑于2022年,星期五首选非参数检验的情况:首选非参数检验的情况:1、等级资料、等级资料2、开口资料、开口资料3、明显偏态分布的资料、明显偏态分布的资料第8页,共80页,编辑于2022年,星期五第一节第一节Wilcoxon符号秩和检验符号秩和检验一、配对设计的两样本比较一、配对设计的两样本比较例例7.1本例应首先考虑用配对本例应首先考虑用配对t检验,这里只是用于说检验,这里只是用于说明配对设计秩和检验的方法。明配对设计秩和检验的方法。首选配对秩和检验的例子:P291,习题 4.2 4.4 4.6。(用于配对设计的计量资料)第9页,
6、共80页,编辑于2022年,星期五第10页,共80页,编辑于2022年,星期五1、建立假设和确定检验水准、建立假设和确定检验水准H0:差值总体中位数差值总体中位数Md=0H1:差值总体中位数差值总体中位数Md00.052、求差值、求差值3、编秩:、编秩:(1)依依差差值值绝绝对对值值从从小小到到大大编编秩秩,再再根根据据差差值值的的正正负负给秩次冠以正负号;给秩次冠以正负号;(2)差值为零时,舍去不计差值为零时,舍去不计(例数相应减例数相应减1);(3)差值相等,符号相同,按顺序编秩;差值相等,符号相同,按顺序编秩;(4)差值相等,符号不同,取平均秩次差值相等,符号不同,取平均秩次 第11页,
7、共80页,编辑于2022年,星期五4、求秩和并确定检验统计量、求秩和并确定检验统计量(1)分别求正、负秩次之和分别求正、负秩次之和T+、T-本例:本例:T+=3.5;T-=41.5(2)以绝对值较小的秩和为检验统计量,本例以绝对值较小的秩和为检验统计量,本例T-3.5注:总秩和注:总秩和=,本例本例T+T_=45,而而第12页,共80页,编辑于2022年,星期五5、确定、确定P值值 (1)查表法查表法当当n50时,查附表时,查附表6:界值表(配界值表(配对比较的符号秩和检验)对比较的符号秩和检验)以例数以例数n确定查哪一行,然后自左向右用确定查哪一行,然后自左向右用T与每一与每一栏界值相比。栏
8、界值相比。在界值范围之内,在界值范围之内,P值大于表上方相应概率值大于表上方相应概率在界值范围之外,在界值范围之外,P值小于表上方相应概率值小于表上方相应概率(往(往右移一栏继续查)右移一栏继续查)本例本例n9,3.5,在双侧,在双侧P0.05的界值范围的界值范围(540)之外,在双侧)之外,在双侧P0.02的界值范围(的界值范围(342)之内,故之内,故0.02P0.05。第13页,共80页,编辑于2022年,星期五6、推断结论、推断结论0.02P50时时,可可计计算算u值值,确确定定P值值(此此时仍属于非参数法)时仍属于非参数法)注:注:n50时,时,T近似呈正态分布近似呈正态分布而当相同
9、秩次较多而当相同秩次较多(超过超过25%)时时,需计算校正需计算校正u值。值。(注:(注:tj为相同差值的个数)为相同差值的个数)第15页,共80页,编辑于2022年,星期五tj为相同差值的个数。如例7.1,有两个5,两个15。故,t1=2,t2=2,第16页,共80页,编辑于2022年,星期五方法与例方法与例7.1类似。类似。将样本各观察值与已知的总体中位数相减,将样本各观察值与已知的总体中位数相减,用差值进行符号的秩和检验。用差值进行符号的秩和检验。二、单一样本与总体中位数比较二、单一样本与总体中位数比较第17页,共80页,编辑于2022年,星期五第18页,共80页,编辑于2022年,星期
10、五1、建立假设和确定检验水准、建立假设和确定检验水准H0:该厂工人尿氟含量总体中位数等于该厂工人尿氟含量总体中位数等于2.15H1:该厂工人尿氟含量总体中位数大于该厂工人尿氟含量总体中位数大于2.15 单侧单侧0.052、求差值、求差值3、编秩:、编秩:(1)依依差差值值绝绝对对值值从从小小到到大大编编秩秩,再再根根据据差差值值的的正正负负给秩次冠以正负号;给秩次冠以正负号;(2)差值为零时,舍去不计差值为零时,舍去不计(例数相应减例数相应减1);(3)差值相等,符号相同,按顺序编秩;差值相等,符号相同,按顺序编秩;(4)差值相等,符号不同,取平均秩次差值相等,符号不同,取平均秩次 第19页,
11、共80页,编辑于2022年,星期五4、求秩和并确定检验统计量、求秩和并确定检验统计量(1)分别求正、负秩次之和分别求正、负秩次之和T+、T-本例:本例:T+=62.5;T-=3.5(2)以绝对值较小的秩和为检验统计量,本例以绝对值较小的秩和为检验统计量,本例T-3.5注:总秩和注:总秩和=,本例本例T+T_=66,而而第20页,共80页,编辑于2022年,星期五5、确定、确定P值值 用查表法查表法以以n=11查附表查附表6:界值表(配对界值表(配对比较的符号秩和检验)得单侧比较的符号秩和检验)得单侧T0.005=(561),),本例本例,3.5,在单侧,在单侧T0.005范围之外,故范围之外,
12、故P0.005。6、推断结论、推断结论0.005,在在0.05水准上,拒绝水准上,拒绝H0,接,接受受H1,差异有统计学意义,故可认为该厂工人差异有统计学意义,故可认为该厂工人尿氟含量高于当地正常人。尿氟含量高于当地正常人。第21页,共80页,编辑于2022年,星期五第二节成组设计两样本比较的秩和检验第二节成组设计两样本比较的秩和检验(Wilcoxon两样本比较法)两样本比较法)一、原始数据的两样本比较一、原始数据的两样本比较第22页,共80页,编辑于2022年,星期五第23页,共80页,编辑于2022年,星期五1、建立假设和确定检验水准、建立假设和确定检验水准H H0 0:两组雌鼠增重的总体
13、分布相同两组雌鼠增重的总体分布相同 H H1 1:两组雌鼠增重的总体分布不同两组雌鼠增重的总体分布不同0.052、编秩:、编秩:(1)先先将将两两组组数数据据分分别别从从小小到到大大排排序序,再再将将两两组组数数据统一由小到大编秩;据统一由小到大编秩;(2)相同数据在同一组里,按顺序编秩;相同数据在同一组里,按顺序编秩;(3)相同数据在不同组里,取平均秩次。相同数据在不同组里,取平均秩次。第24页,共80页,编辑于2022年,星期五3、求秩和并确定检验统计量、求秩和并确定检验统计量(1)分别求各组的秩和分别求各组的秩和T1、T2;(2)以样本例数以样本例数n较小者的秩和较小者的秩和T为检验统计
14、量(若为检验统计量(若n1=n2,任取一组秩和为任取一组秩和为T);本例);本例T244.54、确定、确定P值值(1)查查表表法法:以以n1和和n2n1,查查附附表表7T界界值值表表(两两样样本本比较的秩和检验用)比较的秩和检验用):在界值范围之内,在界值范围之内,P值大于表上方相应概率值大于表上方相应概率在界值范围之外,在界值范围之外,P值小于表上方相应概率值小于表上方相应概率恰好等于界值,恰好等于界值,P值等于表上方相应概率值等于表上方相应概率第25页,共80页,编辑于2022年,星期五本例本例n17,n2n1=12-7=5;查附表;查附表7T界值表界值表得:得:T界值界值P值值49910
15、.1046940.0542980.02401000.0144.5,在双侧,在双侧P0.05的界值范围之外,在双侧的界值范围之外,在双侧P0.02的界值范围之内,则的界值范围之内,则0.02P0.05。本例本例T244.5第26页,共80页,编辑于2022年,星期五6、推断结论、推断结论本本例例0.02P0.55、推断结论、推断结论本例本例P0.5,在在0.05水准上,不拒绝水准上,不拒绝H0,差别无,差别无统计学意义,尚不能认为该药对两种支气管炎的疗效统计学意义,尚不能认为该药对两种支气管炎的疗效不同。不同。查u值表df=的的t值,值,就是就是u值值!第34页,共80页,编辑于2022年,星期
16、五第三节第三节成组设计多个样本比较的秩和检验成组设计多个样本比较的秩和检验(KruskalWallis法)法)H检验检验一、原始数据多个样本比较一、原始数据多个样本比较二、二、等级资料(频数表资料)等级资料(频数表资料)多多样本比较样本比较第35页,共80页,编辑于2022年,星期五一、原始数据多个样本比较的秩和检验一、原始数据多个样本比较的秩和检验 例 7.5第36页,共80页,编辑于2022年,星期五表表9.4教室在不同时间空气中教室在不同时间空气中CO含量(含量(mg/m)课课前前课课中中课课后后含量含量秩次秩次含量含量秩次秩次含量含量秩次秩次(1)(2)(3)(4)(5)(6)0.48
17、14.4512.52.9570.5324.73143.0780.5534.77153.1890.5544.82163.20100.5854.89173.30110.6265.00184.4512.52192.557.5666第37页,共80页,编辑于2022年,星期五1、建立假设和确定检验水准、建立假设和确定检验水准H0:三个总体的位置相同三个总体的位置相同H1:三个总体的位置不同或不全相同三个总体的位置不同或不全相同0.052、编秩:、编秩:(1)各各组组分分别别从从小小到到大大排排列列,再再将将各各组组数数据据由由小小到到大大统统一编秩;一编秩;(2)在同一组里有相同的数据,按顺序排列;在
18、同一组里有相同的数据,按顺序排列;(3)在不同组里有相同数据,取平均秩次。在不同组里有相同数据,取平均秩次。3、求各组的秩和、求各组的秩和 第38页,共80页,编辑于2022年,星期五4、计算检验统计量、计算检验统计量H而当相同秩次较多而当相同秩次较多(超过超过25%)时时,需计算校正需计算校正HC值值.HC=H/CC=1本例:本例:第39页,共80页,编辑于2022年,星期五5、确定、确定P值值(1)若若k3,每组例数,每组例数5,查附表,查附表8,H界值表界值表;(2)若若组组数数k3,每每组组例例数数5,以以自自由由度度v=k-1,查查附附表表5,X2界值表界值表(此时(此时H服从服从X
19、2分布)。分布)。本本例例k3,每每组组例例数数=5,以以n=15,查查附附表表8,H界界值值表得:表得:H0.01=7.98。H9.245H0.01P0.016、推断结论、推断结论本例本例PX20.005,2P0.0056、推断结论、推断结论本例本例PM0.05,P0.05。由由于于PF0.01(5,25),P0.01由于由于P0.01,按按=0.05,拒绝拒绝H0,接受,接受H1,狗服药,狗服药后不同时间血中药物浓度不同或不全相同。后不同时间血中药物浓度不同或不全相同。第63页,共80页,编辑于2022年,星期五第五节第五节 多个样本两两比较的秩和检验多个样本两两比较的秩和检验当多样本比较
20、秩和检验的结论是拒绝当多样本比较秩和检验的结论是拒绝H0时,时,只能得出各组总的有差别,若要了解某两样本只能得出各组总的有差别,若要了解某两样本之间有无差别,需进一步做组间的两两比较。之间有无差别,需进一步做组间的两两比较。第64页,共80页,编辑于2022年,星期五一、成组设计多个样本两两比较的秩和检验一、成组设计多个样本两两比较的秩和检验(Nemenyi法)法)方法和步骤:1、列出两两对比组;2、求各个对比组的平均秩和差值的绝对值;3、计算界值第65页,共80页,编辑于2022年,星期五对对例例9.5进进行行两两两两比比较较:小小白白鼠鼠接接种种三三种种不不同同菌菌型型的的伤寒杆菌后存活天
21、数两两比较。伤寒杆菌后存活天数两两比较。本例:本例:C=0.9893,V=2,n1=n2=n3=5第66页,共80页,编辑于2022年,星期五第67页,共80页,编辑于2022年,星期五1、建立假设、建立假设H0:任两个总体分布相同任两个总体分布相同H1:任两个总体的位置不同任两个总体的位置不同0.052、求求各各个个对对比比组组平平均均秩秩和和的的差差值值D作作为为检检验验统统计计量量;(;(D=)3、计算界值、计算界值式式中中:C为为相相同同秩秩次次校校正正值值,由由界界值值表表查查得得,N为总例数,为总例数,和和为两个对比组的例数。为两个对比组的例数。第68页,共80页,编辑于2022年
22、,星期五第69页,共80页,编辑于2022年,星期五二、二、随机区组设计资料两两比较的秩和检验随机区组设计资料两两比较的秩和检验 方法和步骤:方法和步骤:1、列出两两对比组;、列出两两对比组;2、求两两对比组秩和之差的绝对值、求两两对比组秩和之差的绝对值、|RA-RB|;3、计算检验界值、计算检验界值;v=(b-1)(k-1)第70页,共80页,编辑于2022年,星期五对例对例9.10做两两比较。做两两比较。第71页,共80页,编辑于2022年,星期五表表7.9各组间的两两比较各组间的两两比较对比组对比组|RA-RB|P对比组对比组|RA-RB|P1与与3|15-34|=190.051与与4|
23、15-32|=170.013与与5|34-15|=190.053与与6|34-6|=280.011与与6|15-6|=90.014与与5|32-15|=170.012与与3|24-34|=100.014与与6|32-6|=260.012与与4|24-32|=80.015与与6|15-6|=90.012与与5|24-15|=90.01除服药后除服药后0.5小时与小时与24小时,小时,6和和8小时药物浓度无差别外,小时药物浓度无差别外,其余时间药物浓度不同。其余时间药物浓度不同。第72页,共80页,编辑于2022年,星期五首选秩和检验的资料:首选秩和检验的资料:1、明显偏态分布资料、明显偏态分布资
24、料2、等级资料;、等级资料;3、开口资料、开口资料秩和检验小结秩和检验小结第73页,共80页,编辑于2022年,星期五狗服用阿司匹林后不同时间血中药物浓度(狗服用阿司匹林后不同时间血中药物浓度(r/ml)0.5小时小时1小时小时6小时小时8小时小时24小时小时48小时小时15243432156第74页,共80页,编辑于2022年,星期五某年某地一次伤寒爆发潜伏期频数表 潜伏期 发病人数 3 3 5-24 7-20 9-17 11-14 13-7 15-6 17-2 19-1 21 2明显偏态开口资料第75页,共80页,编辑于2022年,星期五用某药治疗不同病情的老慢支病人,疗效见下表,用某药治疗不同病情的老慢支病人,疗效见下表,试比较两种病情的疗效有无差别?试比较两种病情的疗效有无差别?疗效疗效 单纯性单纯性 单纯性合并肺气肿单纯性合并肺气肿 合计合计 控制控制 65 42 107 65 42 107 显效显效 18 6 24 18 6 24 有效有效 30 23 53 30 23 53 无效无效 13 11 24 13 11 24 合计合计 126 82 208 126 82 208 等级资料等级资料第76页,共80页,编辑于2022年,星期五