十六章logistic回归分析LogisticRegression.ppt

《十六章logistic回归分析LogisticRegression.ppt》由会员分享，可在线阅读，更多相关《十六章logistic回归分析LogisticRegression.ppt（40页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、十六章logistic回归分析LogisticRegression Still waters run deep.流静水深流静水深,人静心深人静心深 Where there is life,there is hope。有生命必有希望。有生命必有希望ContentnLogistic regression nConditional logistic regression nApplication 239讲述内容讲述内容:第一节第一节 logistic logistic回归回归第二节第二节 条件条件logisticlogistic回归回归第三节第三节 logistic logistic回归的应用回归的

2、应用 及其注意事项及其注意事项339目的：目的：作出以多个自变量（危险因素）估计作出以多个自变量（危险因素）估计应变量（应变量（结果因素）结果因素）的的logistic回归方程。回归方程。属于概率型非线性回归。属于概率型非线性回归。资料：资料：1.应变量为反映某现象发生与不发生的应变量为反映某现象发生与不发生的二值变量；二值变量；2.自变量宜全部或大部分为分类自变量宜全部或大部分为分类变量，可有少数数值变量。分类变量要数量变量，可有少数数值变量。分类变量要数量化。化。439用途：用途：研究某种疾病或现象发生和多个危研究某种疾病或现象发生和多个危 险因素（或保护因子）的数量关系。险因素（或保护因

3、子）的数量关系。用用 检验（或检验（或u检验）的检验）的局限性局限性：1.只能研究只能研究1个危险因素；个危险因素；2.只能得出定性结论。只能得出定性结论。539种类种类:1.成组（非条件）成组（非条件）logistic回归方程。回归方程。2.配对（条件）配对（条件）logistic回归方程。回归方程。639第一节第一节 logistic回归回归（非条件（非条件logistic回归回归）739 一、基本概念一、基本概念，在m个自变量的作用下阳性结果发生的概率记作个自变量的作用下阳性结果发生的概率记作:839若令若令:回回归归模模型型概率概率P：01，logitP：。取值范围取值范围 939图图

4、16-1 logistic函数的图形函数的图形1039模模型型参参数数的的意意义义常数常数项项 表示暴露剂量为表示暴露剂量为0时个体发病时个体发病与不发病概率之比的自然对数。与不发病概率之比的自然对数。回归系数回归系数 表示自变量表示自变量 改变一个单位时改变一个单位时logitP 的改变量。的改变量。1139流行病学衡量危险因素作用大小的流行病学衡量危险因素作用大小的比数比例比数比例指标。指标。计算公式为：计算公式为：优势比优势比OR(odds ratio)1239与与 logisticP 的关系的关系：13391439二、logistic回归模型的参数估计回归模型的参数估计 1.参数估计参

5、数估计2.2.原理：最大似然原理：最大似然(likelihood)估计估计 1539 2.优势比估计优势比估计 可反映某一因素两个不同水平（c1，c0）的优势比。1639例16-1 表16-1是一个研究吸烟、饮酒与食道癌关系的病例对照资料，试作logistic回归分析。确确定定各各变变量量编编码码 1739表表16-1 吸烟与食道癌关系的病例对照调查资料吸烟与食道癌关系的病例对照调查资料 1839经logistic回归计算后得:的95可信区间:1939 三、三、logistic回归模型的假设检验回归模型的假设检验 2.1.似然比检验2039方法：方法：前进法、后退法和逐步法前进法、后退法和逐步

6、法。检验统计量：检验统计量：不是不是 F 统计量，而是似然比统计量、统计量，而是似然比统计量、Wald 统计量和计分统计量之一。统计量和计分统计量之一。四、变量筛选四、变量筛选例16-2 为了探讨冠心病发生的有关危险因素，对为了探讨冠心病发生的有关危险因素，对26例冠心病病人和例冠心病病人和28例对照者进行病例例对照者进行病例 对照研究，各对照研究，各因素的说明及资料见表因素的说明及资料见表16-2和表和表16-3。试用。试用logistic 逐逐步回归分析方法筛选危险因素。步回归分析方法筛选危险因素。2139表表16-2 冠心病冠心病8个可能的危险因素与赋值个可能的危险因素与赋值2239表表

7、16-3 冠心病危险因素的病例冠心病危险因素的病例 对照调查资料对照调查资料 2339 表表16-4 例例16-2进入方程中的自变量及有关参数的估计值进入方程中的自变量及有关参数的估计值 学会看结果！学会看结果！24392539一、原理 第二节 条件logistic回归 2639表表16-5 1:M 条件条件logistic回归数据的格式回归数据的格式*t=0 为病例，其他为对照 2739条件 logistic 模型2839二、应用实例2939表表16-7 喉癌喉癌1:2配对病例配对病例 对照调查资料整理表对照调查资料整理表 P344:3039表16-8 例16-3进入方程中的自变量及有关参数

8、的估计值 采用逐步法 6个危险因素 变量筛选 4个进方程，结果见表16-9。3139一、logistic回归的应用1流行病学危险因素分析 logistic回归分析的特点之一是参数意义清楚，即得到某一因素的回归系数后，可以很快估计出这一因素在不同水平下的优势比或近似相对危险度，因此非常适合于流行病学研究。logistic回归既适合于队列研究(cohort study)，也适合于病例-对照研究(case-control study),同样还可以用于断面研究(cross-sectional study)第三节 logistic回归的应用及其注意事项32392临床试验数据分析 临床试验的目的大多是为了

