医学统计学第七章 秩和检验课件

卫生统计学,欢迎学习,卫生统计学,第七章 秩和检验,参数统计和非参数统计的概念,一、参数统计概念 参数统计是指在样本资料满足一定条件（正态分布，方差齐）的基础上，对样本所来自的总体的参数（总体均数或率）进行估计或是否相等进行检验，称为参数统计。 参数统计均有一定的条件要求，若不满足条件则不能使用，可考虑进行变量变换或采用不要求条件的其它统计方法，如秩和检验。,二、非参数统计概念、优点和缺点 1、概念 非参数统计是指不考虑资料的分布形式及其总体参数，而对资料的分布是否相同进行检验，这种统计方法称非参数检验。秩和检验、符号检验、游程检验、Ridit分析、X2检验等均属于非参数检验方法。,2、优点 （1）不受总体分布的限定，适用范围广；可用于各种统计资料，主要用于偏态分布资料、分布不明资料； （2）对数据要求不严格；可用于不能准确定量的资料，主要用于等级资料，开口资料； （3）有些方法在样本例数不多时，尚简便易行。,3、缺点 不能充分利用资料所提供的信息（仅考虑位次大小），故检验效率较参数检验低，犯第二类错误的概率较参数检验大，同一资料要达到相同的检验效能（1-），则非参数检验比参数检验所需的样本例数多。 因此，在进行统计分析时，应首先考虑是否满足参数检验，不满足参数检验时才考虑使用非参数检验。,配对设计差值符号的秩和检验 （Wilcoxon配对法）,表9.1 两种方法测定10名健康人尿汞值( ug/L ) 编号 离子交换法 蒸馏法 差值 秩次 (1) (2) (3) (4)=(2)-(3) (5) 1 0.5 0.0 0.5 2 2 2.2 1.1 1.1 7 3 0.0 0.0 0.0 - 4 2.3 1.3 1.0 6 5 6.2 3.4 2.8 8 6 1.0 4.6 -3.6 -9 7 1.8 1.1 0.7 3.5 8 4.4 4.6 -0.2 -1 9 2.7 3.4 -0.7 -3.5 10 1.3 2.1 -0.8 -5 ,无效假设为差值的总体中位数等于0，备择假设为差值的总体中位数不等于0。假设无效假设成立的话（两种方法的测定结果结果无差别，两个总体分布的位置相同），其配对数值之差应服从以0为中心的对称分布，理论上正秩和与负秩和在理论上应该相等，即使有差别，也只能是随机误差，其差别不会很大。如果正负秩和差别很大的话，超过了允许的界值，我们就有理由拒绝无效假设。反之，就不拒绝。,1、建立假设和确定检验水准 H0: 差值总体中位数Md=0 H1: 差值总体中位数Md0 0.05,2、求差值 3、编秩： (1)依差值绝对值从小到大编秩，再根据差值的正负给秩次冠以正负号； (2)差值为零时，舍去不计(例数相应减1)； (3)差值相等，取其平均秩次；,4、求秩和并确定检验统计量 (1)分别求正、负秩次之和T+、T- 本例：T+=3.5 ；T- =41.5 (2)以绝对值较小的秩和为检验统计量，本例 T+3.5 注：总秩和= , 本例T+ T_=45，而,5、确定P值 (1)查表法 当n50时，查附表6（P268）: 界值表（配对比较的符号秩和检验） 以例数n确定查哪一行，然后自左向右用T与每一 栏界值相比。 在界值范围之内，P值大于表上方相应概率 在界值范围之外，P值小于表上方相应概率 等于界值，P值等于表上方相应概率 本例 n9，3.5，在双侧P0.05的界值范围（837）之外，P0.02。,(2)正态近似法 当n50时，可计算u值，确定P值 （属于非参数法） 注：n50时，T近似呈正态分布 而当相同秩次较多(超过25%)时,需计算校正u值。 （注：为相同差值的,6、推断结论 0.02P25时，T分布较好地近似正态分布。,单一样本与总体中位数比较,例已知某地正常人尿氟含量的中位数为2.15mmol/l。今在该地某厂随机抽取12名工人尿氟含量如表。问该厂工人尿氟是否高于当地正常人？,12名工人的尿氟含量测定结果,1、建立假设和确定检验水准 H0: 该厂工人尿氟含量总体中位数Md=2.15 H1:该厂工人尿氟含量总体中位数Md02.15 0.05,2、求差值 3、编秩： (1)依差值绝对值从小到大编秩，再根据差值的正负给秩次冠以正负号； (2)差值为零时，舍去不计(例数相应减1)； (3)差值相等，取平均秩次；,4、求秩和并确定检验统计量 (1)分别求正、负秩次之和T+、T- 本例：T+=62.5 ；T- =3.5 (2)以绝对值较小的秩和为检验统计量，本例 T-3.5,5、确定P值 (1)查表法 当n50时，查附表6: 界值表（配对比较的符号秩和检验） 以例数n确定查哪一行，然后自左向右用T与每一 栏界值相比。 在界值范围之内，P值大于表上方相应概率 在界值范围之外，P值小于表上方相应概率 等于界值，P值等于表上方相应概率 本例 n11，3.5，在单侧P0.005的界值范围（837）之外，P0.005。 6、下结论,成组设计两样本比较的秩和检验 （Wilcoxon两样本比较法）,原始数据两样本比较的秩和检验,表9.3 两组小鼠发癌后生存日数 实验组 对照组 生存日数 秩次 生存日数 秩次 10 9.5 2 1 12 12.5 3 2 15 15 4 3 15 16 5 4 16 17 6 5 17 18 7 6 18 19 8 7 20 20 9 8 23 21 10 9.5 90以上 22 11 11 12 12.5 13 14 n1=10 T1=170 n2=12 T2=83,1、建立假设和确定检验水准 H0:喂高蛋白与喂低蛋白的雌鼠增重的总体分 布相同 H1:喂高蛋白与喂低蛋白的雌鼠总体分布不同 0.05 2、编秩： (1)先将两组数据分别从小到大排序，再将两组数据统一由小到大编秩； (2)遇相同数据，取平均秩次。,3、求秩和并确定检验统计量 (1)分别求各组的秩和T1、T2 ； (2)以样本例数n较小者的秩和T为检验统计量（若n1=n2,任取一组秩和为T）；本例T244.5 4、确定P值 (1)查表法：以n1和n2n1，查附表7， T界值表（两样本比较的秩和检验用）: 在界值范围之内，P值大于表上方相应概率 在界值范围之外，P值小于表上方相应概率 恰好等于界值， P值等于表上方相应概率,本例 n17，n2n1=12-7=5；查附表7 T界值表 得： T界值 P值 49 91 0.1 46 94 0.05 42 98 0.02 40 100 0.0144.5，在双侧P0.05的界值范围（4694）之外，则 P0.05。,(2)正态近似法： 当n1或/和n2-n1超出附表7范围时，可计算u值，确定P值 。 而当相同秩次较多(超过25%)时,需计算校正u值。,6、推断结论 本例P0.05， 在0.05水准上，不拒绝H0，故尚不能认为该药对两种支气管病人的疗效分布不同。,成组设计多个样本比较的秩和检验 （KruskalWallis法）,原始数据多个样本比较的秩和检验,表9.4 教室在不同时间空气中CO含量（mg/m） 课 前 课 中 课 后 含量 秩次 含量 秩次 含量 秩次 （1） （2） （3） （4） （5） （6） 0.48 1 4.45 12.5 2.95 7 0.53 2 4.73 14 3.07 8 0.55 3 4.77 15 3.18 9 0.55 4 4.82 16 3.20 10 0.58 5 4.89 17 3.30 11 0.62 6 5.00 18 4.45 12.5 21 92.5 57.5 6 6 6,1、建立假设和确定检验水准 H0:三个总体的分布位置相同 H1: 三个总体的分布位置不同或不全相同 0.05 2、编秩： (1)各组分别从小到大排列，再将各组数据由小到大统一编秩； (2) 有相同的数据，取平均秩次。 3、求各组的秩和,4、计算检验统计量H 而当相同秩次较多(超过25%)时,需计算校正HC值. HC =H/C C=1本例：,5、确定P值 (1)若k3，每组例数5，查附表8 ,H界值表 ； (2)若组数k3，每组例数5，以自由度v=k-1，查附表8，X2界值表 （此时H服从X2分布）。 本例 k3，每组例数P0.01, P0.01。 6、推断结论 本例PX20.005,2 ，P0.005 6、推断结论 本例 P0.05一与三 5 5 8.6 6.89 8.54 0.05,5、推断结论 一 (9D) 二（11C）、三(DSC1) =3.7, =8, =12.3,随机区组设计（配伍设计）资料的秩和检验,一、M检验( Friedman法)查表法,5名受试者穿5种不同的防护服测得的脉搏数（次/分） 防护服A 防护服B 防护服C 防护服D 防护服E编号 脉搏 秩次 脉搏 秩次 脉搏 秩次 脉搏 秩次 脉搏 秩次 1 130 1 144 5 143 3.5 133 2 143 3.5 2 111 1.5 116 3 119 5 118 4 111 1.5 3 114 3 106 1 115 4 113 2 116 5 4 123 4 98 1 120 3 104 2 133 5 5 115 5 104 2 111 4 110 3 101 1 b 5 5 5 5 5 14.5 12 19.5 13 16 15 15 15 15 15 ( - ) 0.5 3 4.5 2 1 ( - )2 0.25 9 20.25 4 1,M检验方法： 1、将每个配伍组编秩，有相同的观察值则取平均秩次； 2、求每个处理组的秩和 ； 3、求平均秩和 ； = 4、求M；M= 5、查附表9，M界值表，确定P值。,H0:此药不同剂量时血清中DT值的总体分布相同 H1:此药不同剂量时血清中DT值的总体分布不同 =0.05 用b=7,k=4,查附表9，M界值表得：M0.05=92, M=213 M0.05, P0.05 由于P0.05，拒绝H0，接受H1，差别有统计学意义，故可以认为此药不同剂量时血清中DT值的影响有差异。,二、X2分布近似法,当处理组或配伍组超过附表9时，可以采用近似X2分布法。V=k-1,当各区组的相同秩次较多时，需要校正：,V=4-1=3查X2界值表，p F0.01(5,25),P0.01 由于P0.01, 按 =0.05,拒绝H0，接受H1，狗服药后不同时间血中药物浓度不同或不全相同。,随机区组设计资料两两比较的秩和检验,方法和步骤：1、列出两两对比组；2、求两两对比组秩和之差的绝对值、|RA-RB|;3、计算检验界值 ; v=(b-1)(k-1),对例9.10做两两比较。,表9.11 各组间的两两比较对比组 |RA-RB| P 对比组 |RA-RB| P1与3 |15-34|=19 0.051与4 |15-32|=17 0.05 3与6 |34- 6 |=28 0.011与6 |15- 6