线性分类|机器学习推导系列（四）

一、从线性回归到线性分类

1. 线性回归的特性

$\left\{\begin{matrix} 线性\left\{\begin{matrix} 属性线性\overset{打破}{\rightarrow}属性非线性：特征转换（多项式回归）\\ 全局线性\overset{打破}{\rightarrow}全局非线性：线性分类（**函数是非线性）\\ 系数线性\overset{打破}{\rightarrow}系数非线性：神经网络 \end{matrix}\right.\\ 全局性\overset{打破}{\rightarrow}线性样条回归、决策树\\ 数据未加工\overset{打破}{\rightarrow}PCA、流形 \end{matrix}\right.$⎩⎪⎪⎪⎪⎨⎪⎪⎪⎪⎧​线性⎩⎨⎧​属性线性→打破属性非线性：特征转换（多项式回归）全局线性→打破全局非线性：线性分类（激活函数是非线性）系数线性→打破系数非线性：神经网络​全局性→打破线性样条回归、决策树数据未加工→打破PCA、流形​

线性回归$f(w,b)=w^{T}x+b$f(w,b)=wTx+b具备线性、全局性和数据未加工的特性。

线性包括三个方面，其中属性线性指的是$f(w,b)$f(w,b)关于$x$x是线性的；全局线性指的是$w^{T}x+b$wTx+b只是一个线性组合，然后直接就输出得到$f(w,b)$f(w,b)；系数线性指的是$f(w,b)$f(w,b)关于$w^{T}$wT是线性的。

全局性指的是线性回归是在整个特征空间上学习，并没有将特征空间进行划分然后在每个划分上学习。

数据未加工指的是线性回归直接在给定数据上进行学习没有对数据进行其他的加工。

其他的学习算法跟线性回归比较起来打破了其某些特性，在上面的树状图中也举出了一些例子。

2. 从线性回归到线性分类

$线性回归\overset{**函数}{\underset{降维}{\Rightarrow }}线性分类$线性回归降维⇒​激活函数​线性分类

线性回归经过一个**函数然后根据一个阈值来获得分类的结果，如此就成为线性分类，也可以理解为将$w^{T}x+b$wTx+b降维到一维而获得分类结果。

$**函数y=f(w^{T}x+b)，y\in \left\{\begin{matrix} \left \{0,1\right \}，硬分类\\ [0,1]，软分类 \end{matrix}\right.\\ 函数f叫做**函数（activation\: function）\\ 函数f^{-1}叫做链接函数（link\: function）\\ f:w^{T}x+b \mapsto \left \{0,1\right \}/[0,1]\\ f^{-1}:\left \{0,1\right \}/[0,1]\mapsto w^{T}x+b$激活函数y=f(wTx+b)，y∈{{0,1}，硬分类[0,1]，软分类​函数f叫做激活函数（activationfunction）函数f−1叫做链接函数（linkfunction）f:wTx+b↦{0,1}/[0,1]f−1:{0,1}/[0,1]↦wTx+b

3. 硬分类和软分类

$线性分类\left\{\begin{matrix} 硬分类\left\{\begin{matrix} 线性判别分析\\ 感知机 \end{matrix}\right.\\ 软分类\left\{\begin{matrix} 高斯判别分析\\ 朴素贝叶斯\\ 逻辑回归 \end{matrix}\right. \end{matrix}\right.$线性分类⎩⎪⎪⎪⎪⎨⎪⎪⎪⎪⎧​硬分类{线性判别分析感知机​软分类⎩⎨⎧​高斯判别分析朴素贝叶斯逻辑回归​​

二、感知机

1. 概述

假设有一可以被线性分类的样本集$\left \{(x_{i},y_{i})\right \}_{i=1}^{N}${(xi​,yi​)}i=1N​，其中$x_{i}\in \mathbb{R}^{p},y_{i}\left \{-1,+1\right \}$xi​∈Rp,yi​{−1,+1}。感知机算法使用随机梯度下降法（SGD）来在特征空间$ \mathbb{R}^{{p}$寻找一个超平面$w}{T}x+b=0$来将数据划分为正、负两类，其中$来将数据划分为正、负两类，其中w\in \mathbb{R}^{p}$，是超平面的法向量。

2. 学习策略

感知机的思想是错误驱动。其模型是$f(x)=sign(w^{T}x+b)$f(x)=sign(wTx+b)，$f(x)$f(x)输出该样本点的类别，定义集合M为误分类点的集合。

可以确定对于误分类的数据$(x_{i},y_{i})$(xi​,yi​)来说，满足以下关系：

$\underset{=\left |w^{T}x_{i}+b \right |}{\underline{-y_{i}(w^{T}x_{i}+b)}}>0$=∣wTxi​+b∣−yi​(wTxi​+b)​​>0

损失函数的一个自然选择是误分类点的个数，即$L(w)=\sum_{i=1}^{N}I\left \{-y_{i}(w^{T}x_{i}+b)>0\right \}$L(w)=∑i=1N​I{−yi​(wTxi​+b)>0}，但是这样的损失函数是不可导的，不易优化。因此采用另一种损失函数，即误分类点到超平面的总距离。

在${R}^{p}$Rp空间中任一点$x_{0}$x0​到超平面的距离为:

$\frac{\left | w^{T}x_{0}+b\right |}{\left \| w\right \|}\\ (可以参考初中知识，平面中点到直线的距离：d=\frac{\left | Ax+By+C\right |}{\sqrt{A^{2}+B^{2}}})$∥w∥∣∣​wTx0​+b∣∣​​(可以参考初中知识，平面中点到直线的距离：d=A2+B2​∣Ax+By+C∣​)

因此所有误分类点到超平面$w^{T}x+b=0$wTx+b=0的总距离为：

$\sum _{x_{i}\in M}\frac{\left | w^{T}x_{0}+b\right |}{\left \| w\right \|}=-\frac{1}{{\left \| w\right \|}}\sum _{x_{i}\in M}y_{i}(w^{T}x_{i}+b)$xi​∈M∑​∥w∥∣∣​wTx0​+b∣∣​​=−∥w∥1​xi​∈M∑​yi​(wTxi​+b)

不考虑$\frac{1}{{\left \| w\right \|}}$∥w∥1​，就得到感知机的损失函数：

$L(w,b)=-\sum _{x_{i}\in M}y_{i}(w^{T}x_{i}+b)$L(w,b)=−xi​∈M∑​yi​(wTxi​+b)

3. 学习算法

计算损失函数的梯度：

$\frac{\partial L(w,b)}{\partial w}=-\sum _{x_{i}\in M}y_{i}x_{i}\\ \frac{\partial L(w,b)}{\partial b}=-\sum _{x_{i}\in M}y_{i}$∂w∂L(w,b)​=−xi​∈M∑​yi​xi​∂b∂L(w,b)​=−xi​∈M∑​yi​

感知机的学习算法使用随机梯度下降法（SGD），这里选取$\eta (0<\eta\leq 1)$η(0<η≤1)作为学习率，其学习的步骤如下：
①选取初值$w_{0},b_{0}$w0​,b0​；
②在训练集中选取数据$(x_{i},y_{i})$(xi​,yi​)；
③如果$y_{i}(w^{T}x_{i}+b)\leq 0$yi​(wTxi​+b)≤0，则更新参数：

$w\leftarrow w+\eta y_{i}x_{i}\\ b\leftarrow b+\eta y_{i}$w←w+ηyi​xi​b←b+ηyi​

④转至②，直到训练集中没有误分类点。

截止这里我们都是假设数据是线性可分的，如果线性不可分，可以用口袋算法（pocket algorithm），这里不做过多介绍。

三、线性判别分析

1. 概述

线性判别分析可用于处理二分类问题，其过程是寻找一个最佳的投影方向，使得样本点在该方向上的投影符合类内小、类间大的思想，具体指的是类内的方差之和小，类间的均值之差大。

假设有以下数据：

$X=(x_{1},x_{1},\cdots ,x_{N})^{T}=\begin{pmatrix} x_{1}^{T}\\ x_{2}^{T}\\ \vdots \\ x_{N}^{T} \end{pmatrix}_{N\times p}Y=\begin{pmatrix} y_{1}\\ y_{2}\\ \vdots \\ y_{N} \end{pmatrix}_{N \times 1}\\ \left \{(x_{i},y_{i})\right \}_{i=1}^{N},x_{i}\in \mathbb{R}^{p},y_{i}\in \{\underset{C_{1}}{\underline{+1}},\underset{C_{2}}{\underline{-1}}\}\\ x_{C_{1}}=\left \{x_{i}|y_{i}=+1\right \},x_{C_{2}}=\left \{x_{i}|y_{i}=-1\right \}\\ \left | x_{C_{1}}\right |=N_{1},\left | x_{C_{2}}\right |=N_{2},N_{1}+N_{2}=N$X=(x1​,x1​,⋯,xN​)T=⎝⎜⎜⎜⎛​x1T​x2T​⋮xNT​​⎠⎟⎟⎟⎞​N×p​Y=⎝⎜⎜⎜⎛​y1​y2​⋮yN​​⎠⎟⎟⎟⎞​N×1​{(xi​,yi​)}i=1N​,xi​∈Rp,yi​∈{C1​+1​​,C2​−1​​}xC1​​={xi​∣yi​=+1},xC2​​={xi​∣yi​=−1}∣xC1​​∣=N1​,∣xC2​​∣=N2​,N1​+N2​=N

2. 线性判别分析的损失函数

投影轴的方向向量为$w$w，将样本点往该轴上投影以后的值$z_{i}$zi​为$w^{T}x_{i}$wTxi​，均值和方差按照如下方法计算：

$均值\bar{z}=\frac{1}{N}\sum_{i=1}^{N}z_{i}=\frac{1}{N}\sum_{i=1}^{N}w^{T}x_{i}\\ 方差S_{z}=\frac{1}{N}\sum_{i=1}^{N}(w^{T}x_{i}-\bar{z})(w^{T}x_{i}-\bar{z})^{T}$均值zˉ=N1​i=1∑N​zi​=N1​i=1∑N​wTxi​方差Sz​=N1​i=1∑N​(wTxi​−zˉ)(wTxi​−zˉ)T

接下来计算每一类的均值和方差：

$C_{1}：\\ \bar{z}_{1}=\frac{1}{N_{1}}\sum_{i=1}^{N_{1}}w^{T}x_{i}\\ S_{1}=\frac{1}{N_{1}}\sum_{i=1}^{N_{1}}(w^{T}x_{i}-\bar{z}_{1})(w^{T}x_{i}-\bar{z}_{1})^{T}\\ C_{2}：\\ \bar{z}_{2}=\frac{1}{N_{2}}\sum_{i=1}^{N_{2}}w^{T}x_{i}\\ S_{2}=\frac{1}{N_{2}}\sum_{i=1}^{N_{2}}(w^{T}x_{i}-\bar{z}_{1})(w^{T}x_{i}-\bar{z}_{1})^{T}$C1​：zˉ1​=N1​1​i=1∑N1​​wTxi​S1​=N1​1​i=1∑N1​​(wTxi​−zˉ1​)(wTxi​−zˉ1​)TC2​：zˉ2​=N2​1​i=1∑N2​​wTxi​S2​=N2​1​i=1∑N2​​(wTxi​−zˉ1​)(wTxi​−zˉ1​)T

定义损失函数：

$J(w)=\frac{(\bar{z}_{1}-\bar{z}_{2})^{2}}{S_{1}+S_{2}} (\bar{z}_{1}-\bar{z}_{2})^{2}=(\frac{1}{N_{1}}\sum_{i=1}^{N_{1}}w^{T}x_{i}-\frac{1}{N_{2}}\sum_{i=1}^{N_{2}}w^{T}x_{i})^{2}\\ =[w^{T}(\frac{1}{N_{1}}\sum_{i=1}^{N_{1}}x_{i}-\frac{1}{N_{2}}\sum_{i=1}^{N_{2}}x_{i})]^{2}\\ =[w^{T}(\bar{x}_{C_{1}}-\bar{x}_{C_{2}})]^{2}\\ =w^{T}(\bar{x}_{C_{1}}-\bar{x}_{C_{2}})(\bar{x}_{C_{1}}-\bar{x}_{C_{2}})^{T}w\\ S_{1}=\frac{1}{N_{1}}\sum_{i=1}^{N_{1}}(w^{T}x_{i}-\bar{z}_{1})(w^{T}x_{i}-\bar{z}_{1})^{T}\\ =\frac{1}{N_{1}}\sum_{i=1}^{N_{1}}(w^{T}x_{i}-\frac{1}{N_{1}}\sum_{i=1}^{N_{1}}w^{T}x_{i})(w^{T}x_{i}-\frac{1}{N_{1}}\sum_{i=1}^{N_{1}}w^{T}x_{i})^{T}\\ =\frac{1}{N_{1}}\sum_{i=1}^{N_{1}}w^{T}(x_{i}-\bar{x}_{C_{1}})(x_{i}-\bar{x}_{C_{1}})^{T}w\\ =w^{T}[\frac{1}{N_{1}}\sum_{i=1}^{N_{1}}(x_{i}-\bar{x}_{C_{1}})(x_{i}-\bar{x}_{C_{1}})^{T}]w\\ =w^{T}S_{C_{1}}w\\ S_{1}+S_{2}=w^{T}S_{C_{1}}w+w^{T}S_{C_{2}}w=w^{T}(S_{C_{1}}+S_{C_{2}})w\\ \therefore J(w)=\frac{w^{T}(\bar{x}_{C_{1}}-\bar{x}_{C_{2}})(\bar{x}_{C_{1}}-\bar{x}_{C_{2}})^{T}w}{w^{T}(S_{C_{1}}+S_{C_{2}})w}$J(w)=S1​+S2​(zˉ1​−zˉ2​)2​(zˉ1​−zˉ2​)2=(N1​1​i=1∑N1​​wTxi​−N2​1​i=1∑N2​​wTxi​)2=[wT(N1​1​i=1∑N1​​xi​−N2​1​i=1∑N2​​xi​)]2=[wT(xˉC1​​−xˉC2​​)]2=wT(xˉC1​​−xˉC2​​)(xˉC1​​−xˉC2​​)TwS1​=N1​1​i=1∑N1​​(wTxi​−zˉ1​)(wTxi​−zˉ1​)T=N1​1​i=1∑N1​​(wTxi​−N1​1​i=1∑N1​​wTxi​)(wTxi​−N1​1​i=1∑N1​​wTxi​)T=N1​1​i=1∑N1​​wT(xi​−xˉC1​​)(xi​−xˉC1​​)Tw=wT[N1​1​i=1∑N1​​(xi​−xˉC1​​)(xi​−xˉC1​​)T]w=wTSC1​​wS1​+S2​=wTSC1​​w+wTSC2​​w=wT(SC1​​+SC2​​)w∴J(w)=wT(SC1​​+SC2​​)wwT(xˉC1​​−xˉC2​​)(xˉC1​​−xˉC2​​)Tw​

极大化$J(w)$J(w)就可以使得类内的方差之和小，类间的均值之差大。

3. 线性判别分析的求解

$令\left\{\begin{matrix} S_{b}=(\bar{x}_{C_{1}}-\bar{x}_{C_{2}})(\bar{x}_{C_{1}}-\bar{x}_{C_{2}})^{T}\\ S_{w}=S_{C_{1}}+S_{C_{2}} \end{matrix}\right.\\ 则J(w)=\frac{w^{T}S_{b}w}{w^{T}S_{w}w}=w^{T}S_{b}w(w^{T}S_{w}w)^{-1}\\ \frac{\partial J(w)}{\partial w}=2S_{b}w(w^{T}S_{w}w)^{-1}+w^{T}S_{b}w(-1)(w^{T}S_{w}w)^{-2}2S_{b}w\\ =2[S_{b}w(w^{T}S_{w}w)^{-1}+w^{T}S_{b}w(-1)(w^{T}S_{w}w)^{-2}S_{b}w]=0\\ （w^{T}S_{b}w,w^{T}S_{w}w\in \mathbb{R}）\\ \Rightarrow S_{b}w(w^{T}S_{w}w)-w^{T}S_{b}wS_{w}w=0\\ \Rightarrow w^{T}S_{b}wS_{w}w=S_{b}w(w^{T}S_{w}w)\\ \Rightarrow S_{w}w=\frac{w^{T}S_{w}w}{w^{T}S_{b}w}S_{b}w\\ \Rightarrow w=\frac{w^{T}S_{w}w}{w^{T}S_{b}w}S_{w}^{-1}S_{b}w\\ （要注意对于w我们只关注它的方向，不关心它的大小。）\\ \Rightarrow w\propto S_{w}^{-1}S_{b}w\\ \Rightarrow w\propto S_{w}^{-1}(\bar{x}_{C_{1}}-\bar{x}_{C_{2}})\underset{1维}{\underbrace{(\bar{x}_{C_{1}}-\bar{x}_{C_{2}})^{T}w}}\\ \Rightarrow w\propto S_{w}^{-1}(\bar{x}_{C_{1}}-\bar{x}_{C_{2}})\\ \Rightarrow w\propto (S_{C_{1}}+S_{C_{2}})^{-1}(\bar{x}_{C_{1}}-\bar{x}_{C_{2}})$令{Sb​=(xˉC1​​−xˉC2​​)(xˉC1​​−xˉC2​​)TSw​=SC1​​+SC2​​​则J(w)=wTSw​wwTSb​w​=wTSb​w(wTSw​w)−1∂w∂J(w)​=2Sb​w(wTSw​w)−1+wTSb​w(−1)(wTSw​w)−22Sb​w=2[Sb​w(wTSw​w)−1+wTSb​w(−1)(wTSw​w)−2Sb​w]=0（wTSb​w,wTSw​w∈R）⇒Sb​w(wTSw​w)−wTSb​wSw​w=0⇒wTSb​wSw​w=Sb​w(wTSw​w)⇒Sw​w=wTSb​wwTSw​w​Sb​w⇒w=wTSb​wwTSw​w​Sw−1​Sb​w（要注意对于w我们只关注它的方向，不关心它的大小。）⇒w∝Sw−1​Sb​w⇒w∝Sw−1​(xˉC1​​−xˉC2​​)1维(xˉC1​​−xˉC2​​)Tw​​⇒w∝Sw−1​(xˉC1​​−xˉC2​​)⇒w∝(SC1​​+SC2​​)−1(xˉC1​​−xˉC2​​)

进一步如果$S_{w}^{-1}$Sw−1​是各向同性的对角矩阵的话，$w\propto \bar{x}_{C_{1}}-\bar{x}_{C_{2}}$w∝xˉC1​​−xˉC2​​。

四、逻辑回归

1. 概述

逻辑回归是一种二分类算法，通过$sigmoid$sigmoid**函数将线性组合$w^{T}x$wTx压缩到$0$0和$1$1之间来代表属于某一个分类的概率。

假设有如下数据：

$\left \{(x_{i},y_{i})\right \}_{i=1}^{N},x_{i}\in \mathbb{R}^{p},y_{i}\in \left \{0,1\right \}${(xi​,yi​)}i=1N​,xi​∈Rp,yi​∈{0,1}

2. $sigmoid$sigmoid**函数

$\sigma (z)=\frac{1}{1+e^{-z}}$σ(z)=1+e−z1​

其图像为：

3. 逻辑回归的模型

逻辑回归预测$y=1$y=1的概率，然后根据极大似然估计法来求解。

$\left.\begin{matrix} p_{1}=P(y=1|x)=\sigma (w^{T}x)=\frac{1}{1+e^{-w^{T}x}}=\varphi (x;w)\\ p_{0}=P(y=0|x)=1-P(y=1|x)=\frac{e^{-w^{T}x}}{1+e^{-w^{T}x}}=1-\varphi (x;w) \end{matrix}\right\}p(y|x)=p_{1}^{y}p_{0}^{1-y}$p1​=P(y=1∣x)=σ(wTx)=1+e−wTx1​=φ(x;w)p0​=P(y=0∣x)=1−P(y=1∣x)=1+e−wTxe−wTx​=1−φ(x;w)​}p(y∣x)=p1y​p01−y​

4. 逻辑回归的求解

$\hat{w}=\underset{w}{argmax}\: log\underset{likelihood}{\underbrace{P(Y|X)}}\\ =\underset{w}{argmax}\: log\prod_{i=1}^{N}P(y_{i}|x_{i})\\ =\underset{w}{argmax}\sum _{i=1}^{N}logP(y_{i}|x_{i})\\ =\underset{w}{argmax}\sum _{i=1}^{N}(y_{i}log\: p_{1}+(1-y_{i})log\: p_{0})\\ =\underset{w}{argmax}\underset{-\: cross\: entropy}{\underbrace{\sum _{i=1}^{N}(y_{i}log\: \varphi (x;w)+(1-y_{i})log(1-\varphi (x;w))}}$w^=wargmax​loglikelihoodP(Y∣X)​​=wargmax​logi=1∏N​P(yi​∣xi​)=wargmax​i=1∑N​logP(yi​∣xi​)=wargmax​i=1∑N​(yi​logp1​+(1−yi​)logp0​)=wargmax​−crossentropyi=1∑N​(yi​logφ(x;w)+(1−yi​)log(1−φ(x;w))​​

因此这里的极大似然估计就等价于极小化交叉熵损失函数。

求导的过程较为简单，就不做展示了。

五、高斯判别分析

1. 概述

假设有如下数据：

$\left \{(x_{i},y_{i})\right \}_{i=1}^{N},x_{i}\in \mathbb{R}^{p},y_{i}\in \left \{0,1\right \}${(xi​,yi​)}i=1N​,xi​∈Rp,yi​∈{0,1}

2. 高斯判别分析的模型

在高斯判别分析中样本数据的类别$y$y在给定的情况下服从伯努利分布，另外不同类别中的样本数据分别服从多元高斯分布，因此有以下模型：

$y\sim Bernoulli(\phi )\Leftrightarrow \left.\begin{matrix} \phi ^{y},y=1\\ (1-\phi)^{1-y},y=0 \end{matrix}\right\}P(y)=\phi^{y}\phi^{1-y}\\ \left.\begin{matrix} x|y=1\sim N(\mu _{1},\Sigma )\\ x|y=0\sim N(\mu _{2},\Sigma ) \end{matrix}\right\}P(x|y)=N(\mu _{1},\Sigma )^{y}N(\mu _{2},\Sigma )^{1-y}$y∼Bernoulli(ϕ)⇔ϕy,y=1(1−ϕ)1−y,y=0​}P(y)=ϕyϕ1−yx∣y=1∼N(μ1​,Σ)x∣y=0∼N(μ2​,Σ)​}P(x∣y)=N(μ1​,Σ)yN(μ2​,Σ)1−y

这里假设两个高斯分布具有同样的方差。

3. 高斯判别模型的求解

• 损失函数

高斯判别模型的损失函数为其$log$log似然，要估计的参数$\theta$θ为$(\mu _{1},\mu _{2},\Sigma ,\phi )$(μ1​,μ2​,Σ,ϕ):

$L(\theta )=log\prod_{i=1}^{N}P(x_{i},y_{i})\\ =\sum _{i=1}^{N}logP(x_{i},y_{i})\\ =\sum _{i=1}^{N}logP(x_{i}|y_{i})P(y_{i})\\ =\sum _{i=1}^{N}[logP(x_{i}|y_{i})+logP(y_{i})]\\ =\sum _{i=1}^{N}[logN(\mu _{1},\Sigma )^{y_{i}}N(\mu _{2},\Sigma )^{1-y_{i}}+log\phi^{y_{i}}\phi^{1-y_{i}}]\\ =\sum _{i=1}^{N}[\underset{①}{\underbrace{logN(\mu _{1},\Sigma )^{y_{i}}}}+\underset{②}{\underbrace{logN(\mu _{2},\Sigma )^{1-y_{i}}}}+\underset{③}{\underbrace{log\phi^{y_{i}}\phi^{1-y_{i}}}}]$L(θ)=logi=1∏N​P(xi​,yi​)=i=1∑N​logP(xi​,yi​)=i=1∑N​logP(xi​∣yi​)P(yi​)=i=1∑N​[logP(xi​∣yi​)+logP(yi​)]=i=1∑N​[logN(μ1​,Σ)yi​N(μ2​,Σ)1−yi​+logϕyi​ϕ1−yi​]=i=1∑N​[①logN(μ1​,Σ)yi​​​+②logN(μ2​,Σ)1−yi​​​+③logϕyi​ϕ1−yi​​​]

然后使用极大似然估计法来求解：

$\hat{\theta }=\underset{\theta}{argmax}L(\theta )$θ^=θargmax​L(θ)

定义标签为$1$1的样本个数为$N_{1}$N1​，标签为$0$0的样本个数为$N_{2}$N2​，则有$N_{1}+N_{2} `= N$N1​+N2​‘=N。

• 求解$\phi$ϕ

$\phi$ϕ只存在于③式中，因此求解$\phi$ϕ只需要看③式即可：

$③=\sum _{i=1}^{N}[y_{i}log\phi +(1-y_{i})log(1-\phi )]\\ \frac{\partial ③}{\partial \phi}=\sum _{i=1}^{N}[y_{i}\frac{1}{\phi}-(1-y_{i})\frac{1}{1-\phi }]=0\\ \Rightarrow \sum _{i=1}^{N}[y_{i}(1-\phi)-(1-y_{i})\phi ]=0\\ \Rightarrow \sum _{i=1}^{N}(y_{i}-\phi)=0\\ \Rightarrow \sum _{i=1}^{N}y_{i}-N\phi=0\\ \hat{\phi}=\frac{1}{N}\sum _{i=1}^{N}y_{i}=\frac{N_{1}}{N}$③=i=1∑N​[yi​logϕ+(1−yi​)log(1−ϕ)]∂ϕ∂③​=i=1∑N​[yi​ϕ1​−(1−yi​)1−ϕ1​]=0⇒i=1∑N​[yi​(1−ϕ)−(1−yi​)ϕ]=0⇒i=1∑N​(yi​−ϕ)=0⇒i=1∑N​yi​−Nϕ=0ϕ^​=N1​i=1∑N​yi​=NN1​​

• 求解$\mu_{1}、\mu_{2}$μ1​、μ2​

$\phi$ϕ只存在于①式中，因此求解$\mu_{1}$μ1​只需要看①式即可：

$①=\sum _{i=1}^{N}y_{i}log\frac{1}{(2\pi )^{p/2}|\Sigma |^{1/2}}exp\left \{-\frac{1}{2}(x_{i}-\mu _{1})^{T}\Sigma ^{-1}(x_{i}-\mu _{1})\right \}\\ \mu _{1}=\underset{\mu _{1}}{argmax}①\\ =\underset{\mu _{1}}{argmax}\sum _{i=1}^{N}y_{i}[log\frac{1}{(2\pi )^{p/2}|\Sigma |^{1/2}}-\frac{1}{2}(x_{i}-\mu _{1})^{T}\Sigma ^{-1}(x_{i}-\mu _{1})]\\ =\underset{\mu _{1}}{argmax}\sum _{i=1}^{N}y_{i}[-\frac{1}{2}(x_{i}-\mu _{1})^{T}\Sigma ^{-1}(x_{i}-\mu _{1})]\\ =\underset{\mu _{1}}{argmax}\Delta\\ \Delta=\sum _{i=1}^{N}y_{i}[-\frac{1}{2}(x_{i}-\mu _{1})^{T}\Sigma ^{-1}(x_{i}-\mu _{1})]\\ =-\frac{1}{2}\sum _{i=1}^{N}y_{i}(x_{i}^{T}\Sigma ^{-1}-\mu _{1}^{T}\Sigma ^{-1})(x_{i}-\mu _{1})\\ =-\frac{1}{2}\sum _{i=1}^{N}y_{i}(x_{i}^{T}\Sigma ^{-1}x_{i}-2\mu _{1}^{T}\Sigma ^{-1}x_{i}+\mu _{1}^{T}\Sigma ^{-1}\mu _{1})\\ \frac{\partial \Delta}{\partial \mu _{1}}=-\frac{1}{2}\sum _{i=1}^{N}y_{i}(-2\Sigma ^{-1}x_{i}+2\Sigma ^{-1}\mu _{1})=0\\ \Rightarrow \sum _{i=1}^{N}y_{i}(\Sigma ^{-1}\mu _{1}-\Sigma ^{-1}x_{i})=0\\ \Rightarrow \sum _{i=1}^{N}y_{i}(\mu _{1}-x_{i})=0\\ \Rightarrow \sum _{i=1}^{N}y_{i}\mu _{1}=\sum _{i=1}^{N}y_{i}x_{i}\\ \hat{\mu _{1}}=\frac{\sum _{i=1}^{N}y_{i}x_{i}}{\sum _{i=1}^{N}y_{i}}=\frac{\sum _{i=1}^{N}y_{i}x_{i}}{N_{1}}\\ 同理\hat{\mu _{2}}=\frac{\sum _{i=1}^{N}y_{i}x_{i}}{N_{2}}$①=i=1∑N​yi​log(2π)p/2∣Σ∣1/21​exp{−21​(xi​−μ1​)TΣ−1(xi​−μ1​)}μ1​=μ1​argmax​①=μ1​argmax​i=1∑N​yi​[log(2π)p/2∣Σ∣1/21​−21​(xi​−μ1​)TΣ−1(xi​−μ1​)]=μ1​argmax​i=1∑N​yi​[−21​(xi​−μ1​)TΣ−1(xi​−μ1​)]=μ1​argmax​ΔΔ=i=1∑N​yi​[−21​(xi​−μ1​)TΣ−1(xi​−μ1​)]=−21​i=1∑N​yi​(xiT​Σ−1−μ1T​Σ−1)(xi​−μ1​)=−21​i=1∑N​yi​(xiT​Σ−1xi​−2μ1T​Σ−1xi​+μ1T​Σ−1μ1​)∂μ1​∂Δ​=−21​i=1∑N​yi​(−2Σ−1xi​+2Σ−1μ1​)=0⇒i=1∑N​yi​(Σ−1μ1​−Σ−1xi​)=0⇒i=1∑N​yi​(μ1​−xi​)=0⇒i=1∑N​yi​μ1​=i=1∑N​yi​xi​μ1​^​=∑i=1N​yi​∑i=1N​yi​xi​​=N1​∑i=1N​yi​xi​​同理μ2​^​=N2​∑i=1N​yi​xi​​

• 求解$\Sigma$Σ

以下是求解过程中用到的一些预备知识：

$\frac{\partial AB}{\partial A}=B^{T}\\ \frac{\partial |A|}{\partial A}=|A|A^{-1}\\ tr(AB)=tr(BA)\\ tr(ABC)=tr(CAB)=tr(BCA)$∂A∂AB​=BT∂A∂∣A∣​=∣A∣A−1tr(AB)=tr(BA)tr(ABC)=tr(CAB)=tr(BCA)

两类数据按照以下两个集合来表示：

$C_{1}=\left \{x_{i}|y_{i}=1,i=1,2,\cdots ,N\right \}\\ C_{2}=\left \{x_{i}|y_{i}=0,i=1,2,\cdots ,N\right \}\\ |C_{1}|=N_{1},|C_{2}|=N_{2},N_{1}+N_{2}=N$C1​={xi​∣yi​=1,i=1,2,⋯,N}C2​={xi​∣yi​=0,i=1,2,⋯,N}∣C1​∣=N1​,∣C2​∣=N2​,N1​+N2​=N

然后进行求解：

$\hat{\Sigma }=\underset{\Sigma }{argmax}(①+②)\\ ①+②=\sum _{x_{i}\in C_{1}}logN(\mu _{1},\Sigma )+\sum _{x_{i}\in C_{2}}logN(\mu _{2},\Sigma )$Σ^=Σargmax​(①+②)①+②=xi​∈C1​∑​logN(μ1​,Σ)+xi​∈C2​∑​logN(μ2​,Σ)

然后求解上式中的通项：

$\sum_{i=1}^{N}logN(\mu ,\Sigma )=\sum_{i=1}^{N}log\frac{1}{(2\pi )^{p/2}|\Sigma |^{1/2}}exp\left \{-\frac{1}{2}(x_{i}-\mu )^{T}\Sigma ^{-1}(x_{i}-\mu )\right \}\\ =\sum_{i=1}^{N}[log\frac{1}{(2\pi )^{p/2}}+log|\Sigma |^{-\frac{1}{2}}-\frac{1}{2}(x_{i}-\mu )^{T}\Sigma ^{-1}(x_{i}-\mu )]\\ =\sum_{i=1}^{N}[C-\frac{1}{2}log|\Sigma |-\frac{1}{2}(x_{i}-\mu )^{T}\Sigma ^{-1}(x_{i}-\mu )]\\ =C-\frac{1}{2}Nlog|\Sigma |-\frac{1}{2}\sum_{i=1}^{N}(x_{i}-\mu )^{T}\Sigma ^{-1}(x_{i}-\mu )\\ =C-\frac{1}{2}Nlog|\Sigma |-\frac{1}{2}\sum_{i=1}^{N}tr[(x_{i}-\mu )^{T}\Sigma ^{-1}(x_{i}-\mu )]\\ (实数的迹等于其本身)\\ =C-\frac{1}{2}Nlog|\Sigma |-\frac{1}{2}\sum_{i=1}^{N}tr[(x_{i}-\mu )(x_{i}-\mu )^{T}\Sigma ^{-1}]\\ =C-\frac{1}{2}Nlog|\Sigma |-\frac{1}{2}tr[\sum_{i=1}^{N}(x_{i}-\mu )(x_{i}-\mu )^{T}\Sigma ^{-1}]\\ =C-\frac{1}{2}tr(NS\Sigma ^{-1})\\ (S=\frac{1}{N}\sum_{i=1}^{N}(x_{i}-\mu )(x_{i}-\mu )^{T},为协方差矩阵)\\ =-\frac{1}{2}Nlog|\Sigma |-\frac{1}{2}Ntr(S\Sigma ^{-1})+C\\ 则①+②=-\frac{1}{2}N_{1}log|\Sigma |-\frac{1}{2}N_{1}tr(S_{1}\Sigma ^{-1})-\frac{1}{2}N_{2}log|\Sigma |-\frac{1}{2}N_{2}tr(S_{2}\Sigma ^{-1})+C\\ =-\frac{1}{2}Nlog|\Sigma |-\frac{1}{2}N_{1}tr(S_{1}\Sigma ^{-1})-\frac{1}{2}N_{2}tr(S_{2}\Sigma ^{-1})+C\\ =-\frac{1}{2}[Nlog|\Sigma |+N_{1}tr(S_{1}\Sigma ^{-1})+N_{2}tr(S_{2}\Sigma ^{-1})]+C$i=1∑N​logN(μ,Σ)=i=1∑N​log(2π)p/2∣Σ∣1/21​exp{−21​(xi​−μ)TΣ−1(xi​−μ)}=i=1∑N​[log(2π)p/21​+log∣Σ∣−21​−21​(xi​−μ)TΣ−1(xi​−μ)]=i=1∑N​[C−21​log∣Σ∣−21​(xi​−μ)TΣ−1(xi​−μ)]=C−21​Nlog∣Σ∣−21​i=1∑N​(xi​−μ)TΣ−1(xi​−μ)=C−21​Nlog∣Σ∣−21​i=1∑N​tr[(xi​−μ)TΣ−1(xi​−μ)](实数的迹等于其本身)=C−21​Nlog∣Σ∣−21​i=1∑N​tr[(xi​−μ)(xi​−μ)TΣ−1]=C−21​Nlog∣Σ∣−21​tr[i=1∑N​(xi​−μ)(xi​−μ)TΣ−1]=C−21​tr(NSΣ−1)(S=N1​i=1∑N​(xi​−μ)(xi​−μ)T,为协方差矩阵)=−21​Nlog∣Σ∣−21​Ntr(SΣ−1)+C则①+②=−21​N1​log∣Σ∣−21​N1​tr(S1​Σ−1)−21​N2​log∣Σ∣−21​N2​tr(S2​Σ−1)+C=−21​Nlog∣Σ∣−21​N1​tr(S1​Σ−1)−21​N2​tr(S2​Σ−1)+C=−21​[Nlog∣Σ∣+N1​tr(S1​Σ−1)+N2​tr(S2​Σ−1)]+C

然后对$\Sigma$Σ进行求导：

$\frac{\partial ①+②}{\partial \Sigma}=-\frac{1}{2}[N\frac{1}{|\Sigma |}|\Sigma |\Sigma^{-1}+N_{1}\frac{\partial tr(\Sigma ^{-1}S_{1})}{\partial \Sigma}+N_{2}\frac{\partial tr(\Sigma ^{-1}S_{2})}{\partial \Sigma}]\\ =-\frac{1}{2}[N\Sigma^{-1}+N_{1}S_{1}^{T}(-1)\Sigma ^{-2}+N_{2}S_{2}^{T}(-1)\Sigma ^{-2}]\\ =-\frac{1}{2}(N\Sigma^{-1}-N_{1}S_{1}^{T}\Sigma ^{-2}-N_{2}S_{2}^{T}\Sigma ^{-2})\\ =0\\ N\Sigma-N_{1}S_{1}-N_{2}S_{2}=0\\ \hat{\Sigma}=\frac{N_{1}S_{1}+N_{2}S_{2}}{N}$∂Σ∂①+②​=−21​[N∣Σ∣1​∣Σ∣Σ−1+N1​∂Σ∂tr(Σ−1S1​)​+N2​∂Σ∂tr(Σ−1S2​)​]=−21​[NΣ−1+N1​S1T​(−1)Σ−2+N2​S2T​(−1)Σ−2]=−21​(NΣ−1−N1​S1T​Σ−2−N2​S2T​Σ−2)=0NΣ−N1​S1​−N2​S2​=0Σ^=NN1​S1​+N2​S2​​

六、朴素贝叶斯

1. 概述

假设有如下数据：

$\left \{(x_{i},y_{i})\right \}_{i=1}^{N},x_{i}\in \mathbb{R}^{p},y_{i}\in \left \{0,1\right \}${(xi​,yi​)}i=1N​,xi​∈Rp,yi​∈{0,1}

2. 朴素贝叶斯的模型

朴素贝叶斯分类器可以用来做多分类，其基本思想是条件独立性假设，即假设数据的每个特征之间是相互独立的，其形式化表达为

$对于x|y来说x^{i}|y与x^{j}|y是相互独立的（i\neq j），即\\ P(x|y)=\prod_{i=1}^{p}P(x^{i}|y)$对于x∣y来说xi∣y与xj∣y是相互独立的（i​=j），即P(x∣y)=i=1∏p​P(xi∣y)

朴素贝叶斯分类器是最简单的概率图模型（有向图）:

给定$x$x，判断$y$y的类别可以通过以下方法，即将$x$x归为类别的概率中最大的一类：

$\hat{y}=\underset{y}{argmax}P(y|x)\\ =\underset{y}{argmax}\frac{P(x,y)}{P(x)}\\ =\underset{y}{argmax}\frac{P(y)P(x|y)}{P(x)}\\ =\underset{y}{argmax}P(y)P(x|y)$y^​=yargmax​P(y∣x)=yargmax​P(x)P(x,y)​=yargmax​P(x)P(y)P(x∣y)​=yargmax​P(y)P(x∣y)

$P(y)$P(y)是先验概率，如果有两类则服从伯努利分布（Bernoulli distribution），如果有多类则服从类别分布（Categorical distribution）。$P(x|y)$P(x∣y)则符合条件独立性假设$P(x|y)=\prod_{i=1}^{p}P(x^{i}|y)$P(x∣y)=∏i=1p​P(xi∣y)，其中对于$x^{i}$xi，如果$x^{i}$xi是离散的，则可以认为其服从类别分布（Categorical distribution），如果$x^{i}$xi是连续的，则可以认为其服从高斯分布（Gaussian distribution）。

至于其求解过程则可以根据具体情况使用极大似然估计法即可。对于朴素贝叶斯方法重要的是理解其条件独立性假设，这个假设也是其被称为“朴素（Naive）”的原因。

参考资料

ref:李航《统计学习方法》

公众号同步更新

---

---