第四节–朴素贝叶斯(Naive Bayes)法

第四节–朴素贝叶斯(Naive Bayes)法

朴素贝叶斯(Naive Bayes,NB)法是基于贝叶斯定理与特征条件独立假设的分类方法.对于给定的训练数据集,首先基于特征条件独立假设学习输入/输出的联合概率分布;然后基于此模型,对给定的输入x,利用贝叶斯定理求出后验概率最大的输出y.

NB包括以下算法：

• 高斯朴素贝叶斯(Gaussian Naive Bayes)–适用于正态分布
• 伯努利朴素贝叶斯(Bernoulli Naive Bayes)–适用于二项分布
• 多项式朴素贝叶斯(Multinomial Navie Bayes)

朴素贝叶斯法的优缺点：

• 优点：学习和预测的效率高,且易于实现;在数据较少的情况下仍然有效,可以处理分类问题
• 缺点：分类效果不一定很高,特征独立性假设会是朴素贝叶斯变得简单,但是会牺牲一定的分类准确率

一.朴素贝叶斯法的学习与分类

1.基本方法

设输入空间$\mathcal{X}\subseteq {\mathbf{R}}^n$为n维向量的集合,输出空间为类标记集合$y=\{c_1,c_2,\cdots ,c_x\}$,输入为特征向量$x\in \mathcal{X}$,输出为类标记(class label)$y\in \mathcal{Y}$,X是定义在输入空间$\mathcal{X}$上的随机向量,Y是定义在输出空间$\mathcal{Y}$上的随机变量.$P(X,Y)$是X和Y的联合概率分布.训练数据集：
$T=\left\{\left(x_1,y_1\right),\left(x_2,y_2\right),\cdots ,\left(x_N,y_N\right)\right\}$

由$P(X,Y)$独立分布产生

朴素贝叶斯法通过训练数据集学习联合概率分布$P(X,Y)$.具体地,学习以下先验概率分布及条件概率分布.先验概率分布：
$$P\left(Y=c_k\right),k=1,2,\cdots ,K$$

条件概率分布：
$$P\left(X=x|Y=c_k\right)=P\left(X^{(1)}=x^{(1)},\cdots ,X^{(n)}=x^{(n)}|Y=c_k\right),k=1,2,\cdots ,K$$

于是学习到联合概率分布$P(X,Y)$

条件概率分布$P\left(X=x|Y=c_k\right)$有指数级数量的参数,其估计实际是不可行的.事实上,假设$x^{(j)}$可取值有$S_j$个,$j=1,2,\cdots ,n$,Y可取值有K个,那么参数个数$K{\prod }_{j=1}^n{S}_j$

朴素贝叶斯法对条件概率分布作了条件独立性的假设.由于这是一个较强的假设,朴素贝叶斯法也由此得名.具体地,条件独立性假设是：
$$\begin{array}{rl}P\left(X=x|Y=c_k\right)& =P\left(X^{(1)}=x^{(1)},\cdots ,X^{(n)}=x^{(n)}|Y=c_k\right)\\ & =\prod _{j=1}^{n}P\left(X^{(j)}=x^{(j)}|Y=c_k\right)\end{array}$$

朴素贝叶斯法实际上学习到生成数据的机制,所以属于生成模型.条件独立假设等于是说用于分类的特征在类确定的条件下都是条件独立的.这一假设使朴素贝叶斯法变得简单,但有时会牺牲一定的分类准确率

朴素贝叶斯法分类时,对给定的输入x,通过学习到的模型计算后验概率分布$P\left(Y=c_k|X=x\right)$,将后验概率最大的类作为x的类输出,后验概率计算根据贝叶斯定理进行：
$$P\left(Y=c_k|X=x\right)=\frac{P\left(X=x|Y=c_k\right)P\left(Y=c_k\right)}{{\sum }_{k}P\left(X=x|Y=c_k\right)P\left(Y=c_k\right)}$$

