1、 第十六章Logistic回归分析 Logistic regression1复习多元线性回归(multiple linear regression)在在医医学学实实践践中中,常常会会遇遇到到一一个个应应变变量量与与多多个个自自变变量量数数量量关关系系的的问问题题。如如医医院院住住院院人人数数不不仅仅与与门门诊诊人人数数有有关关,而而且且可可能能与与病病床床周周转转次次数数,床床位位数数等等有有关关;儿儿童童的的身身高高不不仅仅与与遗遗传传有有关关还还与与生生活活质质量量,性性别别,地地区区,国国别别等等有有关关;人人的的体体表表面面积积与与体体重重、身高等有关。身高等有关。2 表 多元线性回归
2、分析的数据结构实验对象 y X1 X2 X3 .XP 1 y1 a11 a12 a13 a1p 2 y2 a21 a22 a23 a2p 3 y3 a31 a32 a33 a3p n yn an1 an2 an3 anp 其中:y取值是服从正态分布3多元线性回归模型多元线性回归模型 通通过过实实验验测测得得含含有有p p个个自自变变量量x1,x2,x3,x1,x2,x3,xpxp及及一一个个因因变变量量y y的的n n个个观观察察对对象象值值,利利用用最最小小二二乘乘法法原理原理,建立多元线性回归模型建立多元线性回归模型:其其中中b b0 0为为截截距距,b1 b1,b2,b2 bpbp称称为
3、为偏偏回回归归系系数数.bibi表表示示当当将将其其它它p-1p-1个个变变量量的的作作用用加加以以固固定定后后,XiXi改变改变1 1个单位时个单位时Y Y将改变将改变bibi个单位个单位.4 logistic regression analysis5(一)基本概念和原理(一)基本概念和原理 1.1.应用背景应用背景 LogisticLogistic回归模型是一种概率模型,适合于病例对照研究、随访研究和横断面研究,且结果发生的变量取值必须是二分的或多项分类。可用影响结果变量发生的因素为自变量与因变量,建立回归方程。6 设资料中有一个因变量y、p个自变量x1,x2,xp,对每个实验对象共有n次
4、观测结果,可将原始资料列成表2形式。2、LogisticLogistic回归模型的数据结构7 表2 LogisticLogistic回归模型的数据结构实验对象 y X1 X2 X3 .XP 1 y1 a11 a12 a13 a1p 2 y2 a21 a22 a23 a2p 3 y3 a31 a32 a33 a3p n yn an1 an2 an3 anp 其中:y取值是二值或多项分类8 表3 肺癌与危险因素的调查分析例号 是否患病 性别 吸烟 年龄 地区 1 1 1 0 30 0 2 1 0 1 46 1 3 0 0 0 35 1 30 0 0 0 26 1 注:是否患病中,0代表否,1代表是
5、。性别中1代表男,0代表女,吸烟中1代表吸烟,0代表不吸烟。地区中,1代表农村,0代表城市。9 表4 配对资料(1:1)对子号 病例 对照 x1 x2 x3 x1 x2 x3 1 1 3 0 1 0 1 2 0 3 1 1 3 0 3 0 1 2 0 2 0 10 2 2 2 0 0 0注:X1蛋白质摄入量,取值:0,1,2,3 X2不良饮食习惯,取值:0,1,2,3 X3精神状况,取值:0,1,2 10LogisticLogistic回归回归-Logistic-Logistic回归与回归与多重多重线性回归联系与区别线性回归联系与区别联系联系:用于分析多个自变量与一个因变量的关用于分析多个自变
6、量与一个因变量的关系,目的是矫正混杂因素、筛选自变量和更系,目的是矫正混杂因素、筛选自变量和更精确地对因变量作预测等。精确地对因变量作预测等。区别区别:线性模型中因变量为连续性随机变量,线性模型中因变量为连续性随机变量,且要求呈正态分布且要求呈正态分布.Logistic.Logistic回归因变量的回归因变量的取值仅有两个,不满足正态分布。取值仅有两个,不满足正态分布。113 3、Logistic回归模型l 令令:y=1 发病(阳性、死亡、治愈等)发病(阳性、死亡、治愈等)l y=0 未发病(阴性、生存、未治愈等)未发病(阴性、生存、未治愈等)l 将发病的概率记为将发病的概率记为P,它与自变量
7、它与自变量x x1 1,x x2 2,x xp p之间的之间的Logistic回归模型为:回归模型为:l可知,不发病的概率为:可知,不发病的概率为:l 12经数学变换得:定义:为Logistic变换,即:13 4、回归系数i的意义 流行病学的常用指标优势比(odds ratio,OR)或称比数比,定义为:暴露人群发病优势与非暴露人群发病优势之比。即Xi的优势比为:14故对于样本资料OR=exp()95%置信区间为:可见 是影响因素Xi增加一个单位所引起的对数优势的增量,反映了其对Y作用大小。如果要比较不同因素对Y作用大小,需要消除变量量纲的影响,为此计算标准化回归系数155.假设检验(1)回归
8、方程的假设检验 H0:所有 H1:某个 计算统计量为:G=-2lnL,服从自由度等于n-p的 分布 (2)回归系数的假设检验 H0:H1:计算统计量为:Wald,自由度等于1。16(二)Logistic回归类型及其实例分析 1、非条件Logistic回归 当研究设计为队列研究、横断面研究或成组病例对照研究时,可以用非条件Logistic回归。17实例1 某研讨究者调查了某研讨究者调查了3030名成年人,记录名成年人,记录了同肺癌发病的有关因素情况了同肺癌发病的有关因素情况,数据见表数据见表4 4。其中是否患病中。其中是否患病中,0,0代表否代表否,11代表是;性别中代表是;性别中 11代表男代
9、表男,00代表女;吸烟中代表女;吸烟中 11代表吸烟代表吸烟,00代表不吸烟;地区中代表不吸烟;地区中,1,1代表农代表农村村,0,0代表城市。试分析各因素与肺代表城市。试分析各因素与肺癌间的关系。癌间的关系。18 表5 肺癌与危险因素的调查分析例号 是否患病 性别 吸烟 年龄 地区 1 1 1 0 30 0 2 1 0 1 46 1 3 0 0 0 35 1 30 0 0 0 26 1 注:是否患病中,0代表否,1代表是。性别中1代表男,0代表女,吸烟中1代表吸烟,0代表不吸烟。地区中,1代表农村,0代表城市。19SPSS操作步骤:Analyze-Regression-Binary Logi
10、stic-Dependent框(y)-Covariates框(x1,x2,)-ok非条件Logistic回归SPSS操作步骤:202122 结果表明,性别,吸烟,年龄三个因素都与肺癌有关.由于在对某一因素进行单因素分析时没有控制其它因素的干扰,因此结果不可靠.单因素分析的结果23多因素分析的结果2425 由上最大似然估计分析知因素X2(吸烟),X3(年龄)对肺癌的发生有影响。所得的回归方程为:Logit(P)=-9.7544+2.5152X1+3.9849X2+0.1884X3-1.3037X4.262.条件logistic回归分析配对设计的类型:1:1、1:m、n:m(可采用分层COX模型来
11、拟合)。例如:某市调查三种生活因素与胃癌的关系,资料见表5。27 表6 配对资料(1:1)对子号 病例 对照 x1 x2 x3 x1 x2 x3 1 1 3 0 1 0 1 2 0 3 1 1 3 0 3 0 1 2 0 2 0 10 2 2 2 0 0 0注:X1蛋白质摄入量,取值:0,1,2,3 X2不良饮食习惯,取值:0,1,2,3 X3精神状况,取值:0,1,228配对Logistic回归SPSS操作步骤:Analyze-Survival-COX Regression-Time框(outcome)-Status框(Status)-Define Event:Single value 1:
12、continue-Covariates框(x1、x2、x3)-Strata框(id)-Options-at last step-ok29303132333435363、逐步Logistic回归分析(1 1)向前法)向前法(forward selection)(forward selection)开始方程中没有变量,自变量由少到多一个一个引入回归方程。按自变量对因变量的贡献(P值的大小)由小到大依次挑选,变量入选的条件是其P值小于规定进入方程的P界值Enter,缺省值 P(0.05)。37(2)后退法(backward selection)开始变量都在方程中,然后按自变量因变量的贡献(P值的大小
13、)由大到小依次剔除,变量剔除的条件是其P值小于规定的剔除标准Remove,缺省值 p(0.10)。38(3)逐步回归法逐步引入-剔除法(stepwise selection)前进前进逐步引入-剔除法 是在前进法的思想下,考虑剔除变量,因此有两个p界值Enter,Remove。39 无论是条件还是非条件Logistic回归,在多变量分析时均可以采用逐步回归方法,实现的方法是:在 model后加选项:/selection=forward、backward或stepwise再给出SLE,SLS的界值。调试法:P从大到小取值0.5,0.1,0.05,一般实际用时,SLE,SLS应多次选取调整40解 释
14、 设第i个因素的回归系数为bi,表示当有多个自变量存在时,其它自变量固定不变的情况下,自变量Xi每增加一个单位时,所得到的优势比的自然对数。也就是其它自变量固定不变的情况下,自变量Xi每增加一个单位时,影响因变量Y=0发生的倍数。当bi0时,对应的优势比(odds ratio,记为ORi):ORi=exp(bi)1,说明该因素是危险因素;当bi0时,对应的优势比ORi=exp(bi)1,说明该因素是保护因素。41 在本例中不良饮食习惯X2,取值:0、1,2、3,b=1.2218,OR=3.393,表示不良饮食习惯是危险因素;精神状况X3,取值:0、1、2,b=2.2947,OR=9.922,表
15、示精神状况是危险因素。42 SPSS无论是条件还是非条件Logistic回归,在多变量分析时均可以采用逐步回归方法,实现的方法是:在 method后加选项:Enter:所有变量一次全部进入方程。Forward:逐步向前法 Backward:后退法 变量移出方程所采取的检验方法:Conditional;LR;Ward(尽量不用)4344二分类二分类LogisticLogistic回归回归 method method中文名称中文名称剔除依据剔除依据EnterEnter全部进入全部进入Forward:conditForward:conditionalional向前逐步向前逐步条件参数估计似然比条件参
16、数估计似然比Forward:LRForward:LR向前逐步向前逐步最大偏似然估计似然比最大偏似然估计似然比Forward:WaldForward:Wald向前逐步向前逐步WaldWald统计量统计量Backward:condiBackward:conditionaltional向后逐步向后逐步条件参数估计似然比条件参数估计似然比Backward:LRBackward:LR向后逐步向后逐步最大偏似然估计似然比最大偏似然估计似然比Backward:WaldBackward:Wald向后逐步向后逐步WaldWald统计量统计量45SPSS 软件实现方法 FileFileOpenOpen相应数据相应数据(已存在已存在)Analyze Regression RegressionBinary Binary LogisticLogistic Dependent(y)Dependent(y)Covariates Covariates (自变量)自变量)MothodMothod(逐步回归方法逐步回归方法)Specify model(Specify model(单因素或多因素)单因素或多因素)save(
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1