统计挖掘那些事——强大的逻辑回归(理论+案例)

浏览: 4249

幻灯片1.jpg

幻灯片2.jpg

在上期,浩彬老撕给大家介绍了非线性回归模型,解决了在现实环境中,非线性形式的问题。但是进一步地,我们的因变量也并不总是数值型变量,有可能也是分类型变量,那么对于这种问题,我们能不能也利用回归分析进行适当的扩展,使其也能够解决分类问题?

答案显然也是肯定的。

附回归文章整理:

1.统计挖掘那些事那些情-回归分析

2.统计挖掘那些事-超详尽回归分析指南(理论+动手案例)

3.统计挖掘那些事-9个相关R先生的故事(理论+动手案例)

4.(理论+案例)如何通俗地理解极大似然估计?

5.统计挖掘那些事——婀娜多姿的非线性回归(理论+案例)

1.Logit回归

本期将会大家介绍逻辑回归,虽然逻辑回归并不复杂,但正是由于其简单,高效,可解释性强的特点,在实际用途中十分的广泛,从购物预测到用户营销响应,从流失分析到信用评价,都能看到其活跃的身影,可以说,逻辑回归占据了分类算法中非常重要的地位。

回想在上一期中,我们谈到当因变量与自变量的关系式不再是线性时,通过引入衍生变量y’,使其转换为线性表达形式。那么很自然地,对于我们现在面临的任务,我们就需要一个转换,使得分类变量0和1转化为可用的形式

先考虑一个二分类的预测变量,正如前面所说的,显然由于分类数据的特点,已经不适合运用传统的线性函数进行分析。但是二分类事件的y的期望值E(y)来说,它等价于事件发生概率,从y到E(y),我们就把事件发生与否与值域在[0,1]区间的事件发生概率相联系,这提示我们可以用事件发生的概率进行代替。

既然使用发生概率代替的话,一个自然而然的选择是把回归函数的值域限制在[0,1]区间内,这样当f(xi)接近负无穷时,将有E( yi)趋近于0,而在f(xi)接近正无穷时,将有E(yi )趋近于1,这样看来,显然相比于研究二元变量y与x的关系,研究y发生的条件概率与x更具适应性。

幻灯片3.jpg

在没有任何先验条件的情况下,这里的阈值一般选择0.5。但当我们有进一步明确需求的时候,阈值也是可以调整的,例如我们希望对正例样本有更高的准确率要求,则可以把阈值适当地调高,例如调高到0.6;相反,假如我们希望对正例样本的召回率要求更高,则可以把阈值适当地降低,例如降低到0.4;

一般地,我们选择Logit函数作为转换函数,Logit函数的形式:



幻灯片4.jpg

Logit函数图像是一个典型的S型的曲线,并且它的值域是在[0,1]之间

幻灯片5.jpg

进一步地,我们利用logit函数,可以把事件发生的条件概率与x表示为

幻灯片6.jpg

同样,我们也可以定义一个事件不发生的概率为:

幻灯片7.jpg

为了更显简洁,不妨作如下转换:

幻灯片8.jpg

上式左边实际上就是表示“事件发生的概率”与“事件不发生的概率”之比,称之为事件的发生比,简称odds

幻灯片9.jpg

实际上,转化为线性函数形式后,我们可以看到逻辑回归实际上就是以线性回归的形式去逼近时间发生优势比的对数,因此也有一些文献称之为“对数几率回归”。

2.参数估计

考虑到在Logit回归的推导中,咱们已经把事件发生的概率公式给出,那么我们就可以借助极大似然估计进行参数估计的工作,设

幻灯片10.jpg

求得对数似然函数并整理:

幻灯片11.jpg

3.实战案例:

样例数据:个人收入水平调查分析.xlsx

链接: http://pan.baidu.com/s/1cmmP1W

密码: raxx

