同济医学院-《SAS》课件-SAS的卡方检验(正式).ppt

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,卡方检验,SAS,应用,蒋红卫,Email:JHWCCC21CN.COM,学习目标,掌握四格表普通卡方检验和配对卡方检验方法以及相应的,SAS,程序；,了解,FREQ,过程语句格式；,熟悉,R,C,表资料的分类类型以及相应的统计检验方法；,掌握双向无序,R,C,表资料检验,以及,SAS,程序；,掌握单向有序,R,C,表资料检验以及,SAS,程序；,掌握趋势卡方检验方法以及,SAS,程序；,掌握分层,R,C,表,的分析以及,SAS,程序,；,概述,前面已介绍了两个率比较的检验，在观察例数不够大或拟对多个率进行比较时，检验就不适宜了，因为直接对多个样本率作两两间的检验有可能增加第一类误差。,2,检验可解决此类问题。,卡方检验是用途很广的一种假设检验方法，这里我们主要学习它在分类资料统计推断中的应用，包括：两个率或两个构成比比较的卡方检验；多个率或多个构成比比较的卡方检验以及分类资料的相关分析等。,四格表资料,定性指标分为有序的（如：疗效分为,“,治愈、显效、好转、无效、死亡,”,）和名义的（如：血型分为,“,O,、,A,、,B,、,AB,”,型）类，对于每个受试者来说，有序指标的观测结果只能是该有序指标若干等级中的级（如某人的疗效为,“,显效,”,）；名义指标的观测结果只能是该名义指标若干标志中的个如某人的血型为型），显然，无法像处理定量指标那样去直接分析定性指标，故这类资料常被整理成列联表的形式后再进行分析。,当表中只有个定性指标时，称为维列联表；有个或个以上定性指标时，称为多维列联表。常用,R,、,C,表示维列联表的行数和列数，并称为,R,C,表；当,R=C=,时，称为,表（或四格表）。,表看起来很简单，但根据资料所具备的条件有许多不同的处理方法。,四格表卡方检验的,SAS,程序,在,SAS/STAT,模块中,FREQ,、,TABULATE,和,SUMMARY,等过程可用于分类资料的统计描述，其中,FREQ,过程兼具统计描述和统计推断的功能，对分类变量计算频数分布，产生从一维到,n维的频数表和列联表；对于二维表，可进行,2,检验，对于三维表，可作,Mentel-Hanszel,分层分析。,FREQ,过程是,SAS,用于分析分类资料的一个常用过程。本节将先向大家介绍,FREQ,过程的语句及其格式。,FREQ,过程的语句基本格式如下：,Proc freq data=order=;,Table,分类变量*分类变量,/;,Weight,变量,;,Run;,四格表卡方检验的,SAS,程序,DATA,数据集：规定,PROC FREQ,语句使用的数据集；,ORDER,FREQ,，按频数递减顺序排列；,ORDER,DATA,，按数据集中出现的顺序排列；,ORDER,INTERNAL,按内部值排列,(,缺省,),；,ORDER,FORMATTED,，按外部格式值排列；,Table,语句指定构成表格的变量和表格结构。表格的结构由变量个数和变量排列顺序决定，一个,table,语句允许列出多个表格结构。,PROC FREQ,过程中可有多条,TABLES,语句，,TABLES,语句后可接多个表格请求式，每个请求式可包含任何数量的变量，从而得到所需的表格。,四格表卡方检验的,SAS,程序,如果,TABLES,语句缺省，则,FREQ,过程对数据集中的所有变量都给出相应的一维频数表。不规定任何选项时，若需某变量的一维频数，,FREQ,给出该变量每一水平的频数、累积频数、频数的百分比和累积百分比；若需二维频数表，,FREQ,产生交叉分组列表，即包括各格的频数、总频数的格百分数、行频数的格百分数和列频数的格百分数。,请求式由一个或多个用,“,*,”,连接起来的变量名组成。几个变量可放在括号中，如：,TABLES A*(B C),；等价于,TABLES A*B A*C,；,TABLES(A-C)*D,；等价于,TABLES A*D B*D C*D,；,四格表卡方检验的,SAS,程序,下列选项可用于,TABLES,语句中,“,/,”,的后面：,OUT,数据集：建立一个包含变量值和频数计数的输出数据集。如果,TABLES,语句中不止一个请求式，数据集的内容相应于,TABLES,语句中最后一个请求。,CHISQ,对每层作,c2,检验，包括,Pearson c2,、似然比,c2,和,Mantel-Haenszel c2,。此外还给出与,c2,检验有关的关联指标包括,Phi,系数、列联系数和,Cramer,s V,。对于,2,2,表，给出,Fisher,精确概率。,AGREE,进行配对,c 2,检验。,EXACT,对大于,2,2,的列联表计算,Fisher,精确概率。同时也给出,CHISQ,选项的全部统计量。,四格表卡方检验的,SAS,程序,MEASURES,对每层的二维表计算一系列关联指标及相应的标准误，包括,Pearson,和,Spearman,相关系数，以及,Gamma,和,Kendall,系数等。对于,2,2,表，还给出常用的危险度指标及其标准误。,CMH,给出,Cochran-Mantel-Haenszel,统计量，可检验在调整了,TABLES,语句中其它变量后，行变量与列变量之间的关联程度。对于,2,2,表，,FREQ,过程给出相对危险度估计及其可信区间，还给出各层关联度指标是否齐性的,Breslow,检验。,ALL,给出,CHISQ,、,MEASURES,、,CMH,所请求的全部统计量。,ALPHA,p,给出检验水准。缺省为,0.05,。,四格表卡方检验的,SAS,程序,EXPECTED,给出期望频数。,DEVIATION,给出每格的实际频数与期望频数的差值。,CELLCHISQ,给出每格对总,c2,的贡献，即计算每格的,(,实际频数,-,期望频数,)2/,期望频数。,CUMCOL,给出累积列百分数。,NOFREQ,不给出列联表中的格频数。,NOPERCENT,不给出列联表中的格百分数。,四格表卡方检验的,SAS,程序,NOROW,不给出列联表中各格的行百分数。,NOCOL,不给出列联表中各格的列百分数。,NOCUM,不给出频数表的累积频数和累积百分数。,NOPRINT,不给出表格，但给出,CHISQ,、,MEASURES,或,CMH,等语句所指定的统计量。,Trend,指令系统对,2,C,频数表的,C,个百分率进行,Cochran-Armitage,趋势检验；,WEIGHT,语句：通常每个观察值提供数值,1,给频数计数，当,WEIGHT,语句出现时，每个观察值提供的是该观察值的加权变量值。该值必须非负，但可不必为整数。只能使用一个,WEIGHT,语句，且该语句作用于所有的表。,四格表卡方检验,理论,频数,，记,为,T,。理论,数的计算公式为：,四格表卡方检验,卡方检验的统计量是,2,值，它是每个格子实际频数,A,与理论频数,T,差值平方与理论频数之比的累计和。每个格子中的理论频数,T,是在假定两组的,发病率相等的,情况下计算出来的，故,2,值越大，说明实际频数与理论频数的差别越明显，两组,发病率,不同的可能性越大。,例题,某医院欲比较异梨醇口服液（试验组）和氢氯噻嗪,+,地塞米松（对照组）降低颅内压的疗效。将,200,例颅内压增高症患者随机分为两组。问两组降低颅内压的总体有效率有无差别？,连续性校正公式,2,分布是正态变量的一种分布。设 是,k,个独立的标准正态变量，则 。,2,界值表就是根据这种连续性分布计算出来的。,2,统计量计算公式实质上是正态近似法。分类资料是间断性的，由此计算的,2,值不连续，尤其自由度为,1,的四格表，求出的概率可能偏小，此时需要对,2,值进行连续性校正，公式为,2,检验的应用条件,连续性校正主要针对四格表资料，尤其理论数较小时，连续性校正不可忽略。,四格表,2,检验的应用条件为：,当,n40,且所有,T5,时，用普通的,2,检验，若所得，改用确切概率法；,当,n40,但有,1T5,时，用校正的,2,检验；,当,n40,或有,T1,时，不能用,2,检验，改用确切概率法。,例题,某医师欲比较胞磷胆碱与神经节苷酯治疗脑血管疾病的疗效，将,78,例脑血管疾病患者随机分为两组。问两种药物治疗脑血管疾病的有效率是否有差别？,配对计数资料的卡方检验,把每一份样本平均分成两份，分别用两种方法进行化验，比较此两种化验方法的结果（两类计数资料）是否有本质的不同；或者分别采用甲、乙两种方法对同一批病人进行检查，比较此两种检查方法的结果（两类计数资料）是否有本质的不同，此时要用配对卡方检验。,配对计数资料的卡方检验,比较两法结果有无差别，要着眼于两法结果不一致的部分。表中观察变量是对子中两法的差值或差别，由,b,和,c,两格数据来反映，总体中与,b,和,c,对应的数据可用,B,和,C,表示（,a,格和,d,格表示两法差值为,0,，不予考虑）。,当,40,时，,当,40,时，需作连续性校正：,例题,某实验室分别用乳胶凝集法和免疫荧光法对,58,名可疑系统红斑狼疮患者血清中抗核抗体进行测定。问两种方法的检测结果有无差别？,行,列表资料的,2,检验,前面介绍了两个样本率比较的,2,检验方法，其基本数据有,2,行,2,列，称为,2,2,表或四格表资料。本节介绍的行,列表资料的,2,检验，用于多个样本率的比较、两个或多个构成比的比较以及双向有序或无序分类资料的检验等。,R,C,表可以分为双向无序、单向有序、双向有序属性相同和双向有序属性不同,4,类。,双向无序,RC,表,R,C,表中两个分类变量皆为无序分类变量，对于该类资料：若研究目的为多个样本率（或构成比）的比较，可用行,列表资料的,2,检验；若研究目的为分析两个分类变量之间有无关联性以及关系的密切程度时，可以用行,列表资料的,2,检验以及,Pearson,列联系数进行分析。,例题,1,某医师研究物理疗法、药物治疗和外用膏药三种疗法治疗周围性面神经麻痹的疗效。问三种疗法的有效率有无差别？,例题,2,测得某地,5801,人的,ABO血型和,MN,血型结果，问两种血型系统之间是否有关联？,单向有序,RC,表,有两种形式。一种是,R,C,表中的分组变量是有序的，而指标变量是无序的。此种单向有序,R,C,表资料可用行,列表资料的,2,检验进行分析。另一种情况是,R,C,表中的分组变量是无序的，而指标变量是有序的，此种单向有序,R,C,表资料宜用秩和检验进行分析。,例题,某地城市与农村高血压患者严重程度情况，试比较该地城市和农村高血压患者高血压严重程度是否有差别？,双向有序,RC,表,双向有序属性相同的,R,C,表,R,C,表中的两分类变量皆为有序且属性相同。实际上是,2,2,配对设计的扩展，此时宜用一致性检验（或称,Kappa,检验）。,双向有序属性不同的,R,C,表,R,C,表中的两分类变量皆为有序且属性不相同。对于该类资料，需要分析两有序分类变量间是否存在线性变化趋势，宜用有序分组资料的线性趋势检验。,行列均为顺序变量的相关检验,变量虽然是有序的，但毕竟还不是定量的，需要给有序变量的各等级赋值方可进行相关分析。最简单的赋值法是按顺序赋给秩次,(,即得分,),，即给行变量的等级赋值,1,，,2,，,，,R,和给列变量的等级赋值,1,，,2,，,，,C,。这样,(X,Y),的不同取值就有,R,C,对，表中的,R,C,个频数就是这,R,C,对取值所对应的频数，然后计算,Spearman,秩相关系数，并作显著性检验，这是比较粗糙的分析方法。,行列均为顺序变量的相关检验,Spearman,秩相关分析比较粗糙，这是因为它给有序变量的等级赋值过于简单，不能最大限度地获得有序变量之间的相关信息。而典型相关分析是在使有序变量的相关达到极大的前提下给有序变量的各等级赋值，就是对于表的边缘（指,“,行合计,”,与,“,列合计,”,）设法产生一双变量正态，从而进行相关分析。因产生各等级的得分值的计算过程中涉及矩阵运算，故手工计算较麻烦，将用,SAS,程序实现统计计算。一旦有了各等级的得分值后，就可运用求维频数资料相关系数的公式计算了。典型相关分析可以得出几个典型相关系数,R,，与每一个典型相关系数对应的检验为卡方检验，其数值为,nR2,，其中,n,为,R,C,表中的总频数，与第,k,个卡方值对应的自由度,dfk=R+C-2k-1,。,行列均为顺序变量的相关检验,与,定量资料,的相关与回归分析类似，也可对,R,C,表资料中有序变量之间是否存在线性趋势作显著性检验，通常这两种检验的结果是基本一致的，即相关分析结果显著，线性趋势检验结果也显著。可以通过检验实现线性趋势检验，即利用回归分析思想产生的与线性回归有关的卡分量和偏离线性回归的卡分量。,行列均为顺序变量的相关检验,首先计算,R,C,表的值，然后将总的值分解成线性回归分量与偏离线性回归分量。若两分量均具有统计学意义，说明两分类变量存在相关关系，但关系不是简单的直线关系；若线性回归分量有统计学意义，偏离线性回归分量无统计学意义时，说明两分类变量不仅存在相关关系，而且是线性关系。,行列均为顺序变量的相关检验,计算 总的值,计算线性回归分量 回归,b,为回归系数,计算偏线性回归分量,将总的 值分解成线性回归分量和偏离线性回归分量。若两分量均有统计学意义，说明两分类变量存在相关关系，但关系不是简单的直线关系；若线性回归分量有统计学意义，偏离线性回归分量无统计学意义时，说明两变量不仅存在相关关系，而且是线性关系。,例题,某研究者欲研究年龄与冠状动脉粥样硬化等级之间的关系，将,278,例尸解资料整理成行乘列表，问年龄与冠状动脉粥样硬化等级之间是否存在线性变化趋势？,分层资料的分析,按一个或多个因素分层后，研究行变量与列变量间的联系。,通过控制分层变量的影响，检验研究行变量与列变量的关系。,例题,某研究者欲研究男女性使用别嘌 呤是否引发皮疹。,本章小节,2,检验是用途很广的一种假设检验方法，这里我们主要学习它在分类资料统计推断中的应用，包括：两个率或两个构成比比较的卡方检验；多个率或多个构成比比较的卡方检验以及分类资料的相关分析等。,本章第一节介绍了四格表的普通,2,检验方法，理论频数与实际频数的定义，,2,检验统计量的计算方法，,2,分布以及,2,界值表确定,P,值的方法。说明了四格表,2,检验的应用条件，以及连续性校正公式和四格表专用公式。通过,SAS,分析实例演示了四格表,2,检验过程以及结果解释。我们需要重点理解,2,检验的含义以及应用。掌握,2,检验的应用条件以及连续性校正公式等。,本章小节,R,C,表可以分为双向无序、单向有序、双向有序属性相同和双向有序属性不同,4,类。,双向无序,R,C,表，,R,C,表中两个分类变量皆为无序分类变量。若研究目的为多个样本率（或构成比）的比较，可用行,列表资料的,2,检验；若研究目的为分析两个分类变量之间有无关联性以及关系的密切程度时，可以用行,列表资料的,2,检验以及,Pearson,列联系数进行分析。,单向有序,R,C,表有两种形式：一种是,R,C,表中的分组变量是有序的，而指标变量是无序的，此种单向有序,R,C,表资料可用行,列表资料的,2,检验进行分析；另一种情况是,R,C,表中的分组变量是无序的，而指标变量是有序的，此种单向有序,R,C,表资料宜用秩和检验进行分析。,本章小节,双向有序属性相同的,R,C,表，,R,C,表中的两分类变量皆为有序且属性相同，实际上是,2,2,配对设计的扩展，此时宜用一致性检验（或称,Kappa,检验）；双向有序属性不同的,R,C,表，,R,C,表中的两分类变量皆为有序且属性不相同，对于该类资料，需要分析两有序分类变量间是否存在线性变化趋势，宜用有序分组资料的线性趋势检验。,本章针对不同类型的,R,C,表资料分别予以介绍，并通过实例进行数据分析，并通过,SAS,程序完成统计分析过程，大家在练习不同类型数据的分析方法时，掌握常用数据分析的,SAS,程序，并能够灵活运用。,