《机器学习》周志华-CH7（贝叶斯分类）

7.1贝叶斯决策论

  对分类任务而言，在所有相关概率已知的理想情形下，贝叶斯决策论考虑如何基于这些概率核误判损失来选择最优的类别标记。

$$\begin{array}{c}R(x_i|x)=\sum _{j=1}^N{\lambda }_{ij}P(c_j|x)\end{array}\text{\hspace{1em}(7.1)}$$

  $h^{\ast }$称为贝叶斯最优分类器

  总体风险$R(h^{\ast })$称为贝叶斯风险。

  $1-R(h^{\ast })$反映了分类器所能达到的最佳性能，即模型精度上限。

  具体而言，若目标是最小化分类错误率，则误判损失${\lambda }_{ij}$可写为：

$${\lambda }_{ij}=\{\begin{array}{ll}0& ifi=j\\ 1& other\end{array}$$

  条件风险$R(c|x)=1-P(c|x)$

  最小化分类错误率的贝叶斯最优分类器为

  对每个样本$x$，选后验概率$P(c|x)$最大的类别标记

  机器学习所要实现的是基于有限训练样本集尽可能准确地估计出后验概率$P(c|x)$，主要有两种策略

  类的先验概率$P(c)$：表达了各类样本所占比例，根据大数定律，训练集包含充分独立样本，$P(c)$可通过各类样本出现频率估计。

  类的条件概率$P(x|c)$：由于涉及关于$x$所有属性的联合概率。直接根据样本估计可能遇到困难，使用频率来估计不可行，因为“未被观测到”与“出现概率为零”通常是不同的。

7.2极大似然估计

  估计类条件概率记关于类别$c$的类条件概率为$P(x|c)$,假设$P(x|c)$具有确定的形式并且被参数向量${\theta }_c$唯一确定，则我们的任务就是利用训练集$D$估计参数${\theta }_c$，将$P(x|c)$记为$P(x|{\theta }_c)$

  概率模型的训练过程就是参数估计的过程

  参数估计两种不同方案：

• 频率主义学派，认为参数虽然未知，但客观存在，可通过优化似然函数等准则确定参数值。
• 贝叶斯派认为参数是未观察到的随机变量，其本身可有分布，因此可假设服从一个先验分布，然后基于观测到的数据来计算参数的后验分布。

  极大似然估计MLE，根据数据采样来估计概率分布

  令$D_c$表示训练集$D$的第$c$类样本集合，假设样本独立同分布。

  参数${\theta }_c$对数据集$D_c$的似然说
$$\begin{array}{c}P(D_c|{\theta }_c)=\prod _{x\in D_c}P(x|{\theta }_c)\end{array}\text{\hspace{1em}(7.9)}$$
  对${\theta }_c$进行极大似然估计，就是去寻找最大化似然$P(D_c|{\theta }_c)$的参数值$\widehat{{\theta }_c}$

  （7.9）连乘造成下溢，通常为对数似然（log-likehood）
$$\begin{gathered} \begin{array}{c}LL({\theta }_c)=logP(D_c|{\theta }_c) \\ =\sum _{x\in D_c}logP(x|{\theta }_c)\end{array}\text{\hspace{1em}(7.10)} \end{gathered}$$
  此时参数${\theta }_c$的极大似然估计$\widehat{{\theta }_c}$为

$$\begin{array}{c}\widehat{{\theta }_c}=\underset{{\theta }_c}{\mathrm{argmax}}LL({\theta }_c)\end{array}\text{\hspace{1em}(7.11)}$$

  估计结果的准确性严重依赖所假设的概率分布形式是否符合潜在的真实数据分布。

7.3朴素贝叶斯分类器

  贝叶斯公式来估计后验概率$P(c|x)$困难在于$P(x|c)$是所有属性的联合概率，难以从有限训练样本直接估计。

  朴素贝叶斯分类器采用“属性条件独立性假设”对已知类别，假设所有属性相互独立，假设每个属性独立地对分类结果发生影响

  朴素贝叶斯分类器的训练过程就是基于训练集$D$来估计类先验概率$P(c)$,并为每个属性来估计条件概率$P(x_i|c)$

  $D_c$表示训练集$D$中第$c$类集合，样本充足，则类先验概率：
$$\begin{array}{c}P(c)=\frac{|D_c|}{|D|}\end{array}\text{\hspace{1em}(7.16)}$$

  离散属性而言，$D_{c,x_i}$表示$D_c$中在第$i$个属性上取值为$x_i$的样本组成的集合，则条件概率$P(x_i|c)$:
$$\begin{array}{c}P(x_i|c)=\frac{|D_{c,x_i}|}{|D_c|}\end{array}\text{\hspace{1em}(7.17)}$$

  下面用西瓜数据集3.0训练一个朴素贝叶斯分类器，对测试例“测1”进行分类

  需注意，若某个属性值在训练中没有与某个类同时出现过，如

  $P_{清脆|是}=P(敲声=清脆|好瓜=是)=\frac{0}{8}=0$

  此时乘积永远是$0$,避免这种情况，在估计概率值时通常进行“平滑”，查用“拉普拉斯修正”

  令$N$表示训练集$D$中可能的类别数，$N_i$表示第$i$个属性可能取值

7.4半朴素贝叶斯分类器

  人们尝试对属性条件独立性假设进行一定程度的放松。

  半朴素贝叶斯分类器基本想法：适当考虑一部分属性间的相互依赖信息，从而既不需要进行完全联合概率计算，又不至于彻底忽略了比较强的属性依赖关系。

  “独依赖估计（One-Dependent Estimator，ODE）” 是半朴素贝叶斯分类器最常用的一种策略。顾名思义，所谓“独依赖”就是假设每个属性在类别之外最多仅依赖一个其他属性，即：

  如何确定每个属性的父属性？不同策略的独依赖分类器又不同

  所有属性都依赖于同一个属性，称为“超父”，$x_1$必是超父属性

  AODE（Averaged One-Dependent Estimator） 尝试将每个属性作为超父构建SPODE，然后将那些具有足够训练数据支持的SPODE集成为最终结果，即：

  其中$D_{xi}$是第$i$个属性取值为$x_i$样本集合，$m^{{}^{\prime }}$为阈值常数

  AODE无需模型选择，既能通过预计算节省预测时间，也能采取懒惰学习方法在预测时在再进行计数，易于实现增量学习。

7.5贝叶斯网

  贝叶斯网（Bayesian network） 亦称“信念网”（belief network），借助有向五环图（Directed Acyclic Graph，简称DAG） 刻画属性依赖关系，并使用条件概率表（Conditional Probability Table，简称CPT）来描述属性联合概率分布。

  一个贝叶斯网$B$ 由结构$G$和参数$\theta$两部分组成，$B=<G,\theta >$

  网络结构$G$是一个有向无环图，其中每个结点对应一个属性，两属性有直接依赖关系则由一条边连接

  参数$\theta$定量描述这种依赖关系，假设属性$x_i$在$G$中的父结点集为${\pi }_i$,则$\theta$包含了每个属性的条件概率表是${\theta }_{x_i|{\pi }_i}=P_B(x_i|{\pi }_i)$

7.5.1结构

  贝叶斯网结构有效地表达了属性间的条件独立性，给定父结点集，贝叶斯网假设每个属性与它的非后裔属性独立。

  $B=<G,\theta >$将属性$x_1,x_2,...,x_d$的联合概率分布定义为
$$\begin{array}{c}{P}_B(x_1,x_2,...,x_d)=\prod _{i=1}^d{P}_B(x_i|{\pi }_i)=\prod _{i=1}^d{\theta }_{x_i|{\pi }_i}\end{array}\text{\hspace{1em}(7.26)}$$
  以图7.2为例，联合概率分布定义为
$$P(x_1,x_2,x_3,x_4,x_5)=P(x_1)P(x_2)P(x_3|x_1)P(x_4|x_1,x_2)P(x_5|x_2)$$
  显然，$x_3$和$x_4$在给定$x_1$的取值时独立，$x_4$和$x_5$在给定$x_2$的取值时独立，

  为了分析有向图中变量间的条件独立性，可使用“有向分离”

  先把有向图转变为一个无向图

• 找出有向图中所有$V$型结构，在$V$型结构两个父结点之间加一条无向边
• 将所有有向边设为无向边

  由此产生的无向图称为“道德图”，令父结点相连的过程称为“道德化”

7.5.2学习

  贝叶斯网学习的首要任务就是根据训练数据集来找出结构最“恰当”的贝叶斯网

  “评分搜索” 先定义一个评分函数，以此来评估贝叶斯网与训练数据的契合程度，基于评分函数在寻找结构最优的贝叶斯网

  “最小描述长度” MDL准则

  每个贝叶斯网描述了一个在训练数据上的概率分布，自有一套编码机制能使那些经常出现的样本有更短的编码。选择综合长度最短的网。

  给定训练集 D = { x 1 , x 2 , . . . , x m } D=\{x_1,x_2,...,x_m\} D={x1​,x2​,...,xm​},贝叶斯网$B=<G,\theta >$在$D$上的评分函数：

$$\begin{array}{c}s(B|D)=f(\theta )|B|-LL(B|D)\end{array}\text{\hspace{1em}(7.28)}$$
  其中$|B|$是贝叶斯网的参数个数；$f(\theta )$表示描述每个参数$\theta$所需的字节数；
$$\begin{array}{c}LL(B|D)=\sum _{i=1}^{m}log{P}_B(x_i)\end{array}\text{\hspace{1em}(7.29)}$$
  $S(B|D)$第1项是描述网$B$字节数，第2项是$B$对应概率分布$P_B$的字节数

• $f(\theta )=1$得AIC评分函数
  $$\begin{array}{c}AIC(B|D)=|B|-LL(B|D)\end{array}\text{\hspace{1em}(7.30)}$$
• $f(\theta )=\frac{1}{2}lo{g}^m$得BIC评分函数
  $$\begin{array}{c}BIC(B|D)=\frac{1}{2}lo{g}^m|B|-LL(B|D)\end{array}\text{\hspace{1em}(7.31)}$$
• $f(\theta )=0$评分函数退化为负对数似然
    若网$B=<G,\theta >$中$G$固定，则$S(B|D)$第1项为常数，参数${\theta }_{x_i|{\pi }_i}$可直接由数据集$D$得到：

$$\begin{array}{c}{\theta }_{x_i|{\pi }_i}={\hat{P}}_D(x_i|{\pi }_i)\end{array}\text{\hspace{1em}(7.32)}$$

7.6EM算法

  前面训练样本都是"完整"的，实际上并不一定都是“完整”。