如何用spss17.0进行二元和多元logistic回归分析

如何用spss17.0 进行二元和多元logistic 回归分析一、二元logistic 回归分析二元logistic 回归分析的前提为因变量是可以转化为0、1 的二分变量，如：死亡或者生存，男性或者女性，有或无，Yes 或No，是或否的情况。

下面以医学中不同类型脑梗塞与年龄和性别之间的相互关系来进行二元logistic 回归分析。

（一）数据准备和SPSS 选项设置第一步，原始数据的转化：如图1-1 所示，其中脑梗塞可以分为ICAS、ECAS 和NCAS 三种，但现在我们仅考虑性别和年龄与ICAS 的关系，因此将分组数据ICAS、ECAS 和NCAS 转化为1、0 分类，是ICAS 赋值为1，否赋值为0。

年龄为数值变量，可直接输入到spss中，而性别需要转化为（1、0）分类变量输入到spss当中，假设男性为1，女性为0，但在后续分析中系统会将1，0 置换（下面还会介绍），因此为方便期间我们这里先将男女赋值置换，即男性为“0”，女性为“1”。

图1-1第二步：打开“二值Logistic 回归分析”对话框：沿着主菜单的“分析（Analyze）→回归（Regression）→二元logistic（Binary Logistic）” 的路径（图1-2）打开二值Logistic 回归分析选项框（图1-3）。

如图1-3左侧对话框中有许多变量，但在单因素方差分析中与ICAS显著相关的为性别、年龄、有无高血压，有无糖尿病等（P<0.05），因此我们这里选择以性别和年龄为例进行分析。

图1-2图1-3在图1-3中，因为我们要分析性别和年龄与ICAS的相关程度，因此将ICAS选入因变量（Dependent）中，而将性别和年龄选入协变量（Covariates）框中，在协变量下方的“方法（Method）”一栏中，共有七个选项。

采用第一种方法，即系统默认的强迫回归方法（进入“Enter”）。

接下来我们将对分类（Categorical），保存（Save），选项（Options）按照如图1-4、1-5、1-6中所示进行设置。

在“分类”对话框中，因为性别为二分类变量，因此将其选入分类协变量中，参考类别为在分析中是以最小数值“0（第一个）”作为参考，还是将最大数值“1（最后一个）”作为参考，这里我们选择第一个“0”作为参考。

在“存放”选项框中是指将不将数据输出到编辑显示区中。

在“选项”对话框中要勾选如图几项，其中“exp（B）的CI（X）”一定要勾选，这个就是输出的OR和CI值，后面的95%为系统默认，不需要更改。

图 1-5 保存另外在“选项”对话框中，“输出”一栏中，系统默认为“在每个步骤中”，这里更改 为“在最后一个步骤中”，即：输出结果将仅仅给出最终结果，而省略每一步的计算过程。

由于我们采用强迫回归，逐步回归概率选项可以不管此外还有一个选项需要说明。

一是分类临界值(Classification cutoff )，默认值为0.5， 即按四舍五入的原则将概率预测值化为0 或者1。

如果将数值改为0.6，则大于等于0.6 的概 率值才表示为1，否则为0。

其情况余依此类推。

二是最大迭代值(Maximum Iterations )， 规定系统运算的迭代次数，默认值为20 次，为安全起见，我们将迭代次数增加到50。

原因 是，有时迭代次数太少，计算结果不能真正收敛。

三是模型中包括常数项(Include constant in model )，即模型中保留截距。

除了迭代次数之外，其余两个选项均采用系统默认值。

完成后，点击各项中“继续(Continue )”按钮。

返回图1-3，单击“确定”按钮。

(二)结果解读其他结果参照文章《利用SPSS 进行Logistic 回归分析》中解读，这里重点将两点： 第一，分类变量编码(图1-7)，由于这里包括性别分类变量，而我们对性别赋值为1 和0，但在spss 中系统会默认把我们的数值进行置换，即1→参数编码0，0→参数编码1，而 最终输出结果是以1来计算的，而0为参考数据。

所以这也就是为什么我么之前要对研究组男 性的赋值进行置换了。

如果男性为1那么spss 中最终输出的将是女性的分析结果。

图 1-7第二，最终输出数据(图1-8)在该结果中，Exp(B)即为文献中提及的OR 值，而EXP(B) 的95%C.I.即为文献中提及的CI 值。

其中Exp(B)表示某因素(自变量)内该类别是其相应参 考类别具有某种倾向性的倍数。

而有的文献中提到的Crode OR 和Adjust OR 则分别为单因素 优势率(Crode odds ratio )和多因素优势率(Adjust odds ratio )，即仅对性别单个变 量的单因素分析或者对性别和年龄等多个变量进行多因素分析后所得到的不同结果。

CI 则为 可信区间(Confidence interval)。

Sig.即我们常说的P 值，P<0.05为显著(无效假说不成 立，具有统计学意义)，P>0.05为不显著(无效假说成立，不具有统计学意义)。

