机器学习全解 - Logistic回归（含算法讲解，公式全解，手写代码实现，调包实现）-CSDN博客

本文链接：https://blog.csdn.net/qq_51222843/article/details/138721407

1. Logistic回归

1.1 概念和适用范围

Logistic回归（逻辑回归）是一种统计方法，主要用于预测二分类（多分类不常用，不做介绍，下同）结果，称为因变量，可以是某疾病是否复发、是否死亡、是否再入院等。逻辑回归的基本思想是使用逻辑函数（通常是Sigmoid函数）将线性回归模型的输出转换成概率。这种转换使得逻辑回归模型能够处理分类问题，尤其是二分类问题。

通俗来讲就是根据数据集进行分类，分类的依据标准是概率。找到一个概率阈值，例如50%：样本属于A类的概率大于50%则为A类，反之为B类。

但实际上，不是所有的数据我们都可以一股脑的将他们填充进模型内部进行训练，在统计学中，是有前置条件：这些特征需要相互独立。

1.2 回归模型

这是一个通用的回归公式，根据不同的自变量x我们得到不同的因变量值，他们之间有着线性关系。因变量是一个连续变量，取值可以任意限度。通过上述模型的定义，我们想把它限制在一个概率中，也就是将因变量的值限制在0到1之间，第一反应就是对数函数。

对数函数将因变量的数值限制在[0,1)之间，我们将函数进行变换：

感觉万事俱备了？

但还有一个问题，作为一个二分类问题，发生了一件事那么必然另一件事没有发生。例如我是长头发，那么说明我一定不是短头发(有点废话文学那意思)。简单来说，也就是说这两个分类的概率相加一定为1。

引出几率(odd)就可以很好的将这个问题解决。

几率（odd)是指事件发生的概率与不发生的概率之比，假设事件 A 发生的概率为P，不发生的概率为1-P，那么事件 A 的几率为:

相反B的几率为1 - odd(A)。

根据上述取极限的情况，这样子有一个很好的相对关系。万事俱备，只欠东风，我们将函数进行汇总一下：

将上述公式进行化简得到我们最后的sigmoid函数：

以下为sigmoid函数图像：

1.3 梯度上升法

所以说，Logistic算法就是找到一个合适的线性系数。

那么问题来了！如何得到合适的参数向量θ?

梯度上升法是一种很好的方式：

实际上就是一种迭代算法，不断去逼近函数最大值， $\alpha$ 代表步长，每次迭代向前走多少步。

梯度上升法求解：

推倒方式在这边文章就不多加细说，下面链接将梯度上升法推导讲解的简单易懂，有兴趣的可以看看。

Logistic回归基础篇之梯度上升算法-CSDN博客

2. 源代码算法实现

def sigmoid(inX):
    return 1.0 / (1+np.exp(-inX))

首先将我们的sigmoid函数实现。

def gradAscnet(dataMatIn,classLabels):
    dataMatrix = np.mat(dataMatIn) #进行矩阵的创建
    labelMat = np.mat(classLabels).transpose() #进行矩阵的转置
    m,n = np.shape(dataMatrix)
    alpha = 0.001 #定义步长
    maxCycles = 500 #迭代次数
    weights = np.ones((n,1)) #系数
    for k in range(maxCycles):
        h = sigmoid(dataMatrix*weights)
        error = (labelMat - h) #误差
        weights = weights + alpha*dataMatrix.transpose() * error #更新系数
    return weights
gradAscnet(dataMat,classMat)

我们实现梯度上升算法进行迭代系数的值，通过迭代maxCycles次数得到最终回归模型的系数值，计算出误差，再代入公式最后得到回归系数。

梯度上升算法在每次更新回归系数的时候都需要遍历整个数据集，如果有成千上万个特征这个方法的计算复杂度就太高了。随机梯度上升算法就可以很好的去解决这个问题：

def stocGradAscent0(dataMatrix,classLabels):
    m,n = np.shape(dataMatrix)
    alpha = 0.01
    weights = np.ones(n)
    for i in range(m):
        h = sigmoid(sum(dataMatrix[i]*weights))
        error = (classLabels[i] - h)
        weights = weights + alpha*error*dataMatrix[i]
    return weights

3. 调包实现

import numpy as np
from sklearn.model_selection import train_test_split
from sklearn.linear_model import LogisticRegression
dataMat = np.array([[99,98,1],[97,100,1],[1,2,0],[2,4,0]])
classMat = np.array([1,1,0,0])
Model = LogisticRegression()
Model.fit(dataMat,classMat)