朴素贝叶斯以及三种常见模型推导

朴素贝叶斯

在机器学习中，朴素贝叶斯分类器是一系列以假设特征之间强（朴素）独立下运用贝叶斯定理为基础的简单概率分类器。

朴素贝叶斯算法Naive Bayes定义中有两个关键定义:特征之间强假设独立和贝叶斯定理.这两个定义就是朴素贝叶斯的关键.接下来先了解一下这两个定义.

贝叶斯定理

贝叶斯定义是概率论中的一个定理,它跟随机变量的条件概率以及边缘概率分布有关.

通常，事件A在事件B（发生）的条件下的概率，与事件B在事件A（发生）的条件下的概率是不一样的，然而,这两者之间是有确定的关系的，贝叶斯定理就是这种关系的陈述。贝叶斯公式的一个用途在于通过已知的三个概率函数推出第四个.

直接给出公式:

P(A∣B)=P(B∣A)P(A)P(B)P(A|B) = \frac{P(B|A)P(A)}{P(B)}P(A∣B)=P(B)P(B∣A)P(A)​

其中,P(A|B)是指事件B发生的情况下事件A发生的概率(条件概率).在贝叶斯定理中,每个名词都有约定俗成的名称:

• P(A|B)是已知B发生后A的条件概率,也由于得知B的取值而被称作A的后验概率;
• P(A)是A的先验概率(或边缘概率).之所以称为"先验"是因为它不考虑任何B方面的因素;
• P(B|A)是已知A发生后B的条件概率,也由于得知A的取值而成称作B的后验概率;
• P(B)是B的先验概率(或边缘概率).

按这些术语,贝叶斯定理可以表述为:

后验概率 = (似然性 * 先验概率)/标准化常量

也就是说,后验概率与先验概率和相似度的乘积成正比.

同时,分母P(B),可以根据全概率公式分解为:

P(B)=∑i=1nP(Ai)P(B∣Ai)P(B)=\sum_{i=1}^nP(Ai)P(B|Ai)P(B)=∑i=1n​P(Ai)P(B∣Ai)

条件独立性假设