将上面两式联合：
$$P\left(Y=c_k|X=x\right)=\frac{P\left(Y=c_k\right){\prod }_{j}P\left(X^{(j)}=x^{(j)}|Y=c_k\right)}{{\sum }_{k}P\left(Y=c_k\right){\prod }_{j}P\left(X^{(j)}=x^{(j)}|Y=c_k\right)},k=1,2,\cdots ,K$$

这是朴素贝叶斯法分类的基本公式,于是朴素贝叶斯分类器可表示为：
$$y=f(x)=\arg \underset{c_k}{\max }\frac{P\left(Y=c_k\right){\prod }_{j}P\left(X^{(j)}=x^{(j)}|Y=c_k\right)}{{\sum }_{k}P\left(Y=c_k\right){\prod }_{j}P\left(X^{(j)}=x^{(j)}|Y=c_k\right)}$$

注意到,在上式中分母对所有$c_k$都是相同的,所以：
$$y=\arg \underset{c_k}{\max }P\left(Y=c_k\right)\prod _{j}P\left(X^{(j)}=x^{(j)}|Y=c_k\right)$$

2.后验概率最大化的含义

朴素贝叶斯法将实例分到后验概率最大的类中,这等价于期望风险最小化.假设选择0-1损失函数：
$$L(Y,f(X))=\{\begin{array}{ll}1,& Y̸=f(X)\\ 0,& Y=f(X)\end{array}$$

式中$f(X)$是分类决策函数.这时期望风险函数为：
$$R_{\mathrm{e}\mathrm{x}\mathrm{p}}(f)=E[L(Y,f(X))]$$

期望是对联合分布$P(X,Y)$取的.由此取条件期望：
$$R_{\mathrm{e}\mathrm{x}\mathrm{p}}(f)=E_X\sum _{k=1}^K\left[L\left(c_k,f(X)\right)\right]P\left(c_k|X\right)$$

为了使期望风险最小化,只需对$X=x$逐个极小化,由此得到：
$$\begin{array}{rl}f(x)& =\arg \underset{y\in y}{\min }\sum _{k=1}^{K}L\left(c_k,y\right)P\left(c_k|X=x\right)\\ & =\arg \underset{y\in y}{\min }\sum _{k=1}^{K}P\left(y̸=c_k|X=x\right)\\ & =\arg \underset{y\in \mathcal{Y}}{\min }\left(1-P\left(y=c_k|X=x\right)\right)\\ & =\arg \underset{y\in y}{\max }P\left(y=c_k|X=x\right)\end{array}$$

这样一来,根据期望风险最小化准则就得到了后验概率最大化准则：
$$f(x)=\arg \underset{c_k}{\max }P\left(c_k|X=x\right)$$

即朴素贝叶斯法所采用的原理

二.朴素贝叶斯法的参数估计

1.极大似然估计

在朴素贝叶斯法中,学习意味着估计$P\left(Y=c_k\right)$和$P\left(X^{(j)}=x^{(j)}|Y=c_k\right)$.可以应用极大似然估计相应的概率.先验概率$P\left(Y=c_k\right)$的极大似然估计是：
$$P\left(Y=c_k\right)=\frac{{\sum }_{i=1}^{N}I\left(y_i=c_k\right)}{N},k=1,2,\cdots ,K$$

设第j个特征$x^{(j)}$可能取值的集合为$\left\{a_{j1},a_{j2},\cdots ,a_{j{S}_j}\right\}$.条件概率$P\left(X^{(j)}=a_{jl}|Y=c_k\right)$的极大似然估计是：

$$P\left(X^{(j)}=a_{jl}|Y=c_k\right)=\frac{{\sum }_{i=1}^{N}I\left(x_i^{(j)}=a_{jl}{y}_i=c_k\right)}{{\sum }_{i=1}^{N}I\left(y_i=c_k\right)}$$

$$j=1,2,\cdots ,n;l=1,2,\cdots ,S_j:k=1,2,\cdots ,K$$

式中,$x_i^{(j)}$是第i个样本的第j个特征;$a_{jl}$是第j个特征可能取的第l个值：I为指示函数

2.学习与分类算法

输入：训练数据$T=\left\{\left(x_1,y_1\right),\left(x_2,y_2\right),\cdots ,\left(x_N,y_N\right)\right\}$,其中$x_i={\left(x_i^{(1)},x_i^{(2)},\cdots ,x_i^{(n)}\right)}^{\mathrm{T}}$,$x_i^{(j)}$是第i个样本的第j个特征,$x_i^{(j)}\in \left\{a_{j1},a_{j2},\cdots ,a_{j{s}_j}\right\}$,$a_{jl}$是第j个特征可能取的第l个值,$j=1,2,\cdots ,n,l=1,2,\cdots ,S_j,y_i\in \left\{c_1,c_2,\cdots ,c_K\right\}$,实例x;

输出：实例x的分类

1. 计算先验概率及条件概率
   $$\begin{array}{l}P\left(Y=c_k\right)=\frac{{\sum }_{i=1}^{N}I\left(y_i=c_k\right)}{N},k=1,2,\cdots ,K\\ P\left(X^{(j)}=a_{jl}|Y=c_k\right)=\frac{{\sum }_{i=1}^{N}I\left(x_i^{(j)}=a_{jl},y_i=c_k\right)}{{\sum }_{i=1}^{N}I\left(y_i=c_k\right)}\\ j=1,2,\cdots ,n;l=1,2,\cdots ,S_j;k=1,2,\cdots ,K\end{array}$$

2. 对于给定的实例$x={\left(x^{(1)},x^{(2)},\cdots ,x^{(n)}\right)}^{\mathrm{T}}$,计算：
   $$P\left(Y=c_k\right)\prod _{j=1}^{n}P\left(X^{(j)}=x^{(j)}|Y=c_k\right),k=1,2,\cdots ,K$$

3. 确定实例x的类：
   $$y=\arg \underset{c_k}{\max }P\left(Y=c_k\right)\prod _{j=1}^{n}P\left(X^{(j)}=x^{(j)}|Y=c_k\right)$$

实例1：通过下表的训练数据学习一个朴素贝叶斯分类器并确定$x=(2,S{)}^T$的类标记y.表中$X^{(1)},X^{(2)}$为特征,取值的集合分别为$A_1=\{1,2,3\},A_2=\{S,M,L\}$,Y为类标记,$Y\in C=\{1,-1\}$

1	1	2	3	4	5	6	7	8	9	10	11	12	13	14	15
$X^{(1)}$	1	1	1	1	1	2	2	2	2	2	3	3	3	3	3
$X^{(2)}$	S	M	M	S	S	S	M	M	L	L	L	M	M	L	L
Y	-1	-1	1	1	-1	-1	-1	1	1	1	1	1	1	1	-1

from IPython.display import Image
Image(filename="./data/4_2.png",width=500)

3.贝叶斯估计

用极大似然估计可能会出现所要估计的概率值为0的情况.这时会影响到后验概率的计算结果.是分类产生偏差.解决这一问题的方法是采用贝叶斯估计,具体地,条件概率的贝叶斯估计是：
$$P_{\lambda }\left(X^{(j)}=a_{jl}|Y=c_k\right)=\frac{{\sum }_{i=1}^{N}I\left(x_i^{(j)}=a_{jl},y_i=c_k\right)+\lambda }{{\sum }_{i=1}^{N}I\left(y_i=c_k\right)+S_j\lambda }$$

式中$\lambda ⩾0$,等价于在随机变量各个取值的频数上赋予一个正数$\lambda >0$.当$\lambda =0$时就是极大似然估计.常取$\lambda =1$,这时称为拉普拉斯平滑(Laplace smoothing).显然对任何$l=1,2,\cdots ,S_j,k=1,2,\cdots ,K$,有：
$$\begin{array}{l}{P}_{\lambda }\left(X^{(j)}=a_{jl}|Y=c_k\right)>0\\ {\sum }_{i=1}^{s_j}P\left(X^{(j)}=a_{jl}|Y=c_k\right)=1\end{array}$$

同样,先验概率的贝叶斯估计是：
$$P_{\lambda }\left(Y=c_k\right)=\frac{{\sum }_{i=1}^{N}I\left(y_i=c_k\right)+\lambda }{N+K\lambda }$$

实例2：对实例1,按照拉普拉斯平滑估计概率,即取$\lambda =1$

Image(filename="./data/4_1.png",width=500)

三.代码实现

%matplotlib inline
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt

from sklearn.datasets import load_iris
from sklearn.model_selection import train_test_split

from collections import Counter
import math

def load_data():
    iris=load_iris()
    df=pd.DataFrame(iris.data,columns=iris.feature_names)
    df["label"]=iris.target
    df.columns=["sepal lenght","sepal width","petal length","petal width","label"]
    data=np.array(df.iloc[:100,:])
    return data[:,:-1],data[:,-1]

X,y=load_data()
X_train,X_test,y_train,y_test=train_test_split(X,y,test_size=0.3)

X_test[0],y_test[0]

(array([4.5, 2.3, 1.3, 0.3]), 0.0)

1.自定义GaussianNB

特征的可能性被假设为高斯

概率密度函数：
$$P\left(x_i|y_k\right)=\frac{1}{\sqrt{2\pi {\sigma }_{yk}^2}}\exp \left(-\frac{{\left(x_i-{\mu }_{yk}\right)}^2}{2{\sigma }_{yk}^2}\right)$$

数学期望(mean)：$\mu$,方差：${\sigma }^2=\frac{\sum (X-\mu {)}^2}{N}$

class NaiveBayes(object):
    def __init__(self):
        self.model=None
        
    # 数学期望
    @staticmethod
    def mean(X):
        return sum(X)/float(len(X))
    
    # 标准差(方差)
    def stdev(self,X):
        avg=self.mean(X)
        return math.sqrt(sum([pow(x-avg,2) for x in X])/float(len(X)))
    
    #概率密度函数
    def gaussian_probability(self,x,mean,stdev):
        exponent=math.exp(-(math.pow(x-mean,2)/(2*math.pow(stdev,2))))
        return (1/(math.sqrt(x*math.pi)*stdev))*exponent
    
    # 处理X_train
    def summarize(self,train_data):
        summaries=[(self.mean(i),self.stdev(i)) for i in zip(*train_data)]
        return summaries
    
    # 分类别求出数学期望和标准差
    def fit(self,X,y):
        labels=list(set(y))
        data={label:[] for label in labels}
        for f,label in zip(X,y):
            data[label].append(f)
        self.model={label:self.summarize(value) for label,value in data.items()}
        return "GaussianNB train done"
    
    # 计算概率
    def calculate_probabilities(self,input_data):
        probabilities={}
        for label,value in self.model.items():
            probabilities[label]=1
            for i in range(len(value)):
                mean,stdev=value[i]
                probabilities[label]*=self.gaussian_probability(input_data[i],mean,stdev)
        return probabilities
    
    # 类别
    def predict(self,X_test):
        label=sorted(self.calculate_probabilities(X_test).items(),key=lambda x:x[-1])[-1][0]
        return label
    
    def score(self,X_test,y_test):
        right=0
        for X,y in zip(X_test,y_test):
            label=self.predict(X)
            if label==y:
                right+=1
        return right/float(len(X_test))
    

model=NaiveBayes()

model.fit(X_train,y_train)

'GaussianNB train done'

print(model.predict([4.4,3.2,1.3,0.2]))

0.0

model.score(X_test,y_test)

1.0

2.sklearn Naive_Bayes

from sklearn.naive_bayes import GaussianNB

clf=GaussianNB()
clf.fit(X_train,y_train)

GaussianNB(priors=None, var_smoothing=1e-09)

clf.score(X_test,y_test)

1.0

clf.predict([[4.4,3.2,1.3,0.2]])

array([0.])