医学统计学(列联表资料分析).ppt

1,【例7-1】某研究人员为研究不同地区的人群血型分类构成是否一样，收集了亚洲，欧洲和北美洲人的A，B，O，AB的血型资料。,2,第七章列联表资料分析,列联表：由两个或两个以上的变量进行交叉分类的频数分布表。,3,列联表（行列表）资料模式,横向变量R类,纵向变量C类,简称RC表,4,教学目的与要求,掌握：2检验的基本思想和用途四格表独立性检验配对四格表独立性检验了解：Fisher精确概率法的计算单变量拟合度检验频数分布的拟合度检验,5,教学内容提要,重点讲解：2检验的基本思想四格表独立性检验配对四格表独立性检验介绍：Fisher精确概率法的计算拟合度检验,6,第一节列联表2检验,2检验及原理四格表资料的2检验双向无序行列表的2检验配对计数资料的2检验四格表确切概率法,7,双向无序列联表资料单向有序列联表资料双向有序且属性相同列联表资料双向有序且属性不同列联表资料,一、列联表（行列表）资料类型,8,三组消化性溃疡住院病人四种证型的构成,1.双向无序列联表资料,9,2.单向有序列联表资料,10,3.双向有序且属性不同列联表资料,11,4.双向有序且属性相同列联表资料,12,四格表资料,13,二、四格表资料的2检验,2检验是一种检验实际频数（Actualfrequecy）与理论频数（Theoreticalfrequecy）之间偏离度的方法。偏离度用2值（统计量）表示。,2=（A-T）2/T,1.2检验的基本思想,14,2.四格表资料的2检验,比较两个样本率、两种治疗方法之间的差异有无显著性。这类资料由4个基本数据构成：两种治疗方法的有效例数和无效例数，统计学称这类资料为四格表资料。,15,2.四格表资料的2检验,四格表资料的模式,16,【例7-2】为了解铅中毒病人是否有尿棕色素增加现象，分别对病人组和对照组的尿液作尿棕色素定性检查，结果如下，问铅中毒病人和对照人群的尿棕色素阳性率有无差别？两组人群尿棕色素阳性率比较,建立检验假设和确定检验水准H0：总体阳性率1=2H1：12=0.05计算检验统计量,其中为Ai实际频数，Ti为理论频数,nR是行合计，nC是列合计，n是四格之和。,【例7-2】假设检验步骤：,理论频数是在H0成立的条件下每个格子理论上的数据。（1=2=两样本合计率）若H0成立，则合计的阳性率为38/73=52.05%则铅中毒病人呈阳性的理论上应有T11=3638/73=18.74人,2值的大小决定于：（1）A-T的差值；（2）格子数（严格地说是df）的多少。,T12=17.26T21=19.26T22=17.74=23.12确定p值=(R-1)(C-1)=1,查2界值表（附表6）得P40，T值2=4.125确定P值作推论,【例7-4】某矿石粉厂当生产一种矿石粉石时，在数天内即有部分工人患职业性皮肤炎，在生产季节开始，随机抽取15名车间工人穿上新防护服，其余仍穿原用的防护服，生产进行一个月后，检查两组工人的皮肤炎患病率，结果如下：,穿新旧两种防护服工人的皮肤炎患病率比较,问两组工人的皮肤炎患病率有无差别？,H0：1=2；H1：12。=0.05计算最小的理论频数T11=1511/43=3.84,140，所以宜用2值的校正公式=2.94查2界值表，得0.05p0.10按=0.05水平不拒绝H0，尚不能认为穿不同防护服的皮肤炎患病率有差别。若不校正，2=4.33p0.05,【例7-5】研究复方哌唑嗪对高血压病治疗效果的临床试验并与复方降压片和安慰剂作对照，结果如下表，问三种药物效果有无差别？,三种药物治疗高血压病的有效率比较,27,【例7-6】某省为了解花生黄曲霉素污染，随机观察了三个地区的花生受黄曲霉素B1污染的情况，结果见下表。请问这三个地区花生的黄曲霉素B1污染率是否不同？某省三个地区花生的黄曲霉素B1污染率,28,三、双向无序列联表的2检验,1.列联表的2值计算2=(R-1)(C-1),29,【例7-5】三种药物治疗高血压病有效率比较,2.多个率的比较,H0：三种药物的总体有效率相等，即1=2=3H1：1、2、3不相等或不全相等=0.05,2=(2-1)(3-1)=2查2界值表得p0.005,按=0.05水平拒绝H0，认为三种药物的治疗效果不相同或不全相同。,31,【例7-7】1979年某地爆发松毛虫病，333例患者按年龄以14岁为界分成2组，资料如表，试考察两组病人病变类型的构成比有无差别？,某地两组松毛虫病患者型别构成比较,3.构成比的比较,32,a.H0：两组患者中患各型松毛虫病的构成比相同H1：两组患者中患各型松毛虫病的构成比不同=0.05b.2=70.14c.查表得p0.005d.按=0.05水平拒绝H0，认为两组患者中患各型松毛虫病的构成比不同。,3.构成比的比较,行列表的2检验要求理论频数T不宜太小，要求不宜有1/5以上的格子的T5，或有一个T1。如果以上条件不能满足，可采用：增加样本含量删去某行或某列合理地合并部分行或列用精确概率法多个率或构成比比较的2检验，结论为拒绝H0时，仅表示几组有差别，并非任2组之间都有差别。若要了解之，可进行多个率之间的两两比较(如2表的分割法)。,4.行列表的2检验注意事项,34,5.多个率(构成比)两两比较行列表的分割（略）,【例7-5】三种药物治疗高血压病的有效率比较,2=,35,分割表1分割表2,21=0.58322=21.873,【例7-8】（续）哌唑嗪组与降压片组的有效率接近，将此2组合并，再与另一个样本比较，得到2个分割表：,20.05,1=3.84,36,可以得出21+22=22.456约等于原来的2=22.36，两者不相等是由于分割表中理论频数的数值与总表中理论频数不同所致。结论：哌唑嗪与降压片总体有效率相同,但与安慰剂的有效率不相同。,37,例题：考察一种新型透析疗法的效果,四、22列联表资料的2检验,38,【例7-8】用两种检验方法对某食品作沙门氏菌检验，结果如下，试比较两种方法的阳性结果是否有差别？配对设计两种方法结果比较,39,配对设计的计量资料（观察值）比较，比较两总体均数的差异，用配对t检验配对设计的计数资料，我们用22列联表资料的2检验,40,1.22列联表资料的2检验,配对四格表资料分析(略)（）独立性检验（是否有关联）用四格表2检验法（）优势性检验（是否相同）用McNemar检验法（）一致性检验用Kappa检验法,41,【例7-8】用两种检验方法对某食品作沙门氏菌检验，结果如下，试比较两种方法的阳性结果是否有差别？配对设计两种方法结果比较,着眼于两法结果不一致的部分。,42,两法结果不一致的部分：H0：总体B=C；H：总体BC检验统计量,当b+c40时2=当b+c40时2=1样本例数较少时还可以用精确概率法,原理同配对设计的符号检验。,43,H0：甲法阳性乙法阴性例数与甲法阴性乙法阳性例数相同，即总体B=C；H1：总体BC=0.052=12.90查2界值表p0.005，拒绝H0，认为两法检验结果不一样，荧光抗体法阳性结果高于常规培养法。,44,配对（列联表）资料的2检验,【例7-9】有93份咽喉涂抹标本，每份标本分别接种在甲乙两种白喉杆菌培养基上，观察其生长情况，结果如下表，试分析两种培养基的效果是否不同？表5甲乙两种培养基的生长情况,45,1.建立假设H0：B=C；H1：BC确定显著性水平=0.052.计算统计量:b+c=10+19=2940，校正公式：2=（|b-c|-1）2/（b+c）=2.207=（2-1）（2-1）=13.确定P值4.作推断,46,五、四格表的确切概率法,【例7-10】在某牧区观察慢性布鲁氏病患者植物血凝素皮肤试验反应，资料如下，问活动型与稳定型布氏患者反应阳性率有无差别？,47,五、四格表的确切概率法,适用条件四格表资料n40时；四格表中有A=0时；四格表中有T1时；用其它检验方法所得概率接近检验水准时。,48,五、四格表的确切概率法,【例7-10】在某牧区观察慢性布鲁氏病患者植物血凝素皮肤试验反应，资料如下，问活动型与稳定型布氏患者反应阳性率有无差别？,49,五、四格表的确切概率法,基本思想,50,H0：两型阳性率相等，即1=2H1：12在将原四格表周边数字（行和与列和）不变的情况下，变动a值，以得到不同情况：,51,p为|A-T|大于等于原来表中的|A-T|值（0.14）的情况下的p(i)的和，即p=p(0)+p(1)+p(4)=0.2668按=0.05水准接受H0,认为两型阳性率相等。,计算每种情况下的|A-T|和|A-T|大于等于原来表中的|A-T|值（0.14）的p(i)值，p(i)值按下列公式来计算：,52,【例7-11】某医师为研究乙肝免疫球蛋白预防胎儿宫内感染HBV的效果，将33例HBsAg阳性孕妇随机分为预防注射组和非预防组，结果见下表。问两组新生儿的HBV总体感染率有无差别？,53,2211,T11=6A-T=2.0,54,第二节分类资料的拟合优度检验,此例只有两种分类。因此应有理论次数=总数分类项数，即=700.5=35（人）,=,=,【例7-12】随机抽取70名学生，调查他们对高中分文理科的意见，回答赞成的有42人，反对的有28人。问对分科的意见有无差异？,一、单变量拟合优度检验,55,二、频数分布的拟合优度检验,1.建立假设，确定检验水准：H0：评定结果服从正态分布，H1：评定结果不服从正态分布=0.05,【例7-13】今对某校100名学生进行操行评定，分优、良、中、差四等，评定结果为：优19人、良39人、中35人、差7人。试检验其分布的形式是否属于正态分布？,2.计算值：,第一等分（优）：1.53之间的面积比率为：0.49870.4332=0.0655，即7%。第二等分（良）：01.5之间面积比率为0.4332，即43%。第三等分（中）：01.5之间面