逻辑斯蒂回归（LR）原理详解及公式推导

说白了，逻辑斯蒂回归实际上就是label y = 1的对数几率是输入x的线性函数的模型。

• 对数几率
  一个事件发生的概率为P$P$, 则该事件发生的几率为P1−P$\frac{P}{1-P}$, 即为该事件发生的概率除以该事件不发生的概率。
  该事件的对数几率即为：

  logit(p)=logp1−p$$logit(p)=log\frac{p}{1-p}$$

• 逻辑斯蒂回归模型
  事件y=1$y=1$ 的对数几率为输入x$x$ 的线性函数，即

  logP(y=1|x)1−P(y=1|x)=ωx+b　　　　　　(１)$$log\frac{P(y=1|x)}{1-P(y=1|x)}=\omega x+b (１)$$
  
  上述模型即为逻辑斯蒂回归模型，因此逻辑斯蒂回归是一种对数线性模型（PS：有可能数据集实际上是线性不可分的，这时候若要用逻辑斯蒂回归，可引入核函数，利用核技巧来实现（PPS：核技巧会在后续SVM算法中提到），利用核函数可以将线性不可分数据集映射到高维空间，在高维空间中，数据集变为线性可分的了）。通过定义合理的损失函数（下文会分析推导），学习得到ω,b$\omega ,b$ 的最优值。对于一个新的未知样本x$x$， 代入即可得到P(y=1|x)$P(y=1|x)$ 的值，若大于0.5，即可认为y=1$y=1$ 的概率大于y=0$y=0$ 的概率，判定y=1$y=1$是合理的。

• 逻辑斯蒂回归损失函数与最大似然估计
  首先说结论，最小化逻辑斯蒂回归的损失函数与最大化对数似然函数是等价。下面开始推导：
  由逻辑斯蒂回归模型(1)容易得到

  π(x)=P(y=1|x)=eωx+beωx+b+1$$\pi (x)=P(y=1|x)=\frac{e^{\omega x+b}}{e^{\omega x+b}+1}$$
  
  则P(y=0|x)=1−π(x)$P(y=0|x)=1-\pi (x)$ 于是样本(x,y)$(x,y)$ 的似然值为L(x,y)=π(x)y×(1−π(x))(1−y)$L(x,y)=\pi (x{)}^y\times (1-\pi (x){)}^{(1-y)}$ （y=1时L(x,y)=π(x);当y=0时，L(x,y)=1−π(x)$y=1时L(x,y)=\pi (x);当y=0时，L(x,y)=1-\pi (x)$ 综合两种情况即可得出结论）于是参数ω,b$\omega ,b$的似然函数为：
  L(ω,b)=∏xi,yi{π(xi)yi×(1−π(xi))1−yi}$$L(\omega ,b)=\prod _{x_i,y_i}\{\pi {(x_i)}^{y_i}\times (1-\pi (x_i){)}^{1-y_i}\}$$
  
  对上式取对数将连乘变为累加：
  L(ω,b)=∑i=1N{yilog(π(xi))+(1−yi)log(1−π(xi))}$$L(\omega ,b)=\sum _{i=1}^N\{y_{i}log(\pi (x_i))+(1-y_i)log(1-\pi (x_i))\}$$
  
  而逻辑斯蒂回归的损失函数定义为：
  J(ω,b)=−1N∑i=1N{yilog(π(xi))+(1−yi)log(1−π(xi))}$$J(\omega ,b)=-\frac{1}{N}\sum _{i=1}^N\{y_{i}log(\pi (x_i))+(1-y_i)log(1-\pi (x_i))\}$$
  
  对比以上两式可知最小化损失函数与极大化似然函数是等价的。

• 逻辑斯蒂回归优化方法
  在对数据进行训练时，即是对损失函数最小化的过程，在训练过程中可以利用随机梯度下降法进行训练更新权值。即可以没进来一个样本或者一批（batch）样本后计算梯度，对权值以及偏置进行更新。更新公式如下
  ωt=ωt−η∂L∂ω$${\omega }_t={\omega }_t-\eta \frac{\mathrm{\partial }L}{\mathrm{\partial }\omega }$$
  
  令zi=ωxi+b$z_i=\omega x_i+b$ ，对于单个样本的代价函数为：
  L(xi,yi)=−(yilog(π(xi))+(1−yi)log(1−π(xi)))$$L(x_i,y_i)=-(y_{i}log(\pi (x_i))+(1-y_i)log(1-\pi (x_i)))$$
  利用链式求导法则：
  ∂L∂ω=∂L∂π(xi)×∂π(xi)∂zi×∂zi∂ω$$\frac{\mathrm{\partial }L}{\mathrm{\partial }\omega }=\frac{\mathrm{\partial }L}{\mathrm{\partial }\pi (x_i)}\times \frac{\mathrm{\partial }\pi (x_i)}{\mathrm{\partial }{z}_i}\times \frac{\mathrm{\partial }{z}_i}{\mathrm{\partial }\omega }$$
  
  具体梯度更新公式推导敲公式太麻烦了，下面是我的推导过程。
  
  权值向量ω$\omega$ 和偏置b的梯度即更新公式如上图所示，如果每次更新时都用到了所有的样本，则这种优化算法为普通的梯度下降法GD，若每次更新权值都只利用部分样本甚至是单个样本则为随机梯度下降法SGD。