统计学回归分析

统计学回归分析第1页，共56页，2023年，2月20日，星期二问题提出：

医学研究中常研究某因素存在条件下某结果是否发生？以及之间的关系如何？因素（X）疾病结果（Y）

x1，x2，x3…XK

发生Y=1

不发生Y=0例：暴露因素冠心病结果高血压史(x1)：有或无有或无高血脂史(x2)：有或无吸烟(x3)：有或无第2页，共56页，2023年，2月20日，星期二研究问题可否用多元线性回归方法？

1.多元线性回归方法要求Y的取值为计量的连续性随机变量。2.多元线性回归方程要求Y与X间关系为线性关系。3.多元线性回归结果不能回答“发生与否”logistic回归方法补充多元线性回归的不足第3页，共56页，2023年，2月20日，星期二Logistic回归方法

该法研究是当y取某值（如y=1）发生的概率（p）与某暴露因素（x）的关系。

P（概率）的取值波动0～1范围。基本原理：用一组观察数据拟合Logistic模型，揭示若干个x与一个因变量取值的关系，反映y对x的依存关系。第4页，共56页，2023年，2月20日，星期二

一、基本概念

1.变量的取值

logistic回归要求应变量（Y）取值为分类变量（两分类或多个分类）

自变量（Xi）称为危险因素或暴露因素，可为连续变量、等级变量、分类变量。可有m个自变量X1，X2，…Xm

第5页，共56页，2023年，2月20日，星期二2.两值因变量的logistic回归模型方程一个自变量与Y关系的回归模型如：y：发生=1,未发生=0x有=1无=0，记为p（y=1/x）表示某暴露因素状态下，结果y=1的概率（P）模型。或模型描述了应变量p与x的关系第6页，共56页，2023年，2月20日，星期二P概率10.5Z值0123-1-2-3图16-1Logistic回归函数的几何图形Β为正值，x越大，结果y=1发生的可能性（p）越大。第7页，共56页，2023年，2月20日，星期二几个logistic回归模型方程第8页，共56页，2023年，2月20日，星期二logistic回归模型方程的线性表达对logistic回归模型的概率（p）做logit变换，截距（常数）回归系数Y～（-∞至+∞）线性关系方程如下：第9页，共56页，2023年，2月20日，星期二在有多个危险因素（Xi）时多个变量的logistic回归模型方程的线性表达：或第10页，共56页，2023年，2月20日，星期二2.模型中参数的意义Β0（常数项）：暴露因素Xi=0时，个体发病概率与不发病概率之比的自然对数比值。第11页，共56页，2023年，2月20日，星期二

的含义：某危险因素，暴露水平变化时，即Xi=1与Xi=0相比，发生某结果（如发病）优势比的对数值。

P1（y=1/x=1）的概率P0（y=1/x=0）的概率第12页，共56页，2023年，2月20日，星期二

危险因素

Yx=1x=0发病=130（a）10（b）不发病=070（c）90（d）

a+cb+d

危险因素

Yx=1x=0发病=1p1p0

不发病=01-p11-p0

有暴露因素人群中发病的比例第13页，共56页，2023年，2月20日，星期二

反映了在其他变量固定后，X=1与x=0相比发生Y事件的对数优势比。回归系数β与ORX与Y的关联

β=0，OR=1，无关

β＞0，OR＞1，有关，危险因素

β＜0，OR＜1，有关，保护因子事件发生率很小，OR≈RR。多元回归模型的的概念第14页，共56页，2023年，2月20日，星期二二、Logistic回归模型Logistic回归的分类

二分类多分类条件Logistic回归非条件Logistic回归第15页，共56页，2023年，2月20日，星期二Logit变换

也称对数单位转换

logitP=第16页，共56页，2023年，2月20日，星期二流行病学概念：设P表示暴露因素X时个体发病的概率，则发病的概率P与未发病的概率1-P

之比为优势（odds），logitP就是odds的对数值。第17页，共56页，2023年，2月20日，星期二Logistic回归模型

Logistic回归的logit模型

Logistic回归模型第18页，共56页，2023年，2月20日，星期二三、参数估计最大似然估计法（Maximumlikehoodestimate）似然函数：L=∏Pi

对数似然函数：

lnL=∑(lnP)=lnP1+lnP2+…+lnPn

非线性迭代方法——

Newton-Raphson法第19页，共56页，2023年，2月20日，星期二四、参数检验似然比检验（likehoodratiotest）

通过比较包含与不包含某一个或几个待检验观察因素的两个模型的对数似然函数变化来进行，其统计量为G

（又称Deviance）。

G=-2(lnLp-lnLk)

样本量较大时，G近似服从自由度为待检验因素个数的２分布。第20页，共56页，2023年，2月20日，星期二比分检验（scoretest）以未包含某个或几个变量的模型为基础，保留模型中参数的估计值，并假设新增加的参数为零，计算似然函数的一价偏导数（又称有效比分）及信息距阵，两者相乘便得比分检验的统计量S

