Logistic回归分析

单击此处编辑母版标题样式,*,*,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,Logistic回归分析,1,LOGISTIC回归分析,因变量,Y,是两分类变量：,0,，,1,(,也可以是多分类或等级资料,),自变量可以是连续性和离散性定量资料或定性变量。,可以分析两变量的交互作用,Z=X1*X2,回归系数与比数比的关系,OR=,e,b,2,1.非条件Logistic回归分析（成组设计）,PROC,LOGISTIC,;,MODEL Y=X1-X5;,RUN,;,WEIGHT,指定变量为频数,Logistic过程后可加选项：,DES 指定Y=0为病例组水平,SIMPLE 输出各变量的简单统计描述结果,OUTEST= 可把参数估计值输出到指定数据集,NOPRINT 禁止输出任何结果,3,MODEL 语句后可加的选项：,/,STB,计算标准回归系数,CL,计算总体回归系数,95%,可信区间,NOINT,指定方程中无常数项,(a=0),RSQ,输出,R,2,SELECTION,指定逐步回归时筛选变量的方法,常用有：,STEPWISE,逐步法,SLE=SLS=0.05,FORWARD,前进法,SLE=0.05,BACKWARD,后退法,SLS=0.05,4,练习1：肾细胞癌转移有关因素分析,数据：unit2logistic1.xls,NO 病人编号,X1 确诊时病人年龄（岁）,X2 肾细胞癌血管内皮生长因子(分3级),X3 肾细胞癌组织内微血管数,X4 肾细胞癌组织学分级(由低到高共4级),X5 肾细胞癌分期(由低到高共4期),Y 肾细胞癌转移情况（有=1，无=0）,程序： unit2logistic.sas,5,练习2：Logistic回归预测,为分析新生儿低出生体重与支气管发育不良的关系，调查了233名低出生体重和新生儿，试作分析。,出生体重,组中值,患BPD人数,未患BPD人数,合计人数,550,750,49,19,68,950,1150,18,62,80,13501750,1550,9,66,75,程序： unit2example3.sas,6,练习3： 无序分类变量设,哑变量的方法,为研究冠心病发病与种族的关系，某研究所调查了,100个,样本，数据及变量、亚变量赋值方法如下表，试估计各种族间患冠心病的相对危险度。,程序： unit2example4.sas,7,练习4：不考虑交互作用的混杂因素控制,某研究所调查饮酒对食管癌的影响，定义因变量Y=1为食管癌，Y0为对照；数据如表4，试分析饮酒与食管癌的关系，并将年龄作为混杂因素。,模型1：仅包含自变量drink；,模型2：包含drink和age两个自变量，其中age取1，2，3；,模型3：仅包含age，定义3个（0，1）亚变量age1age3；,模型4：包含drink和age的3个亚变量age1age3 ；,程序： unit2example5.sas,8,练习5：,考虑交互作用的混杂因素控制,在一项回顾性研究中，调查了818个人的吸烟饮酒等生活方式以及各种慢性疾病的患病情况，表5是部分调查结果。试分析吸烟饮酒对患胃病的影响程度以及它们的交互作用。,程序： unit2example6.sas,9,练习6：多重共线性诊断,复旦大学金力实验室研究吸毒飘感（y)与某基因代谢通路10个基因（x1x10）的关系，X11-x13为性别、年龄和首吸年龄。,自变量之间不独立，引起分析结果无法解释,借用REG过程作共线性诊断的选项。,可以人为解决具有共线性自变量，或采用逐步回归分析的方法筛选变量。,数据： unit2example7.xls,程序： unit2example7.sas,10,练习2：Logistic回归预测,分析新生儿低出生体重与支气管发育不良的关系,程序： unit2example3.sas,练习3： 无序分类变量设哑变量的方法,研究冠心病发病与种族的关系,程序： unit2example4.sas,练习4：不考虑交互作用的混杂因素控制,调查饮酒对食管癌的影响，将年龄作为混杂因素,程序： unit2example5.sas,练习5：,考虑交互作用的混杂因素控制,分析吸烟饮酒对患胃病的影响以及它们的交互作用,程序： unit2example6.sas,练习6：多重共线性诊断,研究吸毒飘感（y)与某基因代谢通路10个基因（x1x10）的关系，X11-x13为性别、年龄和首吸年龄。,数据： unit2example7.xls,程序： unit2example7.sas,11,2. 条件Logistic回归分析（配对设计）, LOGISTIC 过程,仅用于1：1配对,先计算出每个变量病例与对照的差值,用差值作无截距的Logistic回归模型,DATA AA；SET B；Y=0;,Z1=X1-Y1；Z2=X2-Y2；,PROC,LOGISTIC,;,MODEL Y=Z1-Z7/NOINT;,RUN,;,12,练习7,胃癌危险因素1:1病例对照研究资料,X1X7病例的数据，Y1Y7对照的数据,X1/Y1 接触有害物质 （0，1）,X2/Y2 蛋白质摄入量（0，1，2，3，4，5）,X3/Y3 蔬菜水果摄入量（0，1，2）,X4/Y4 盐摄入量 （0，1，2，3，4，5）,X5/Y5 不良饮食习惯 （0，1，2，3，4）,X6/Y6 吸烟指数 （吸烟支数/日*吸烟年数/100）,X7/Y7 精神因素（0，1，2，3）,程序： unit2配对logistic1.sas,数据： unit2logistic2.xls,13,2. 条件Logistic回归分析（匹配设计）, PHREG 过程,可用于1：M配对（Y：病例=0，对照=1）,增加1个变量配对号 I,PROC,PHREG,;,MODEL Y=X1-X7/RL,TIES=DISCRETE;,指定作配对Logistic回归,STRATA I；,指定配对号作分层变量,RUN,;,14,练习7,胃癌危险因素1:1病例对照研究资料,X1 接触有害物质 （0，1）,X2 蛋白质摄入量（0，1，2，3，4，5）,X3 蔬菜水果摄入量（0，1，2）,X4 盐摄入量 （0，1，2，3，4，5）,X5 不良饮食习惯 （0，1，2，3，4）,X6 吸烟指数 （吸烟支数/日*吸烟年数/100）,X7 精神因素（0，1，2，3）,ID 配对号,Y 病例=0，对照=1,程序： unit2配对logistic2.sas,15,