如果P(X,Y|Z)=P(X|Z)P(Y|Z)，或等价地P(X|Y,Z）=P(X|Z），则称事件X,Y对于给定事件Z是条件独立的，也就是说，当Z发生时，X发生与否与Y发生与否是无关的。

应用在自然语言处理中,就是说在文章类别确定的条件下,文章的各个特征(单词)在类别确定的条件下是独立的,并不相关,用通俗的话说,在文章类别确定的条件下,文章各个词之间出现与否没有相关性(事实上,并不成立).这是一个非常强的假设,但对问题的求解来说变得更加简单.

朴素贝叶斯的概率模型

设输入空间X⊆RnX \subseteq R^nX⊆Rn为n为向量的集合,输出空间为类标记集合Y=c1,c2,...,ckY = {c_1,c_2,...,c_k}Y=c1​,c2​,...,ck​.输入为特征向量x∈Xx \in Xx∈X,输出为类标记y∈Yy \in Yy∈Y. X是定义在输入空间X上的随机变量,Y是定义在输出空间Y上的随机变量.P(X,Y)是X和Y的联合概率分布.训练数据集:

T={(x1,y1),(x2,y2),...,(xN,yN)}T = \{(x_1,y_1), (x_2,y_2),...,(x_N,y_N)\}T={(x1​,y1​),(x2​,y2​),...,(xN​,yN​)}

由P(X,Y)独立同分布产生.因此朴素贝叶斯模型也是一个生成模型.

朴素贝叶斯算法通过训练集学习联合概率分布P(X,Y),具体地,学习先验概率分布以及条件概率分布.其中先验概率分布

P(Y=ck),k=1,2,...,KP(Y=c_k), k=1,2,...,KP(Y=ck​),k=1,2,...,K

条件概率分布

P(X=x∣Y=ck)=P(X(1)=x(1),X(2)=x(2),...,X(n)=x(n)∣Y=ck)P(X=x|Y=c_k)=P(X^{(1)}=x^{(1)},X^{(2)}=x^{(2)},...,X^{(n)}=x^{(n)}|Y=c_k)P(X=x∣Y=ck​)=P(X(1)=x(1),X(2)=x(2),...,X(n)=x(n)∣Y=ck​), k=1,2,…,K

通过两个概率得到联合概率分布P(X,Y) = P(X|Y)P(Y).

条件概率分布P(X=x|Y=c_k)有指数级数量的参数,其估计实际上不可行的.假设x(j)x^{(j)}x(j)有SjS_jSj​个取值,j=1,2,…,n,Y有K个取值,那么参数个数为K∏j=1nSjK\prod_{j=1}^n S_jK∏j=1n​Sj​.

指数级的参数估计事实上并不可行,因此朴素贝叶斯算法针对各个特征之间做了假设,也就是对条件概率分布作了条件独立性假设,这是一个很强的假设,通过这个假设,我们的参数求解变得可行,这也就是朴素贝叶斯的"朴素"的由来.这种情况下,我们同样假设x(j)x^{(j)}x(j)有SjS_jSj​个取值,j=1,2,…,n,Y有K个取值,那么参数个数为K∑i=1nSjK\sum_{i=1}^nS_jK∑i=1n​Sj​,需要估计的参数量大大减少.条件独立性假设是

P(X=x∣Y=ck)=P(X(1)=x(1),X(2)=x(2),...,X(n)=x(n)∣Y=ck)=∏j=1nP(X(j)=x(j)∣Y=ck)P(X=x|Y=c_k)=P(X^{(1)}=x^{(1)},X^{(2)}=x^{(2)},...,X^{(n)}=x^{(n)}|Y=c_k) \\=\prod_{j=1}^n P(X^{(j)}=x^{(j)}|Y=c_k)P(X=x∣Y=ck​)=P(X(1)=x(1),X(2)=x(2),...,X(n)=x(n)∣Y=ck​)=∏j=1n​P(X(j)=x(j)∣Y=ck​)

朴素贝叶斯算法分类时,对给定输入x,通过学习到的模型计算后验概率分布P(Y=ck∣X=x)P(Y=c_k|X=x)P(Y=ck​∣X=x),将后验概率最大的类作为输入x的类输出.后验概率根据贝叶斯定理计算:

P(Y=ck∣X=x)=P(Y=ck)P(X=x∣Y=ck)∑kP(X=x∣Y=ck)P(Y=ck)=P(Y=ck)∏j=1nP(X(j)=x(j)∣Y=ck)∑kP(Y=ck)∏j=1nP(X(j)=x(j)∣Y=ck)P(Y=c_k|X=x)=\frac{P(Y=c_k)P(X=x|Y=c_k)}{\sum_k P(X=x|Y=c_k)P(Y=c_k)}\\ =\frac{P(Y=c_k)\prod_{j=1}^n P(X^{(j)}=x^{(j)}|Y=c_k)}{\sum_k P(Y=c_k)\prod_{j=1}^n P(X^{(j)}=x^{(j)}|Y=c_k)}P(Y=ck​∣X=x)=∑k​P(X=x∣Y=ck​)P(Y=ck​)P(Y=ck​)P(X=x∣Y=ck​)​=∑k​P(Y=ck​)∏j=1n​P(X(j)=x(j)∣Y=ck​)P(Y=ck​)∏j=1n​P(X(j)=x(j)∣Y=ck​)​

上面的公式是后验概率分布中的一项,由于对于相同输入x下不同类别的后验概率的分母都相同,而最终的类输出是后验概率分布中概率最大对应的类别,所以我们可以简化为只比较分子的大小就可以确定最终的结果,也就是说,最终类输出为:

y=argmaxckP(Y=ck)∏j=1nP(X(j)=x(j)∣Y=ck)y=arg max_{c_k}P(Y=c_k)\prod_{j=1}^n P(X^{(j)}=x^{(j)}|Y=c_k)y=argmaxck​​P(Y=ck​)∏j=1n​P(X(j)=x(j)∣Y=ck​).

如果我们对右边的乘积概率取log,连乘积就可以转换成为和,计算更简单(加法总比乘法简单),上诉公式存在一种变种:

y=argmaxcklogP(Y=ck)+∑j=1nlogP(X(j)=x(j)∣Y=ck)y=arg max_{c_k} logP(Y=c_k) + \sum_{j=1}^n logP(X^{(j)}=x^{(j)}|Y=c_k)y=argmaxck​​logP(Y=ck​)+∑j=1n​logP(X(j)=x(j)∣Y=ck​).

同时这种形式,也可以看做是一种线性回归,权重系数为1.

介绍完,朴素贝叶斯的概率模型之后,我们目前的主要问题就集中在如何估计这个模型的K∑i=1nSjK\sum_{i=1}^nS_jK∑i=1n​Sj​个参数:P(Y=ck),P(X(j)=x(j)∣Y=ck)P(Y=c_k),P(X^{(j)}=x^{(j)}|Y=c_k)P(Y=ck​),P(X(j)=x(j)∣Y=ck​),估算出参数,我们就可以对输入向量x做预测.针对这些参数的求解方法不同,存在不同的朴素贝叶斯类型,具体介绍三种:伯努利朴素贝叶斯,多项式朴素贝叶斯和高斯朴素贝叶斯.不同类型的朴素贝叶斯对参数的求法不同,而根源在于对P条件概率(X=x|Y=c_k)的假设分布不同,也就是说在给定类别的情况下,对X假设的分布不同:伯努利假设是伯努利分布(其实应该是多变量伯努利分布),多项式假设是多项式分布,而高斯也就是假设是高斯分布(其实是多变量高斯分布).然后,我们细化到三种不同类型的朴素贝叶斯理论中.

Bernoulli Naive Bayes 伯努利朴素贝叶斯

伯努利朴素贝叶斯,其实应该叫"Multi-variate Naive Bayes",假设P(X=x|Y=c_k)是多变量伯努利分布.在了解多变量伯努利分布之前,先介绍一下什么是(单变量的)伯努利分布.

伯努利分布

伯努利分布,又叫做两点分布或0-1分布,是一个离散型概率分布.称随机变量X有伯努利分布,参数为p(0< p <1),它分别以概率p和1-p取1和0为值.

最简单的例子就是抛硬币,硬币结果为正或反.

P(X=x)=px(1−p)(1−x)=px+(1−p)(1−x)P(X=x)=p^x (1-p)^{(1-x)}=px+(1-p)(1-x)P(X=x)=px(1−p)(1−x)=px+(1−p)(1−x)

幂次运算变成乘法运算,更简单.当x=1时,概率是P(X=1)=p,当x=0时,概率P(X=0)=1-p,这样就可以将两种情况合在一起.

了解了什么是伯努利分布之后,我们再看什么是多元伯努利分布(多变量 multi-variate Bernoulli).

多元伯努利分布

多元伯努利分布,通俗的讲就是同时进行多个不同的伯努利实验,P(X=x)=θP(X=x)=\thetaP(X=x)=θ,其中x是一个向量,θ\thetaθ也是一个向量,表示不同伯努利实验的参数.

伯努利多项式将文档的生成模型P(X=x|Y=c_k)假设服从为多元伯努利分布,由于我们之前做的特征独立性假设,x=(x(1),x(2)=x(2),...,(n))x=(x^{(1)},x^{(2)}=x^{(2)},...,^{(n)})x=(x(1),x(2)=x(2),...,(n))是一个向量形式,而其中x(j)∈{0,1}x^{(j)} \in \{0, 1\}x(j)∈{0,1},也就是说x向量是onehot形式的向量(每个维度值是0或1),表示这个维度的特征是否出现.特征集F={f1,f2,...,fn}F = \{f_1,f_2,...,f_n \}F={f1​,f2​,...,fn​}有n个特征,特征集的维度决定了输入空间X的维度,而且特征集的维度可以对应到输入空间的各个维度上.

P(X=x∣Y=ck)=P(X(1)=x(1),X(2)=x(2),…,X(n)=x(n)∣Y=ck)=∏i=1np(ti∣Y=ck)xi+(1−p(ti∣Y=ck))(1−xi)P(X=x|Y=c_k)=P(X^{(1)}=x^{(1)},X^{(2)}=x^{(2)},…,X^{(n)}=x^{(n)}|Y=c_k)\\=\prod_{i=1}^{n} p(t_i|Y=c_k)x^i+(1-p(t_i|Y=c_k))(1-x^i)P(X=x∣Y=ck​)=P(X(1)=x(1),X(2)=x(2),…,X(n)=x(n)∣Y=ck​)=∏i=1n​p(ti​∣Y=ck​)xi+(1−p(ti​∣Y=ck​))(1−xi)

因为特征之间的独立性,所以多元伯努利变成各个伯努利分布的连乘积,需要注意的一点是因为是伯努利分布,0-1,特征出现有一个概率p,特征不出现也有一个概率1-p.最终模型的参数估计完成之后,对新样本进行预测时,如果某个特征不出现,需要乘上这个特征不出现的概率,不能只计算特征出现的概率!!!两个向量直接相乘,并不能得到最终的结果.

对应的伯努利朴素贝叶斯模型为:

P(Y=ck∣X=x)=P(Y=ck)P(X=x∣Y=ck)∑kP(X=x∣Y=ck)P(Y=ck)=P(Y=ck)∏j=1np(tj∣Y=ck)xj+(1−p(tj∣Y=ck))(1−xj)∑kP(Y=ck)∏j=1nP(X(j)=x(j)∣Y=ck)P(Y=c_k|X=x)=\frac{P(Y=c_k)P(X=x|Y=c_k)}{\sum_kP(X=x|Y=c_k)P(Y=c_k)}\\=\frac{P(Y=c_k)\prod_{j=1}^n p(t_j|Y=c_k)x^j+(1-p(t_j|Y=c_k))(1-x^j)}{\sum_k P(Y=c_k)\prod_{j=1}^n P(X^{(j)}=x^{(j)}|Y=c_k)}P(Y=ck​∣X=x)=∑k​P(X=x∣Y=ck​)P(Y=ck​)P(Y=ck​)P(X=x∣Y=ck​)​=∑k​P(Y=ck​)∏j=1n​P(X(j)=x(j)∣Y=ck​)P(Y=ck​)∏j=1n​p(tj​∣Y=ck​)xj+(1−p(tj​∣Y=ck​))(1−xj)​

为了简化运算,我们可以将分母忽略,虽然对应的结果不是真正的概率,但是相同样本的各个后验概率之间的大小关系保持不变,同时如果两边同时做log运算,后验概率之间的大小关系同样保持不变.因此,

y=argmaxckP(Y=ck)∏j=1nP(X(j)=x(j)∣Y=ck)=argmaxcklogP(Y=ck)+∑j=1nlogP(X(j)=x(j)∣Y=ck)y=arg max_{c_k} P(Y=c_k)\prod_{j=1}^n P(X^{(j)}=x^{(j)}|Y=c_k) \\ =arg max_{c_k} logP(Y=c_k)+\sum_{j=1}^nlogP(X^{(j)}=x^{(j)}|Y=c_k)y=argmaxck​​P(Y=ck​)∏j=1n​P(X(j)=x(j)∣Y=ck​)=argmaxck​​logP(Y=ck​)+∑j=1n​logP(X(j)=x(j)∣Y=ck​).

了解了多元伯努利分布之后,接下来的工作就是对参数进行估计,计算P(xi∣Y=ck)P(x^i|Y=c_k)P(xi∣Y=ck​),P(Y=ck)P(Y=c_k)P(Y=ck​).

参数估计

参数估计的过程也是朴素贝叶斯分类器学习的过程.而参数估计可以使用极大似然估计.先验概率P(Y=ck)P(Y=c_k)P(Y=ck​)的极大似然估计是

P(Y=ck)=∑i=1NI(yi=ck)NP(Y=c_k)=\frac{\sum_{i=1}^N I(y_i = c_k)}{N}P(Y=ck​)=N∑i=1N​I(yi​=ck​)​, k=1,2,…,K

其中,I(x)是一个指示函数,如果x为真,I(x)结果为1,如果x为假,I(x)=0.用语言描述来说,P(Y=ck)P(Y=c_k)P(Y=ck​)这个概率等于在N个样本的数据集中,类别为ckc_kck​的样本所占的比例.

条件概率P(X(i)=x(i)∣Y=ck)P(X^{(i)}=x^{(i)}|Y=c_k)P(X(i)=x(i)∣Y=ck​)的极大似然估计是:

P(Xi=xi∣Y=ck)=∑i=1NI(Xi=xi,yi=ck)∑i=1NI(yi=ck)P(X^i=x^i|Y=c_k) = \frac{\sum_{i=1}^N I(X^i=x^i, y_i=c_k)}{\sum_{i=1}^N I(y_i=c_k)}P(Xi=xi∣Y=ck​)=∑i=1N​I(yi​=ck​)∑i=1N​I(Xi=xi,yi​=ck​)​

用语言描述来说,条件概率P(Xi=xi∣Y=ck)P(X^i=x^i|Y=c_k)P(Xi=xi∣Y=ck​)等于在类别为ckc_kck​的样本集合(数据集的子集)中,第i个特征等于xix_ixi​的概率,xix^ixi是0或1,而且P(Xi=xi∣Y=ck)P(X^i=x^i|Y=c_k)P(Xi=xi∣Y=ck​)服从伯努利分布,所以只需要计算一个,比如P(Xi=1∣Y=ck)(X^i=1|Y=c_k)(Xi=1∣Y=ck​),P(Xi=0∣Y=ck)=1−P(Xi=1∣Y=ck)P(X^i=0|Y=c_k) = 1 - P(X^i=1|Y=c_k)P(Xi=0∣Y=ck​)=1−P(Xi=1∣Y=ck​),因为两个概率和为1(这是同一个变量).

这些参数估计完之后,朴素贝叶斯就完成了学习过程,接下来就可以使用它去进行预测了(应用才是最终的目的).

0概率处理

由于P(Xi=xi∣Y=ck)P(X^i=x^i|Y=c_k)P(Xi=xi∣Y=ck​)是伯努利分布,参数p在[0,1]之间,有可能存在P(Xi=xi∣Y=ck)=0P(X^i=x^i|Y=c_k) = 0P(Xi=xi∣Y=ck​)=0,也就是存在0概率.

举例来说,在当前类别下的所有样本中特征i都出现了(=1),根据上面的条件概率极大似然估计,可以知道P(Xi=1∣Y=ck)=1P(X^i=1|Y=c_k) = 1P(Xi=1∣Y=ck​)=1,那么对应的,P(Xi=0∣Y=ck)=0P(X^i=0|Y=c_k) = 0P(Xi=0∣Y=ck​)=0,那么当新样本来的时候,加入存在一条记录x,它很巧地没有第i个特征(这不巧了吗?不是),由于0概率的存在,那么使用上面的贝叶斯公式,就会出现属于某个列别的概率为0,P(Y=ck∣X=x)=0P(Y=c_k|X=x)=0P(Y=ck​∣X=x)=0.但这种情况是应该避免的,那么如何避免呢?

我们在对条件概率进行极大似然估计时,针对分子和分母做一些小变动,

P(Xi=xi∣Y=ck)=α+∑i=1NI(Xi=xi,yi=ck)Siα+∑i=1NI(yi=ck)=α+∑i=1NI(Xi=xi,yi=ck)2α+∑i=1NI(yi=ck)P(X^i=x^i|Y=c_k)=\frac{\alpha+\sum_{i=1}^NI(X^i=x^i,y_i=c_k)}{S_i\alpha+\sum_{i=1}^NI(y_i=c_k)} \\ =\frac{\alpha+\sum_{i=1}^NI(X^i=x^i,y_i=c_k)}{2\alpha+\sum_{i=1}^NI(y_i=c_k)}P(Xi=xi∣Y=ck​)=Si​α+∑i=1N​I(yi​=ck​)α+∑i=1N​I(Xi=xi,yi​=ck​)​=2α+∑i=1N​I(yi​=ck​)α+∑i=1N​I(Xi=xi,yi​=ck​)​

其中,SiS_iSi​表示第i个特征不同取值的个数,由于这里是one-hot,取值为2,所以Si=2S_i = 2Si​=2,α\alphaα乘以SiS_iSi​是保证SiS_iSi​个不同取值对应的条件概率之和为1,不偏袒任意一种情况,一视同仁.

代码实现

To Be Continued.

数学推导—矩阵形式,批量计算,高效快捷.

Multinomial Naive Bayes 多项式朴素贝叶斯

多项式朴素贝叶斯,假设P(X=x|Y=c_k)是多项式分布.在了解多项式朴素贝叶斯之前,先介绍一下什么是多项式分布?

多项式分布

将伯努利分布的单变量扩展到d维向量x⃗\vec xx,其中xi∈0,1x_i \in {0,1}xi​∈0,1,且∑i=1dxi=1\sum_{i=1}^d x_i =1∑i=1d​xi​=1,假设xi=1x_i=1xi​=1的概率是μ∈[0,1]\mu \in [0,1]μ∈[0,1],并且∑i=1dμi=1\sum_{i=1}^d \mu_i=1∑i=1d​μi​=1,则将得到离散分布:

P(x∣μ)=∏i=1dμixiP(x|\mu)= \prod_{i=1}^d \mu_i^{x_i}P(x∣μ)=∏i=1d​μixi​​.

其中x,μ\muμ都是d维向量形式.在此的基础上扩展二项分布到多项式分布(Multinomial distribution),该分布描述的是在n次独立实验中有mim_imi​词xi=1x_i=1xi​=1的概率,其密度函数可以表达为如下形式:

p(m1,m2,...,md∣n,μ)=n!m1!m2!...md!∏i=1dμimip(m_1,m_2,...,m_d|n, \mu)=\frac{n!}{m_1!m_2!...m_d!} \prod_{i=1}^d \mu_i^{m_i}p(m1​,m2​,...,md​∣n,μ)=m1​!m2​!...md​!n!​∏i=1d​μimi​​

多项式分布的期望,方差如下: E(x)=nμiE(x)=n \mu_iE(x)=nμi​

var(x)=nμi(1−μi)var(x)=n \mu_i(1-\mu_i)var(x)=nμi​(1−μi​)

多项式朴素贝叶斯

多项式分布应用到朴素贝叶斯上,对于文档分类问题来说,假设给定文档类型的基础上文档生成模型P(X=x∣Y=ck)P(X=x|Y=c_k)P(X=x∣Y=ck​)是一个多项式分布.这样对应关系就是:

• 文档分类中的d维字典(d个特征)对应于多项式分布中的向量的d个维度;
• 文档分类中,词wiw_iwi​出现与否,对应于d维向量中xi∈0,1x_i \in {0,1}xi​∈0,1,两种取值情况,且∑i=1dxi=1\sum_{i=1}^d x_i =1∑i=1d​xi​=1;
• 文档分类中,词wiw_iwi​出现的概率,对应于离散分布中xi=1x_i=1xi​=1的概率是μ∈[0,1]\mu \in [0,1]μ∈[0,1],并且∑i=1dμi=1\sum_{i=1}^d \mu_i=1∑i=1d​μi​=1;
• 文档分类中,给定类别下,对应一次抽样结果x(d维向量)的概率为:P(x∣μ)=∏i=1dμixiP(x|\mu)= \prod_{i=1}^d \mu_i^{x_i}P(x∣μ)=∏i=1d​μixi​​,因为如果一个词不出现,即xi=0x_i=0xi​=0,那么对应μi0=1\mu_i^0=1μi0​=1,所以P(x∣μ)P(x|\mu)P(x∣μ)可以简写为P(xi∣μ)=μiP(x_i|\mu)=\mu_iP(xi​∣μ)=μi​;
• n次独立实验中有nin_ini​词xi=1x_i=1xi​=1的概率,对应到文档模型中特征i(第i个词)出现了nin_ini​次,n次实验对应到文档模型中表示这篇文档的长度为n(一共有n个词),对应概率密度函数为:p(m1,m2,...,md∣n,μ)=n!m1!m2!...md!∏i=1dμimip(m_1,m_2,...,m_d|n, \mu)=\frac{n!}{m_1!m_2!...m_d!} \prod_{i=1}^d \mu_i^{m_i}p(m1​,m2​,...,md​∣n,μ)=m1​!m2​!...md​!n!​∏i=1d​μimi​​

需要注意的是,应用在文本分类的多项式朴素贝叶斯模型之前,一般多项式条件概率如下:

P(X=x∣Y=ck)=P(X1=x1,X2=x2,...,Xd=xd∣Y=ck)=P(∣x∣)n!x1!x2!...xd!∏i=1dP(wi∣Y=ck)xiP(X=x|Y=c_k)=P(X^1=x^1,X^2=x^2,...,X^d=x^d|Y=c_k) \\=P(|x|)\frac{n!}{x^1!x^2!...x^d!} \prod_{i=1}^d P(w_i|Y=c_k)^{x^i}P(X=x∣Y=ck​)=P(X1=x1,X2=x2,...,Xd=xd∣Y=ck​)=P(∣x∣)x1!x2!...xd!n!​∏i=1d​P(wi​∣Y=ck​)xi

我们的多项式朴素贝叶斯概率模型为:

P(Y=ck∣X=x)=P(Y=ck)n!x1!x2!...xd!∏i=1dP(wi∣Y=ck)xiP(X=x)P(Y=c_k|X=x)=\frac{P(Y=c_k)\frac{n!}{x^1!x^2!...x^d!} \prod_{i=1}^d P(w_i|Y=c_k)^{x^i}}{P(X=x)}P(Y=ck​∣X=x)=P(X=x)P(Y=ck​)x1!x2!...xd!n!​∏i=1d​P(wi​∣Y=ck​)xi​

这里为了方便,首先我们假设文章长度和文章的类别之间没有相关性(事实上也并不成立,比如说相对较长的邮件,相比垃圾邮件而言,正常邮件的概率更大),也就是说P(|x|)的分布与文章所属类别无关.另一方面,由于最终所属类别是后验概率最大对应的类别,所以,我们可以将文章长度P(|x|)建模的概率,忽略掉,就像我们之前忽略伯努利分布中的分母一样.

P(X=x∣Y=ck)=P(X1=x1,X2=x2,...,Xd=xd∣Y=ck)=n!x1!x2!...xd!∏i=1dP(wi∣Y=ck)xiP(X=x|Y=c_k) = P(X^1=x^1,X^2=x^2,...,X^d=x^d|Y=c_k)=\frac{n!}{x^1!x^2!...x^d!} \prod_{i=1}^d P(w_i|Y=c_k)^{x^i}P(X=x∣Y=ck​)=P(X1=x1,X2=x2,...,Xd=xd∣Y=ck​)=x1!x2!...xd!n!​∏i=1d​P(wi​∣Y=ck​)xi

再者,为了更加方便,我们通常对两边取log运算,将幂次运算转换成线性运算:

logP(X=x∣Y=ck)=logP(X1=x1,X2=x2,...,Xd=xd∣Y=ck)=n!x1!x2!...xd!∏i=1dP(wi∣Y=ck)xilogP(X=x|Y=c_k) = logP(X^1=x^1,X^2=x^2,...,X^d=x^d|Y=c_k)=\frac{n!}{x^1!x^2!...x^d!} \prod_{i=1}^d P(w_i|Y=c_k)^{x^i}logP(X=x∣Y=ck​)=logP(X1=x1,X2=x2,...,Xd=xd∣Y=ck​)=x1!x2!...xd!n!​∏i=1d​P(wi​∣Y=ck​)xi

我们也可以将文章长度阶乘省略,然后变成:

y=argmaxckP(Y=ck)P(X=x∣Y=ck))=argmaxcklogP(Y=ck)+∑j=1dx(j)logP(wj∣Y=ck)y=arg max_{c_k} P(Y=c_k)P(X=x|Y=c_k))\\=arg max_{c_k} logP(Y=c_k)+\sum_{j=1}^d x^{(j)}logP(w_j|Y=c_k)y=argmaxck​​P(Y=ck​)P(X=x∣Y=ck​))=argmaxck​​logP(Y=ck​)+∑j=1d​x(j)logP(wj​∣Y=ck​).