。样本量较大时，S近似服从自由度为待检验因素个数的２分布。第21页，共56页，2023年，2月20日，星期二Wald检验（waldtest）即广义的t检验，统计量为u

u服从正态分布，即为标准正态离差。

Logistic回归系数的区间估计

第22页，共56页，2023年，2月20日，星期二

上述三种方法中，似然比检验最可靠，比分检验一般与它相一致，但两者均要求较大的计算量；而Wald检验未考虑各因素间的综合作用，在因素间有共线性时结果不如其它两者可靠。第23页，共56页，2023年，2月20日，星期二五、回归系数的意义

单纯从数学上讲，与多元线性回归分析中回归系数的解释并无不同，亦即bi表示xi改变一个单位时，logitP的平均变化量。第24页，共56页，2023年，2月20日，星期二流行病学中的一些基本概念：相对危险度（relativerisk）:RR=P1/P2比数

Odds=P/(1-P)比数比OR=[P１/(1-P１)]/[P２/(1-P２)]在患病率较小情况下，OR≈RR第25页，共56页，2023年，2月20日，星期二Logistic回归中的常数项（b0）表示，在不接触任何潜在危险／保护因素条件下，效应指标发生与不发生事件的概率之比的对数值。Logistic回归中的回归系数（bi

）表示，某一因素改变一个单位时，效应指标发生与不发生事件的概率之比的对数变化值，即OR的对数值。第26页，共56页，2023年，2月20日，星期二Logistic回归系数的意义分析因素xi为二分类变量时，存在（暴露）xi

＝１，不存在（未暴露）xi

＝０，则Logistic回归中xi的系数bi就是暴露与非暴露优势比的对数值．即OR=exp(bi)=e(bi)第27页，共56页，2023年，2月20日，星期二分析因素xi为多分类变量时，为方便起见，常用1，2，…，k分别表示k个不同的类别。进行Logistic回归分析前需将该变量转换成k-1个指示变量或哑变量（design/dummyvariable），这样指示变量都是一个二分变量，每一个指示变量均有一个估计系数，即回归系数，其解释同前。第28页，共56页，2023年，2月20日，星期二分析因素xi为等级变量时，如果每个等级的作用相同，可按计量资料处理：如以最小或最大等级作参考组，并按等级顺序依次取为0，1，2，…。此时，e(bi)

表示xi增加一个等级时的优势比，e(k*bi)表示xi增加k个等级时的优势比。如果每个等级的作用不相同，则应按多分类资料处理。分析因素xi为连续性变量时，e(bi)表示xi增加一个计量单位时的优势比。第29页，共56页，2023年，2月20日，星期二

多因素Logistic回归分析时，对回归系数的解释都是指在其它所有自变量固定的情况下的优势比。存在因素间交互作用时，Logistic回归系数的解释变得更为复杂，应特别小心。第30页，共56页，2023年，2月20日，星期二

根据Wald检验，可知Logistic回归系数bi服从u分布。因此其可信区间为进而，优势比e(bi)的可信区间为第31页，共56页，2023年，2月20日，星期二六、Logistic回归分析方法基本思想同线性回归分析。从所用的方法看，有强迫法、前进法、后退法和逐步法。在这些方法中，筛选变量的过程与线性回归过程的完全一样。但其中所用的统计量不再是线性回归分析中的F统计量，而是以上介绍的参数检验方法中的三种统计量之一。第32页，共56页，2023年，2月20日，星期二

为计算方便，通常向前选取变量用似然比或比分检验，而向后剔除变量常用Wald检验。第33页，共56页，2023年，2月20日，星期二七、条件Logistic回归对配对/比调查资料，应该用条件Logistic回归分析。

对于配比资料，第i个配比组可以建立一个Logistic回归：第34页，共56页，2023年，2月20日，星期二假设自变量在各配比组中对结果变量的作用是相同的，即自变量的回归系数与配比组无关。配比设计的Logistic回归模型其中不含常数项。第35页，共56页，2023年，2月20日，星期二可以看出此回归模型与非条件Logistic回归模型十分相似，只不过这里的参数估计是根据条件概率得到的，因此称为条件Logistic回归模型。条件Logistic回归的回归系数检验与分析，和非条件Logistic回归完全相同。第36页，共56页，2023年，2月20日，星期二

八、logistic回归的应用

1.疾病（某结果）的危险因素分析和筛选

用回归模型中的回归系数（βi）和OR说明危险因素与疾病的关系。适用的资料：

前瞻性研究设计、病例对照研究设计、横断面研究设计的资料。三类研究计算的logistic回归模型的β意义是一致。仅常数项不同。（证明略）第37页，共56页，2023年，2月20日，星期二2.校正混杂因素，对疗效做评价在临床研究和疗效的评价，组间某些因素构成不一致干扰疗效分析，通过该法可控制非处理因素，正确评价疗效。3.预测与判别预测个体在某因素存在条件下，发生某事件（发病）的概率，为进一步治疗提供依据。第38页，共56页，2023年，2月20日，星期二输精管切除术与动脉粥样硬化疾病的研究1．问题的描述（1）输精管切除术是否与动脉粥样硬化疾病有关？（2）如果存在联系，与其他已知的危险因素相比，输精管切除术的相对重要性有多大？（3）哪些男性亚群在输精管切除术以后发生动脉粥样硬化疾病的可能性特别大？

