网站设计与建设实训,免费注册公司名字大全,沈阳做网站哪个好,互联网是指哪些工作本文内容来自《中华流行病学杂志》2019年第40卷第8期#xff0c;作者为冯国双教授#xff0c;原题目为《观察性研究中的logistic回归分析思路》。将这篇文章分享给医咖会的伙伴们#xff0c;希望大家能从领域大咖的见解中有所收获#xff0c;指导医学研究之路。#xff08…本文内容来自《中华流行病学杂志》2019年第40卷第8期作者为冯国双教授原题目为《观察性研究中的logistic回归分析思路》。将这篇文章分享给医咖会的伙伴们希望大家能从领域大咖的见解中有所收获指导医学研究之路。感谢冯国双教授的授权观察性研究在研究设计中占有非常重要的地位实际应用中比较常见的是病例对照研究和队列研究。尽管其应用广泛但在数据分析中却存在不少问题。在分析时往往只考虑数据本身而未能结合研究类型从而导致结果的偏倚。甚至在已发表的文章中也存在一些不严谨用语。本文从观察性研究的类型出发基于不同研究类型的研究目的以logistic回归分析为例探讨观察性研究的不同分析思路希望为医学科研工作者提供一定的参考和借鉴。logistic回归假定有m个自变量X1X2…Xmlogistic回归模型的基本形式可表达为只从数据本身考虑的话logistic回归模型都是包括一个分类因变量及若干自变量(可以是分类变量也可以是连续变量)反映了m个自变量对因变量的线性影响。无论对于病例对照研究还是队列研究这种形式都是不变的。部分研究在数据分析时忽略了前期的设计思路只是简单地把因变量和所有自变量纳入统计软件中相应位置点击运行直接给出结果。从数据上来看病例对照研究和队列研究的数据形式完全一样软件操作过程也并无不同都是指定因变量和自变量然后给出参数估计值及统计检验结果。统计软件无法判断研究者采用的是病例对照研究还是队列研究也并不清楚作者的主要研究目的是什么只是对指定的变量进行参数估计。而统计分析的思路需要根据研究目的和研究类型而定对于病例对照研究或队列研究而言它们的分析思路显然不同。一味依靠统计软件不仅容易出现一些错误分析思路也会导致错误的分析结果。病例对照研究中的logistic回归从数据分析的角度来看病例对照研究大致有两大类目的一是探索危险因素二是验证危险因素。1、以探索危险因素为目的的分析思路危险因素的探索常见于临床研究中通常用于研究初期此时研究者并不清楚哪些因素可能会影响结局的发生因此先进行初步探索。根据专业知识和经验收集一些可能的指标然后从中寻找可能对结局影响较大的因素。例如探索儿童打鼾的危险因素研究者并无太多的前期基础只是为了发现可能与儿童打鼾有关的因素这种情况下会根据文献报道、专业经验等收集一些可能有关的指标并从中找出与儿童打鼾有关的部分因素。危险因素探索的文章中最常见的表述错误是“校正其他混杂因素”后发现共K个变量对结局有影响。混杂因素是相对主要研究因素而言而危险因素探索的研究中并无明确的主要研究因素所有变量都是待研究的因素目的是从这些变量中找出哪些有影响。此类研究中“校正其他混杂因素”是一种不严谨的表达方式。对于这种分析思路需要有一定的分析经验和技巧。实际分析中需要考虑几个问题① 线性问题由于logistic回归本质上仍属于“线性模型”因此一定要确认自变量与因变量logit P之间是否存在线性关系如果不是需要考虑进行相应的变换否则可能会产生错误结果。例1某研究分析老年人高血压二分类变量是或否的危险因素研究因素包括gender、age、ox-LDL、Adiponectin、ox-LDL IgG 和 ox-LDL IgM 共6个指标。其中gender 为二分类变量其余变量均为连续变量。如果把6个自变量直接纳入统计软件分析所得结果见表 1。可以看出6个变量均差异无统计学意义。然而对数据重新分析后发现并不是这些变量对结局均无影响只是未能发现它们之间的真实关系而已。经仔细观察发现 age 和 ox-LDL IgM 对结局的影响是有统计学意义的但不是线性影响而是二次项关系表 2。② 共线性问题共线性即自变量之间存在高度相关从而导致结果不可靠 [1]。共线性是大多数回归模型都需要考虑的一个问题一旦发现该问题需要采取不同措施来解决。常见的解决方案包括删除某一自变量、主成分分析、Lasso回归等。例2某研究分析乳腺增生的危险因素自变量同时包括妊娠次数三分类变量用1、2、3表示相应次数和流产次数三分类变量用0、1、2表示相应次数。在单因素分析中妊娠次数差异有统计学意义2 vs. 1P0.0263 vs. 1P0.035然而多因素分析中则差异无统计学意义P值分别为0.635、0.594。分析原因发现主要是由于妊娠次数和流产次数有较强的共线性二者相关系数高达0.55从而导致妊娠次数变得无统计学意义。解决方案采用了删除法删除妊娠次数变量保留了流产次数变量。③ 单因素和多因素的问题目前危险因素筛选的一种分析思路先进行单因素分析将单因素分析中差异有统计学意义P0.05的变量再纳入多因素分析选出最终有统计学意义的变量作为危险因素。然而这一思路并非十分可靠有些情况下可能会出现单因素分析无统计学意义而多因素分析有统计学意义的情况此时就容易漏掉某些重要的因素。例3某研究分析两个血清学指标分别用阳性和阴性表示对胃癌的影响数据结果见表 3。该数据采用单因素分析的话可以发现 x1 差异无统计学意义P0.114而在多因素分析中却变得有统计学意义P0.018。如果只将单因素分析中有统计学意义的变量纳入多因素分析的话就会漏掉x1变量。为什么会出现这种情况主要是因为 x1 和 x2 之间存在负相关而 x1、x2 与结局之间均为正相关。因此数据分析过程中不要盲目套用所谓的“分析套路”而应结合实际情况具体问题具体分析。总之在筛选危险因素时建议不要仅将单因素分析有统计学意义的变量纳入多因素分析一定要厘清变量之间的关系否则容易遗漏重要的变量或纳入无意义的变量。2、以验证危险因素为目的的分析思路验证危险因素说明研究者在研究开始时已经有明确的主要研究因素主要目的是为了验证该因素是不是真正的影响因素。基于这种目的研究者在设计时会突出主要因素但同时也会收集其他可能的混杂因素。例如探索肺癌与吸烟的关系吸烟是主要研究因素因此问卷调查中会详细设置各种与吸烟有关的问题。考虑到其他因素可能也会影响肺癌发生因此调查时也会加入其他有关因素的调查但这些因素不是研究者关心的只是为了校正这些因素以便真正明确吸烟与肺癌的关系。因此对于这种研究目的关键的问题是如何控制混杂因素以便真正明确主要研究因素与结局的关系。混杂因素在流行病学中已有详细定义[2]不再赘述。从数据分析的角度来看要判断一个因素是否为混杂因素可以从两个方面来考虑第一分析该因素是否对结局有较大影响通常可采用χ2检验或单因素logistic回归来实现第二分析该因素在主要研究因素中的分布情况通常采用χ2检验来实现。例4某研究分析性别与幽门螺杆菌Hp的关系现在考虑吸烟是否为影响二者关系的混杂因素。具体数据见表 4。首先分析吸烟对结局的影响采用χ2检验或单因素logistic回归不难发现吸烟人群与不吸烟人群相比Hp阳性的风险更高OR1.8495%CI:1.44~2.35。其次分析吸烟在性别中的分布χ2检验结果显示男性和女性中吸烟的比例差异有统计学意义χ2396.97P0.001男性的吸烟比例远高于女性。由此看出以性别作为主要分析变量在分析性别与Hp感染时吸烟可能是影响二者关系的混杂因素必须加以校正。校正前结果显示性别对Hp的影响有统计学意义男性有更高的Hp阳性风险OR1.6295%CI:1.26~2.07校正后发现性别对Hp的影响无统计学意义(OR1.2695%CI:0.94~1.68)。因此对于以验证危险因素为目的的logistic回归分析分析思路主要是明确哪些因素可能是混杂因素并加以校正以发现主要研究因素与结局的真实关系。建议尽量避免的两种思路①把所有变量都进行校正。除非样本量足够大否则这种方式不可取。因为纳入的自变量越多所消耗的自由度越大用于估计主要研究因素的样本量相对越小结果的精确度也越低。②采用逐步回归筛选变量。作为主要研究变量一定要保留在模型中同时要纳入混杂因素。逐步回归筛选适用于探索危险因素不适用于验证危险因素。队列研究中的logistic回归队列研究绝大多数都是为了验证某一危险因素这是由研究性质决定的。因为队列研究在一开始就需要指定暴露和非暴露也就相当于确定了主要研究因素。因此从数据分析角度来讲队列研究主要是为了排除混杂因素与前文介绍的思路并无不同。但队列研究在时间顺序上可以证明研究因素发生在前结局发生在后因此其验证能力更强。由于队列研究有明确的时间先后顺序此时在说明主要研究因素与结局的关联强度时可采用RRrisk ratio而非ORodds ratio。队列研究中当结局发生率很低时10%OR是RR的一个很好的替代指标此时用logistic回归可直接求得OR值用来说明暴露的危险度。但如果结局发生率不是很低OR与RR差别较大此时用OR来说明危险度可能会有一定的偏倚 [3]。例5某研究分析Hp感染与胃黏膜病变进展的关系观察数据见表 5。本研究如果计算OR值则OR2.4495%CI:1.05~5.70如果计算RR值则RR1.7795%CI:1.01~3.12。由于病变进展的发生率较高两个指标差别较大。队列研究中RR值的计算通常可采用对数二项分布回归log-binomial regression。通常需要借助软件实现如SAS的proc genmod过程[4]。小结本文介绍了病例对照研究和队列研究中logistic回归分析的不同思路以及常见的一些应用错误。然而本文的思路并不仅限于logistic回归分析完全可以推广到其他广义线性模型。例如队列研究的观察结局如果是计数资料则可考虑Poisson回归或负二项回归此时仍需考虑混杂因素的校正问题。因此本文思路对各种常见的回归模型均有一定借鉴意义至于模型的选择主要取决于研究结局类型及其分布。在各种常见的回归分析中一定要分清研究类型及其目的到底是探索危险因素还是验证危险因素。危险因素的筛选过程较为复杂需要考虑较多问题包括变量筛选方式等验证危险因素相对较为简单不需要考虑变量筛选但要明确混杂因素并加以校正。一定要避免“把数据完全交给软件”这种分析方式软件主要用来解决计算问题分析思路必须由研究者来确定。统计分析不是简单的参数估计而应结合研究类型明确研究思路才能给出合理的结果。利益冲突所有作者均声明不存在利益冲突参考文献1.Mennard S. Applied logistic regression analysis[M]. Newbury Park, California: SAGE Publications, Inc, 2001.2.徐飙. 流行病学原理[M]. 上海: 复旦大学出版社, 2007.Xu B. Epidemic theory[M]. Shanghai: Fudan University Press, 2007.3.Stokes ME, Davis CS, Koch GG. Categorical data analysis using the SASsystem[M]. 2nd ed. Cary, NC: John Willy Sons, Inc, 2000.4.冯国双, 刘德平. 医学研究中的logistic回归分析及SAS实现[M]. 2版. 北京: 北京大学医学出版社,, 2015.Feng GS, Liu DP. Logistic regression analysis and SAS application in medical research[M]. 2nd ed. Beijing: Peking University Medical Press, 2015.扫码关注“医咖会”公众号及时获取最新统计教程