这样就变成线性运算,就和线性回归一样,运算高效而简单.

将文档模型对应到多项式分布上得到多项式朴素贝叶斯,在我们对其做出假设分布之后,剩下的工作就是对假设分布下每个类别下的d个条件概率以及先验分布进行估计.此外,还需要说明的一点是:多项式朴素贝叶斯模型采用词袋模型,每个xix_ixi​表示第i个特征出现的次数,也就是词频term-frequency,有时候也可以使用tf-idf作为值.

参数估计

参数估计的过程也是朴素贝叶斯分类器学习的过程.而参数估计可以使用极大似然估计.先验概率P(Y=ck)P(Y=c_k)P(Y=ck​)的极大似然估计是

P(Y=ck)=∑i=1NI(yi=ck)NP(Y=c_k)=\frac{\sum_{i=1}^NI(y_i=c_k)}{N}P(Y=ck​)=N∑i=1N​I(yi​=ck​)​, k=1,2,…,K

其中,I(x)是一个指示函数,如果x为真,I(x)结果为1,如果x为假,I(x)=0.用语言描述来说,P(Y=ck)P(Y=c_k)P(Y=ck​)这个概率等于在N个样本的数据集中,类别为ckc_kck​的样本所占的比例.

条件概率P(wt∣Y=ck)P(w_t|Y=c_k)P(wt​∣Y=ck​)的极大似然估计是:

P(wt∣Y=ck)=∑i=1NI(wt=1,yi=ck)xi(t)∑i=1N∑s=1dI(ws=1,yi=ck)xi(s)P(w_t|Y=c_k)=\frac{\sum_{i=1}^N I(w_t=1,y_i=c_k)x_i^{(t)}}{\sum_{i=1}^N \sum_{s=1}^d I(w_s=1,y_i=c_k)x_i^{(s)}}P(wt​∣Y=ck​)=∑i=1N​∑s=1d​I(ws​=1,yi​=ck​)xi(s)​∑i=1N​I(wt​=1,yi​=ck​)xi(t)​​

用语言描述来说,条件概率P(wt∣Y=ck)P(w_t|Y=c_k)P(wt​∣Y=ck​)等于在类别为ckc_kck​的样本集合(数据集的子集)中,第t个特征出现的概率等于ckc_kck​类样本第t个特征出现的总次数(考虑词频,不再是0,1)占ckc_kck​类样本的总词数(文章长度之和,文章单词特征是固定的,考虑了词频)的比例.

为了方便理解,将Nt,kN_{t,k}Nt,k​表示为第k类样本集合中第t个特征出现的总次数,NkN_kNk​表示为在所有样本中第k类样本的总词数(第k类样本长度之和,考虑频数),简写成:

P(wt∣Y=ck)=Nt,kNkP(w_t|Y=c_k)=\frac{N_{t,k}}{N_k}P(wt​∣Y=ck​)=Nk​Nt,k​​

需要注意的是,∑tP(wt∣Y=ck)=1\sum_t P(w_t|Y=c_k)=1∑t​P(wt​∣Y=ck​)=1,意思是给定分类下,各个维度的概率之和为1,在文章分类中,就是给定文章分类的情况下,各个词出现的条件概率之和为1,和每个词出现多少词没有关系.

0概率处理

和伯努利朴素贝叶斯模型类似,有可能存在某一维度,数据集在这一维度上都是0,对应到文档分类上,就是这个词在所有文章中都没有出现过(词典选的不好,特征选择不好),这种情况就会出现0概率.所以我们需要对条件概率做一点小改动:

P(wt∣Y=ck)=α+Nt,kd∗α+NkP(w_t|Y=c_k)=\frac{\alpha + N_{t,k}}{d*\alpha + N_k}P(wt​∣Y=ck​)=d∗α+Nk​α+Nt,k​​

其中,d表示数据维度为d(有d个特征,每个特征都加α\alphaα,保证概率和为1,α\alphaα需要乘d).当α=1\alpha=1α=1时,叫做Laplace Smoonthing拉普拉斯平滑,当然α\alphaα也可以小于1.