图 1-6 选项图 1-4 分类图1-8二、多项（多元、多分类、Multinomial）logistic 回归分析前面讲的二元logistic回归分析仅适合因变量Y 只有两种取值（二分类）的情况，当Y 具有两种以上的取值时，就要用多项logistic 回归（Mutinomial Logistic Regression）分析了。

这种分析不仅可以用于医疗领域，也可以用于社会学、经济学、农业研究等多个领域。

如不同阶段（初一、初二、初三）学生视力下降程度，不同龋齿情况（轻度、中度、重度）下与刷牙、饮食、年龄的关系等。

下面我们以图1-2 中，对apoba1（ApoB/AI）项中数值做四分位数后，将病人的ApoB/AI 的比值划分为低、较低、中、高四个分位后利用多项logistic 回归分析其与ICAS 之间的相互关系。

首先来做四分位数，很多人在做四分位数的时候都是自己算出来的，其实在SPSS 里面给出了做四分位数的程度即分析（Aanlyze）→描述统计（Descriptive Statistics）→频率（Frequencies）。

打如图2-1 开频率对话框。

将我们要分析的数值变量Apoba1 选入到变量对话框中。

图2-1选择统计量，按照图2-2 中勾选四分位数选项，其他选项按照自己需要勾选，然后点击图2-1 中的确定按钮，开始运算。

在图2-3 中可以读取我们的四分位数值。

图中百分数表示的是对该变量做的四分位数的百分比，25表示前25%的， 50表示前 50%的，75表示前 75%的。

每一项对应的后面数值即为相应的四分 位数，如 0.5904，即为前 25%的个体与后 75%个体的分位数。

图 2-3按照如上方法得出 ApoB/AI 的比率后我们可以把该比值划分为四个区 间，即当 ApoB/AI 的比率＜0.5904 为低、当0.5904≤ApoB/AI 的比率≤0.88 时为较低、当 0.89≤ApoB/AI 的比率≤1.0886 时为中，当 ApoB/AI 的比 率＞1.0886 时为高。

然后将这一划分如图 1-1 中“四分位数”一项用分类数值 表示即 1 代表低，2 代表较低，3 代表中，4 代表高。

这里还要强调的是我们要 研究其与 ICAS 之间的相互关系，那么我们需要将其设为二分类变量，即是 ICAS 的情况为 1，否则为 0，但多项 logistic 回归分析也会将 1，0 置换，所以我们 需要在这里将我们需要研究的情况置换为 0，然后将其他置换为 1。

下面就可以 进行多项 logistic 回归分析了。

如图2-4打开多项 logistic 回归分析对话框（图 2-5）。

如图 2-5 所示，在”因变量”中选入刚才我们输入的四分位数分类变量，在因 子中输入分类变量 ICAS （这里一定是分类变量，可以是一个也可以是多个），在 “协变量”中输入数值变量如年龄（这里一定是数值变量，可以是一个也可以是 多个），但因本次没有对年龄进行分析，仅对 ICAS 进行了单因素分析，所以我 们把年龄移出协变量选项。

在 SPSS 中对因变量的定义是，如果因变量 Y 有 J 个值（即 Y 有 J 类），以 其中一个类别作为参考类别，其他类别都同他相比较生成J-1 个冗余的Logit 变 换模型，而作为参考类别的其模型中所有系数均为 0。

在SPSS 中可以对所选因 变量的参考类别进行设置，如图 2-5 在因变量对话框下有一“参考类别”选项。

点击后会弹出图 2-6对话框。

在该对话框中图 2-2我们选中设定，输入数值1，这代表我们以分类数值 1 所代表的类别作为参考类别，即最低数值作为参考类别。

单击 继续。

当然也可以选择“第一类别”和“最后类别”，入选中分别表示以最低数 值或最高数值作为参考类别。

其他设置与二元 Logistic 分析相似，将我们要输 出的项勾选即可，点击图 2-5 中确定，输出数据。

输出数据基本与二元 Logistic 分析相似，我们重点讲下最后一项“参考估 计”，如图 2-7 所示，其中参考类别为 ICAS=1 的分类情况，而其中的 ICAS=0 分为 2、3、4 三种，分别给出了 ICAS=0 时的数值。

而其中Exp(B)(即OR 值)表示某因素（自变量）内该类别是其相应参考类别具有图 2-4图 2-5 图 2-6某种倾向性的倍数。

如Exp（B）=2.235 时，即表示在较轻这一类别下ICAS 患者数为其他类别（ECAS 和NCAS ）的 2.235 倍。

这里面的显著水平即为P 值。

这里要强调的是，一些文献中在输出数据的时候经常会给出“Referent （参考）”项，这里的Referent，即为我们这里所选的参考类别1，因为1 作为参考类别，所以其所有数值为0，即无数据输出。

因此在文中需标注其为Referent。