医学统计学列联表资料分析

1、关于医学统计学列联表资料分析1第一张，PPT共五十八页，创作于2022年6月2第七章 列联表资料分析列联表：由两个或两个以上的变量进行交叉分类的频数分布表。第二张，PPT共五十八页，创作于2022年6月3列联表（行列表）资料模式横向变量R类, 纵向变量C类, 简称RC表第三张，PPT共五十八页，创作于2022年6月4教学目的与要求 掌握：2检验的基本思想和用途四格表独立性检验配对四格表独立性检验了解：Fisher精确概率法的计算单变量拟合度检验频数分布的拟合度检验 第四张，PPT共五十八页，创作于2022年6月5教学内容提要 重点讲解：2检验的基本思想四格表独立性检验配对四格表独立性检验介绍：

2、Fisher精确概率法的计算拟合度检验 第五张，PPT共五十八页，创作于2022年6月6第一节 列联表 2检验2检验及原理四格表资料的2检验双向无序行列表的2检验配对计数资料的2检验四格表确切概率法第六张，PPT共五十八页，创作于2022年6月7双向无序列联表资料单向有序列联表资料双向有序且属性相同列联表资料双向有序且属性不同列联表资料一、列联表（行列表）资料类型第七张，PPT共五十八页，创作于2022年6月8三组消化性溃疡住院病人四种证型的构成 1. 双向无序列联表资料第八张，PPT共五十八页，创作于2022年6月92. 单向有序列联表资料第九张，PPT共五十八页，创作于2022年6月103

3、. 双向有序且属性不同列联表资料第十张，PPT共五十八页，创作于2022年6月114. 双向有序且属性相同列联表资料第十一张，PPT共五十八页，创作于2022年6月12四 格 表 资 料第十二张，PPT共五十八页，创作于2022年6月13二、四格表资料的2 检验 2检验是一种检验实际频数（Actual frequecy）与理论频数（Theoretical frequecy ）之间偏离度的方法。 偏离度用2值（统计量）表示。2 =（A-T）2/T1. 2 检验的基本思想第十三张，PPT共五十八页，创作于2022年6月142. 四格表资料的2检验 比较两个样本率、两种治疗方法之间的差异有无显著性。

4、这类资料由4个基本数据构成：两种治疗方法的有效例数和无效例数，统计学称这类资料为四格表资料。第十四张，PPT共五十八页，创作于2022年6月152. 四格表资料的2检验 四格表资料的模式 组 别阳 性阴 性合计甲aba+b乙cdc+d合计a+cb+dnA11 (T11)第十五张，PPT共五十八页，创作于2022年6月16【例7-2】为了解铅中毒病人是否有尿棕色素增加现象，分别对病人组和对照组的尿液作尿棕色素定性检查，结果如下，问铅中毒病人和对照人群的尿棕色素阳性率有无差别？两组人群尿棕色素阳性率比较第十六张，PPT共五十八页，创作于2022年6月 建立检验假设和确定检验水准 H0：总体阳性率1

5、=2 H1：12 =0.05 计算检验统计量其中为Ai实际频数，Ti为理论频数, nR是行合计，nC是列合计，n是四格之和。【例7-2】假设检验步骤：第十七张，PPT共五十八页，创作于2022年6月T11T12 理论频数是在H0成立的条件下每个格子理论上的数据。 （1=2= 两样本合计率） 若H0成立，则合计的阳性率为 38/73=52.05% 则铅中毒病人呈阳性的理论上应有 T11=3638/73=18.74人第十八张，PPT共五十八页，创作于2022年6月2值的大小决定于： （1）A-T的差值； （2）格子数（严格地说是df）的多少。 T12=17.26 T21=19.26 T22=17.

6、74 =23.12 确定p值 =(R-1)(C-1)=1 , 查2界值表（附表6）得 P40，T值 2=4.125确定P值作推论1486100309012044176220第二十三张，PPT共五十八页，创作于2022年6月【例7-4】 某矿石粉厂当生产一种矿石粉石时，在数天内即有部分工人患职业性皮肤炎，在生产季节开始，随机抽取15名车间工人穿上新防护服，其余仍穿原用的防护服，生产进行一个月后，检查两组工人的皮肤炎患病率，结果如下：穿新旧两种防护服工人的皮肤炎患病率比较第二十四张，PPT共五十八页，创作于2022年6月问两组工人的皮肤炎患病率有无差别？ H0：1=2；H1：12。=0.05 计算