九、logistic回归的应用举例第39页，共56页，2023年，2月20日，星期二２．研究目的（1）一般目的确定输精管切除术与其后因动脉粥样硬化疾病而住院之间是否有因果联系；如果有，则确定输精管切除术与患动脉粥样硬化性疾病的危险性是否因观察对象具有其他冠心病危险因素，如吸烟、高血压和高胆固醇等而增强。第40页，共56页，2023年，2月20日，星期二（2）特殊目的A．估计输精管切除术及其男性动脉粥样硬化性疾病危险因素的相对危险度（使用单变量分析方法）；B．估计输精管切除术对动脉粥样硬化性疾病独立的作用（利用条件Logistic回归模型）；C．检测输精管切除术导致动脉粥样硬化性疾病的可能期限；D．检测在输精管切除术、吸烟和高血压之间可能的协同作用。第41页，共56页，2023年，2月20日，星期二３．研究变量（1）按照世界卫生组织标准来诊断动脉粥样硬化性疾病；（2）病人特征：年龄，出生日期，宗教，教育，职业，家庭史，婚姻状况。（3）生育史：存活子女的数量与性别，妻子的生育情况。（4）生活方式：吸烟习惯，饮酒量，饮食习惯，食盐量，喝咖啡量，体力活动。（5）医学史：有无可能导致不育的疾病或手术、高血压、糖尿病或高胆固醇血症。第42页，共56页，2023年，2月20日，星期二４．研究设计

①抽样：以完全随机抽样方式确定样本。样本量大小的计算是基于事先确定显著性水平、把握度和打算测知的相对危险度水平。假定样本量确定为病例组500例，对照组1000例。本设计力图避免或减少在病例一对照研究中常见的偏倚与误差。回忆偏倚、选择偏倚和输精管切除术报告偏倚对研究结果的真实性是极其重要的。第43页，共56页，2023年，2月20日，星期二②资料收集：病例组的选取：调查员审阅每日住院病人情况→如果诊断适合研究的范围，将病例转给心脏病主任医师作评估→由他做出病例诊断是否合格的决定→调查人员核对病人背景资料是否合格→如果病人满足诊断标椎和背景资料合格→调查人员开始询问并填写调查表→每完成5个病例和10个配对对照以后，请研究中心的工作人员对调查表进行评估→重复以上步骤。第44页，共56页，2023年，2月20日，星期二对照组的选取：调查人员评估入院单并选择符合配对标椎和入院诊断合适的对照病例→核对病人背景资料是否符合入院标准→如对照组成员合格，则进行调查。第45页，共56页，2023年，2月20日，星期二③统计分析：A．计算未校正的相对危险比（OR）与95％可信限范围，以估计输精管切除术和其他危险因素的相对危险度。B．利用配对的条件Logistic回归模型计算校正相对危险比，用此比来评价输精管切除术对动脉粥样硬化的独立影响。C．利用动脉粥样硬化作为因变量，输精管切除术后间隔时间作为自变量，用条件Logistic回归模型来检验输精管切除术后时间长短对发生动脉粥样硬化的时间效应。第46页，共56页，2023年，2月20日，星期二二、Logistic回归应用的注意事项

1.模型中自变量的取值自变量（X）可为计量数据、分类数据和等级数据。计量数据常重新划为有序组段，OR的实际意义较大。例１：年龄（岁，x1）

十、logistic回归的注意事项第47页，共56页，2023年，2月20日，星期二数据的几种赋值形式1）两分类变量，赋值为：有=1，无=02）有序变量，赋值；无=0，少=1，中=2，多=3

例；年龄＜45=145-54=255-64=3≥65=43.）多分类无序变量：赋值为：哑变量（dummyvariable）形式见例：注：变量取值不同，方程的系数和符号将发生变化。第48页，共56页，2023年，2月20日，星期二

冠心病可能的危险因素与赋值因素变量名赋值说明年龄(岁)X1<45=1,4554=2,5564=3,65=4高血压史X2无=0,有=1高血压家族史X3无=0,有=1吸烟X4不吸=0,吸=1第49页，共56页，2023年，2月20日，星期二年龄（X）化为哑变量的赋值年龄（岁)有序变量哑变量（方法一）X水平D1D2D3<40100040～

210050～301060～4001方程1：有序变量方程Β含义：x每增加1个单位（10岁），发病的lnOR平均增加β1第50页，共56页，2023年，2月20日，星期二方程2：哑变量方程（哑变量个数=分类数－1）方程系数的解释：

表示40-岁/＜40岁相比的对数优势比表示50-岁/＜40岁相比的对数优势比表示60-岁/＜40岁相比的对数优势比第51页，共56页，2023年，2月20