秩和检验1113课件

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,第十章 基于秩次的非参数检验,单样本和配对设计资料的符号秩和检验,两组独立样本比较的秩和检验,多组独立样本比较的秩和检验,本章内容,例10-3 用两种药物杀灭钉螺，每批用200300只活钉螺，用药后清点钉螺的死亡数，并计算死亡率，结果如表10-3，问两种药物杀死钉螺的效果有无差异？,一、定量变量两组独立样本的秩和检验,甲药,乙药,死亡率（）,秩次,死亡率（）,秩次,32.5,5.5,16.0,1,35.5,7,22.5,2,40.5,10,26.0,3,40.5,10,28.5,4,49.0,12,32.5,5.5,49.0,13,38.0,8,51.5,14,40.5,10,n,1,=7,T,1,=71.5,n,2,=7,T,2,=33.5,表10-3 两种药物杀灭钉螺死亡率的比较,1,、,建立检验假设,H,0,：两种药物杀灭钉螺死亡率的总体中位,数相同,H,1,：两种药物杀灭钉螺死亡率的总体中位,数不同,=0.05,该资料为百分率数据，不服从正态分布，采用秩和检验。,2,、,计算检验统计量,T,值,1编秩,两组数据从小到大统一混合编秩。,相同数据，在同一组内，顺次编秩；在不同组内，取平均秩次。,甲药,乙药,死亡率（）,秩次,死亡率（）,秩次,32.5,5.5,16.0,1,35.5,7,22.5,2,40.5,10,26.0,3,40.5,10,28.5,4,49.0,12,32.5,5.5,49.0,13,38.0,8,51.5,14,40.5,10,n,1,=7,T,1,=71.5,n,2,=7,T,2,=33.5,表10-3 两种药物杀灭钉螺死亡率的比较,2求各组秩和：n 较小为n1，较大为n2，相应秩和T1、T2。,T,1,71.5,，,T,2,33.5,T,1,T,2,N(N,1)/2=14(14,1)/2=105,3确定统计量T：TT1 71.5,N=,n,1,n,2,假设n1 n2,TT1或TT2,3,、,确定,P,值，做出推断结论,T,0.05,=36,69,，,T,0.01,=32,73,，,0.01,P,0.05,。,按a=0.05水准，拒绝H0，可认为两种药物杀死钉螺的效果有差异。,1查表法：当n110，n1-n2 10时，查附表10。,谁高？,1正态近似法：当n110，n1-n2 10时，T 逼近正态分布，故计算Z 值。,当相持较多25%时，校正Z值。,t,j,：第,j,个相同差值的个数。,根本思想,H,0,成立,，两样本来自分布相同总体，两组数据相间排列，理论上：,n,1,n,2,，,T,1,T,2,T,1,/n,1,T,2,/n,2,n,1,n,2,T,1,T,2,如,T,1,超出范围，,H,0,成立可能性小，按,水准，拒绝,H,0,，接受,H,1,。,T1/n1N1/2，T1n1N1/2。,T1/n1N1/2，T1n1N1/2。,实际上，由于抽样误差的存在。,例10-4 某医科大学营养教研室为了解居民体内核黄素营养状况，于某年冬夏两季收集成年居民口服5mg核黄素后4小时的负荷尿，测定体内核黄素含量，结果见表10-4，试比较该地居民冬夏两季体内核黄素含量有无差异？,二、有序分类变量两组独立样本的秩和检验,状况,例数,合计,秩次范围,平均秩次,夏季,冬季,缺乏,10,22,32,1,32,16.5,不足,14,18,32,33,64,48.5,适宜,16,4,20,65,84,74.5,合计,40,44,84,表,10-4,某地居民冬夏两季体内核黄素营养状况,1,、,建立检验假设,H,0,：冬夏两季居民体内核黄素含量的中位,数相同,H,1,：冬夏两季居民体内核黄素含量的中位,数不同,=0.05,2,、,计算检验统计量,T,值,1编秩：两组数据按等级从小到大统一编秩。,计算各等级合计数,确定各等级秩次范围。,求各等级平均秩次：等于秩次范围上、下,限之和除以,2,。,2求各组秩和：n 较小为n1=40。,T,1,16.510+48.514+74.516=2036,T,2,16.522+48.518+74.54=1534,T,N(N,1)/2=84(84,1)/2=3570,3确定统计量T：TT1 2036,状况,例数,合计,秩次范围,平均秩次,夏季,冬季,缺乏,10,22,32,1,32,16.5,不足,14,18,32,33,64,48.5,适宜,16,4,20,65,84,74.5,合计,40,44,84,表,10-4,某地居民冬夏两季体内核黄素营养状况,3,、,确定,P,值，做出推断结论,n,1,=40,超出,T,值范围，计算,Z,。,3.0048,查t界值表v=，得 P0.001。按a=0.05水准，拒绝H0，可认为冬夏两季居民体内核黄素含量有差异。,8.5365,第三节 多组独立样本比较的秩和检验,又称KruskalWallis H 检验，用于推断计量资料或等级资料的多个独立样本所来自的多个总体分布有无差异。,一、定量变量多组独立样本的秩和检验,例10-5 某医院用三种不同方法治疗15例胰腺癌患者，每种方法各治疗5例。治疗后生存月数见表10-5，问这三种方法对胰腺癌患者的疗效有无差异？,甲法,乙法,丙法,生存月数,秩次,生存月数,秩次,生存月数,秩次,3,2.5,6,6,2,1,4,4,9,12,3,2.5,7,7.5,10,13,5,5,8,10,12,14,7,7.5,8,10,13,15,8,10,R,i,34,60,26,n,i,5,5,5,表10-5 3种方法治疗胰腺癌患者的生存月数比较,1,、,建立检验假设,H,0,：,3,法治疗后患者生存时间总体中位数相同,H,1,：,3,法治疗后患者生存时间总体中位数不同,=0.05,该资料首选单因素,F,分析，但生存时间不服从正态分布，故采用,Kruskal,Wallis,H,检验。,2,、,计算检验统计量,H,值,1编秩,三组数据从小到大统一混合编秩。,相同数据，在同一组内，顺次编秩；在不同组内，取平均秩次。,甲法,乙法,丙法,生存月数,秩次,生存月数,秩次,生存月数,秩次,3,2.5,6,6,2,1,4,4,9,12,3,2.5,7,7.5,10,13,5,5,8,10,12,14,7,7.5,8,10,13,15,8,10,R,i,34,60,26,n,i,5,5,5,表10-5 3种方法治疗胰腺癌患者的生存月数比较,2求各组秩和Ri：,R,1,34,，,R,2,60,，,R,3,26,3计算统计量H：,R,N(N,1)/2=15(15,1)/2=120,当相持较多,(25%),时,需计算校正,H,C,值：,H,C,=,H/C,t,j,：第,j,个相同值的个数。,3,、,确定,P,值，做出推断结论,H,0.05,=5.78,，,H,0.01,=7.98,，,0.01,P,0.05,。,按a=0.05水准，拒绝H0，可认为三种方法治疗胰腺癌患者的疗效有差异。,1查表H界值表法：当组数k=3，ni 5时，查附表11。,2查2界值表法：当k 和ni 超出H界值表时，由于H0 成立时H 值近似服从=k-1的分布2，查2界值表确定P 值。,二、有序变量多组独立样本的秩和检验,例10-6 某医院用3种方法治疗慢性喉炎，结果见表10-6。问3种方法的疗效有无差异？,疗效等级,例数,秩次范围,平均,秩次,甲法,乙法,丙法,合计,无效,24,20,20,64,1,64,32.5,好转,26,16,22,64,65,128,96.5,显效,72,24,14,110,129,238,183.5,治愈,186,32,22,240,239,478,358.5,合计,308,92,78,478,表10-6 3种方法治疗慢性喉炎的疗效比较,1,、,建立检验假设,H,0,：,3,种疗效的总体分布位置相同,H,1,：,3,种疗效的总体分布位置不同,=0.05,2,、,计算检验统计量,H,值,1编秩：方法同例10-4。,2求各组秩和Ri：,R,1,32.524+96.526+183.572+358.5186=83182,R,2,32.520+96.516+183.524+358.532=18070,R,3,32.520+96.522+183.514+358.522=13229,R,N(N,1)/2=478(478,1)/2=114481,疗效等级,例数,秩次范围,平均,秩次,甲法,乙法,丙法,合计,无效,24,20,20,64,1,64,32.5,好转,26,16,22,64,65,128,96.5,显效,72,24,14,110,129,238,183.5,治愈,186,32,22,240,239,478,358.5,合计,308,92,78,478,表10-6 3种方法治疗慢性喉炎的疗效比较,3计算统计量H：,3,、,确定,P,值，做出推断结论,按a=0.05水准，拒绝H0，可认为三种方法治疗慢性喉炎的效果有差异。,=,k,-1=3-1=2,查,2,界值表得：,P,0.005,。,三、多个独立样本间的多重比较,经过多个独立样本比较的KruskalWallis H 检验拒绝H0，接受H1，认为总体的位置不同或不全相同时，假设要进一步推断是哪两个总体分布位置不同，需进行两两比较。,两两比较的方法很多，现介绍t 检验。,=,N,-,k,H,为,Kruskal,Wallis,H,检验的统计量,H,或,H,C,。,例10-7 对例10-6资料做两两比较。,1,、,建立检验假设,H,0,：第,i,组与第,j,组疗效的总体分布位置相同,H,1,：第,i,组与第,j,组疗效的总体分布位置不同,=0.05,2,、,计算检验统计量,t,值,（,1,）求各组平均秩次,R,i,甲组：,R,1,83182/308=270.07,乙组：,R,2,18070/92=196.41,丙组：,R,3,13229/78=169.60,2列两两比较计算表,对比组,t,P,甲与乙,73.66,4.742,0.001,甲与丙,100.47,6.062,0.001,乙与丙,26.81,1.322,0.10,P,0.20,表10-7 例10-7资料的两两比较,3,、,确定,P,值，做出推断结论,按a=0.05水准，甲组与乙组，甲组与丙组，均拒绝H0；乙组与丙组不拒绝H0。故可认为三种方法治疗慢性喉炎疗效的差异在于甲法与其它两法之间，而乙法与丙法间无差异。,=,N,-,k,=478-3=475,查,t,界值表得,P,值。,