7、 最小的理论频数T11=1511/43=3.84, 1T1140，所以宜用2值的校正公式 =2.94 查2界值表，得0.05p 0.10 按=0.05水平不拒绝H0，尚不能认为穿不同防护服的皮肤炎患病率有差别。 若不校正，2=4.33 p0.05第二十五张，PPT共五十八页，创作于2022年6月 【例7-5】 研究复方哌唑嗪对高血压病治疗效果的临床试验并与复方降压片和安慰剂作对照，结果如下表，问三种药物效果有无差别？组别 有效 无效 合计 有效率% 复方哌唑嗪30104075.00复方降压片20103066.67安慰剂7253221.88合计574510255.88三种药物治疗高血压病的有效率

8、比较第二十六张，PPT共五十八页，创作于2022年6月27 【例7-6】某省为了解花生黄曲霉素污染，随机观察了三个地区的花生受黄曲霉素B1污染的情况，结果见下表。请问这三个地区花生的黄曲霉素B1污染率是否不同？某省三个地区花生的黄曲霉素B1污染率调查地区受检样品合计污染率（）未污染污染甲 6232979.3乙30144431.8丙 8 31127.3合计44408447.6第二十七张，PPT共五十八页，创作于2022年6月28三、双向无序列联表的2检验 1. 列联表的2值计算 2 = = =(R-1)(C-1) 第二十八张，PPT共五十八页，创作于2022年6月29【例7-5】三种药物治疗高血

9、压病有效率比较 2. 多个率的比较组别 有效 无效 合计 有效率% 复方哌唑嗪30104075.00复方降压片20103066.67安慰剂7253221.88合计574510255.88第二十九张，PPT共五十八页，创作于2022年6月H0：三种药物的总体有效率相等，即1=2=3 H1：1、2、3不相等或不全相等 =0.05 2=(2-1)(3-1)=2 查2界值表得 p0.005 按=0.05水平拒绝H0，认为三种药物的治疗效果不相同或不全相同。 第三十张，PPT共五十八页，创作于2022年6月31【例7-7】 1979年某地爆发松毛虫病，333例患者按年龄以14岁为界分成2组，资料如表，试

10、考察两组病人病变类型的构成比有无差别？ 某地两组松毛虫病患者型别构成比较 3. 构成比的比较第三十一张，PPT共五十八页，创作于2022年6月32a. H0：两组患者中患各型松毛虫病的构成比相同 H1：两组患者中患各型松毛虫病的构成比不同 =0.05b. 2=70.14c. 查表得p0.005d.按=0.05水平拒绝H0，认为两组患者中患各型松毛虫病的构成比不同。3. 构成比的比较第三十二张，PPT共五十八页，创作于2022年6月行列表的2检验要求理论频数T不宜太小，要求不宜有1/5以上的格子的T5，或有一个T1。如果以上条件不能满足，可采用： 增加样本含量 删去某行或某列 合理地合并部分行或

11、列 用精确概率法多个率或构成比比较的2检验，结论为拒绝H0时，仅表示几组有差别，并非任2组之间都有差别。若要了解之，可进行多个率之间的两两比较(如2表的分割法)。4. 行列表的2检验注意事项 第三十三张，PPT共五十八页，创作于2022年6月345. 多个率(构成比)两两比较行列表的分割（略） 【例7-5】 三种药物治疗高血压病的有效率比较 2=组别 有效 无效 合计 有效率% 复方哌唑嗪30104075.00复方降压片20103066.67安慰剂7253221.88合计574510255.88第三十四张，PPT共五十八页，创作于2022年6月35分割表1 分割表2组别有效无效合计组别有效无效

12、合计哌唑嗪301040用药组502070降压片201030安慰剂72532合计502070合计574510221=0.583 22=21.873【例7-8】（续） 哌唑嗪组与降压片组的有效率接近，将此2组合并，再与另一个样本比较，得到2个分割表：20.05,1=3.84第三十五张，PPT共五十八页，创作于2022年6月36 可以得出21+22=22.456约等于原来的2=22.36，两者不相等是由于分割表中理论频数的数值与总表中理论频数不同所致。 结论：哌唑嗪与降压片总体有效率相同,但与安慰剂的有效率不相同。第三十六张，PPT共五十八页，创作于2022年6月37例题： 考察一种新型透析疗法的效

13、果病人序号透析前透析后 131.618.2 220.7 7.3 336.426.5 433.123.7 529.522.6 620.710.7 750.325.1 831.220.9 936.623.71028.116.5四、22列联表资料的2检验 第三十七张，PPT共五十八页，创作于2022年6月38【例7-8】用两种检验方法对某食品作沙门氏菌检验，结果如下，试比较两种方法的阳性结果是否有差别？ 配对设计两种方法结果比较荧光抗体法（甲法）常规培养法（乙法）合计+-+160(a)26(b)186-5(c)48(d)53合计16574239第三十八张，PPT共五十八页，创作于2022年6月39配

