逻辑回归
一、什么是逻辑回归
逻辑回归的原理是用逻辑函数把线性回归的结果(-∞,∞)映射到(0,1)
线性回归函数
(资料图片)
线性回归函数的数学表达式:
其中xi是自变量,y是因变量,y的值域为(-∞,∞),θ0是常数项,θi(i=1,2,...,n)是待求系数,不同的权重θi反映了自变量对因变量不同的贡献程度。
我们初中学过的一元一次方程:y=a+bx,这种只包括一个自变量和一个因变量的回归分析称为一元线性回归分析。
初中学过的二元一次方程:y = a+b1x1+b2x2,三元一次方程:y = a+b1x1+b2x2+b3x3,这种回归分析中包括两个或两个以上自变量的回归分析,称为多元线性回归分析。
不管是一元线性回归分析还是多元线性回归分析,都是线性回归分析。
逻辑函数(Sigmoid函数)
二、如何求解逻辑回归中的参数
极大似然函数
先来看一个小例子:如果小华这次考试考了90分以上,妈妈99%会奖励小华一个手机,如果没有考到90分以上,妈妈99%不会奖励小华手机。现在小华没有得到手机,问小华这次有没有考到90分。
可能我们的第一反应是小华大概率没有考到90分以上。这种利用已知样本结果,反推最有可能导致这样结果的参数值,就是极大似然估计。
结合逻辑回归函数,如果我们已经积累了大量的违约客户和正常客户的样本数据,利用极大似然函数由果溯因,估计出使得目前结果的可能性最大参数(系数)θ,有了参数我们就可以求任何一个客户违约的概率了。
我们上文提到过客户违约的后验概率
相应的可以得到客户不违约的概率:
如果令
违约的后验概率可以写成:
不违约的后验概率可以写成:
对于某一个客户,我们采集到了样本数据(x,y)。对于这个样本,他的标签是y的概率可以定义成:
其中y∈{0,1}。当y=0时,上式为不违约的后验概率,当y=1时,上式为违约的后验概率。
现在我们有m个客户的观测样本
将每一个样本发生的概率相乘,就是这个合成在一起得到的合事件发生的总概率(利用概率中的乘法公式),即为似然函数,可以写成:
其中θ为待求参数。
注:我们总是希望出现目前结果的可能性最大,所以想要得到极大化似然函数对应的参数θ。
为便于求解,我们引入不改变函数单调性的对数函数ln,把连乘变成加法,得到对数似然函数:
至此,可以用梯度上升法求解对数似然函数,求出使得目前结果的可能性最大的参数θ。也可以由对数似然函数构造损失函数,用梯度下降法求出使得损失最小对应的参数θ,接下来看下逻辑回归中的损失函数。
注:使用对数似然函数,不仅仅把连乘变成加法,便于求解,而且对数似然函对应的损失函数是关于未知参数的高阶连续可导的凸函数,便于求其全局最优解。
构造损失函数
在机器学习中有损失函数的概念,我们知道损失函数一般定义为预测值和真实值的差,比如我们预测小华在这次考试中能考98分,成绩出来了小华实际考了97分,小华的成绩预测值和真实值差为1,这个1通俗理解就是损失函数的值。
从上面的案例知,如果损失函数越小,说明模型预测越准。所以在函数比较复杂没有确定解(解析解)或很难求出确定解的情况下,一般求的是数值解(近似解)。一般模型求数值解可以求出使得损失函数最小对应的参数θ。
结合逻辑回归中的极大似然函数,如果取整个数据集上的平均对数似然损失,我们可以得到:
其中J(θ)为损失函数,由对数似然函数前面添加负号取平均得到。
即在逻辑回归模型中,最大化似然函数和最小化损失函数实际上是等价的(求最大化对数似然函数对应的参数θ和求最小化平均对数似然损失对应的参数θ是一致的),即:
那如何求得损失函数最小对应的参数呢?可以用下节讲到的方法:梯度下降法。
用梯度下降法求解参数
先以一个人下山为例讲解梯度下降法的步骤:
step1:明确自己现在所处的位置;step2:找到现在所处位置下降最快的方向;step3: 沿着第二步找到的方向走一个步长,到达新的位置,且新位置低于刚才的位置;step4:判断是否下山,如果还没有到最低点继续步骤一,如果已经到最低点,则停止。从上面的分析知,用梯度下降法求解参数最重要的是找到下降最快的方向和确定要走的步长。
那么什么是函数下降最快的方向?
如果学过一元函数的导数,应该知道导数的几何意义是某点切线的斜率。除此之外导数还可以表示函数在该点的变化率,导数越大,表示函数在该点的变化越大。
可以发现p2点的斜率大于p1点的斜率,即p2点的导数大于p1点的导数。
对于多维向量
它的导数叫做梯度(偏导数),当求某个变量的导数时,把其它变量视为常量,对整个函数求导,也就是分别对于它的每个分量求导数,即
对于函数的某个特定点,它的梯度就表示从该点出发,函数值变化最为迅猛的方向。至此梯度下降法求解参数的方向已经找到,那就是函数的梯度方向。
接下来推导损失函数的梯度(偏导数):
由损失函数的公式知:
对损失函数求偏导:
损失函数 J(θ)中,是 θ和 x 都是一个向量形式,也就是:
当我们仅仅针对θ进行求导的时候:
可以看到上面需要利用sigmoid函数的求导方法:
至此,LR的损失函数的求导形式已经结束,这里面主要用到了sigmoid函数的求导,推导起来其实是比较简单的。最后写出梯度下降的更新公式:
关键词:
推荐阅读
日全食过程 日全食对人类和地球的影响有哪些?
日全食过程可以包括以下五个阶段:初亏、食既、食甚、生光、复原1 初亏因为月亮从西向东绕着地球转,所以日总是从太阳圆形表面的西边开始。 【详细】
海兔的名字怎么来的?海兔是鱿鱼吗?
海兔不是长在陆地上的兔子,而海兔是海里的软体动物,可以在海兔吃。沿海的渔民一般在海兔,晒干,可以在海兔直接晒干,或者在海兔蒸熟。海 【详细】
苹果CEO去年薪酬超1400万美元 苹果高级管理人员的工资是多少?
苹果CEO去年薪酬超1400万美元根据美国科技博客MacBrothers的报道,苹果公司提交给美国证券交易委员会(SEC)的2020年委托书显示,苹果公司CEO 【详细】
毅力号录下了来自火星的声音 科学家如何采集宇宙中的声音?
毅力号录下了来自火星的声音当地时间2月22日,NASA(美国国家航空航天局)发布了毅力号火星车录下的首段来自火星的音频。此外,NASA还公布了 【详细】
将禁用不可降解塑料袋 可以使用那些塑料袋替代品呢?
将禁用不可降解塑料袋国家发改委、生态环境部等九部门联合印发《关于扎实推进塑料污染治理工作的通知》,提出自明年1月1日起,在直辖市、省 【详细】
相关新闻
- 逻辑回归
- 全球快看:别傻了,增程式电动汽车是省油不省钱!
- 世界观天下!碳纳米管薄膜电学输运性能与其手性结构的依存关系 - 进展
- 全球消息!四水是哪四水_中国 ldquo 四水 rdquo 是那些
- 三星W2019经典升级:翻盖双摄真机或将上线!-天天实时
- 华硕发布ATX规格B365主板:ROG Strix B365-F Gaming
- 索尼发布A7M3 升级换代更新快-速看料
- 小白智能摄像头电池版:充电一次,待机100天,这续航能力,跪了~_全球微资讯
- 热议:影驰名人堂hof pro m.2 SSD 让玩家体验“快”人一步
- 新能源还是建议买插混,15-20万就看这两款,销量高,口碑也好 每日快播
- 一招教你识别新能源车牌不同字母的含义-全球热头条
- 核显性能对标RTX 4070!AMD Zen5 APU全线曝光 环球关注
- 【全球独家】特斯拉的“充电革命”?10分钟续航250km,或将淘汰油车和换电站
- MSI正式开战,BLG成最大BOSS,Bin连放三句狠话,但状态令人担心
- 租户装修私拆承重墙致200户居民被疏散,律师:当事人涉嫌危害公共安全罪,施工队也难逃责任
- Mac Mini 全身大改造,速来围观!
- 全球快资讯:雷神笔记本具有哪几个功能较好的系列
- 尼康D6:全旗舰单反,你值得拥有
- 丰田首席科学家:单一纯电模式是错误的
- 江西佳时特数控技术有限公司深耕机床领域-“数字大脑”指挥生产制造-每日头条