logistic回归模型 SPSS例析

合集下载

多因素logistic回归分析spss

多因素logistic回归分析spssLogistic回归分析是一种用来研究影响离散变量的因素的方法，该方法的输出是一个logistic模型，这一模型可以用于预测变量的值，即预测该变量的值有多高的概率会取各种可能的取值。

简言之，logistic回归分析的主要目的是把客观的结果（例如，是否改变某个政策，是否感染某种疾病等）变成可预测的离散变量，以便分析影响客观结果的各种因素。

Spss可以提供多因素logistic回归分析，这种分析可用于识别影响离散变量（例如，是否改变某个政策，是否感染某种疾病等）的多个因素之间的关联。

该分析需要有一个组合变量作为自变量，以及一个离散变量作为因变量。

例如，如果您要研究性别和年龄两个因素如何影响某种疾病的发生率，那么性别和年龄两个因素就是组合变量，而疾病的发生率则是因变量。

1.建立变量和分类（上述示例中需要建立性别和年龄两个变量，以及分类变量的可能的取值）。

2.执行logistic回归分析。

打开spss，并在“分析”菜单中打开多元分析，然后点击“逻辑回归”，并选择您要研究的变量和分类。

3.生成回归模型和检验其统计学意义。

在spss中，您可以使用类似“回归系数”之类的描述性统计学方法来估算回归模型，并可以使用“p-值”来判断回归模型中各变量的统计学意义。

4.Interpret模型。

根据p值判断各变量的统计学意义，进而分析影响离散变量的多个因素之间的关联。

四、总结Logistic回归分析是一种用来研究影响离散变量的因素的方法，spss可以提供多因素logistic回归分析，这种分析可用于识别影响离散变量的多个因素之间的关联，spss中步骤：建立变量和分类，执行logistic回归分析，生成回归模型和检验其统计学意义，Interpret模型。

SPSS数据分析—配对Logistic回归模型

Logistic回归模型也可以用于配对资料，但是其分析方法和操作方法均与之前介绍的不同，具体表现在以下几个方面1.每个配对组共有同一个回归参数，也就是说协变量在不同配对组中的作用相同2.常数项随着配对组变化而变化，反映了非实验因素在配对组中的作用，但是我们并不关心其大小，因此在拟合时采用条件似然函数代替了一般似然函数，从而在拟合中消去了反映层因素的参数。

SPSS中没有直接拟合配对Logistic回归模型的过程，需要对数据进行一些处理，采用其他方法进行拟合，拟合方法有变量差值拟合和COX模型一、变量差值拟合只适用于1:1配对，通过求出同一对中案例组与对照组多有变量的差值，对差值进行不含常数项的无序多分类Logistic回归模型拟合来达到目的例：收集了一组数据，希望分析服用雌激素与子宫内膜癌之间的关系，除了研究因素之外，还额外收集了两个变量，数据为配对数据，1为病例，0为对照，case为是否患病，也就是因变量采用变量差值进行拟合，首先求出所有变量之间的差值，可以使用计算变量过程，但是该过程每次只能处理一个变量，比较麻烦，我们使用语法编辑器进行程序编写，如下全部选中之后运行，在原数据中就会依次出现新生成的差值变量，接下来，我们对这些差值变量进行无序多分类Logistic回归分析—回归—多项Logistic==================================================二、分层Cox模型该方法最常用来进行生存函数估计，但是由于在拟合方法上和配对Logistic模型一致，因此也可以用来拟合配对Logistic回归模型，它不仅可以拟合1:1配对，还可以拟合1:r，n:m配对，应用范围比较广。

在数据组成上，和变量差值拟合有所不同，需要给每个个案一个虚拟生存时间，默认案例组比对照组生存时间段，具体值不限，两个时间差距大小也不限。

案例发生算为失效事件，对照组为删失，并且对子号作为分层因素，我们还是通过变量差值拟合的案例来进行说明我们将原数据调整如下我们将原来并排排列的案例组和对照组合并为竖列，并且增加虚拟生存时间变量time，案例组为1，对照组为2，接下来按照Cox回归模型进行拟合分析—生存函数—Cox回归。

SPSS实验8-二项Logistic回归分析

SPSS作业8：二项Logistic回归分析为研究和预测某商品消费特点和趋势，收集到以往胡消费数据.数据项包括是否购买，性别,年龄和收入水平。

这里采用Logistic回归的方法，是否购买作为被解释变量（0/1二值变量）,其余各变量为解释变量，且其中性别和收入水平为品质变量，年龄为定距变量。

变量选择采用Enter方法，性别以男为参照类，收入以低收入为参照类。

（一）基本操作：（1)选择菜单Analyz e－Regression－Binary Logistic；（2)选择是否购买作为被解释变量到Dependent框中，选其余各变量为解释变量到Covariates框中，采用Enter方法，结果如下：消费的二项Logistic分析结果（一）（强制进入策略）Categorical Variables CodingsFrequency Parameter coding （1) (2)收入低收入132 .000 .000中收入144 1.000 。

000高收入155 。

000 1。

000性别男191 。