14、对设计的计量资料（观察值）比较，比较两总体均数的差异，用配对 t检验配对设计的计数资料，我们用22列联表资料的 2检验 第三十九张，PPT共五十八页，创作于2022年6月401. 22列联表资料的2检验 配对四格表资料分析(略)（）独立性检验（是否有关联）用四格表2检验法 （）优势性检验（是否相同）用McNemar检验法 （）一致性检验用Kappa检验法第四十张，PPT共五十八页，创作于2022年6月41【例7-8】用两种检验方法对某食品作沙门氏菌检验，结果如下，试比较两种方法的阳性结果是否有差别？ 配对设计两种方法结果比较荧光抗体法（甲法）常规培养法（乙法）合计+-+160(a)26(b)1

15、86-5(c)48(d)53合计16574239着眼于两法结果不一致的部分。第四十一张，PPT共五十八页，创作于2022年6月42两法结果不一致的部分： H0：总体B=C； H：总体BC检验统计量当b+c40时 2=当b+c40时 2= =1 样本例数较少时还可以用精确概率法, 原理同配对设计的符号检验。第四十二张，PPT共五十八页，创作于2022年6月43H0：甲法阳性乙法阴性例数与甲法阴性乙法阳性例数相同，即总体B=C；H1：总体BC =0.05 2= =12.90查2界值表 p0.005，拒绝H0，认为两法检验结果不一样，荧光抗体法阳性结果高于常规培养法。 荧光抗体法常规培养法合计+-+

16、160(a)26(b)186-5(c)48(d)53合计16574239第四十三张，PPT共五十八页，创作于2022年6月44配对（列联表）资料的2检验 【例7-9】 有93份咽喉涂抹标本，每份标本分别接种在甲乙两种白喉杆菌培养基上，观察其生长情况，结果如下表，试分析两种培养基的效果是否不同？表5 甲乙两种培养基的生长情况乙培养法甲培养法合计阳性阴性阳性37 1047阴性19 2746合计563793第四十四张，PPT共五十八页，创作于2022年6月451. 建立假设 H0 ： B=C ；H1 ： BC 确定显著性水平 =0.052. 计算统计量: b+c=10+19=2940， 校正公式：

17、2 =（|b-c|- 1）2/ （b+c）=2.207 =（2-1）（2-1）=13. 确定P值4. 作推断第四十五张，PPT共五十八页，创作于2022年6月46五、四格表的确切概率法【例7-10】 在某牧区观察慢性布鲁氏病患者植物血凝素皮肤试验反应，资料如下，问活动型与稳定型布氏患者反应阳性率有无差别？第四十六张，PPT共五十八页，创作于2022年6月47五、四格表的确切概率法适用条件 四格表资料 n40 时；四格表中有A= 0时； 四格表中有T1时；用其它检验方法所得概率接近检验水准时。第四十七张，PPT共五十八页，创作于2022年6月48五、四格表的确切概率法【例7-10】 在某牧区观察

18、慢性布鲁氏病患者植物血凝素皮肤试验反应，资料如下，问活动型与稳定型布氏患者反应阳性率有无差别？第四十八张，PPT共五十八页，创作于2022年6月49五、四格表的确切概率法基本思想 第四十九张，PPT共五十八页，创作于2022年6月50H0：两型阳性率相等，即1=2 H1：12在将原四格表周边数字（行和与列和）不变的情况下，变动a值，以得到不同情况：第五十张，PPT共五十八页，创作于2022年6月51p为| A- T|大于等于原来表中的| A- T|值（0.14）的情况下的p(i)的和，即 p=p(0)+p(1)+p(4)=0.2668按=0.05水准接受H0,认为两型阳性率相等。 计算每种情况

19、下的|A- T|和|A- T|大于等于原来表中的| A- T|值（0.14）的p(i)值，p(i)值按下列公式来计算：第五十一张，PPT共五十八页，创作于2022年6月52【例7-11】 某医师为研究乙肝免疫球蛋白预防胎儿宫内感染HBV的效果，将33例HBsAg阳性孕妇随机分为预防注射组和非预防组，结果见下表。问两组新生儿的HBV总体感染率有无差别？第五十二张，PPT共五十八页，创作于2022年6月532211T11=6 A-T=2.0第五十三张，PPT共五十八页，创作于2022年6月54第二节 分类资料的拟合优度检验 此例只有两种分类。因此应有理论次数= 总数分类项数，即=700.5=35（人）= = 【例7-12】随机抽取70名学生，调查他们对高中分文理科