卡方检验-研究生

2检验（Chi-squaretest）（3学时）吴成秋公共卫生学院卫生统计学教研室,2检验的应用,1.推断两个(多个)总体率或构成比之间是否有差别2.多个样本率间的多重比较3.两个分类变量之间有无关联性4.频数分布拟合优度检验,一、2检验的基本思想,例7-1某医生用两种疗法治疗前列腺癌，出院后随访3年，甲疗法治疗86例。乙疗法治疗95例，结果见下表。问两种疗法治疗前列腺癌病人的3年生存率是否有差别？,表10-7甲乙两种疗法治疗前列腺癌的3年生存率处理生存数死亡数合计生存率(%)甲疗法a52(46.09)b34(39.91)(a+b)8660.47乙疗法c45(50.91)d50(44.09)(c+d)9547.37合计(a+c)97(b+d)84(n)18153.59表中为基本数据，称为实际频数，用AR.C表示,表中是根据A计算得到的，称为理论频数，用TR.C表示。该资料称为四格表资料。,若假设两种疗法治疗前列腺癌病人的3年生存率相同，即H0:1=2=53.59%，则T1.1=8697/181=46.09T1.2=8684/181=39.91T2.1=9597/181=50.91T2.2=9584/181=44.09,TR.C=nRnc,n,2=(AT)2=(R1)(C1)=(21)(21)=1,T,假如H0:1=2成立，则A、T不会相差太大，那么2会小；(A与T的差别是由抽样误差所致)反之若检验假设H0不成立，则A、T会相差很大，那么2会大;(抽样误差不能引起A与T如此大的差别)判断结果，根据自由度的大小查2界值表当22,时，P，拒绝H0，接受H1当22,时，P，不拒绝H0.2反映了实际频数与理论频数的吻合程度,（一）、四格表的一般公式,二、一般四格表的2检验,2=(AT)2=(R1)(C1),T,以例10-7为例做2检验H0:1=2即两种疗法治疗前列腺癌的3年生存率相同H1:12即两种疗法治疗前列腺癌的3年生存率不相同=0.05T1.1=8697/181=46.09T1.2=8684/181=39.91T2.1=9597/181=50.91T2.2=9584/181=44.092=(AT)2=(52-46.09)2+(34-39.91)2+(45-50.91)2+(50-44.09)2=3.11=(21)(21)=1查2界值表得20.05,1=3.84,20.05。按=0.05的水准，不拒绝H0，还不能认为两种疗法治疗前列腺癌的3年生存率不相同。,T,44.09,50.91,39.91,46.09,(二)、四格表的专用公式,当n40，T5时，用2=,(adbc)2n,(ab)(cd)(ac)(bd),a、b、c、d为四格表的实际频数(ab)，(cd)，(ac)，(bd)是四格表周边的合计数n=abcd,为总的样本含量对例10-7的计算2=3.11,(5250-3445)2181,(86959784),显然，计算结果与前相同，但该法避免了计算理论数的麻繁,(三)、四格表的校正公式,1.2检验的应用条件(1)当n40且T5时，用2检验的基本公式或四格表的专用公式；当P时，改用四格表的确切概率法。(2)当n40时但有1T5时，用四格表2的校正公式或用四格表的确切概率法(3)当n40时，或T1时，用四格表的确切概率法2.四格表2的校正公式,2=(AT0.5)2,T,2=,(adbcn/2)2n,(ab)(cd)(ac)(bd),例7-2某医师用两种疗法治疗心绞痛，结果如表10-10，问两种疗法的疗效有无差别？表10-10两种疗法治疗心绞痛的效果率,H0:1=2即两种疗法的疗效相等H1:12即两种疗法的疗效不等=0.05T=299/59=4.42(计算行合计和列合计均为最小的理论数)c2=0.61=1,查2界值表得,P0.05。按=0.05的水准不拒绝H0，尚不能认为两种两种疗法的疗效不等.,(23362759/2)259,2930509,三、配对四格表资料的2检验,配对四格表资料的2检验是对配对设计研究所获得的资料为计数资料进行比较。,配对设计包括：同一批样品用两种不同的处理方法观察对象根据配对条件配成对子，同一对子内的不同个体分别接受不同的处理在病因或危险因素的研究中，将病人和对照按配对条件配成对子，研究是否存在某种病因或危险因素。,观察的结果只有阳性、阴性两种可能，清点成对资料时只有四种情况：a：甲+乙+b：甲+乙-c：甲-乙+d：甲-乙-列成下面四格表,适用于配对设计的计数资料配对资料的结果有甲乙结果abcd,配对设计的计数资料,2=(bc)2/(bc),=1,(bc)40c2=(bc1)2/(bc),=1(bc)40,例对120名确诊的乳腺癌病人，分别用甲乙两种方法来检测，甲方法检出阳性率为60%，乙方法检出阳性率为50%，甲乙一致检出阳性率为35%，问甲乙两方法检出阳性率有无差别？,表：两种方法的检测结果,H0:B=C(1=2)H1:BC(12)=0.05,bc=3018=48402=(3018)2/(3018)=3.0=1，查2界值表得P0.05。按=0.05水准，不拒绝H0,还不能认为两种方法的检验结果不同。,配对四格表资料的关联性检验(P761),H0：两种方法无关(无一致性)H1：两种方法有关(有一致性)=0.052两种方法检测结果的相关系数rr必定等于kappa检验的kappa值,K0.75一致性好0.40k0.75一致性较好K0.40一致性差,例73某实验室分别用乳胶凝集法和免疫荧光法对58名可疑系统红斑狼疮患者血清中抗体进行测定，结果见表73。问两种方法的检验结果有无差别?表73两种方法的检测结果,H0:B=CH1:B=C=0.05,bc=124=1440c2=(1221)2/(122)=5.79=1，查2界值表得0.01P0.025。按=0.05水准拒绝H0,接受H1,认为两种方法的检验结果不同，免疫荧光法的阳性检出率较高。,四、行列表资料的2检验,可以用于多个率及多个构成比资料的比较和双向无序分类资料关联性检验。其基本数据有三种情况：1、多个样本率比较。有R行2列，称为R2表2、两个样本的构成比比较，有2行C列，称为2C表3、多个样本的构成比比较，以及双向无序分类资料关联性检验，有R行C列，称为RC表,行列表检验的专用公式2=n(A2/nRnC1)=(R1)(C1),(一)多个样本率的比较,例7-6某医师研究物理疗法、药物治疗和外用膏药三种疗法治疗周围性面神经麻痹的疗效，资料见表7-8。问三种疗法的有效率有无差别？表7-8三种疗法有效率的比较,H0：1=2=3，即三种疗法的有效率相等H1：三种疗法的有效率不全相等=0.052=532(1992/20648172/20651262/144511)=21.04=(31)(21)=2查2界值表:得P0.005.按=0.05拒绝H0，接受H1，可以认为三种疗法治疗周围性面神经麻痹的有效率有差别。,(二)样本构成比的比较例7-7某医师在研究血管紧张素转化酶(ACE)基因I/D多态与2型糖尿病肾病(DN)的关系时，将249例型糖尿病患肾病分为两组，资料见表79。问两组型糖尿病患者的ACE基因型总体分布有无差别？表7-9DN组与无DN组型糖尿病患者ACE基因型分布,H0：两组2型糖尿病患者ACE基因型的总体构成比相同H1：两组2型糖尿病患者ACE基因型的总体构成比不同=0.05按公式（710）计算2=249(422/11172482/111120362/138571)=7.91=(21)(31)=2查2界值表得0.01P0.025。按=0.05拒绝H0,接受H1，认为DN与无DN的2型糖尿病患者的ACE基因型分布不同。,(三)双向无序分类资料的关联性检验,例7-8测得某地5801人的ABO血型和MN血型结果如表7-10，问两种血型系统之间是否有关联？表7-10某地5801人的血型,H0：两种血型系统间无关联H1：两种血型系统间有关联=0.052=58014322/(18231451)4902/(18231666)322/(3482684)1=213.16=(41)(31)=6查2界值表得P0.005。按=0.05拒绝H0，接受H1，认为两种血型系统间有关联。,若需进一步分析关系的密切程度时，可计算Pearson列联系数CC取值范围在01之间。0表示独立：1表示完全相关；愈接近于0，关系愈不密切，愈接近于1，关系愈密切。,(四)、行列表资料检验的注意事项,1、理论频数T不应太小；不宜有1/5以上的格子1T0.05按=0.05检验水准不拒绝H0，尚不能认为胆囊线癌与胆囊线瘤的P53基因表达阳性率不等。,频数分布拟合优度的2检验,1、频数分布拟合优度检验：是一种推断某种现象的频数分布是否符合某一理论分布的假设检验方法2、2值反映了实际频数和理论频数的吻合程度，因此，2检验可用于频数分布拟合优度检验，可推断一个资料的频数分布是否符合正态分布、二项分布、Poisson分布、负二项分布等。3、步骤,例7-12观察某克山病区克山病患者的空间分布情况，调查者将该区划分为279个取样单位，统计各取样单位历年的累计病例数，资料如表7-15，问此资料是否符合Poisson分布？,表7-15Poisson分布的拟合与检验,解：n=279,fx=686,fx