1.1、Logistics Regression模型
1、线性可分VS线性不可分
对于一个分类问题，通常可以分为线性可分与线性不可分两种 。如果一个分类问题可以使用线性判别函数正确的分类，则称该问题为线性可分。如图所示为线性可分，否则为线性不可分：
下图为线性不可分：
1.2、Logistics Regression模型
Logistics Regression模型为广义的线性模型的一种，属于线性的分类模型。对于线性可分问题，需要找到一条直线，能够将两个不同的类分开，这条直线也称为超平面。对于上述超平面，可以使用如下的线性函数表示：
其中W为权重，b为偏置。若在多维的情况下，权重W和偏置b均为向量。在Logistic Regression算法中，通过对训练样本的学习，最终得到该超平面，将数据分成正负两个类别。此时可以使用阈值函数，将样本映射到不同的类别中，常见的阈值函数有Sigmoid函数，其形式如下：
Sigmoid函数的图像如图所示。
从Sigmoid函数的图像可以看出，其函数值域为（0,1），在0附近的变化比较明显。其导函数
为：
以下为Python实现Sigmoid函数，为了能够使用numpy中的函数，我们首先导入numpy
# -*- coding: UTF-8 -*-
# date:2018/5/24
# User:WangHong
import numpy as np
def sig(x):'''name:Sigmoid函数:param x::return:Sigmoid函数的值'''return 1.0 / (1 + np.exp(-x))
下列代码绘制自变量取值在[-7,7]的Sigmoid函数图像：
# -*- coding: UTF-8 -*-
# date:2018/5/24
# User:WangHong
import numpy as np
import matplotlib.pyplot as plt
def sig(x):'''name:Sigmoid函数:param x::return:Sigmoid函数的值'''return 1.0 / (1 + np.exp(-x))z = np.arange(-7,7,0.1)
phi_z = sig(z)
plt.plot(z,phi_z)
plt.axvline(0.0,color='k')
plt.axhspan(0.0,1.0,facecolor='1.0',alpha=1.0,ls='dotted')
plt.axhline(y=0.0,ls='dotted',color='k')
plt.axhline(y=0.5,ls='dotted',color='k')
plt.axhline(y=1.0,ls='dotted',color='k')
plt.yticks([0.0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.0])
plt.ylim(-0.1,1.1)
plt.xlabel('z')
plt.ylabel('$\phi (z)$')
plt.show() 
可以看出，sigmoid函数连续，光滑，严格单调，以(0,0.5)中心对称，是一个非常良好的阈值函数。
当x趋近负无穷时，y趋近于0；趋近于正无穷时，y趋近于1；x=0时，y=0.5。当然，在x超出[-6,6]的范围后，函数值基本上没有变化，值非常接近，在应用中一般不考虑。
Sigmoid函数的值域范围限制在(0,1)之间，我们知道[0,1]与概率值的范围是相对应的，这样sigmoid函数就能与一个概率分布联系起来了。
Sigmoid函数的导数是其本身的函数，即f ′ (x)=f(x)(1?f(x)) ，
计算非常方便，也非常节省计算时间。推导过程如下：
在程序中，Sigmoid函数的输出为Simoid值，对于输入向量X，其属于正例的概率为：
其中，
表示的是Sigmoid函数。那么，对于输入向量X，其属于负例的概率为：
对于Logistics Regression算法来说，需要求解的分隔超平面中的参数，即为权重矩阵W和偏置向量b，那么这些参数该如何求解呐？为了求解模型的两个参数，首先必须定义损失函数。
1.3、损失函数
对于上述的Logistics Regression算法，其属于类别y的概率为：
要求上述问题的参数W和b，可以使用极大似然法对其进行估计，假设训练数据集有m个训练样本{（X1，y1），（X2，y2）,....,(Xm,ym)}，则其似然函数为：
其中，假设函数hw,b(X(i))为：
对于似然函数的极大值的求解，通常使用Log似然函数，在Logistics Regression算法中，通常是将负的Log似然函数作为其损失函数，即the negative log-likelihoo（NLL）作为其损失函数，此时，需要计算的NLL的最小值。损失函数lW,b为：
此时，我们需要求解的问题为：
为了求解德损失函数lW,b的最小值，可以使用基于梯度的方法进行求解。
1.4、梯度下降法：
在机器学习算法中，对于很多监督学习模型，需要对原始的模型构建损失函数，接下来便是通过优化算法对损失函数进行优化，以便寻找到最优的参数W，字求解机器学习参数W的优化算法时，使用较多的是基于梯度下降的优化算法（Gradient Descent，GD）。
梯度下降法的优点，只需要求解损失函数的一阶导数，计算的成本比较小，这使得梯度下降法能在很多大规模数据集上得到应用。梯度下降法的含义是通过当前点出的梯度方向寻找新的迭代点，并从当前点移动到新的迭代点继续寻找新的迭代点，直到找到最优解。
1.4.1、梯度下降法的流程
梯度下降法是一种迭代性的优化算法，根据初始点在每一次迭代的过程中选择下降的方向，静儿改变需要修改的参数，对于优化问题min f（W），梯度下降法的详细过程如下所示：
-
-
- 随机选择一个初始点W0
- 重复以下过程
- 决定梯度下降方向：
- 选择步长α
- 更新：Wi+1=Wi+α*di
- 直到满足终止条件。
-
在初始时，在 W0处，选择下降的方向d0，选择步长α，更新W的值，此时到达W1处，判断是否满是终止条件，发现并未到达最优解W*，重复上述的过程，直至到达W*
1.5、凸优化与非凸优化
简单来讲，凸优化问题是指只存在一个最优解的优化问题，即在任何一个局部最优解即全局最优解，如图所示：
非凸优化是指在解空间存在多个局部最优解，而全局最优解是其中的某一个局部最优解，如图
最小二程（Least Sqares）、岭回归（Ridge Regression）和Logistics回归（Logistics Regression）的损失函数都是凸优化问题。
1.6、利用梯度下降法训练Logistics Regression模型
对于上述的Logistics Regression算法的损失函数可以通过梯度下降法对其进行求解，其梯度为：
其中xj（i）表示的是样本X（i）的第j个分量。取w0=b，且将偏置项的变量x0设置为1，则可以将上述的梯度合并为：
根据梯度下降法，得到如下更新公式：
利用上述的Logistics Regression中的权重的更新公式，可以实现Logistics Regression模型的训练，下面是利用梯度下降法训练的具体过程，程序如下：
def lr_train_bgd(feature,label,maxCycle,alpha):''':param feature: 特征:param label: 标签:param maxCycle: 最大迭代次数:param alpha: 学习率:return w:权重'''n = np.shape(feature)[1] #特征个数w = np.mat(np.ones((n,1)))#初始化权重i = 0while i<=maxCycle:#在最大迭代次数的范围内i += 1#当前的跌打次数h = sig(feature*w) #计算Sigmoid值err = label - hif i%100==0:print("\t----------iter="+str(i)+\",train error rate="+str(error_rate(h,label)))w = w+alpha*feature.T*err#修正权重return w
误差函数为：
def error_rate(h,label):'''计算当前的损失函数值:param h: 预测值:param label: 实际值:returnsum_ree / m: 错误率'''m = np.shape(h)(0)sum_err = 0.0for i in range(m):if h[i,0]>0 and (1 - h[i,0])>0:sum_err -= (label[i,0]*np.log(h[i,0]) + \(1-label[i,0]) * np.log(1-h[i,0]))else:sum_err -= 0return sum_err / m
1.7梯度下降法的若干问题
1.为了求解优化问题f（w）的最小值，我们希望每次迭代的结果能够接近最优值W*，对于一维的情况，如图所示：
若当前的梯度为负，则最小值在当前点的右侧，若当前点的梯度为正，则最小值在当前点的左侧，负的梯度即为下降的方向。对于上述的一维的情况，有下述的更新规则：
其中，αi为步长。对于二维的情况，此时更新的规则如下：
▽代表的含义是矢量求偏导
2.步长的选择：
对于步长的选择，若选择太小，会导致收敛速度比较慢，若选择较大，则会出现震荡的现象，即跳过最优解，在最优解附近徘徊，如图：
步长选择过大 步长选择过小
因此，选择合适的步长对于梯度下降的收敛效果显得尤为重要，如图：