代码实现

To Be Continued

Gaussian Naive Bayes 高斯朴素贝叶斯

高斯朴素贝叶斯,假设P(X=x|Y=c_k)是多元高斯分布.在了解高斯朴素贝叶斯之前,先介绍一下什么是高斯分布,什么是多元高斯分布?

高斯分布

高斯分布又称正态分布，在实际应用中最为广泛。对于单变量x∈(−∞,+∞)x \in(−\infty,+\infty)x∈(−∞,+∞),高斯分布的参数有两个，分别是均值μ∈(−∞,+∞)\mu \in (−\infty,+\infty)μ∈(−∞,+∞)和方差σ2&gt;0\sigma^2&gt;0σ2>0，其概率密度函数为

N(x∣μ,σ2)=12πσexp{−(x−μ)22σ2}N(x|\mu,\sigma^2)=\frac1{\sqrt{2\pi}\sigma}exp\{-\frac{(x-\mu)^2}{2\sigma^2}\}N(x∣μ,σ2)=2π​σ1​exp{−2σ2(x−μ)2​}

多元高斯分布

image

其中,μ\muμ是D维均值向量,∑\sum∑是DxD的协方差矩阵,∣∑∣|\sum|∣∑∣是∑\sum∑的行列式.多元高斯分布的期望是μ\muμ,方差是∑\sum∑

