如何用spss进行二元和多元logistic回归分析.doc

如何用spss进行二元和多元logistic回归分析

———————————————————————————————— 作者：
———————————————————————————————— 日期：

9

个人收集整理 勿做商业用途
如何用spss17.0进展二元和多元ｌｏgisｔic回归分析
一、二元loｇisｔic回归分析
二元ｌoｇisｔic回归分析的前提为因变量是可以转化为０、1的二分变量,如：死亡或者生存，男性或者女性，有或无，Yes或No，是或否的情况。
下面以医学中不同类型脑堵塞与年龄和性别之间的相互关系来进展二元ｌogｉｓtiｃ回归分析。
(一〕数据准备和ＳPSS选项设置
图 1-1
第一步，原始数据的转化:如图１－１所示,其中脑堵塞可以分为IＣAＳ、EＣAS和ＮCＡS三种,但现在我们仅考虑性别和年龄与ICAS的关系，因此将分组数据ICＡＳ、ＥＣAＳ和ＮCAS转化为1、0分类,是IＣAS赋值为1,否赋值为0。年龄为数值变量,可直接输入到spsｓ中，而性别需要转化为〔1、０)分类变量输入到ｓpｓs当中，假设男性为1，女性为０,但在后续分析中系统会将1,0置换〔下面还会介绍)，因此为方便期间我们这里先将男女赋值置换,即男性为“0〞，女性为“1〞。
第二步:翻开“二值Ｌoｇistｉc 回归分析〞对话框:
沿着主菜单的“分析(Ａnａlyze)→回归〔Regreｓsion〕→二元logistic(Bｉnaｒy　Logistic)〞的路径(图1-2)翻开二值Lｏgistic 回归分析选项框〔图1-3)。
如图1-３左侧对话框中有许多变量,但在单因素方差分析中与ICAS显著相关的为性别、年龄、有无高血压,有无糖尿病等〔P＜０.05)，因此我们这里选择以性别和年龄为例进展分析。图1－3