9、评价某种药物或治疗方法的效果，如果有其他影响效果的非处理因素(如年龄、病情等)在试验组和对照组中分布不均衡，就有可能夸大或掩盖试验组的治疗效果。尽管在分组时要求按随机化原则分配，但由于样本含量有限，非处理因素在试验组和对照组内的分布仍有可能不均衡，需要在分析阶段对构成混杂的非处理因素进行调整。当评价指标为二值变量时(如有效和无效)，可以利用logistic回归分析得到调整后的药物评价结果。对于按分层设计的临床试验可以用相同的方法对分层因素进行调整和分析。33393分析药物或毒物的剂量反应 在一些药物或毒物效价的剂量-反应实验研究中，每一只动物药物耐受量可能有很大的不同，不同剂量使动物发生“阳性

10、反应”的概率分布常呈正偏态，将剂量取对数后则概率分布接近正态分布。由于正态分布函数与logistic分布函数十分接近，如果用P表示在剂量为X时的阳性率，可用下述模型表示它们之间的关系用这一模型可以求出任一剂量的阳性反应率用这一模型可以求出任一剂量的阳性反应率传传统统的的一一些些方方法法往往往往对对实实验验设设计计有有严严格格的的要要求求，如如剂剂量量按按等等比比级级数数排排列列，各各剂剂量量组组的的例例数数必必须须相相同同等等,采采用用logistic回回归归的的方方法法则则没有这些限制。没有这些限制。34394预测与判别 logistic回归是一个概率型模型，因此可以利用它预测某事件发生的概

11、率。例如在临床上可以根据患者的一些检查指标，判断患某种疾病的概率有多大。关于判别问题见第十八章。3539二、logistic回归应用的注意事项3639本章小结：本章小结：本章小结：本章小结：目的：目的：作出以多个自变量（危险因素）估计应变量（结果作出以多个自变量（危险因素）估计应变量（结果因素）的因素）的logistic回归方程。属于概率型非线性回归。回归方程。属于概率型非线性回归。资料：资料：1.应变量为反映某现象发生与不发生的二值变量；应变量为反映某现象发生与不发生的二值变量；2.自变量宜全部或大部分为分类变量，可有少数数自变量宜全部或大部分为分类变量，可有少数数值值 变量。分类变量要数量

12、化。变量。分类变量要数量化。用途：用途：研究某种疾病或现象发生和多个危险因素（或保护研究某种疾病或现象发生和多个危险因素（或保护因子）的数量关系。因子）的数量关系。种类种类:1.成组（非条件）成组（非条件）logistic回归方程。回归方程。2.配对（条件）配对（条件）logistic回归方程。回归方程。3739课后应用思考题：课后应用思考题：课后应用思考题：课后应用思考题：为了分析影响为了分析影响为了分析影响为了分析影响医院抢救急性心肌梗死（医院抢救急性心肌梗死（医院抢救急性心肌梗死（医院抢救急性心肌梗死（AMIAMIAMIAMI）患者能否成）患者能否成）患者能否成）患者能否成功的因素，某医

13、院收集了功的因素，某医院收集了功的因素，某医院收集了功的因素，某医院收集了5 5 5 5年中所有的年中所有的年中所有的年中所有的AMIAMIAMIAMI患者的抢救病史患者的抢救病史患者的抢救病史患者的抢救病史（有关因素很多，由于篇幅有限，本例仅列出（有关因素很多，由于篇幅有限，本例仅列出（有关因素很多，由于篇幅有限，本例仅列出（有关因素很多，由于篇幅有限，本例仅列出3 3 3 3个），共个），共个），共个），共200200200200例见下表。其中例见下表。其中例见下表。其中例见下表。其中P=0P=0P=0P=0表示抢救成功，表示抢救成功，表示抢救成功，表示抢救成功，P=1P=1P=1P=1表

14、示抢救未成功而死表示抢救未成功而死表示抢救未成功而死表示抢救未成功而死亡；亡；亡；亡；X1=1X1=1X1=1X1=1表示抢救前已发生休克，表示抢救前已发生休克，表示抢救前已发生休克，表示抢救前已发生休克，X1=0 X1=0 X1=0 X1=0表示抢救前未发生休表示抢救前未发生休表示抢救前未发生休表示抢救前未发生休克；克；克；克；X2=1 X2=1 X2=1 X2=1表示抢救前已发生心衰，表示抢救前已发生心衰，表示抢救前已发生心衰，表示抢救前已发生心衰，X2=0 X2=0 X2=0 X2=0表示抢救前未发生表示抢救前未发生表示抢救前未发生表示抢救前未发生心衰；心衰；心衰；心衰；X3=1 X3=

15、1 X3=1 X3=1表示患者从开始表示患者从开始表示患者从开始表示患者从开始AMIAMIAMIAMI症状到抢救时已超过症状到抢救时已超过症状到抢救时已超过症状到抢救时已超过12121212小时，小时，小时，小时，X3=0X3=0X3=0X3=0表示患者从开始表示患者从开始表示患者从开始表示患者从开始AMIAMIAMIAMI症状到抢救时未超过症状到抢救时未超过症状到抢救时未超过症状到抢救时未超过12121212小时。小时。小时。小时。请问最好采用哪种分析方法？为什么？请问最好采用哪种分析方法？为什么？请问最好采用哪种分析方法？为什么？请问最好采用哪种分析方法？为什么？分析结果有哪些？分析结果有哪些？分析结果有哪些？分析结果有哪些？3839AMIAMI患者的抢救危险因素资料患者的抢救危险因素资料患者的抢救危险因素资料患者的抢救危险因素资料P=0(在医院抢救成功)P=1(在医院抢救未成功而死亡)X1X2X3NX1X2X3N0003500040013400110010170104011190111510017100610161019110611061116111639394039