流行病学常用多因素回归统计分析.pptx

EpidemiologyEpidemiology&Statistics&Statistics流行病学资料的多因素回归(hugu)分析郑卫军郑卫军 流行病学流行病学(li xn bn xu)(li xn bn xu)与卫生统计学教研室与卫生统计学教研室第一页，共七十八页。disease disease source source gene gene environment environment health health habit habit SES SES gender gender age age?Example of complexity of medical studyExample of complexity of medical study第二页，共七十八页。Example of complexity in medical studyExample of complexity in medical studyphysicalactivity eatingdrinkingsmokingsocial economic statusenduredglucosebloodpressureweightcholesteroldiabetescoronary heart diseaseenvironmentlife styleendogeneousfactorsdiseases第三页，共七十八页。流行病学多变量分析(fnx)工具回归多因素、多元线性回归多因素、多元logistic回归COX 比例(bl)风险模型主成分分析因子分析聚类判别分析路径分析结构方程模型第四页，共七十八页。本章(bn zhn)我们来交流点什么？多因素(yn s)线性回归多因素logistic回归多因素COX 回归第五页，共七十八页。本讲学习(xux)目标掌握(zhngw)多因素回归的技巧及注意点第六页，共七十八页。一、多因一果的回归(hugu)分析方法比较研究实验因素(yn s)与疗效实验组VS对照组：t检验、卡方、秩和、方差分析相关性研究相关分析回归分析第七页，共七十八页。多因素问题常见(chn jin)情况常见于观察性研究探讨多个因素作用模式(msh)分析由于分组不均衡造成存在着混杂因素实验性研究也可以遇到。基本的处理方法：多因素回归第八页，共七十八页。举例(j l)关于胰脏癌病人生存时间的研究。该研究的终点为死亡，包括很多可能影响生存的因素第九页，共七十八页。本例目标分析影响生存的多个因素，因素间地位(dwi)平等第十页，共七十八页。举例 关于(guny)胰脏癌病人生存时间的研究。该研究的终点为死亡，探讨术中放疗对生存的影响第十一页，共七十八页。本例目标探讨术中放疗(fn lio)较未放疗(fn lio)是否影响生存，其它为混杂因素，需要排除混杂干扰第十二页，共七十八页。COXCOX回归分析回归分析(fnx)(fnx)结果结果第十三页，共七十八页。多因素(yn s)回归多因素线性回归：因变量定量资料多因素logistic回归：因变量是分类资料（常见二分类）多因素Poisson回归：因变量是泊松分布资料多因素COX比例风险模型：因变量是时间生存(shngcn)资料多因素回归的自变量：可以是定量，也可以是分类资料第十四页，共七十八页。不同方法的应用(yngyng)场合病例对照研究(ynji)？队列研究？实验性研究？第十五页，共七十八页。不同方法的应用(yngyng)场合线性回归 队列/实验性研究(ynji)Logistic回归 病例对照/队列/实验性研究Poisson回归 队列COX模型 队列/实验性研究第十六页，共七十八页。思考思考(sko)(sko)：什么样的设计？采用什么样的方法：什么样的设计？采用什么样的方法？例1 某研究者跟踪随访已发病冠心病病人，研究病人的年龄age、性别sex（0为女性、1为男性）、心电图检验(jinyn)是否异常ecg（ST段压低、0为正常、1为轻度异常、2为重度异常）与冠心病ca复发是否有关。例1 某研究者希望研究病人的年龄age、性别sex（0为女性、1为男性）、心电图检验是否(sh fu)异常ecg（ST段压低、0为正常、1为轻度异常、2为重度异常）与冠心病ca发病是否有关。第十七页，共七十八页。横截面研究可以(ky)做多因素回归吗？横截面调查的数据疾病：心血管病、高血压相关因素：性别、年龄、红肉日均摄入量、吸烟、饮酒(yn ji)、高血糖、血管弹性指数第十八页，共七十八页。横截面研究可以(ky)做多因素回归吗？横截面调查的数据疾病：心血管病、高血压相关因素：性别(xngbi)、年龄、红肉日均摄入量、吸烟、饮酒、高血糖、血管弹性指数哪些因素可以做呢？第十九页，共七十八页。横截面研究可以做多因素(yn s)回归吗？横截面调查的数据疾病：心血管病、高血压相关因素(yn s)：性别、年龄、红肉日均摄入量、吸烟、饮酒、高血糖、血管弹性指数固有的因素，理论上和结局时间逻辑明确的因素可以开展多因素回归分析第二十页，共七十八页。多因素回归(hugu)的注意点样本量回归方程中，因素的筛选(shixun)多重共线性问题第二十一页，共七十八页。样本量线性回归：纳入回归模型(mxng)变量数的5-10倍Logistic回归：阳性数纳入模型变量数的5-10倍，总样本量至少乘以2样本量过少，则导致抽样误差增加，回归模型拟合效果变差第二十二页，共七十八页。因素(yn s)筛选1.常见方法单因素加多因素方法 即单因素回归P值0.05,0.1,0.2的纳入(nr)模型分析。此举是为了减少变量个数，提高拟合效果最优子集方程法（逐步回归法）向前法向后法逐步法单因素加多因素逐步回归法第二十三页，共七十八页。24（1）前进(qinjn)法（2）后退法（3）逐步法自变量从无到有、从少到多自变量从无到有、从少到多 先将全部自变量放入方程先将全部自变量放入方程(fngchng)(fngchng)，然后逐步剔除，然后逐步剔除 双向筛选：引入有意义的变量双向筛选：引入有意义的变量(binling)(binling)（前进法），剔除无意义（前进法），剔除无意义变量变量(binling)(binling)（后退法）（后退法）最优方程筛选第二十四页，共七十八页。25逐步回归分析的基本思想-拟合(n h)效果线性回归 保证(bozhng)纳入因素后，模型的R2值最大Logistic等回归方法 AIC值、SC 值或-2 log likelihood值越小好第二十五页，共七十八页。26因素(yn s)筛选注意事项如果(rgu)有变量的引入导致具有关键作用的因子其统计学意义发生改变（从无到有，或从有到无），必须发现原因！是否由于样本量过少，变量过多导致是否存在多重共线性是否属于混杂因素是否属于中间变量:回归模型不能同时纳入影响因素及其中间变量第二十六页，共七十八页。27多重共线性是指在进行多因素回归分析时，自变量间存在较强的线性相较强的线性相关关系。关关系。共线关系的存在，可使得估计系数方差加大，系数估计不稳，结果分析困难。因此在回归分析时，特别是当回归结果难以用专业知识解释时，要进行共线性诊断，找出存在共线性且不重要的那些自变量，剔出方程，另行回归分析。对于存在共线性的资料，可以利用共线性诊断有选择的保留自变量以消除共线性；或者采用岭回归岭回归、主成分主成分(chng fn)(chng fn)回归回归等回归分析方法以避免共线性指标对结果的影响。多重共线性问题(wnt)第二十七页，共七十八页。28多重共线性多重共线性的表现的表现(bioxin)(bioxin)在实际应用中主要表现为：在实际应用中主要表现为：（1）模型拟合效果很好，但偏回归系数几乎都无统计学意义；（2）偏回归系数估计值的方差(fn ch)很大；（3）偏回归系数估计值不稳定，随着样本含量的增减各偏回归系数发生较大变化或当一个自变量被引入或剔除时其余变量偏回归系数有很大变化；（4）偏回归系数估计值的大小与符号可能与事先期望的不一致或与经验相悖，结果难以解释 出现以上情况，提示存在出现以上情况，提示存在多重共线性多重共线性问题，应进行问题，应进行多重共线性多重共线性诊诊断。断。第二十八页，共七十八页。回归回归(hugu)(hugu)案例分析案例分析（logisticlogistic）某妇幼保健院采用病例对照研究探索乳腺增生的影响因素。随机抽取某年在本院就诊(ji zhn)的200名乳腺增生患者及200名非乳腺增生患者，采用问卷调查法对各种可能的影响因素进行调查，以探索乳腺增生的危险因素。第二十九页，共七十八页。回归回归(hugu)(hugu)案例分析案例分析第三十页，共七十八页。1 1：单因素：单因素(yn s)(yn s)分析分析回归回归(hugu)(hugu)案例分析案例分析第三十一页，共七十八页。单因素分析初步结论：由于年龄的影响在单因素分析中无统计学意义，且OR值也较为(jio wi)接近1，在多因素分析中不再考虑年龄因素，仅对其它4个因素进行分析。回归案例回归案例(n l)(n l)分析分析第三十二页，共七十八页。多因素(yn s)分析结果：回归案例回归案例(n l)(n l)分析分析第三十三页，共七十八页。多因素分析小结：为什么rs变得无统计学意义了？相关分析结果显示，妊娠(rnshn)次数与流产次数具有较强的相关性（r=0.55，P=25与初产年龄=25与初产年龄25岁的人相比，其效应（参数估计值）为0.2717+1.4587=1.7304。回归案例回归案例(n l)(n l)分析分析第三十八页，共七十八页。本例结论：流产次数多是乳腺增生的危险因素，初产年龄和是否(sh fu)母乳喂养的交互作用对乳腺增生的影响也有统计学意义。初产年龄25且非母乳喂养的人，其乳腺增生的危险明显增高。而单纯初产年龄25的人，或单纯非母乳喂养的人，发生乳腺增生的危险并不高。这提示可针对不同人群开展不同的干预措施。回归案例回归案例(n l)(n l)分析分析第三十九页，共七十八页。当混杂混入(hn r)流行病学匹配(ppi)分层多因素回归多因素回归匹配多因素回归倾向得分多因素回归倾向得分第四十页，共七十八页。回归中的混杂(hnz)因素混杂因素的识别：（1）该因素对结局有影响，可利用单因素回归或2检验来验证。（2）该因素在分析因素中的分布不均衡，可用2检验验证一下。（3）从专业角度来判断，即该因素不能是分析因素与结局关系的一个中间环节。也就是说，不能是分析因素引起该因素，通过该因素再引起结局。这一点主要根据(gnj)专业知识来确定。如果同时满足这三个条件，基本可以断定是混杂因素。如果有一条不满足，该因素就不是混杂因素。第四十一页，共七十八页。回归中的混杂回归中的混杂(hnz)(hnz)因素因素例2：分析性别、吸烟对幽门螺杆菌（HP）的影响，判断(pndun)吸烟是否是混杂因素性别性别吸烟吸烟状况状况HP阴性阴性HP阳性阳性男男不不吸烟吸烟146343吸烟吸烟64178女女不不吸烟吸烟211吸烟吸烟110530第四十二页，共七十八页。回归(hugu)中的混杂因素（1）判断吸烟是否对）判断吸烟是否对Hp有影响：单因素有影响：单因素logistic回归回归提示：吸烟对Hp的发生(fshng)有一定影响，吸烟发生(fshng)的危险更高一些第四十三页，共七十八页。回归中的混杂(hnz)因素（2）判断吸烟在不同性别）判断吸烟在不同性别(xngbi)中分布是否不同：卡方检验中分布是否不同：卡方检验提示：吸烟在性别中分布不同，男性中比例更高第四十四页，共七十八页。回归(hugu)中的混杂因素未校正吸烟因素未校正吸烟因素(yn s)(yn s)，性别分析结果，性别分析结果校正吸烟校正吸烟(x yn)(x yn)因素后，性别分析结果因素后，性别分析结果第四十五页，共七十八页。回归中的混杂(hnz)因素不吸烟的人之中，性别的分析不吸烟的人之中，性别的分析(fnx)(fnx)结果结果吸烟的人之中，性别的分析吸烟的人之中，性别的分析(fnx)(fnx)结果结果第四十六页，共七十八页。胰岛素治疗的胰岛素治疗的2 2型糖尿病患者型糖尿病患者(hunzh)(hunzh)中胰岛素和胰岛中胰岛素和胰岛素类似物的剂量与癌症之间的关系素类似物的剂量与癌症之间的关系混杂因素分析(fnx)案例第四十七页，共七十八页。研究(ynji)方法及设计l本研究(ynji)克服了以往相关研究(ynji)的不足1-4 -巢式病例对照研究 -多因素 logistic回归分析 -纳入每种胰岛素平均日用剂量作为变量l入选标准 -2型糖尿病;入选前5年未用过胰岛素;无恶性肿瘤史l患者平均随访75.9月，在此期间患者暴露于甘精胰岛素、地特胰岛素、赖脯胰岛素、门冬胰岛素、人胰岛素（基础、餐时）第四十八页，共七十八页。第四十九页，共七十八页。预先计划的比较(bjio)在肿瘤组和对照组之间比较(bjio)使用各种胰岛素的人数所占比例使用各种胰岛素的平均日剂量(Mean daily dose,MDD)使用各种胰岛素的平均日剂量MDD 0.3 IU/kg 的人数所占比例Mannucci et al.Diabetes Care.2010 Jun 14.Epub ahead of print第五十页，共七十八页。肿瘤风险(fngxin)与胰岛素日用剂量的关联性病例组甘精胰岛素平均日用病例组甘精胰岛素平均日用(ryng)(ryng)剂量较对照组高剂量较对照组高 （0.24 vs 0.16 IU/kg/day,p=0.036）人胰岛素与其他胰岛素类似物肿瘤发生与平均日用剂量无关人胰岛素与其他胰岛素类似物肿瘤发生与平均日用剂量无关 Mannucci et al.Diabetes Care.2010 Jun 14.Epub ahead of print第五十一页，共七十八页。肿瘤风险与胰岛素日用剂量(jling)的关联性每种胰岛素“平均日剂量0.3IU/day*kg”：“大剂量”经Chalrson同病指数(zhsh)、二甲双胍暴露、总胰岛素平均日剂量校正后纳入所有新发肿瘤病例及对照纳入所有新发肿瘤病例及对照除外随访中首除外随访中首12月内新发肿瘤病例及其对照月内新发肿瘤病例及其对照 Mannucci et al.Diabetes Care.2010 Jun 14.Epub ahead of print第五十二页，共七十八页。结论结论(jiln)(jiln)使用每种胰岛素的人数比例，肿瘤组和对照组之间差异无统计学意义。校正混杂因素之后，高平均日剂量(0.3 IU/kg)的甘精胰岛素与肿瘤风险有联系，其它(qt)胰岛素无此现象。当评价胰岛素及其类似物与肿瘤风险的可能联系时，必须考虑剂量。第五十三页，共七十八页。多因素(yn s)回归匹配常见logistic回归(hugu)非条件logistic：非匹配条件logistic回归：匹配第五十四页，共七十八页。条件(tiojin)logistic回归：1:1匹配，可以通过SPSS软件COX回归模块实施1:n匹配，则无法采用SPSS软件，而要SAS或者(huzh)其它软件实施第五十五页，共七十八页。倾向的分法是Rosenbaum和Rubin(1984)年首次提出的，其主要目的是均衡各对比组间各个特征变量的可比性。倾向得分：是指在一定协变量条件下，一个观察对象可能接受某种处理（或暴露）因素的可能性。倾向的分值相同的两个个体，其协变量的分布也趋于一致，即具有(jyu)相同的接受处理/对照的概率。倾向的分主要用来降低选择偏倚。多因素倾向(qngxing)得分第五十六页，共七十八页。倾向(qngxing)得分应用方法倾向得分倾向得分是协变量的一个综合指标，可以校正已测量的混杂因素(yn s)。可以作为协变量进行多调整可以作为匹配的因素开展或不开展多因素调整也可以用于分层。第五十七页，共七十八页。不同胰岛素促泌剂与二甲双胍在有或无心肌梗不同胰岛素促泌剂与二甲双胍在有或无心肌梗 塞史二型糖尿病患者中死亡率和心血管风险塞史二型糖尿病患者中死亡率和心血管风险 的比较的比较:丹麦丹麦(dn mi)(dn mi)研究研究多因素倾向(qngxing)得分第五十八页，共七十八页。数据来源数据来源:丹麦国家注册丹麦国家注册(zhc)(zhc)登记研究登记研究所有丹麦居民各有一个唯一、永久的注册号，与国内各种注册相连接.“丹麦国家处方注册”，记录了1995年来所有的处方.“全国患者注册”，1978年以来，所有的住院患者出院时记录了主要诊断（以及必要时,次要(cyo)诊断）.“全国死因注册”，记录死因信息。第五十九页，共七十八页。研究(ynji)人群丹麦居民年龄大于20岁，在1997年至2006年之间应用胰岛素促泌剂或者二甲双胍单药治疗者 总人数为107,806；随访时间：中位数3.3年 (最长9年)剔除(tch)了接受胰岛素单一治疗的8183例和多个药物联合治疗的3434例。入组人群分为有或无心肌梗死病史研究期间允许更改治疗药物 77%的患者在研究期间仅使用了一种药物治疗MI,myocardial infarction第六十页，共七十八页。研究(ynji)设计研究目的：评估与二甲双胍相比，常用的胰岛素促泌研究目的：评估与二甲双胍相比，常用的胰岛素促泌剂（包括瑞格列奈）的死亡率和心血管风险剂（包括瑞格列奈）的死亡率和心血管风险(fngxin)(fngxin)。大型、回顾性、非随机分组的观察性研究。大型、回顾性、非随机分组的观察性研究。对居住在丹麦、接受单一口服降糖药治疗的对居住在丹麦、接受单一口服降糖药治疗的T2DM患者患者的临床注册资料进行统计分析。的临床注册资料进行统计分析。Slide No 61Date第六十一页，共七十八页。无心梗史患者的各组基线无心梗史患者的各组基线(jxin)(jxin)数据数据二甲双胍二甲双胍格列美脲格列美脲格列格列齐特齐特格列本格列本脲脲格列格列吡嗪吡嗪甲苯甲苯磺磺丁脲丁脲瑞格瑞格列奈列奈N(%)43,340(54.3)36,313(37.0)5,926(6.0)12,495(12.7)6,965(6.1)5,335(5.4)2,513(2.6)年龄(年)52.5 14.060.9 13.360.0 13.263.2 13.763.013.564.413.557.912.6男性(%)50.955.356.554.454.153.856.0治疗时间(年)1.76 1.582.11 1.752.10 1.752.35 2.082.35 2.082.36 2.131.97 1.76各组基线各组基线(jxin)水平不均衡水平不均衡,不具可比性不具可比性.第六十二页，共七十八页。既往心梗史患者的各组基线既往心梗史患者的各组基线(jxin)(jxin)数据数据二甲双胍二甲双胍格列美脲格列美脲格列格列齐特齐特格列格列本脲本脲格列格列吡嗪吡嗪甲苯甲苯磺磺丁脲丁脲瑞格列瑞格列奈奈N(%)2,906(30.2)3,894(43.3)517(6.9)1168(12.2)660(7.3)501(5.6)186(2.1)年龄(年)65.8 10.770.9 11.070.5 10.970.9 11.070.5 10.471.2 11.468.2 10.3男性(%)73.170.369.370.069.764.970.4治疗时间(年)1.67 1.481.98 1.601.96 1.782.28 1.962.19 1.902.12 1.982.04 1.76各组基线水平不均衡各组基线水平不均衡(jnhng),不具可比不具可比性性.第六十三页，共七十八页。统计学方法统计学方法：倾向：倾向(qngxing)(qngxing)匹配匹配第六十四页，共七十八页。倾向得分倾向得分(d fn)(d fn)计算计算第六十五页，共七十八页。倾向性匹配倾向性匹配(ppi)(ppi)(propensity-matching)(propensity-matching)3.每一个个体都有一个每一个个体都有一个倾向性得分倾向性得分,与该个体出现在两与该个体出现在两组之中的概率组之中的概率P和和1-P有关。有关。4.根据促泌剂组中某患者根据促泌剂组中某患者倾向性得分倾向性得分,到服用二甲双胍到服用二甲双胍的人群中找一个的人群中找一个“对象对象”,两者两者“倾向性得分倾向性得分“最接近最接近,若最接近的候选者很多若最接近的候选者很多,可随机抽取可随机抽取(chu q)(chu q)一位。一位。5.事先规定得分相差多少算事先规定得分相差多少算接近接近；若二甲双胍的人群；若二甲双胍的人群中没有得分中没有得分“接近接近”的对象的对象,促泌剂组中的该成员就找促泌剂组中的该成员就找不到不到对象对象。第六十六页，共七十八页。第六十七页，共七十八页。基线基线(jxin)水平的均衡性有所改善。水平的均衡性有所改善。倾向性匹配倾向性匹配(ppi)后，各组人后，各组人数数第六十八页，共七十八页。第六十九页，共七十八页。EpidemiologyEpidemiology&Statistics&Statistics二甲二甲(r ji)双胍双胍格列美脲格列美脲格列齐特格列齐特格列苯脲格列苯脲格列吡嗪格列吡嗪甲苯磺丁脲甲苯磺丁脲瑞格列奈瑞格列奈1.7结结果果(ji(ji gu)gu)：无无心心梗梗患患者者各各组组总总死死亡亡率率第七十页，共七十八页。EpidemiologyEpidemiology&Statistics&Statistics1.7结果结果(ji gu)(ji gu)：既往心梗患者各组：既往心梗患者各组总死亡率总死亡率二甲二甲(r ji)双胍双胍格列美脲格列美脲格列齐特格列齐特格列苯脲格列苯脲格列吡嗪格列吡嗪甲苯磺丁脲甲苯磺丁脲瑞格列奈瑞格列奈第七十一页，共七十八页。结论结论(jiln)(jiln)第七十二页，共七十八页。结论结论(jiln)(jiln)第七十三页，共七十八页。其它(qt)要点：实验性研究的多因素方法严谨的实验性研究，RCT，分组均衡，基线资料经统计学比较，无统计学差异，混杂偏倚小实验因素和疗效之间的关系，t、卡方、秩和、K-M分析即可一些情况(qngkung)下，可能会存在：1.研究者设计、采集信息不严谨2.样本量不多3.失访导致分组不均衡，需要采用线性回归开展处理第七十四页，共七十八页。其他要点(yodin)：非独立性的数据无论是实验性研究、还是(hi shi)观察性研究，数据非独立性是导致结果失真的重要因素样本不是完全随机抽取，而是整群性多中心临床试验整群随机抽样研究数据存在着聚集性，个体观察值Y，存在关联性（不是X）第七十五页，共七十八页。对于(duy)流行病学家来说，是个挑战对于统计学家来说，是个机遇重复测量方差分析广义估计方程广义混合线性模型随机效益模型多水平模型。其他(qt)要点：非独立性的数据第七十六页，共七十八页。总结(zngji)：流行病学的分析方法首先，统计描述其次，初步的分析：t、F、卡方、秩和、相关再次，单因素回归(hugu)，掌握混杂因素的信息开展各式多因素回归第七十七页，共七十八页。内容(nirng)总结流行病学资料(zlio)的多因素回归分析。多因素logistic回归：因变量是分类资料(zlio)（常见二分类）。病例对照/队列/实验性研究。Logistic回归：阳性数纳入模型变量数的5-10倍，总样本量至少乘以2。（3）偏回归系数估计值不稳定，随着样本含量的增减各偏回归系数发生较大变化或当一个自变量被引入或剔除时其余变量偏回归系数有很大变化。这提示可针对不同人群开展不同的干预措施。（2）判断吸烟在不同性别中分布是否不同：卡方检验第七十八页，共七十八页。