《2_2列联表资料的统计分析与SAS软件实现.pdf》由会员分享，可在线阅读，更多相关《2_2列联表资料的统计分析与SAS软件实现.pdf（5页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、Correspondence:Prof.Liang2ping HU;Tel:010266931130;E2mail:lphu812 Medical Statistics 医学统计学22列联表资料的统计分析与SAS软件实现胡良平,王琪军事医学科学院生物医学统计学咨询中心,北京100850关键词:统计学;医学;数据分析,统计;定性资料;SAS软件Hu LP,Wang Q.J Chin Integr Med.2009;7(7):6782682.Received May 31,2009;accepted June 10,2009;published online July 15,2009.Indexe
2、d/abstracted in and full text link2out at PubMed.Journal title in PubMed:Zhong Xi Yi J ie He Xue Bao.Free full text(HTML and PDF)is available at .Forward linking and reference linking via CrossRef.DOI:10.3736/jcim20090715Open AccessStatistical analysis using SAS soft ware package for dataof22conting
3、ency tableLiang2ping HU,Qi WANGConsulting Center of Biomedical Statistics,Academy of Military Medical Sciences,Beijing 100850,ChinaKeywords:statistics;medicine;data analysis,statistical;qualitative data;SAS software 定性资料常以列联表形式呈现出来,通常分为22表、2C表和R2表、RC表和高维列联表。本文结合科研中的部分实例,介绍4种22列联表资料并用SAS实现其资料的统计分析。1
4、横断面研究设计的22表资料 横断面研究设计的22表资料是指在某个时间断面(时点或很短时间内)进行调查或实验研究时,对一组受试对象同时按两个定性变量来划分,每个定性变量都只有两个水平,通常一个为原因变量,另一个为结果变量,其目的是了解某个时间点的现状1。例1 多发伤的诊断对临床医师来说是最棘手的问题之一,容易发生漏诊。将医师按资历分为两组,比较两组医师的骨折漏诊情况,资料见表12。表1 接诊医师的资历与骨折漏诊情况医师资历例数漏诊情况:漏诊未漏诊合计3年以下342512853年以上31381412合计65632697 分析资料所需要的SAS程序:data y;do A=1 to 2;do B=1
5、 to 2;inputF;output;end;end;cards;34 25131 381;proc freq;weightF;tables A3B/CHISQ;run;运行结果:AB表的统计量统计量自由度值概率卡方13.866 7 0.049 3似然比卡方13.800 5 0.051 2连续校正卡方13.363 3 0.066 7Mantel2Haenszel卡方13.861 2 0.049 4Phi系数0.074 5列联系数0.074 3Cramer V统计量0.074 5876中西医结合学报2009年7月第7卷第7期 Journal of Chinese Integrative Med
6、icine,July 2009,Vol.7,No.7Fisher精确检验单元格(1,1)频数(F)34左侧Pr=F0.034 1表概率(P)0.015 6双侧Pr=P0.063 0 由于2=3.866 7,P=0.049 3,似乎表明两组不同资历医师的骨折漏诊情况不同,即3年以下医师资历的接诊医师的漏诊率(11.9%)大于3年以上医师资历的接诊医师(7.5%)。但这样下结论是不够科学严谨的!因为2检验所依赖的理论依据是2分布,而此分布原本是用来描述连续性随机变量变化规律的,用来分析定性资料纯属一种近似计算,对资料有一些特殊要求。本例虽然总样本量大于40,也没有小于5的理论频数,但计算所得到的检
7、验统计量值刚超过最低临界值3.841,此时,应以校正2检验结果或Fisher精确检验结果为准,即校正 2=3.363 3,P=0.066 7;或采用Fisher精确检验结果P=0.063。对本例而言,正确的专业结论是:虽然3年以下医师资历的接诊医师的漏诊率(11.9%)大于3年以上医师资历的接诊医师(7.5%),但二者之间的差异没有统计学意义,可认为其漏诊率基本相同。2 队列研究设计的22表资料 队列研究设计的22表资料是依据专业知识,通过对不同暴露水平的对象进行追踪观察,确定其疾病发生情况,从而分析暴露因素与疾病发生之间的因果关系。它是由因溯果的分析性研究设计,是论证疾病的可疑原因是否的确为
8、致病原因的一种重要研究方法。例2为探讨保守治疗掌侧Barton骨折的指征和预后,保守治疗23例掌侧Barton骨折患者,按复位满意与否分为两组,追踪观察7个月,观察术后并发症即创伤性关节炎的发生情况,资料见表23。表2 创伤性关节炎与复位结果的比较组别例数创伤性关节炎 无创伤性关节炎合计复位满意3912复位不满意10111合计131023 分析资料所需要的SAS程序:data y;do A=1 to 2;do B=1 to 2;inputF;output;end;end;cards;3 910 1;proc freq;weightF;tables A3B/CHISQ CMH;run;运行结果:
9、AB表的统计量统计量自由度值概率卡方110.144 90.001 4似然比卡方111.294 30.000 8连续校正卡方17.640 20.005 7Mantel2Haenszel卡方19.703 80.001 8Phi系数-0.664 1列联系数0.555 3Cramer V统计量-0.664 125%的单元格的期望计数比5小,卡方可能不是有效检验。Fisher精确检验单元格(1,1)频数(F)3左侧Pr=F0.999 9表概率(P)0.002 1双侧Pr=P0.002 8 因n 34170111 3151122273合计192192384 分析资料所需要的SAS程序:data y;do
10、A=1 to 2;do B=1 to 2;inputF;output;end;end;cards;41 70151 122;proc freq;weightF;tables A3B/CHISQ CMH;run;运行结果:AB表的统计量统计量自由度值概率似然比卡方110.751 60.001 0连续校正卡方19.934 90.001 6Mantel2Haenszel卡方110.629 40.001 1Phi系数-0.166 6列联系数0.164 3Cramer V统计量-0.166 6Fisher精确检验单元格(1,1)频数(F)41左侧Pr=F0.999 7表概率(P)4.31810-4双侧P
11、r=P0.001 6 2=10.657 2,P=0.001 1,说明两行的频数分布不同,也就是说看电视距离远近与是否患原发性闭角型青光眼可能有关(因为是回顾性调查研究,因果关系并非十分肯定)。采用Fisher精确检验得到的结果,P=0.001 6,结论同上。AB表的汇总统计量如下。Cochran2Mantel2H aenszel统计量(基于表得分)统计量对立假设自由度值概率1非零相关110.629 40.001 12行均值得分差值110.629 40.001 13一般关联110.629 40.001 1普通相对风险的估计值(行1/行2)研究类型方法值95%置信限案例对照Mantel2Haens
12、zel0.473 2 0.300 7 0.744 6(优比)Logit0.473 2 0.300 7 0.744 6CohortMantel2Haenszel0.667 8 0.512 1 0.870 8(第1列风险)Logit0.667 8 0.512 1 0.870 8CohortMantel2Haenszel1.411 2 1.162 2 1.713 5(第2列风险)Logit1.411 2 1.162 2 1.713 5 关于优势比(odds rate,OR)是否等于1的假设检验结果为2=10.629 4,P=0.001 1,说明总体OR1。因为本例属于病例对照研究设计,看上表“普通相
13、对风险的估计值(行1/行2)”前两行的结果,OR=0.473 2,其95%置信区间为0.300 7,0.744 6。专业结论:因OR=0.473 2,其95%置信区间为0.300 7,0.744 6,且对于“H0:总体OR=1”的假设检验结果为P 3 m对患原发性闭角型青光眼具有一定的保护作用;换句话说,看电视距离 3 m086中西医结合学报2009年7月第7卷第7期 Journal of Chinese Integrative Medicine,July 2009,Vol.7,No.7的OR=1/0.473 2=2.113 3,即看电视距离 3 m大2.113 3倍。4 配对研究设计的22表
14、资料 配对研究设计的22表资料是指按照配对原则分别接受两种不同的处理方法,每种处理方法的结果都可分为阳性和阴性两种,数出两种处理方法同时判定为阳性和阴性的频数以及它们结果不一致的频数,将结果列成配对设计的22表形式。例4 为探讨脑钠肽(brain natriuretic peptide,BNP)在脓毒症血瘀证诊断中的临床价值,以1988年北京血瘀证国际会议制订的血瘀证诊断参考标准为金标准,将174例脓毒症患者分为血瘀证组和非血瘀证组。分别检测患者血清BNP水平,BNP 100 ng/L为血瘀证组,结果见表45。表4 脓毒症BNP与血瘀证金标准诊断结果BNP法人数金标准:血瘀证非血瘀证合计血瘀证
15、8635121非血瘀证94453合计9579174 因为1988年北京血瘀证国际会议制订的血瘀证诊断参考标准是金标准,故该资料为特设金标准的配对设计四格表资料。若原文的目的是为了比较两种检测方法不一致部分的频数之间的差别是否具有统计学意义,应选用配对设计定性资料的2检验,即McNemar2检验;若原文的目的是为了检验两种方法的诊断结果是否一致,则应采用一致性检验,即Kappa检验。分析资料所需要的SAS程序:data y;do A=1 to 2;do B=1 to 2;inputF;output;end;end;cards;86 359 44;proc freq;weightF;tables
16、A3B;test Kappa;run;变量A、B分别代表BNP法检测的结果和根据金标准检测的结果,1代表阳性,2代表阴性。运行结果:McNemar检验统计量(S)15.363 6自由度1Pr SZ|Z|0.000 1 这是用McNemar2检验计算的结果,2=15.363 6(此处使用的是未校正的计算公式),P0.000 1,说明两种方法在检验脓毒症血瘀证方面的差别有统计学意义。Kappa=0.475 4(反映两种方法检测结果一致程度大小的统计量),其渐近标准误EKappa=0.064 7,总体Kappa的95%置信区间为0.348 6,0.602 2。“H0:总体Kappa=0”的假设检验结
17、果为Z=6.595 9,P 0.000 1,说明两种方法之间的检验一致率与期望一致率之间的差别有统计学意义,即两种方法的检验结果具有较高的一致性。专业结论之一:经差异性检验可知,两种检测方法之间不一致部分的差别有统计学意义,说明BNP法检测结果假阳性率很高。专业结论之二:经一致性检验可知,BNP法与金标准在检测脓毒症血瘀证方面具有较高的一致性。至于能否用BNP检测法取代金标准,需要提出专业要求,即BNP法检测结果与金标准的正确率应达到多大数值,本例中的正确率为74.71%(86+44)/174100%,相对来说,还是比较低的,可能不具有多大的实际意义。REFERENCES1Hu LP.Desi
18、gn of scientific researches and statisticalanalysis in dentistry.Beijing:Peoples Military MedicalPress.2007:1522163.Chinese.胡良平.口腔医学科研设计与统计分析.北京:人民军医出版社.2007:1522163.186中西医结合学报2009年7月第7卷第7期 Journal of Chinese Integrative Medicine,July 2009,Vol.7,No.72Li Q,Wang CQ,Liu G,Li K,Deng J,LuoKJ.Analysis for
19、 causes ofmissed diagnosis of fractureduring treatment ofmultiple trauma.ZhongguoGuShang.2008;21(2):89290.Chinese.李青,王春庆,刘钢,李昆,邓进,罗开俭.多发伤救治中骨折漏诊原因分析.中国骨伤.2008;21(2):89290.3Tang WJ,Wang MY,Gong XY,An GS.Clinicalinvestigation of conservative treatment for volar Bartonfracture.Zhongguo Gu Shang.2008;21
20、(5):3832385.Chinese with abstract in English.汤文杰,王满宜,贡小英,安贵生.保守治疗掌侧Barton骨折的临床探讨.中国骨伤.2008;21(5):3832385.4Zhuang X,Zhu RR,Guan HJ,Huang CH,Shi WP,Jiang SY.A case2control study of factors associatedwith primary angle2closure glaucoma.Zhonghua Yan KeZa Zhi.2008;44(6):5032506.Chinese with abstract inEn
21、glish.庄勋,朱蓉嵘,管怀进,黄春华,施文平,姜声扬.原发性闭角型青光眼相关因素的病例对照研究.中华眼科杂志.2008;44(6):5032506.5Chen Y,Zhu L,Qian YM.Diagnostic value of assay ofbrain natriuretic peptide in patients with sepsis of blood2stasis syndrome.J Chin Integr Med.2008;6(3):2662269.Chinese with abstract in English.陈越,朱亮,钱义明.脑钠肽检测对脓毒症患者血瘀证的临床诊断价
22、值.中西医结合学报.2008;6(3):2662269.2009全国中西医结合周围血管疾病学术交流会报到通知“2009全国中西医结合周围血管疾病学术交流会”拟定于2009年10月在郑州召开,同时召开中国中西医结合学会周围血管疾病专业委员会第七次换届会议。会议由中国中西医结合学会主办,周围血管疾病专业委员会承办,河南中医学院第一附属医院协办。会议主题为中西医结合周围血管疾病的新进展。现将会议有关事项通知如下。1 会议安排(1)会议日程:2009年10月30日报到,2009年10月31日11月1日开会;(2)报到地点:河南明珠大酒店(河南省郑州市金水区人民路6号,河南中医学院第一附属医院对面);(
23、3)会务费:每位参会者800元(含资料费、餐饮费等),住宿费回单位报销;(4)每位参会者授予国家级继续教育学分6分。2注意事项(1)由于住宿紧张,务必在10月15日前寄回会议回执表(以邮戳为准),以便安排。(2)由飞机场到会址乘车路线:乘机场大巴到民航大酒店下车,换乘52路公交车,在工人新村下车,马路对面向北100 m即到。需要接机者,请务必提前与会务组联系。(3)由火车站到会址乘车路线:由火车站乘坐1路电车(每位车费1元),到工人新村下车,向北100 m即到,会议不在车站设接待处。(4)有会议发言者,请提前做好发言准备。(5)第七届周围血管疾病专业委员会的委员候选人请务必参会,否则视为自动放弃委员资格。3联系方式 联系人:崔公让(手机:13323819068;E2mail:zhzhcgr ),崔炎(手机:037166598866;办公室电话:0371266212427),马海涛(手机:13903717598;E2mail:mahaitao6968 );联系地址:郑州市人民路19号河南中医学院第一附属医院;邮政编码:450000。中国中西医结合学会286中西医结合学报2009年7月第7卷第7期 Journal of Chinese Integrative Medicine,July 2009,Vol.7,No.7