000女240 1.000分析：上表显示了对品质变量产生虚拟变量的情况，产生的虚拟变量命名为原变量名（编码)。

可以看到，对收入生成了两个虚拟变量名为Income（1）和Income（2），分别表示是否中收入和是否高收入，两变量均为0时表示低收入；对性别生成了一个虚拟变量名为Gedder（1），表示是否女，取值为0时表示为男。

消费的二项Logistic 分析结果（二）（强制进入策略）Block 0: Beginning BlockClassification Table a,bObserved Predicted是否购买 Percentage Correct不购买购买Step 0是否购买不购买 269 0 100。

购买162。

0 Overall Percentage62。

4a 。

Constant is included in the model 。

SPSS—回归—二元Logistic回归案例分析

SPSS—回归—⼆元Logistic回归案例分析数据分析真不是⼀门省油的灯，搞的⼈晕头转向，⽽且涉及到很多复杂的计算，还是书读少了，⼩学毕业的我，真是死了不少脑细胞，学习⼆元Logistic回归有⼀段时间了，今天跟⼤家分享⼀下学习⼼得，希望多指教！⼆元Logistic，从字⾯上其实就可以理解⼤概是什么意思，Logistic中⽂意思为“逻辑”但是这⾥，并不是逻辑的意思，⽽是通过logit变换来命名的，⼆元⼀般指“两种可能性”就好⽐逻辑中的“是”或者“否”⼀样，Logistic 回归模型的假设检验——常⽤的检验⽅法有似然⽐检验（likelihood ratio test）和 Wald检验）似然⽐检验的具体步骤如下：1：先拟合不包含待检验因素的Logistic模型，求对数似然函数值INL02：再拟合包含待检验因素的Logistic模型，求新的对数似然函数值InL13:最后⽐较两个对数似然函数值的差异，若两个模型分别包含l个⾃变量和P个⾃变量，记似然⽐统计量G的计算公式为 G=2(InLP - InLl). 在零假设成⽴的条件下，当样本含量n较⼤时，G统计量近似服从⾃由度为 V = P-l 的 x平⽅分布，如果只是对⼀个回归系数（或⼀个⾃变量）进⾏检验，则 v=1.wald 检验，⽤u检验或者X平⽅检验，推断各参数βj是否为0，其中u= bj / Sbj, X的平⽅=（bj / Sbj), Sbj 为回归系数的标准误这⾥的“⼆元”主要针对“因变量”所以跟“曲线估计”⾥⾯的Logistic曲线模型不⼀样，⼆元logistic回归是指因变量为⼆分类变量是的回归分析，对于这种回归模型，⽬标概率的取值会在（0-1），但是回归⽅程的因变量取值却落在实数集当中，这个是不能够接受的，所以，可以先将⽬标概率做Logit变换，这样它的取值区间变成了整个实数集，再做回归分析就不会有问题了，采⽤这种处理⽅法的回归分析，就是Logistic 回归设因变量为y, 其中“1” 代表事件发⽣， “0”代表事件未发⽣，影响y的 n个⾃变量分别为 x1, x2 ,x3 xn等等记事件发⽣的条件概率为 P那么P= 事件未发⽣的概理为 1-P事件发⽣跟”未发⽣的概率⽐为( p / 1-p ) 事件发⽣⽐，记住Odds将Odds做对数转换，即可得到Logistic回归模型的线性模型：还是以教程“blankloan.sav"数据为例，研究银⾏客户贷款是否违约（拖⽋）的问题，数据如下所⽰：上⾯的数据是⼤约700个申请贷款的客户，我们需要进⾏随机抽样，来进⾏⼆元Logistic回归分析，上图中的“0”表⽰没有拖⽋贷款，“1”表⽰拖⽋贷款，接下来，步骤如下：1：设置随机抽样的随机种⼦，如下图所⽰：选择“设置起点”选择“固定值”即可，本⼈感觉200万的容量已经⾜够了，就采⽤的默认值，点击确定，返回原界⾯、2：进⾏“转换”—计算变量“⽣成⼀个变量（validate)，进⼊如下界⾯：在数字表达式中，输⼊公式：rv.bernoulli（0.7），这个表达式的意思为：返回概率为0.7的bernoulli分布随机值如果在0.7的概率下能够成功，那么就为1，失败的话，就为"0"为了保持数据分析的有效性，对于样本中“违约”变量取缺失值的部分，validate变量也取缺失值，所以，需要设置⼀个“选择条件”点击“如果”按钮，进⼊如下界⾯：如果“违约”变量中，确实存在缺失值，那么当使⽤"missing”函数的时候，它的返回值应该为“1”或者为“true"，为了剔除”缺失值“所以，结果必须等于“0“ 也就是不存在缺失值的现象点击 ”继续“按钮，返回原界⾯，如下所⽰：将是“是否曾经违约”作为“因变量”拖⼊因变量选框，分别将其他8个变量拖⼊“协变量”选框内，在⽅法中，选择：forward.LR⽅法将⽣成的新变量“validate" 拖⼊"选择变量“框内，并点击”规则“设置相应的规则内容，如下所⽰：设置validate 值为1，此处我们只将取值为1的记录纳⼊模型建⽴过程，其它值（例如：0）将⽤来做结论的验证或者预测分析，当然你可以反推，采⽤0作为取值记录点击继续，返回，再点击“分类”按钮，进⼊如下页⾯在所有的8个⾃变量中，只有“教育⽔平”这个变量能够作为“分类协变量” 因为其它变量都没有做分类，本例中，教育⽔平分为：初中，⾼中，⼤专，本科，研究⽣等等, 参考类别选择：“最后⼀个” 在对⽐中选择“指⽰符” 点击继续按钮，返回再点击—“保存”按钮，进⼊界⾯：在“预测值"中选择”概率，在“影响”中选择“Cook距离” 在“残差”中选择“学⽣化”点击继续，返回，再点击“选项”按钮，进⼊如下界⾯：分析结果如下：1：在“案例处理汇总”中可以看出：选定的案例489个，未选定的案例361个，这个结果是根据设定的validate = 1得到的，在“因变量编码”中可以看出“违约”的两种结果“是”或者“否” 分别⽤值“1“和“0”代替，在“分类变量编码”中教育⽔平分为5类，如果选中“为完成⾼中，⾼中，⼤专，⼤学等，其中的任何⼀个，那么就取值为 1，未选中的为0，如果四个都未被选中，那么就是”研究⽣“ 频率分别代表了处在某个教育⽔平的个数，总和应该为 489个1：在“分类表”中可以看出：预测有360个是“否”（未违约）有129个是“是”（违约）2：在“⽅程中的变量”表中可以看出：最初是对“常数项”记性赋值，B为-1.026，标准误差为：0.103那么wald =( B/S.E)²=(-1.026/0.103)² = 99.2248, 跟表中的“100.029⼏乎接近，是因为我对数据进⾏的向下舍⼊的关系，所以数据会稍微偏⼩，B和Exp(B) 是对数关系，将B进⾏对数抓换后，可以得到：Exp(B) = e^-1.026 = 0.358, 其中⾃由度为1， sig为0.000，⾮常显著1：从“不在⽅程中的变量”可以看出，最初模型，只有“常数项”被纳⼊了模型，其它变量都不在最初模型内表中分别给出了，得分，df , Sig三个值, ⽽其中得分（Score)计算公式如下：（公式中（Xi- X¯) 少了⼀个平⽅）下⾯来举例说明这个计算过程：(“年龄”⾃变量的得分为例）从“分类表”中可以看出：有129⼈违约，违约记为“1” 则违约总和为 129，选定案例总和为489那么： y¯ = 129/489 = 0.2638036809816x¯ = 16951 / 489 = 34.664621676892所以：∑(Xi-x¯)² = 30074.9979y¯（1-y¯）=0.2638036809816 *（1-0.2638036809816 ）=0.19421129888216则：y¯（1-y¯）* ∑(Xi-x¯)² =0.19421129888216 * 30074.9979 = 5 840.9044060372则：[∑Xi(yi - y¯）]^2 = 43570.8所以：=43570.8 / 5 840.9044060372 = 7.4595982010876 = 7.46 （四舍五⼊）计算过程采⽤的是在 EXCEL ⾥⾯计算出来的，截图如下所⽰：从“不在⽅程的变量中”可以看出，年龄的“得分”为7.46，刚好跟计算结果吻合！！答案得到验证~1:从“块1” 中可以看出：采⽤的是：向前步进的⽅法，在“模型系数的综合检验”表中可以看出：所有的SIG ⼏乎都为“0” ⽽且随着模型的逐渐步进，卡⽅值越来越⼤，说明模型越来越显著，在第4步后，终⽌，根据设定的显著性值和⾃由度，可以算出卡⽅临界值，公式为：=CHIINV(显著性值,⾃由度) ，放⼊excel就可以得到结果2：在“模型汇总“中可以看出：Cox&SnellR⽅和 Nagelkerke R⽅拟合效果都不太理想，最终理想模型也才：0.305 和 0.446，最⼤似然平⽅的对数值都⽐较⼤，明显是显著的似然数对数计算公式为：计算过程太费时间了，我就不举例说明计算过程了Cox&SnellR⽅的计算值是根据：1：先拟合不包含待检验因素的Logistic模型，求对数似然函数值INL0 （指只包含“常数项”的检验）2：再拟合包含待检验因素的Logistic模型，求新的对数似然函数值InLB （包含⾃变量的检验）再根据公式：即可算出：Cox&SnellR⽅的值！提⽰：将Hosmer 和 Lemeshow 检验和“随机性表” 结合⼀起来分析1：从 Hosmer 和 Lemeshow 检验表中，可以看出：经过4次迭代后，最终的卡⽅统计量为：11.919，⽽临界值为：CHINV(0.05,8) =15.507卡⽅统计量< 临界值，从SIG ⾓度来看： 0.155 > 0.05 , 说明模型能够很好的拟合整体，不存在显著的差异。

SPSS数据分析—二分类Logistic回归模型

对于分类变量，我们知道通常使用卡方检验，但卡方检验仅能分析因素的作用，无法继续分析其作用大小和方向，并且当因素水平过多时，单元格被划分的越来越细，频数有可能为0，导致结果不准确，最重要的是卡方检验不能对连续变量进行分析。

使用线性回归模型可以解决上述的部分问题，但是传统的线性模型默认因变量为连续变量，当因变量为分类变量时，传统线性回归模型的拟合方法会出现问题，因此人们继续发展出了专门针对分类变量的回归模型。

此类模型采用的基本方法是采用变量变换，使其符合传统回归模型的要求。

根据变换的方法不同也就衍生出不同的回归模型，例如采用Logit变换的Logistic回归模型，采用Probit变换的Probit回归模型等，相比之下，Logistic是使用最为广泛的针对分类数据的回归模型。

Logistic回归模型的适用条件1.因变量为二分类变量或是某事件的发生率2.自变量与Logit变换后的因变量呈线性关系3.残差合计为0，且服从二项分布4.各观测值之间独立由于Logistic回归模型的残差项服从二项分布而不是正态分布，因此不能使用最小二乘法进行参数估计，而是要使用最大似然法。

和其他回归分析一样，Logistic回归也放在分析—回归过程下面，下面我们通过一个例子来说明具体操作收集了一组数据，希望通过这些数据分析出低出生体重儿的影响因素，数据如下可见，数据集中变量比较多，且数据类型丰富，因变量为二分类变量Low，有两个水平：0-正常体重，1-低出生体重，我们先做一个最简单的单变量Logistic 回归，只考虑smoke这个因素分析—回归—二元Logistic回归前面我们只引入了一个自变量，可以看到模型的效果并不理想，而且Logistic 回归和传统回归模型一样，也可以引入多个自变量并且可以对自变量进行筛选，尽量引入对因变量存在强影响的自变量，下面我们继续加入自变量并进行筛选。

spsslogistic回归分析结果解读

spsslogistic回归分析结果解读
本文分析了使用SPSS Logistic回归分析的结果，以了解不同变量之间
是否存在潜在关系。

Logistic回归是一种用于预测调查中的变量组合能够预测调查的结果的
机器学习技术。

在这种情况下，我们使用Logistic回归来预测一个变量
（假设为购买行为）和其他变量（价格，品牌认知度等）之间的关系。

特别是，我们可以评估价格是否是客户决定购买商品的重要影响因素。

SPSS Logistic回归分析的结果表明，在本例中，我们发现价格是一个
重要的影响因素。

我们看到，价格的变化程度会影响客户购买商品的可能性：客户可能更愿意购买相对较低的价格，而对于较高的价格则更不可能购买。

此外，品牌认知度也会影响客户是否愿意购买：客户对品牌认知度越高，购
买概率越高。

这可能是因为客户更倾向于信任已经熟悉的品牌而忽略未熟悉
的品牌，或者可能是因为客户更了解该品牌的商品及其优缺点，因此可以作
出的更明智的购买决策。

因此，本次分析表明，价格和品牌认知度在客户决定购买商品时都有重
要的影响。

商家应考虑这些因素，以确保它们的产品在客户面前具有足够的
吸引力和优势，使其愿意购买。

详解利用SPSS进行Logistic_回归分析

第8 章利用SPSS 进行Logistic 回归分析现实中的很多现象可以划分为两种可能，或者归结为两种状态，这两种状态分别用0和1 表示。

如果我们采用多个因素对0－1 表示的某种现象进行因果关系解释，就可能应用到logistic 回归。

Logistic 回归分为二值logistic 回归和多值logistic 回归两类。

首先用实例讲述二值logistic 回归，然后进一步说明多值logistic 回归。

在阅读这部分内容之前，最好先看看有关SPSS 软件操作技术的教科书。

§8.1 二值logistic 回归8.1.1 数据准备和选项设置我们研究2005 年影响中国各地区城市化水平的经济地理因素。

城市化水平用城镇人口比重表征，影响因素包括人均GDP、第二产业产值比重、第三产业产值比重以及地理位置。

地理位置为名义变量，中国各地区被分别划分到三大地带：东部地带、中部地带和西部地带。

我们用各地区的地带分类代表地理位置。

第一步：整理原始数据。

这些数据不妨录入Excel 中。

数据整理内容包括两个方面：一是对各地区按照三大地带的分类结果赋值，用0、1 表示，二是将城镇人口比重转换逻辑值，变量名称为“城市化”。

以各地区2005 年城镇人口比重的平均值45.41%为临界值，凡是城镇人口比重大于等于45.41%的地区，逻辑值用Yes 表示，否则用No 表示（图8-1-1）图8-1-1 原始数据（Excel 中，局部）将数据拷贝或者导入SPSS 的数据窗口（Data View）中（图8-1-2）。

图8-1-2 中国31 个地区的数据（SPSS 中，局部）第二步：打开“聚类分析”对话框。

沿着主菜单的“Analyze→Regression→Binary Logistic K”的路径（图8-1-3）打开二值Logistic 回归分析选项框（图8-1-4）。

图8-1-3 打开二值Logistic 回归分析对话框的路径对数据进行多次拟合试验，结果表明，像二产比重、三产比重等对城市化水平影响不显著。

利用SPSS进行logistic回归分析(二元、多项)

线性回归是很重要的一种回归方法，但是线性回归只适用于因变量为连续型变量的情况，那如果因变量为分类变量呢？比方说我们想预测某个病人会不会痊愈，顾客会不会购买产品，等等，这时候我们就要用到logistic回归分析了。

Logistic回归主要分为三类，一种是因变量为二分类得logistic回归，这种回归叫做二项logistic回归，一种是因变量为无序多分类得logistic回归，比如倾向于选择哪种产品，这种回归叫做多项logistic回归。

还有一种是因变量为有序多分类的logistic回归，比如病重的程度是高，中，低呀等等，这种回归也叫累积logistic回归，或者序次logistic回归。

二值logistic回归：选择分析——回归——二元logistic，打开主面板，因变量勾选你的二分类变量，这个没有什么疑问，然后看下边写着一个协变量。

有没有很奇怪什么叫做协变量？在二元logistic回归里边可以认为协变量类似于自变量，或者就是自变量。

把你的自变量选到协变量的框框里边。

细心的朋友会发现，在指向协变量的那个箭头下边，还有一个小小的按钮，标着a*b，这个按钮的作用是用来选择交互项的。

我们知道，有时候两个变量合在一起会产生新的效应，比如年龄和结婚次数综合在一起，会对健康程度有一个新的影响，这时候，我们就认为两者有交互效应。

那么我们为了模型的准确，就把这个交互效应也选到模型里去。

我们在右边的那个框框里选择变量a，按住ctrl，在选择变量b，那么我们就同时选住这两个变量了，然后点那个a*b的按钮，这样，一个新的名字很长的变量就出现在协变量的框框里了，就是我们的交互作用的变量。

然后在下边有一个方法的下拉菜单。

默认的是进入，就是强迫所有选择的变量都进入到模型里边。

除去进入法以外，还有三种向前法，三种向后法。

一般默认进入就可以了，如果做出来的模型有变量的p值不合格，就用其他方法在做。

再下边的选择变量则是用来选择你的个案的。

用SPSS做logistic回归分析解读

如何用SPSS做logistic回归分析解读————————————————————————————————作者：————————————————————————————————日期：如何用进行二元和多元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<），因此我们这里选择以性别和年龄为例进行分析。

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

Spss软件之logistic回归分析

logit(P)=0 1X1 2 X 2 m X m
Logistic regression analysis
（二）模型参数的意义如果把logistic模型中的 P看作是在某一暴露状态下发
病的概率，则 β0：表示所有暴露剂量为0时发病与不发病概率之比的自然对数，反映了疾病的基准状态。 βj ：表示当因素 Xj 改变一个单位时logit（P）的改变量。
G 2（ln L1 ln L0）
当样本含量较大时，在零假设下得到的G统计量
近似服从自由度为d（d=p-l）的
2
分布。
由例13-1可以算得
lnL(X1 ) 585.326
•对于 H0：β1＝0和 H0：β2＝0
lnL(X1 , X2 ) 579.711
Hypothesis test
lnL(X2 ) 597.436
G1 2[lnL(X1 , X2 ) lnL(X2 )]=35.45>3.84 G2 2[lnL(X1 , X2 ) lnL(X1 )]=11.23>3.84
Hypothesis test
上面计算结果说明：在α=0.05检验水准上拒绝H0，接受H1，说明平衡了饮酒因素的影响后，食管癌与吸烟有显著性关系；同理，平衡了吸烟因素的影响后，食管癌与饮酒有显著性关系。
Hypothesis test
2.Wald检验
z bj , Sbj
2
bj Sbj
2
对于大样本资料，在零假设下z 近似
服从标准正态分布，而则近似服从
自由度=1的分布。
2
2
Abraham Wald
Hypothesis test
似然比检验可以对自变量增减时所得到的不同回归模型进行比较，既适合单个自变量的假设检验，又适合多个自变量的同时检验。Wald检验比较适合单个自变量的检验，但结果略为保守。

利用SPSS进行Logistic回归分析

Classification Tablea,b
Predicted
Observed
Step 0 城市化
Yes
No
Overall Percentage
a. Constant is included in the model.
b. The cut value is .500
城市化
Yes 0
No 11
0
20
Sig. .111
Exp(B) 1.818
6. Variable not in the Equation（不在初始方程中的变量）。人均 GDP 和代表地理位置的中部地带的系数初始值设为 0，这相当于，在初始模型中不考虑这两个变量（图 8-1-14）。表中给出了 Score 检验值及其对应的自由度 df 和 P 值，即 Sig.值。Score 检验是一种初始检验，在建模之初根据变量之间的结构关系判断自变量与因变量之间的密切程度。Score 检验值的计算公式为
图 8-1-1 原始数据（Excel 中，局部）将数据拷贝或者导入 SPSS 的数据窗口（Data View）中（图 8-1-2）。
1
研究生地理数学方法（实习）
Part 2 统计分析软件 SPSS
图 8-1-2 中国 31 个地区的数据（SPSS 中，局部）第二步：打开“聚类分析”对话框。
沿着主菜单的“Analyze→Regression→Binary LogisticK ”的路径（图 8-1-3）打开二值
研究生地理数学方法（实习）
Part 2 统计分析软件 SPSS
第 8 章利用 SPSS 进行 Logistic 回归分析
现实中的很多现象可以划分为两种可能，或者归结为两种状态，这两种状态分别用 0 和 1 表示。如果我们采用多个因素对 0－1 表示的某种现象进行因果关系解释，就可能应用到 logistic 回归。Logistic 回归分为二值 logistic 回归和多值 logistic 回归两类。首先用实例讲述二值 logistic 回归，然后进一步说明多值 logistic 回归。在阅读这部分内容之前，最好先看看有关 SPSS 软件操作技术的教科书。

手把手教你SPSS二分类Logistic回归分析

手把手教你SPSS二分类Logistic回归分析本教程手把手教您用SPSS做Logistic回归分析，目录如下：一、数据格式二、对数据的分析理解三、SPSS做Logistic回归分析操作步骤3.1 线性关系检验假设3.2 多重共线检验假设3.3 离群值、杠杆点和强影响点的识别3.4 Logistic回归分析四、SPSS计算结果的解释五、结果结论的撰写一、数据格式某研究者想了解年龄、性别、BMI和总胆固醇（TC）预测患心脏病（CVD）的能力，招募了100例研究对象，记录了年龄（age）、性别（gender）、BMI，测量血中总胆固醇水平（TC），并评估研究对象目前是否患有心脏病（CVD）。

部分数据如图1。

二、对问题分析使用Logistic模型前，需判断是否满足以下7项假设。

假设1：因变量（结局）是二分类变量。

假设2：有至少1个自变量，自变量可以是连续变量，也可以是分类变量。

假设3：每条观测间相互独立。

分类变量（包括因变量和自变量）的分类必须全面且每一个分类间互斥。

假设4：最小样本量要求为自变量数目的15倍，但一些研究者认为样本量应达到自变量数目的50倍。

假设5：连续的自变量与因变量的logit转换值之间存在线性关系。

假设6：自变量之间无多重共线性。

假设7：没有明显的离群点、杠杆点和强影响点。

假设1-4取决于研究设计和数据类型，本研究数据满足假设1-4。

那么应该如何检验假设5-7，并进行Logistic回归呢？三、SPSS操作3.1 检验假设5：连续的自变量与因变量的logit转换值之间存在线性关系。

连续的自变量与因变量的logit转换值之间是否存在线性关系，可以通过多种方法检验。

这里主要介绍Box-Tidwell方法，即将连续自变量与其自然对数值的交互项纳入回归方程。

本研究中，连续的自变量包括age、BMI、TC。

使用Box-Tidwell 方法时，需要先计算age、BMI、TC的自然对数值，并命名为ln_age、ln_BMI、ln_TC。

应用SPSS软件进行多分类Logistic回归分析

应用SPSS软件进行多分类Logistic回归分析应用SPSS软件进行多分类Logistic回归分析一、简介Logistic回归是一种常用的统计分析方法，在很多领域中都有广泛的应用。

它主要用于预测一个分类变量的可能性或概率，例如判断一个疾病的患病风险、判断学生成绩的优劣、预测金融市场的涨跌等。

本文将介绍如何使用SPSS软件进行多分类Logistic回归分析，并以一个具体案例来说明其应用。

二、SPSS软件介绍SPSS软件是统计分析的常用工具之一，它具有友好的用户界面和丰富的分析功能。

在进行Logistic回归分析时，SPSS可以帮助我们进行数据处理、模型建立、模型拟合、模型评估等步骤，并输出详细的分析结果。

三、案例描述我们假设有一份数据集，包含了500个样本和5个自变量，要根据这些自变量对样本进行多分类。

自变量包括性别、年龄、教育水平、收入和职业。

而多分类的目标变量是购买冰淇淋的偏好，包括三个分类：喜欢巧克力口味、喜欢草莓口味和喜欢香草口味。

四、数据处理首先，我们需要对数据进行处理。

SPSS可以读取各种文件格式，如Excel、CSV等。

我们将数据导入SPSS后，可以进行缺失值处理、异常值处理等预处理步骤。

这些步骤是为了保证后续的分析结果的准确性和可靠性。

五、模型建立在SPSS中，我们可以使用多分类Logistic回归模型进行建模。

它采用最大似然估计方法来估计模型参数，以便进行分类预测。

我们需要将自变量和目标变量进行指定，SPSS会自动计算出各个自变量对目标变量的系数和统计学意义。

六、模型拟合在模型拟合阶段，SPSS会对模型进行拟合优度的检验，包括卡方拟合优度检验、Hosmer-Lemeshow检验等。

这些检验可以帮助我们评估模型的拟合程度和可靠性。

如果模型的拟合程度不好，我们可以对模型进行进一步调整和改进。

七、模型评估在模型评估阶段，SPSS提供了一系列的统计指标和图表，用于评估多分类Logistic回归模型的性能。

如何用SPSS做logistic回归分析

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

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

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

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

图1-1第二步：打开“二值Logi‎s tic 回归分析”对话框：沿着主菜单的‎“分析（Analyz‎e）→回归（Regres‎s ion）→二元logi‎s tic （Binary‎Logist‎i c）”的路径（图1-2）打开二值Lo‎g istic‎回归分析选项‎框（图1-3）。

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

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

SPSS-配对条件 Logistic 回归分析

SPSS配对调查资料的条件 Logistic 回归分析（1:1或1:n）1. 1:1 病例对照研究的基本概念在管理工作中，我们也经常要开展对照调查。

例如为什么有的人患了胃癌，有的人却不会患胃癌？如果在同一居住地选取同性别、年龄相差仅±2 岁的健康人作对照调查，调查他们与患胃癌有关的各种影响因素，这就是医学上很常用的所谓“1:1 病例对照研究”。

病例对照研究资料常用条件Logistic 回归分析。

条件Logistic 回归模型（conditional logistic regression model,CLRM)，下称CLRM 模型。

2. 条件Logistic 回归模型的一个实例某地在肿瘤防治健康教育、社区干预工作中做了一项调查，内容是三种生活因素与胃癌发病的关系。

调查的三种生活因素取值见表 11-6。

请拟合条件Logistic 回归模型，说明胃癌发病的主要危险因素。

表 11-6 三种生活因素与胃癌发病关系的取值------------------------------------------------------------------------------------------ 变量名取值范围------------------------------------------------------------------------------------------ X1 (不良生活习惯) 0，1，2，3，4 表示程度（0 表示无，4 表示很多）X2 (喜吃卤食和盐腌食物) 0，1，2，3，4 表示程度（0 表示不吃，4 表示喜欢吃、吃很多） X3 (精神状况) 0 表示差，1 表示好------------------------------------------------------------------------------------------表 11-7 50 对胃癌病例（S=1）与对照（S=0）三种生活习惯调查结果------------------------------------------------------------------------------------------ 病例对照病例对照-----------------------------------------------------------------------------No S X1 X2 X3 No S X1 X2 X3 No S X1 X2 X3 No S X1 X2 X3------------------------------------------------------------------------------------------1 12 4 0 1 03 1 0 26 1 2 2 0 26 0 1 1 02 13 2 1 2 0 0 1 0 27 1 2 0 1 27 0 0 2 13 1 3 0 0 3 0 2 0 1 28 1 1 1 1 28 0 3 0 14 1 3 0 0 4 0 2 0 1 29 1 2 0 1 29 0 4 0 05 1 3 0 1 5 0 0 0 0 30 1 3 1 0 30 0 0 2 16 1 2 2 0 6 0 0 1 0 31 1 1 0 1 31 0 0 0 07 1 3 1 0 7 0 2 1 0 32 1 4 2 1 32 0 1 0 18 1 3 0 0 8 0 2 0 0 33 1 4 0 1 33 0 2 0 19 1 2 2 0 9 0 1 0 1 34 1 2 0 1 34 0 0 0 110 1 1 0 0 10 0 2 0 0 35 1 1 2 0 35 0 2 0 111 1 3 0 0 11 0 0 1 1 36 1 2 0 0 36 0 2 0 112 1 3 4 0 12 0 3 2 0 37 1 0 1 1 37 0 1 1 013 1 1 1 1 13 0 2 0 0 38 1 0 0 1 38 0 4 0 014 1 2 2 1 14 0 0 2 1 39 1 3 0 1 39 0 0 1 015 1 2 3 0 15 0 2 0 0 40 1 2 0 1 40 0 3 0 116 1 2 4 1 16 0 0 0 1 41 1 2 0 0 41 0 1 0 117 1 1 1 0 17 0 0 1 1 42 1 3 0 1 42 0 0 0 118 1 1 3 1 18 0 0 0 1 43 1 2 1 1 43 0 0 0 019 1 3 4 1 19 0 2 0 0 44 1 2 0 1 44 0 1 0 020 1 0 2 0 20 0 0 0 0 45 1 1 1 1 45 0 0 0 121 1 3 2 1 21 0 3 1 0 46 1 0 1 1 46 0 0 0 022 1 1 0 0 22 0 2 0 1 47 1 2 1 0 47 0 0 0 023 1 3 0 0 23 0 2 2 0 48 1 2 0 1 48 0 1 1 024 1 1 1 1 24 0 0 1 1 49 1 1 2 1 49 0 0 0 125 1 1 2 0 25 0 2 0 0 50 1 2 0 1 50 0 0 3 1------------------------------------------------------------------------------------------- 3. 条件Logistic 回归模型的拟合原理与方法本例以 SPSS 软件包来拟合 CLRM 模型。

无序多分类Logistic回归-SPSS教程

无序多分类Logistic回归-SPSS教程一、问题与数据为了探讨基因X突变与恶性肿瘤Y不同组织类型发生风险的关系，某医生设计了一项病例对照研究。

该医生纳入所在科室一年收治的145名该恶性肿瘤患者，并从医院体检数据库中随机选择了100名未患该肿瘤的体检者作为对照。

相关信息整理成表1：表1 各病例组织类型与突变情况变量赋值情况如表2：表2 变量及变量赋值情况二、对问题分析该研究中，“病例”与“对照”的关系不再是简单的“患病”与“不患病”，而是病例分为四类（本例中包含对照组共四类），且各类别无次序关系。

或者说，因变量Y不再是二分类的，而是无序多分类的。

通过无序多分类的Logistic回归分析可以将三种不同组织类型的病例分别与对照组进行对比，分别得到基因X 突变与三种肿瘤组织类型的暴露-风险关系。

三、SPSS操作A. 数据录入SPSS若数据格式如表1所示，则首先在SPSS变量视图（Variable View）中新建三个变量：ID代表患者编号，Y代表组织类型，X代表是否突变，赋值参考表2.然后在数据视图（Data View）中录入数据。

B. 选择Analyze →Regression →Multinomial LogisticC. 选项设置将变量Y选入因变量（Dependent）位置，变量X选入因子（Factors）位置。

如果自变量中还有连续型变量，则需要放入协变量（Covariate）位置。

由于因变量Y有多个分类，而无序多分类Logistic回归的原理是先指定一个类别为参考类别，然后将其他类别分别与参考类别对比。

故需点击Reference Category 设置参考类别（本例中作为参考类别的为对照组）。

SPSS默认选择因变量赋值中按升序排列后最后类别（即赋值最大者）为参考类别（即对照组），而本研究中参考类别Y赋值为0，故可以点击First Category 或直接在Custom中输入0，点击Continue。

SPSS专题2回归分析线性回归Logistic回归对数线性模型

• 建立这种关系的过程就叫做回归。
2
回归分析 • 一旦建立了回归模型 • 可以对各种变量的关系有了进一步的定量理解 • 还可以利用该模型（函数）通过自变量对因变量做
预测。 • 这里所说的预测，是用已知的自变量的值通过模型
对未知的因变量值进行估计；它并不一定涉及时间先后的概念。
3
例1 有50个从初中升到高中的学生.为了比较初三的成绩是否和高中的成绩相关,得到了他们在初三和高一的各科平均成绩(数据:highschool.sav)
Si g. .0 00 a
Coeffi c ie ntsa
Unstandardi zed Coeffi ci ents
M odel
1
(Con st a nt )
B
Std. Error
2 6. 4 4 4
5 . 39 6
j3
. 6 51
. 0 72
a. Dependent Vari abl e: s1
M odel Summary
ANOVAc
M odel 1
2
R .9 30 a
.9 38 b
R Square . 8 66 . 8 79
Adj usted R Square
. 8 63
. 8 75
Std. Error of the Esti mate
1 2. 4 8 44 1
1 1. 9 5 60 2
1 . 00 0 .
N
1 64
1 92
Spearman's rho cl eanwateraccess_ rural (%)
Correl ati on Coeffi ci ent
1 . 00 0
.6 76 **

SPSS的Logistics回归

SPSS的Logistics回归实验⽬的学会使⽤SPSS的简单操作，Logistic回归。

实验要求使⽤SPSS。

实验内容实验步骤（1）⼆项分类Logistic回归SPSS分析，使⽤Hosmer和Lemeshow于1989年研究低出⽣体重婴⼉的影响因素作为演⽰例⼦。

结果变量为“是否娩出低出⽣体重⼉”，考虑影响因素有8个，详见Logistics_step.sav⽂件。

本例题主要演⽰“⾃变量的筛选与逐步回归”。

操作如下：点击【分析】→【回归】→【⼆元Logistics回归】，在打开的对话框中，把待结果变量LOW选⼊【因变量】中，将变量LWT,AGE,SMOKE,PTL,HT,UI,FTV,RACE选⼊【协变量】中。

点击【分类】，把RACE选⼊【分类协变量】→【第⼀个】→【变化量】→【继续】，【块】⾥的【⽅法（M）】选【向前：LP】，【选项】→【Exp（B）的置信区间】→【继续】，单击【运⾏】。

主要分析结果如下：分类变量编码频率参数编码(1)(2)种族⽩⼈96.000.000⿊⼈26 1.000.000其他种族67.000 1.000 上表输出race在产⽣哑变量时的编码情况，以⽩⼈为参照⽔平。

未包括在⽅程中的变量得分⾃由度显著性步骤 0变量产妇体重 4.6161.032产妇年龄 2.4071.121产妇在妊娠期间是否吸烟 4.9241.026本次妊娠前早产次数7.2671.007是否患有⾼⾎压 4.3881.036应激性 4.2051.040随访次数.9341.334种族 5.0052.082种族(1) 1.7271.189种族(2) 1.7971.180总体统计29.1409.001 输出的是拟合包含常数项和任⼀⾃变量的Logistics回归模型检验统计量、⾃由度及P值。

其中race产⽣两个哑变量，因此其总⾃由度为2。

由上表可以发现，本次妊娠前早产次数（ptl）的score统计量最⼤，P=0.007，⼩于SPSS默认选⼊变量的标准（0.05）因此下⼀步将它⾸先选⼊模型。