该数据集是某地区的个人收入调查分析,包含32561条记录,其中目标变量是收入水平(分别是<=50k以及>50k),其他自变量包括年龄,受教育时间,性别,资产净增,资产损失,一周工作时间。

幻灯片12.jpg

模型流如下所示:

幻灯片13.jpg

利用Excel源节点读取数据后,接入类型节点,在类型节点中:

(1)把收入水平的测量设为‘标记’,把角色设为‘目标’

(2)把年龄,受教育时间,性别,资产净增,资产损失,一周工作时间设为‘输入’

幻灯片14.jpg

上述介绍的内容在Modeler中是在Logistic节点中实现,因此我们在下方建模选项板中,选中Logistic节点,并将其添加到流。

在回Logistic节点中,模型选项卡下,因为我们的目标变量‘收入水平’属于二分类变量,因此我们选择‘二项式’

在二项式过程中中,我们选择‘向前步进法’建立Logistic回归模型

选择好后,点击运行。

幻灯片15.jpg

运行模型后,点开模型块查看模型结果。

首先看到的是CaseProcessing Summary,我们知道我们一共使用了32561条记录构建模型,其中所有记录无缺失;由于自变量与因变量都含有分类变量,因此需要进行编码。其中因变量,我们把收入水平>50k设为1,收入水平<=50k设为0;另外在自变量部分,只有性别属于分类变量,我们看到其中女性有10771个记录,男性则有21790个记录,其中我们把女性设为‘1’,男性设为‘0’。

幻灯片16.jpg

之后直接看到模型结果,可以看到一共经历6补构建了最终模型,纳入了6个自变量,即我们所有的自变量都被纳入了方程,并且检查系数显著性检验结果,发觉所有显著性检验结果都小于0.05。额外地,看到性别(1)自变量,这是因为性别属于分类变量,我们将其设置为哑变量。即性别为女性的样本进入方程,将减去1.175*1,而性别为男性的样本进入方程则是默认该项取0;

幻灯片17.jpg

根据结果,我们可以写出最终的回归方程有:


幻灯片18.jpg

更重要地,对于分类问题,我们同样十分关心模型的预测准确率,为了进行比较,我们在模型节点后添加‘分析’节点(在输出选型卡下),其中勾选‘重合矩阵(用于字符型目标字段)’,点击运行。

幻灯片19.jpg

通过分析结果,我们可以看到逻辑回归分析的结果还是比较准确地,准确率有82.55%。

幻灯片20.jpg

幻灯片21.jpg


近期热门文章精选(点击标题即可阅读):

1.干货教程|可能是最方便好用的文字云工具

2.(理论+案例)如何通俗地理解极大似然估计?

3.从前,有一位统计学家,他在趟过一条平均水位不足1米深的河流时......

4.数据收集宝鉴|教你找到合适的数据源(1)——国家数据

5.重磅|李飞飞担任谷歌云机器学习最新负责人(附最全报道)

注:由于最近老撕比较忙,所以有些信息没有及时查看回复(微信规定48小时没有回复的话,则不能再次回复),所以如果浩彬老撕没有及时回复,大家可以再次留言。

幻灯片21.jpg

作者简介:浩彬老撕

好玩的IBM数据工程师,

立志做数据科学界的段子手,

致力知识分享,每月至少一次送书活动

推荐 3
本文由 浩彬老撕 创作,采用 知识共享署名-相同方式共享 3.0 中国大陆许可协议 进行许可。
转载、引用前需联系作者,并署名作者且注明文章出处。
本站文章版权归原作者及原出处所有 。内容为作者个人观点, 并不代表本站赞同其观点和对其真实性负责。本站是一个个人学习交流的平台,并不用于任何商业目的,如果有任何问题,请及时联系我们,我们将根据著作权人的要求,立即更正或者删除有关内容。本站拥有对此声明的最终解释权。

2 个评论

有个疑问:为什么逻辑回归的表达式中, 很少有交互项或者自变量的高次项?
我个人觉得是交互项放进去没有问题

要回复文章请先登录注册