图1-3
图1-2
在图1-3中，因为我们要分析性别和年龄与ＩCAS的相关程度，因此将ICＡS选入因变量(Dｅpendｅnｔ)中,而将性别和年龄选入协变量〔Covariatｅs)框中,在协变量下方的“方法(Ｍｅtｈｏd〕〞一栏中，共有七个选项。采用第一种方法,即系统默认的强迫回归方法(进入“Eｎter〞)。
接下来我们将对分类(Categoｒicａl〕,保存〔Save),选项(Optionｓ〕按照如图１-4、1-5、1-6中所示进展设置。在“分类〞对话框中，因为性别为二分类变量,因此将其选入分类协变量中，参考类别为在分析中是以最小数值“0〔第一个)〞作为参考，还是将最大数值“１〔最后一个)〞作为参考，这里我们选择第一个“0〞作为参考。在“存放〞选项框中是指将不将数据输出到编辑显示区中。在“选项〞对话框中要勾选如图几项，其中“exp〔B)的CI(X)〞一定要勾选,这个就是输出的OR和CI值,后面的95％为系统默认,不需要更改。
图1-4 分类
图1-5 保存
图1-6 选项
另外在“选项〞对话框中,“输出〞一栏中,系统默认为“在每个步骤中〞，这里更改为“在最后一个步骤中〞，即：输出结果将仅仅给出最终结果,而省略每一步的计算过程。由于我们采用强迫回归，逐步回归概率选项可以不管
此外还有一个选项需要说明。一是分类临界值〔Ｃlａsｓｉfｉcatiｏn　cutoff〕，默认值为０.5，即按四舍五入的原那么将概率预测值化为0 或者1。如果将数值改为0．６,那么大于等于0.6 的概率值才表示为1，否那么为0。其情况余依此类推。二是最大迭代值〔Maｘimum Iteｒａｔｉｏnｓ〕,规定系统运算的迭代次数,默认值为2０ 次,为平安起见，我们将迭代次数增加到50。原因是，有时迭代次数太少，计算结果不能真正收敛。三是模型中包括常数项〔Ｉncｌudｅ coｎstant in ｍｏdel),即模型中保存截距。除了迭代次数之外，其余两个选项均采用系统默认值。
完成后,点击各项中“继续(Contｉnue〕〞按钮。返回图1-3，单击“确定〞按钮。
〔二〕结果解读
其他结果参照文章?利用ＳPSS进展Logistｉc回归分析?中解读，这里重点将两点：
第一,分类变量编码〔图１-7〕,由于这里包括性别分类变量，而我们对性别赋值为1和0，但在spｓs中系统会默认把我们的数值进展置换，即１→参数编码0，0→参数编码1，而最终输出结果是以1来计算的,而0为参考数据。所以这也就是为什么我么之前要对研究组男性的赋值进展置换了。如果男性为１那么spss中最终输出的将是女性的分析结果。
图1-7
第二，最终输出数据〔图1-８〕在该结果中,Exp〔Ｂ)即为文献中提及的OＲ值,而EXP(B)的9５%Ｃ.I.即为文献中提及的ＣI值。其中Exp(B〕表示某因素〔自变量〕内该类别是其相应参考类别具有某种倾向性的倍数。而有的文献中提到的Ｃroｄe　ＯＲ和Adｊｕst　ＯR那么分别为单因素优势率(Cｒｏdｅ ｏddｓ rａtio〕和多因素优势率〔Ａdjｕｓｔ odｄs rａtｉo〕,即仅对性别单个变量的单因素分析或者对性别和年龄等多个变量进展多因素分析后所得到的不同结果。ＣＩ那么为可信区间〔Cｏnfiｄencｅ inｔｅｒvａl)。Sｉg.即我们常说的Ｐ值，P<0.0５为显著〔无效假说不成立，具有统计学意义〕，P>0.0５为不显著(无效假说成立,不具有统计学意义〕。
图1-8

二、多项(多元、多分类、Multinomial)ｌoｇｉｓtｉｃ回归分析
前面讲的二元lｏgistｉｃ回归分析仅适合因变量Y只有两种取值〔二分类)的情况,当Y具有两种以上的取值时,就要用多项loｇistiｃ回归(Ｍuｔinｏmial Logistic Regrｅsｓiｏｎ〕分析了。这种分析不仅可以用于医疗领域，也可以用于社会学、经济学、农业研究等多个领域。如不同阶段(初一、初二、初三)学生视力下降程度，不同龋齿情况〔轻度、中度、重度)下与刷牙、饮食、年龄的关系等。
下面我们以图1-２中,对aｐoｂａ1〔ApoＢ／AI〕项中数值做四分位数后，将病人的ApoB/AI的比值划分为低、较低、中、高四个分位后利用多项logｉstic回归分析其与ICAＳ之间的相互关系。
图2-1
首先来做四分位数,很多人在做四分位数的时候都是自己算出来的，其实在SPSS里面给出了做四分位数的程度即分析〔Aanlyｚe〕→描述统计〔Descrｉptivｅ Staｔｉstics〕→频率〔Freｑｕeｎcies)。打如图2-1开频率对话框。将我们要分析的数值变量Aｐoba1选入到变量对话框中。
选择统计量,按照图2-2中勾选四分位数选项，其他选项按照自己需要勾选，然后点击图２－１中确实定按钮,开场运算。在图２-3中可以读取我们的四分位数值。图中百分数表示的是对该变量做的四分位数的百分比,25表示前2５%的,５0表示前５0%的,75表示前７5%的。每一项对应的后面数值即为相应的四分位数,图2-2
图2-3
如0.5904,即为前２5%的个体与后７5％个体的分位数。
　　 按照如上方法得出ApｏB/ＡI的比率后我们可以把该比值划分为四个区间，即当ApｏB/AI的比率<0.5９０4为低、当０.5９04≤AｐoB/ＡＩ的比率≤０.88时为较低、当0.８9≤ApoB/AＩ的比率≤1.０88６时为中,当AｐoＢ/ＡI的比率>1．0886时为高。然后将这一划分如图1－１中“四分位数〞一项用分类数值表示即１代表低,2代表较低，3代表中，4代表高。这里还要强调的是我们要研究其与IＣＡS之间的相互关系,那么我们需要将其设为二分类变量,即是ICAＳ的情况为1，否那么为0,但多项ｌogistiｃ回归分析也会将1,0置换，所以我们需要在这里将我们需要研究的情况置换为０,然后将其他置换为１。下面就可以进展多项logistiｃ回归分析了。如图２-4翻开多项ｌogiｓtｉc回归分析对话框〔图2-5〕。
如图2－５所示,在〞因变量〞中选入刚刚我们输入的四分位数分类变量，在因子中输入分类变量ＩCＡS(这里一定是分类变量，可以是一个也可以是多个)，在“协变量〞中输入数值变量如年龄(这里一定是数值变量，可以是一个也可以是多个),但因本次没有对年龄进展分析，仅对ＩＣAＳ进展了单因素分析,所以我们把年龄移出协变量选项。
在SＰＳS中对因变量的定义是,如果因变量Ｙ有J个值〔即Y有Ｊ类〕,以其中一个类别作为参考类别，其他类别都同他相比拟生成J－1个冗余的Logiｔ变换模型，而作为参考类别的其模型中所有系数均为0。在SPSS中可以对所选因变量的参考类别进展设置,如图2－5在因变量对话框下有一“参考类别〞选项。点击后会弹出图2-６对话框。在该对话框中我们选中设定，输入数值1，这代表我们以分类数值1所代表的类别作为参考类别，即最低数值作为参考类别。单击继续。当然也可以选择“第一类别〞和“最后类别〞，入选中分别表示以最低数值或最高数值作为参考类别。其他设置与二元Lｏgｉstｉc分析相似,将我们要输出的项勾选即可，点击图２-5中确定,输出数据。
图2-4
图2-5
图2-6
输出数据根本与二元Logistiｃ分析相似，我们重点讲下最后一项“参考估计〞,如图2－7所示,其中参考类别为ICAS＝1的分类情况,而其中的ICAS=0分为2、３、4三种，分别给出了ＩCAＳ=0时的数值。而其中Eｘp(B)〔即OR值〕表示某因素〔自变量〕内该类别是其相应参考类别具有某种倾向性的倍数。如Exp〔B)=2．２35时，即表示在较轻这一类别下ICＡS患者数为其他类别(ＥCＡS和NCＡＳ〕的2．2３5倍。这里面的显著水平即为Ｐ值。
这里要强调的是,一些文献中在输出数据的时候经常会给出“Ｒeferenｔ〔参考)〞项,这里的Reｆerenｔ,即为我们这里所选的参考类别1，因为１作为参考类别，所以其所有数值为0，即无数据输出。因此在文中需标注其为Refeｒｅnt。
图2-7