(精品)同济医学院-《SAS》课件-SAS的logistic回归(正式).ppt

华中科技大学公卫学院流行病与卫生统计系Logistic回归回归SAS应用应用蒋红卫蒋红卫Email:JHWCCC21CN.COM学习目标学习目标了解了解LogisticLogistic回归模型的建立和假设检验；回归模型的建立和假设检验；了解了解LogisticLogistic回归模型的应用领域；回归模型的应用领域；掌握掌握LogisticLogistic回归模型系数的解释，及回归系数回归模型系数的解释，及回归系数与与OROR值之间的关系；值之间的关系；掌握掌握LogisticLogistic回归过程步；回归过程步；掌握哑变量的设置和结果的解释；掌握哑变量的设置和结果的解释；掌握多元掌握多元LogisticLogistic回归模型的逐步过程法和系数回归模型的逐步过程法和系数的解释；的解释；了解条件了解条件LogisticLogistic回归的应用；回归的应用；掌握条件掌握条件LogisticLogistic回归的回归的SASSAS程序；程序；概述概述线性回归模型和广义线性回归模型要求因变量是线性回归模型和广义线性回归模型要求因变量是连续的正态分布变量，且自变量和因变量呈线性连续的正态分布变量，且自变量和因变量呈线性关系。当因变量是分类型变量时，且自变量与因关系。当因变量是分类型变量时，且自变量与因变量没有线性关系时，线性回归模型的假设条件变量没有线性关系时，线性回归模型的假设条件遭到破坏。这时，最好的回归模型是遭到破坏。这时，最好的回归模型是LogisticLogistic回回归模型，它对因变量的分布没有要求，从数学角归模型，它对因变量的分布没有要求，从数学角度看，度看，LogisticLogistic回归模型非常巧妙地避开了分类回归模型非常巧妙地避开了分类型变量的分布问题，补充完善了线性回归模型和型变量的分布问题，补充完善了线性回归模型和广义线性回归模型的缺陷。从医学研究角度看，广义线性回归模型的缺陷。从医学研究角度看，LogisticLogistic回归模型解决了一大批实际应用问题，回归模型解决了一大批实际应用问题，对医学的发展起着举足轻重的作用。对医学的发展起着举足轻重的作用。非条件非条件LogisticLogistic回归回归 LogisticLogistic回归分析在医学回归分析在医学研究研究中应用广泛。中应用广泛。目前目前主要是用于流行病学研究中危险因素主要是用于流行病学研究中危险因素的筛选，但它同时具有良好的判别和预测的筛选，但它同时具有良好的判别和预测功能，尤其是在资料类型不能满足功能，尤其是在资料类型不能满足FisherFisher判别和判别和BayesBayes判别的条件时，更显示出判别的条件时，更显示出LogisticLogistic回归判别的优势和效能。本研究对回归判别的优势和效能。本研究对LogisticLogistic回归方程的判别分析进行了回归方程的判别分析进行了探讨。探讨。非条件非条件LogisticLogistic回归回归医学研究中经常需要分析分类型变量的问题。比如，生存医学研究中经常需要分析分类型变量的问题。比如，生存与死亡、有病与无病、有效与无效、感染与未感染等二分与死亡、有病与无病、有效与无效、感染与未感染等二分类变量。研究者关心的问题是，哪些因素导致了人群中有类变量。研究者关心的问题是，哪些因素导致了人群中有些人患某种病而有些人不患某种病，哪些因素导致了某种些人患某种病而有些人不患某种病，哪些因素导致了某种治疗方法出现治愈、显效、好转和无效等不同的效果等。治疗方法出现治愈、显效、好转和无效等不同的效果等。这类问题，实质上是一个回归问题，因变量就是上述提到这类问题，实质上是一个回归问题，因变量就是上述提到的这些分类型变量，自变量的这些分类型变量，自变量x x是与之有关的一些因素。但是与之有关的一些因素。但是，这样的问题却不能直接用线性回归分析方法解决，其是，这样的问题却不能直接用线性回归分析方法解决，其根本原因在于因变量是分类型变量，严重违背了线性回归根本原因在于因变量是分类型变量，严重违背了线性回归分析对数据的假设条件。那么应该怎样解决这个问题呢？分析对数据的假设条件。那么应该怎样解决这个问题呢？非条件非条件LogisticLogistic回归回归研究者将所研究的问题转换一个角度，不研究者将所研究的问题转换一个角度，不是直接分析是直接分析y y与与x x的关系，而是分析的关系，而是分析y y取某个取某个值的概率值的概率P P与与x x的关系。例如，令的关系。例如，令y y为为1 1，0 0变变量，量，y=1y=1表示有病，表示有病，y=0y=0表示未患病；表示未患病；x x是是与患病有关的危险因素。如果与患病有关的危险因素。如果P P表示患病的表示患病的概率，即概率，即P P=prob=prob（y=1y=1），那么研究患病），那么研究患病的概率的概率P P与危险因素与危险因素x x的关系就不是很困难的关系就不是很困难的事情了。的事情了。非条件非条件LogisticLogistic回归回归分析因变量分析因变量y y取某个值的概率取某个值的概率P P与自变量与自变量x x的关系，就是寻的关系，就是寻找一个连续函数，使得当找一个连续函数，使得当x x变化时，它对应的函数值变化时，它对应的函数值P P不超不超出出00，11范围。数学上这样的函数是存在且不唯一的，范围。数学上这样的函数是存在且不唯一的，LogisticLogistic回归模型就是满足这种要求的函数之一。与线性回归模型就是满足这种要求的函数之一。与线性回归分析相似，回归分析相似，LogisticLogistic回归分析的基本原理就是利用一回归分析的基本原理就是利用一组数据拟合一个组数据拟合一个LogisticLogistic回归模型，然后借助这个模型揭回归模型，然后借助这个模型揭示总体中若干个自变量与一个因变量取某个值的概率之间示总体中若干个自变量与一个因变量取某个值的概率之间的关系。具体地说，的关系。具体地说，LogisticLogistic回归分析可以从统计意义上回归分析可以从统计意义上估计出在其它自变量固定不变的情况下，每个自变量对因估计出在其它自变量固定不变的情况下，每个自变量对因变量取某个值的概率的数值影响大小。变量取某个值的概率的数值影响大小。LogisticLogistic回归模型有条件与非条件之分，前者适用于配对回归模型有条件与非条件之分，前者适用于配对病例对照资料的分析，后者适用于队列研究或非配对的病病例对照资料的分析，后者适用于队列研究或非配对的病例例-对照研究成组资料的分析。对照研究成组资料的分析。问题的提出问题的提出 在流行病学研究中，经常遇到因变量为离散型分在流行病学研究中，经常遇到因变量为离散型分类变量的情况。如治疗效果的无效好转、显效、类变量的情况。如治疗效果的无效好转、显效、痊愈；不同染毒剂量下小白鼠的存活或死亡；在痊愈；不同染毒剂量下小白鼠的存活或死亡；在某种暴露下的发病与不发病等。最常见的情况是某种暴露下的发病与不发病等。最常见的情况是因变量为二分变量的问题。因变量为二分变量的问题。多元线性回归的局限性多元线性回归的局限性 经典流行病学统计分析方法经典流行病学统计分析方法分层分析的局限性分层分析的局限性1.1.两种主要的流行病学设计两种主要的流行病学设计1 1）病历对照研究）病历对照研究2 2）队列研究）队列研究2.2.判断结局（疾病）和暴露（因素）联系强弱的指标判断结局（疾病）和暴露（因素）联系强弱的指标1 1）相对危险度：相对危险度：RR=p1/p0 RR=p1/p0 p1:p1:暴露于某个危险因素下发病的概率暴露于某个危险因素下发病的概率 p0:p0:不暴露于某个危险因素下发病的概率（对照）不暴露于某个危险因素下发病的概率（对照）2 2）比值比：）比值比：OR=P(D=1|E=1)/P(D=0|E=1)/OR=P(D=1|E=1)/P(D=0|E=1)/P(D=1|E=0)/P(D=0|E=0)P(D=1|E=0)/P(D=0|E=0)D=1:D=1:患某种疾病，患某种疾病，D=0D=0：不患某种疾病：不患某种疾病 E=1:E=1:暴露于某个危险因素，暴露于某个危险因素，E=0:E=0:不暴露于某个危险因素不暴露于某个危险因素 可以简单地表述成：可以简单地表述成：OR=(p1/q1)/(p0/q0)OR=(p1/q1)/(p0/q0)p1:p1:暴露于某个危险因素下发病的概率暴露于某个危险因素下发病的概率 q1:q1:暴露于某个危险因素下不发病的概率暴露于某个危险因素下不发病的概率 p0:p0:不暴露于某个危险因素下发病的概率不暴露于某个危险因素下发病的概率 q0:q0:不暴露于某个危险因素下不发病的概率不暴露于某个危险因素下不发病的概率能够描述病因学意义的指标是能够描述病因学意义的指标是RRRROROR本身没有病因学意义本身没有病因学意义OROR的病因学意义主要体现在发病率的病因学意义主要体现在发病率/患病率比较患病率比较低的情形，此时，低的情形，此时，OROR近似于近似于RRRROROR值是描述暴露与疾病关联强度的指标值是描述暴露与疾病关联强度的指标l分分层层分分析析中中，可可以以分分别别计计算算出出分分层层后后的的各各层层OROR值值，如如果果发发现现与与总总的的OROR有有较较大大的的差差异异，则则可可以以认认为为该该风风层层因因素素是是混混杂杂因因素素。必必须须对对该该因因素素进进行行MHMH调调整整，调调整整后后的的OROR值值才才能能真真正正反反映映因因素素和和结结局局间的关系。间的关系。l如如果果当当分分层层后后各各层层的的OROR值值经经过过一一致致性性检检验验发发现现：各各层层间间的的OROR值值有有统统计计学学差差异异，这这时时说说明明分分析析因因素素在在分分层层因因素素的的不不同同水水平平上上与与结结局局变变量量的的联联系系强强度度是是不不同同的的，这这时时分分层层因因素素和和研研究究因因素素存存在在这这交交互互作作用用（效效应应修修饰饰作作用用）。这这时时应应该该分分层报告层报告OROR值，而不能计算调整值，而不能计算调整OROR值。值。l 分层分析的局限性分层分析的局限性l只只能能控控制制少少数数因因素素（分分层层因因素素过过多多，每个格子中的样本例数太少）每个格子中的样本例数太少）l定量资料需要分组，信息丢失定量资料需要分组，信息丢失l不不能能对对因因素素作作用用大大小小进进行行定定量量分分析析（交互作用）（交互作用）yy=logy=log2 2x x01二、二、Logistic Logistic 回归原理回归原理 Logistic 模型中系数的意义：模型中系数的意义：回归系数的流行病学意义是：在其它自变量都回归系数的流行病学意义是：在其它自变量都不变的条件下，当因素不变的条件下，当因素X X变化一个测量单位时所引起的变化一个测量单位时所引起的OROR值自然对数的改变量。值自然对数的改变量。三、三、Logistic Logistic 回归和回归和OROR值间的关系值间的关系假设建立了如下的假设建立了如下的logisticlogistic回归方程：回归方程：Logit P=Logit P=+x+xx x 为二分变量，当暴露时，取值为为二分变量，当暴露时，取值为1 1；不暴露时，取值为不暴露时，取值为0 0。所以所以暴露暴露时时,Logit(PLogit(P1 1)=)=+，比值比值(odds)=exp(odds)=exp(+)所以所以不暴露不暴露时时,Logit(PLogit(P0 0)=)=，比值比值(odds)=exp(odds)=exp()则，暴露对于不暴露的比值比则，暴露对于不暴露的比值比(odds ratio)(odds ratio)为：为：OR=OR=exp(exp(+)/exp()/exp()=exp()=exp()以四格表为例来说明最大似然求解的意义及过程。以四格表为例来说明最大似然求解的意义及过程。四格表的一般表达形式四格表的一般表达形式发病发病(y=1)(y=1)不发病不发病(y=0)(y=0)暴露暴露(x=1)(x=1)a ab b 不暴露不暴露(x=0)(x=0)c cd d 合计合计a+ca+cb+db+d 暴露者发病概率暴露者发病概率 p p1 1=a/(a+b);=a/(a+b);不暴露者发病概率不暴露者发病概率 p p0 0=c/(c+d)=c/(c+d)OR=ad/(bc)OR=ad/(bc)四、四、Logistic Logistic 回归最大似然建模回归最大似然建模用发病概率来表示四格表，可以得到四格表的另外一种表示形式：用发病概率来表示四格表，可以得到四格表的另外一种表示形式：四格表的另外一种表达形式四格表的另外一种表达形式(1)(1)发病发病(y=1)(y=1)不发病不发病(y=0)(y=0)暴露暴露(x=1)(x=1)p p1 1 1-1-p p1 1 不暴露不暴露(x=0)(x=0)p p0 0 1-1-p p0 0 暴露者发病概率暴露者发病概率:p:p1 1=exp(=exp(+x)/1+x)/1+exp(exp(+x)+x)暴露者不发病概率暴露者不发病概率:q:q0 0=1-p=1-p1 1=1/=1/1+1+exp(exp(+x)+x);不暴露者发病概率不暴露者发病概率:p:p0 0=exp(=exp()/1+)/1+exp(exp()不暴露者不发病概率不暴露者不发病概率:q:q0 0=1-p=1-p0 0=1/=1/1+1+exp(exp();用发病概率来表示四格表，可以得到四格表的另外一种表示形式：用发病概率来表示四格表，可以得到四格表的另外一种表示形式：四格表的另外一种表达形式四格表的另外一种表达形式(2)(2)发病发病(y=1)(y=1)不发病不发病(y=0)(y=0)暴露暴露(x=1)e(x=1)e(+)+)/1+/1+e e(+)+)1/1/1+1+e e(+)+)不暴露不暴露(x=0)(x=0)e e /1+/1+e e 1/1/1+1+e e 因为四格表的四个实际数为因为四格表的四个实际数为a,b,ca,b,c及及d d，故可构造似然函数为：故可构造似然函数为：L=L=e e(+)+)/1+/1+e e(+)+)a a 1/1/1+1+e e(+)+)b be e /1+/1+e e c c 1/1/1+1+e e d d 取对数，有取对数，有Ln(L)=a(Ln(L)=a(+)+)a ln1+e a ln1+e(+)(+)b ln1+e b ln1+e(+)(+)+c +c c ln 1+e c ln 1+e d ln1+e d ln1+e 对以上似然函数分别求对对以上似然函数分别求对 和和 的一阶偏导数，再令两个偏导数为的一阶偏导数，再令两个偏导数为零，就可以解得零，就可以解得 和和 的估计值。的估计值。v在经典流行病学分析里对因素的统计学意义检验在经典流行病学分析里对因素的统计学意义检验是检验是检验OROR值是否为值是否为1 1v对于对于LogisticLogistic回归来说，对模型中变量的统计学检回归来说，对模型中变量的统计学检验就是对检验其回归系数是否为验就是对检验其回归系数是否为0 0v一般常用的方法有三种：似然比检验、一般常用的方法有三种：似然比检验、WaldWald检验、检验、比分检验比分检验五、五、Logistic Logistic 回归模型的统计学回归模型的统计学检验检验1 1）似然比检验）似然比检验似然比检验的原理是通过分析模型中变量似然比检验的原理是通过分析模型中变量变化对似然比的影响，来检验增加变化对似然比的影响，来检验增加/减少的自减少的自变量是否对应变量有统计学上的显著意义。变量是否对应变量有统计学上的显著意义。检验统计量：检验统计量：G=-2 Ln(LkG=-2 Ln(Lk1)-Ln(Lk)1)-Ln(Lk)Ln(Lk Ln(Lk1)1)为不包含检验变量时模型的对数为不包含检验变量时模型的对数似然值似然值 Ln(Lk)Ln(Lk)为包含检验变量时模型的对数似然值为包含检验变量时模型的对数似然值v当检验一个变量时则当检验一个变量时则G G服从自由度为服从自由度为1 1的的 2 2 分分布布v当对模型整个进行检验时就是相当于模型中所当对模型整个进行检验时就是相当于模型中所有有n n个变量的回归系数都为个变量的回归系数都为0 0，G G服从自由度服从自由度为为n n的的 2 2 分布分布2 2）Wald Wald 统计量检验回归系数是否为统计量检验回归系数是否为0 0对于对于方程内单个自变量的检验，采用方程内单个自变量的检验，采用Wald Wald 统计量。统计量。W=/SW=/S S S 为为 的标准误。的标准误。其临界值可取标准正态分布的临界值。同时，其临界值可取标准正态分布的临界值。同时，WW2 2遵循遵循自由度自由度为为1 1的卡的卡方分布。方分布。的的9595可信区间为：可信区间为：（1.961.96 S S ，1.961.96 S S ）OROR的的9595可信区间为：可信区间为：(exp(exp(1.961.96 S S ),exp(),exp(+1.96+1.96 S S)3 3）比分检验）比分检验 在原有模型的基础上增加一个变量，并假设该在原有模型的基础上增加一个变量，并假设该变量的回归系数为变量的回归系数为0 0，似然函数一阶偏导和信息矩阵，似然函数一阶偏导和信息矩阵的乘积为比分检验统计量的乘积为比分检验统计量S S。S S在样本量较大时服从自在样本量较大时服从自由度为检验参数个数的由度为检验参数个数的 2 2 分布分布l上述方法中，似然比检验最可靠，上述方法中，似然比检验最可靠，WaldWald检验和比检验和比分检验一致。分检验一致。lWaldWald检验没有考虑变量间的相互影响，当因素间检验没有考虑变量间的相互影响，当因素间相关时（存在共线性）不可靠相关时（存在共线性）不可靠l参数的可信区间是以参数的可信区间是以WaldWald法计算的，在应用可信法计算的，在应用可信区间时要注意区间时要注意 PROC LOGISTICPROC LOGISTIC过程过程 SASSAS系统中进行系统中进行LogisticLogistic回归分析的过程步是回归分析的过程步是PROC LOGISTICPROC LOGISTIC过程。下面我们先了解以下这过程。下面我们先了解以下这个过程步。个过程步。LOGISTICLOGISTIC过程的语句基本格式如下：过程的语句基本格式如下：PROC LOGISTIC data=descending PROC LOGISTIC data=descending simple order=;simple order=;Model Model 因变量自变量因变量自变量/link=noint rsq cl/link=noint rsq cl rl selection=sle=sls=;rl selection=sle=sls=;freq freq 频数变量频数变量;output out=pred=output out=pred=变量名；变量名；run;run;PROC LOGISTICPROC LOGISTIC过程过程 DESCENDINGDESCENDING：指令系统输出因变量取值由大到小的概：指令系统输出因变量取值由大到小的概率，缺省时，系统输出因变量取值由小到大的概率。使用率，缺省时，系统输出因变量取值由小到大的概率。使用和不使用和不使用DESCENDINGDESCENDING语句的回归模型的回归系数绝对语句的回归模型的回归系数绝对值完全相同，但是符号完全相反。值完全相同，但是符号完全相反。SIMPLESIMPLE：指令系统输出每一个自变量的基本统计量。：指令系统输出每一个自变量的基本统计量。MODELMODEL：给出模型的因变量和自变量。：给出模型的因变量和自变量。LINKLINK：指定因变量概率函数的形式：指定因变量概率函数的形式NOINTNOINT：指令模型不含常数项。：指令模型不含常数项。RSQRSQ：指定输出确定系数。：指定输出确定系数。PROC LOGISTICPROC LOGISTIC过程过程CLCL：指定输出参数估计值的可信区间。：指定输出参数估计值的可信区间。RLRL：指定输出相对危险度（或比数比）的可信区间。：指定输出相对危险度（或比数比）的可信区间。SELECTIONSELECTION、SLESLE和和SLSSLS：和：和REGREG过程步的对应选项一过程步的对应选项一样，用于逐步回归分析。样，用于逐步回归分析。FREQFREQ：是否使用，取决于数据的输入方法。当：是否使用，取决于数据的输入方法。当x x和和y y变量变量的观察值被直接输入到的观察值被直接输入到SASSAS数据集时，不需要使用数据集时，不需要使用FREQFREQ语句。但是，当语句。但是，当x x是分类型变量，且是分类型变量，且x x和和y y变量的观察值不变量的观察值不是直接输入到是直接输入到SASSAS数据集时，而是将数据集时，而是将x x和和y y变量的频数输入变量的频数输入到到SASSAS数据集中时，需要使用数据集中时，需要使用FREQFREQ语句。语句。OUTPUTOUTPUT语句：用来保存指定的分析结果。语句：用来保存指定的分析结果。OutOut选择项用选择项用来给出输出数据集的名字。来给出输出数据集的名字。PredPred选择项用来给出一个预测选择项用来给出一个预测概率的变量名。概率的变量名。例题条件二分类例题条件二分类某北方城市研究喉癌发病的危险因素，用某北方城市研究喉癌发病的危险因素，用1 12 2配对的病例配对的病例 对照研究方法进行了调查。对照研究方法进行了调查。现选取了现选取了6 6个可能的危险因素，试做条件个可能的危险因素，试做条件logisticlogistic逐步回归分析逐步回归分析 。27因素因素变变量名量名赋值说赋值说明明咽炎咽炎X1无无=1，偶，偶尔尔=2，经经常常=3吸烟量吸烟量(支支/日日)X20=1，1=2，5=3，10=4，20=5声嘶史声嘶史X3无无=1，偶，偶尔尔=2，经经常常=3摄摄食新食新鲜鲜蔬菜蔬菜X4少少=1，经经常常=2，每天，每天=3摄摄食水果食水果X5很少很少=1，少量，少量=2，经经常常=3癌症家族史癌症家族史X6无无=0，有，有=1是否患喉癌是否患喉癌对对照照=0，病例，病例=1例题例题2 2多个自变量的二分类多个自变量的二分类为了探讨冠心病发生的有关危险因素，对为了探讨冠心病发生的有关危险因素，对2626例冠心病病人和例冠心病病人和2828例对照者进行病例例对照者进行病例-对对照研究，试用照研究，试用logisticlogistic逐步回归分析方法筛逐步回归分析方法筛选危险因素选危险因素 。28因素变量名赋值说明年龄(岁)X145=1，45=2，55=3，65=4高血压史X2无=0，有=1高血压家族史X3无=0，有=1吸烟X4不吸=0，吸=1高血脂史X5无=0，有=1动物脂肪摄入X6低=0，高=1体重指数(BMI)X724=1，24=2，26=3A型性格X8否=0，是=1冠心病Y对照=0，病例=1例题例题3 3有序分类有序分类某研究人员随机选择某研究人员随机选择8484例患某病的病人做例患某病的病人做临床试验，以探讨性别和治疗方法对该病临床试验，以探讨性别和治疗方法对该病疗效的影响。变量赋值为：性别（疗效的影响。变量赋值为：性别（：男：男=0=0，女，女=1=1）、治疗方法（）、治疗方法（：传统疗法：传统疗法=0=0，新，新型疗法型疗法=1=1）、疗效（）、疗效（：无效：无效=1=1，有效，有效=2=2，痊愈，痊愈=3=3）。请拟合性别、治疗方法对疗）。请拟合性别、治疗方法对疗效的有序效的有序logisticlogistic回归模型。回归模型。29例题例题4 4无序多分类无序多分类某研究人员欲了解不同社区和性别之间成某研究人员欲了解不同社区和性别之间成年居民获取健康知识途径是否不同，对年居民获取健康知识途径是否不同，对2 2个个社区的社区的314314名成人进行了调查，结果见表名成人进行了调查，结果见表16-1116-11。变量赋值为：社区（。变量赋值为：社区（：社区：社区1=01=0，社区，社区2=12=1）、性别（）、性别（：男：男=0=0，女，女=1=1）、）、获取健康知识途径（获取健康知识途径（：传统大众媒介：传统大众媒介=1=1，网络网络=2=2，社区宣传，社区宣传=3=3）。请拟合社区和性）。请拟合社区和性别对居民获取健康知识途径的多分类别对居民获取健康知识途径的多分类logisticlogistic回归模型。回归模型。30条件条件LogisticLogistic回归回归 医学研究中的配对病例对照研究资料可以用条件医学研究中的配对病例对照研究资料可以用条件LogisticLogistic回归模型来分析。条件回归模型来分析。条件LogisticLogistic回归模型回归模型和非条件和非条件LogisticLogistic回归模型的区别在于参数的估回归模型的区别在于参数的估计是否用到了条件概率。所谓的配对病例对照研计是否用到了条件概率。所谓的配对病例对照研究指的是在病例对照研究中，对每一个病例配以究指的是在病例对照研究中，对每一个病例配以性别、年龄或其它条件相似的一个（性别、年龄或其它条件相似的一个（1 1：1 1）或几）或几个（个（1 1：MM）对照，然后分析比较病例组与对照组）对照，然后分析比较病例组与对照组以往暴露于致病因素的经历。分析配对病例对照以往暴露于致病因素的经历。分析配对病例对照研究资料所用的条件研究资料所用的条件LogisticLogistic回归模型中参数的回归模型中参数的估计方法也是采用最大似然估计法，参数和模型估计方法也是采用最大似然估计法，参数和模型的检验方法和非条件的检验方法和非条件LogisticLogistic回归模型也一样。回归模型也一样。条件条件LogisticLogistic回归的回归的SASSAS程序程序 用于条件用于条件LogisticLogistic回归分析的过程步一般借用生回归分析的过程步一般借用生存数据风险回归分析的存数据风险回归分析的PHREGPHREG过程步。实用过程步。实用PHREGPHREG过程步时要注意数据的输入方法和过程步过程步时要注意数据的输入方法和过程步语句的写法。语句的写法。PHREGPHREG过程步的一般格式如下：过程步的一般格式如下：Proc phreg data=;Proc phreg data=;Model Model 因变量因变量*截尾变量截尾变量=协变量协变量/ties=risklimits selection=sle=sls details;ties=risklimits selection=sle=sls details;Strata Strata 分组变量分组变量;Run;Run;条件条件LogisticLogistic回归的回归的SASSAS程序程序MODELMODEL语句：用于指定模型的结构，适用于生语句：用于指定模型的结构，适用于生存时间有右截尾的情况，且生存时间变量作为模存时间有右截尾的情况，且生存时间变量作为模型的因变量，协变量作为模型的自变量。型的因变量，协变量作为模型的自变量。TIESTIES：用来选择处理生存时间结点的方法，每一：用来选择处理生存时间结点的方法，每一种方法使用了不同的公式来计算最大似然值。种方法使用了不同的公式来计算最大似然值。RISKLIMITSRISKLIMITS：指令输出危险比的：指令输出危险比的9595置信区间，置信区间，供选择的有：供选择的有：backwardbackward、forwardforward、stepwisestepwise、scorescore。条件条件LogisticLogistic回归的回归的SASSAS程序程序SLESLE：指定协变量进入模型的显著水平，缺省值：指定协变量进入模型的显著水平，缺省值是是0.050.05。SLSSLS：指定协变量停留在模型中的显著水平，缺：指定协变量停留在模型中的显著水平，缺省值是省值是0.050.05。DETAILSDETAILS：指令输出逐步回归过程中每一步的详：指令输出逐步回归过程中每一步的详细分析结果。细分析结果。STRATASTRATA：指定用于分组计算的分组变量。如果：指定用于分组计算的分组变量。如果分组变量的数值不符合分组要求，则在变量后面分组变量的数值不符合分组要求，则在变量后面的括号内列出分组的端点值。的括号内列出分组的端点值。例题例题1 1二分类二分类研究吸烟、饮酒与食管癌关系的病例研究吸烟、饮酒与食管癌关系的病例-对照对照资料，试作资料，试作logisticlogistic回归分析。回归分析。35本章小节本章小节 LogisticLogistic回归分析在医学回归分析在医学研究研究中应用广泛。中应用广泛。Logistic Logistic 回归分析在流行病学的病因研究回归分析在流行病学的病因研究中，是分析疾病与危险因素间联系的一种中，是分析疾病与危险因素间联系的一种统计方法。在这类研究中，所观察的项目统计方法。在这类研究中，所观察的项目的值，常以二项反应变量取值，即生存与的值，常以二项反应变量取值，即生存与死亡，是否发病，是否接触危险因素等的死亡，是否发病，是否接触危险因素等的反应变量反应变量y y的取值是的取值是0 0或或1 1。本章小节本章小节LogisticLogistic回归模型有条件与非条件之分，前回归模型有条件与非条件之分，前者适用于配对病例对照资料的分析，后者者适用于配对病例对照资料的分析，后者适用于队列研究或非配对的病例适用于队列研究或非配对的病例-对照研究对照研究成组资料的分析。成组资料的分析。在使用分析流行病学的在使用分析流行病学的方法研究疾病病因时，非条件方法研究疾病病因时，非条件LogisticLogistic模型模型是用于分析队列或病例是用于分析队列或病例-对照研究成组资料对照研究成组资料的统计方法，既可以进行因素筛选，也可的统计方法，既可以进行因素筛选，也可以用于混杂因素的控制。以用于混杂因素的控制。SASSAS系统中采用系统中采用LOGISTICLOGISTIC过程步分析非条件过程步分析非条件LogisticLogistic回归。回归。本章小节本章小节在医学研究中经常遇到多分类无序变量，此类变量一般不在医学研究中经常遇到多分类无序变量，此类变量一般不能直接进入能直接进入LogisticLogistic回归模型进行分析，这种变量一般采回归模型进行分析，这种变量一般采用设置哑变量的方法来进行分析。哑变量，就是一组取值用设置哑变量的方法来进行分析。哑变量，就是一组取值1 1和和0 0的二值分类变量，用来表示一个分类变量。的二值分类变量，用来表示一个分类变量。非条件非条件LogisticLogistic回归适合于队列研究、病例对照研究，同回归适合于队列研究、病例对照研究，同样也适合于现况调查中的病因学研究。按样也适合于现况调查中的病因学研究。按3 3种不同抽样方种不同抽样方式作式作LogisticLogistic回归，除病例对照研究资料的常数项与另外回归，除病例对照研究资料的常数项与另外两种不同外，回归系数的意义相同。而配对病例对照研究两种不同外，回归系数的意义相同。而配对病例对照研究资料就需要条件资料就需要条件LogisticLogistic回归模型来分析了。回归模型来分析了。SASSAS系统中系统中采用采用PHREGPHREG过程步进行数据分析。过程步进行数据分析。