特殊的,如果D个维度之间相互独立,那么多元高斯分布就可以表示成单元高斯分布概率密度函数的连乘积.

高斯朴素贝叶斯

高斯朴素贝叶斯模型是假设条件概率P(X=x|Y=c_k)是多元高斯分布,另一方面,由之前的特征的条件独立性假设,我们就可以通过对每个特征的条件概率建模,每个特征的条件概率N(μt,σt2)N(\mu_t,\sigma_t^2)N(μt​,σt2​)也服从高斯分布.

在ccc类下第i个词对应的高斯分布为:

g(xi;μi,c,σi,c)=1σi,c2πexp{−(xi−μi,c)22σi,c2}g(x_i;\mu_{i,c},\sigma_{i,c})=\frac1{\sigma_{i,c}\sqrt{2\pi}}exp\{-\frac{(x_i-\mu_{i,c})^2}{2\sigma_{i,c}^2}\}g(xi​;μi,c​,σi,c​)=σi,c​2π​1​exp{−2σi,c2​(xi​−μi,c​)2​}

其中,μi,c\mu_{i,c}μi,c​,σi,c2\sigma_{i,c}^2σi,c2​表示c类下第i个特征的均值和方差.

由于特征之间的独立性假设,我们可以得到条件概率:

P(X=x∣Y=c)=∏i=1dg(xi;μi,c,σi,c)P(X=x|Y=c)=\prod_{i=1}^d g(x_i;\mu_{i,c},\sigma_{i,c})P(X=x∣Y=c)=∏i=1d​g(xi​;μi,c​,σi,c​)

一共有d个特征.

高斯朴素贝叶斯变成:

P(Y=ck∣X=x)=P(Y=ck)P(X=x∣Y=ck)∑kP(X=x∣Y=ck)P(Y=ck)∝P(Y=ck)P(X=x∣Y=ck)P(Y=c_k|X=x)=\frac{P(Y=c_k)P(X=x|Y=c_k)}{\sum_k P(X=x|Y=c_k)P(Y=c_k)}\propto P(Y=c_k)P(X=x|Y=c_k)P(Y=ck​∣X=x)=∑k​P(X=x∣Y=ck​)P(Y=ck​)P(Y=ck​)P(X=x∣Y=ck​)​∝P(Y=ck​)P(X=x∣Y=ck​)

y=argmaxckP(Y=ck)P(X=x∣Y=ck)y=arg max_{c_k} P(Y=c_k)P(X=x|Y=c_k)y=argmaxck​​P(Y=ck​)P(X=x∣Y=ck​).

了解了多元高斯分布分布之后,接下来的工作就是对参数进行估计,计算μi,c\mu_{i,c}μi,c​和σi,c\sigma_{i,c}σi,c​.

参数估计

先验概率和之前的估算方法相同,不再描述.主要是对高斯分布的均值和方差的估计,采用的方法仍然是极大似然估计.

均值的估计μi,c\mu_{i,c}μi,c​是在样本类别ccc中,所有XiX_iXi​的平均值;

方差的估计σi,k2\sigma_{i,k}^2σi,k2​为样本类别ccc中所有XiX_iXi​的方差.

对于一个连续的样本值,带入高斯分布,就可以求出它的概率分布了.

所有参数估计完成之后,就可以计算给定样本的条件概率P(X=x∣Y=ck)P(X=x|Y=c_k)P(X=x∣Y=ck​),进而可以计算P(Y=ck)P(X=x∣Y=ck)P(Y=c_k)P(X=x|Y=c_k)P(Y=ck​)P(X=x∣Y=ck​),之后就可以确定样本类别,完成模型预测.

Reference

• 概率分布
• sklearn贝叶斯文档
• 《统计学习方法》，李航
• 朴素贝叶斯理论推导与三种常见模型
• 朴素贝叶斯分类器
• Spam Filtering with Naive Bayes – Which Naive Bayes?
• Text Classification using Naive Bayes
• A Comparison of Event Models for Naive Bayes Text Classification