数学建模社区-数学中国

标题: 机器学习之朴素贝叶斯模型及代码示例 [打印本页]

---

作者: zhangtt123    时间: 2020-5-21 10:21
标题: 机器学习之朴素贝叶斯模型及代码示例
一、朴素贝叶斯的推导
朴素贝叶斯学习（naive Bayes）是一种有监督的学习，训练时不仅要提供训练样本的特征向量X，而且还需提供训练样本的实际标记Y，是一种基于贝叶斯定理和特征条件独立假设的分类方法。
/ ]0 Z) |4 z$ z
1. 贝叶斯定理：
' g6 G/ E/ n3 M$ Q  X
- S# Q1 N8 [, `( J; M5 v贝叶斯定理：


0 V0 \% q& V" \2 V8 c对于分类问题，其中        可看作 在样本的特征为X的条件下，样本的类别为Y的条件概率，这正是分类问题中我们想求的；

而右边中的 可看作 训练集中样本类别为Y的概率 ， 可看作 在训练集中特征为X时，样本类别为Y的条件概率， 可看作 训练集中特征为X的样本概率。可看出，右边的值我们是可以通过计算得到的。其中  、 分别为 X和Y的先验概率，其值与训练集的选择有一定的关系。

2. 特征条件独立假设
7 F7 F3 U) Y4 L1 R* R由 1 可知，特征向量 X 的维度不一定是一维的，可能是多维的： 。

2 l3 X' F7 D% K9 T; Y: @因此 样本类别Y 的取值 是跟 样本的每一个维度取值有关的。因此可由贝叶斯定理得出下式：


0 u! c" Z9 e9 P) k( q0 u
+ k1 H  m( B  Y) ^又因为 朴素贝叶斯 对条件概率分布进行了条件独立性的假设，即同一类别中，样本的每一维度的特征都是 独立的 。朴素贝叶斯之所以“朴素”正因为这一假设。因此可有：

- z; F- U' x$ s6 T
- N" Q& a, B1 H: {! c8 [3 `
于是可得：


& [; F  D. J  Z% ~  X
% p8 ?- \5 w8 t" ]6 a# j又由 全概率分布公式 可得：

9 g, b- ~: k" R1 S2 Y5 d
  ~, ~3 Y8 T; U; P7 `1 s
因此我们想要得到的 样本类别y 为：


! x3 H; E' H  `& t& R4 u0 J( q
: C2 ~$ w# n! }3 B1 @即 将概率最高的那个标记 作为预测样本的标签。
2 i. `+ F8 @6 P/ Q6 E, q
+ [* M+ Q( F; \: d又因对于 每一个类别标记 来说：

; \* E; K& I6 t* Y$ W% s 为一常数。

因此 概率最高的样本类别y 可简化表示为：
4 x' \4 X" i8 ?- n7 Z5 s, q% E

1 E+ |8 E. F; X7 n, @
, d$ ~! |+ M( g* U因此 、 的求值是关键。

二、朴素贝叶斯常用模型
2 h+ q  `* Y* K2 P. Y在不同的朴素贝叶斯模型中，  的求值也不同。下列为朴素贝叶斯常见的三种模型。

$ v: C3 n# l4 s$ \  ^* S& w1. 高斯朴素贝叶斯模型
+ p/ \. d) R5 f7 A- w! ~在高斯朴素贝叶斯模型中，特征向量 X 的特征 通常为 连续型变量，并且假定所有特征的取值是符合高斯分布的，即：
" h# R7 Y% i$ G
$ f. }' c3 J9 L+ }7 z1 v( {$ g 。

2 A2 g# W) s9 V其中 参数 、 可通过 极大似然估计（MLE） 求得：

1 M- v% b7 E% R$ M# M
$ ^! `/ ^& k$ M8 C+ C


其中  为 样本类别为 的样本数量， 为 类别为 的第 j 个样本的第 i 个特征的值。
: y- v, G$ n/ ?
由上可取得 概率最高的样本类别y ：



2. 多项分布朴素贝叶斯模型
! P  o  @; }3 D- B) j在多项分布朴素贝叶斯模型中，特征向量 X 的特征 通常为 离散型变量，并且假定所有特征的取值是符合多项分布的，可用于文本分类。对每一个样本类别 ，有对应的参数向量

& M8 b6 }; P2 d- ^  |; S 。
7 J, J! Y9 j6 T- `* ]' @
其中 n 为样本特征的维度，且令 。
2 R& Y0 u: Q: [: ~* T
% C0 b" L& D" W. k其中 参数  同样可以通过 极大似然估计 （MLE）求得：
; ]$ x( v! y0 {+ ?( ^; k4 F& P
( R3 B0 |' s3 B3 p, T% X
7 v! `: u% X  z) t2 g8 f/ m$ K
其中 为 类别为的样本的所有特征出现的总次数，可看作 类别为的样本的所有特征值的总和。

其中  为 类别为的样本的第 i 个特征出现的总次数，可看作 类别为的样本的第 i 个特征值的总和。

又因为为了避免训练集样本对一些特征的缺失，即某一些特征出现的次数为0（特征值的总和为0），从而导致 为 0，进而导致预测的概率为 ，而且甚至会出现 0 除以 0 这种情况。
* D0 j! u6 [/ @( y# ]$ ?
因为我们需要对 参数 进行平滑处理，可有：

+ ]6 A- h: b  i4 R9 A
* Z2 O* l" E( T  j
其中 n 为训练集样本特征的维度， 为平滑系数 ，当 = 1时，称为 Laplace平滑；当  < 1时，成为 Lidstone平滑。
: O7 p9 I* `+ d0 t. F% r
2 J! P- T, p$ d7 n. R- }我们可发现，平滑后的公式 。
% D. L, y- z8 v/ ~
* ~; m1 S- _9 t, J7 R对于其对应的类别  来说，还是保持着 。
! V) I) R9 B+ W- \7 N9 o
而 = 类别为的样本的所有特征的特征值总和 / 训练集所有样本的所有特征的特征值总和 。
0 _* l8 a' c  y: }7 g3 R" ]9 C
& f& S% v6 u2 I- l7 O. i以下为 使用带有Laplace平滑的多项朴素贝叶斯进行文本分类的一个例子：

5 a& f% r4 z) K3 u6 A* {1 ~如下所示：已知在训练集中有两个文本及其标记，下面进行多项朴素贝叶斯模型预测第三个文本的标记。
' E' ?0 L9 f+ b% r
index     X                        Y
  1       Yes，I like it.          1
# a# u& h8 O; f  2       No，I dislike it.        0
$ I, |- k! g: d9 q  V
) w3 @% J1 j0 r7 s2 t4 Z  3       No,No.dislike it.       ?
3 L: E" k, n! Q$ {" W3 C( ~" W————————————————
7 Y  S2 a- R% ?1 C! v: P
2 Y7 v+ l" g' @9 M  Z$ q若使用词袋法，且 以 训练集中的文本 为词汇表，即将训练集中的文本中出现的单词都统计出来作为词典，那么记单词的数目为 n，这代表了文本的n个维度。通过统计，可得 n = 6 。
; C/ w4 T/ G% J2 e: Y* N! v8 ?/ H0 |
以上三个文本在这 6 个特征维度上的表示为：
: T2 Y! ~# `: e1 `. k
[backcolor=rgba(0, 0, 0, 0.05)]    Yes  No  I  like  dislike  it
8 a) j( a5 J1 j1 b[backcolor=rgba(0, 0, 0, 0.05)]X1   1   0   1    1      0     1X2   0   1   1    0      1     1 X3   0   2   0    0      1     1下面需要计算：
4 O2 f4 D) m" V/ p4 z: k
4 u* L5 H2 u* f2 q$ J多项分布朴素贝叶斯是以文本的单词为粒度进行计算的，因此通过训练集的学习可有：

P(Y=1) = 4/8 = 1/2 、 P(Y=0) = 4/8 = 1/2
0 C, Z! P% O$ v) K' }5 Q( _4 S
- G6 n) d9 i- Z5 e4 ?% ~8 v由Laplace平滑公式计算可得：

P(No|Y=1) = (0+1)/(4+6) = 1/10
3 x: _/ b4 b+ p0 @1 JP(dislike|Y=1) = (0+1)/(4+6) = 1/10
* c# E9 g4 I* J: A- G- f& B/ jP(it|Y=1) = (1+1)/(4+6) = 1/5
------------------------------------
P(No|Y=0) = (1+1)/(4+6) = 1/5
P(dislike|Y=0) = (1+1)/(4+6) = 1/5
) ?2 p  \% z5 iP(it|Y=0) = (1+1)/(4+6) = 1/5

6 y  q# u( ^( b# a所以：可计算出X3的样本标记Y=1和Y=0的相对概率：
* y. {2 N: U/ {  B
P(Y=1|X3) = P(Y=1)*P(No|Y=1)*P(No|Y=1)*P(dislike|Y=1)*P(it|Y=1)
; H7 w; i" `/ R' p+ w* d9 d= 1/2*1/10*1/10*1/10*1/5 = 1/10000  
------------------------------------------
. E% M/ v% ^% d4 b& @' qP(Y=0|X3) = P(Y=0)*P(No|Y=0)*P(No|Y=0)*P(dislike|Y=0)*P(it|Y=0)
    = 1/2*1/5*1/5*1/5*1/5 = 1/1250  

0 M- B/ ~4 \; D  S  为什么 P(Y=1|X3) + P(Y=0|X3)不等于1 ?
" |7 @2 U+ k6 g- L6 X' V: ?9 Q2 Q9 {
4 C( F4 J4 V  I2 ~* F* p2 g> 因为计算时，分母P(x1,x2,...,xn)已经被省略掉.
  q- i: L( r4 n* L& O) v  S
由公式  ：
5 R! ]+ y; N5 `
, Q: a( m/ I: f1 ?8 i$ s0 c1 l可知，X3的类别标记应该为 Y=0 。以上即完成了对X3的类别预测。
! I  ^1 m; K/ X& U! i
! |  s- ~! W3 q" K' t词汇表的选择：当预测样本的文本中出现训练集文本从未出现的新词时，这时会出现特征维度的缺失，这时可以以 现成的单词词典 作为词汇表，当然这时样本的维度会提高很多。

停用词的过滤：在文本当中，可以发现一些单词进行词袋法处理后，对文本的预测来说没有太多意义，比如上述的 it 、I等一些代词。我们可以 在词汇表和预测样本中 过滤掉这些单词来简化模型的计算过程。
. ]6 e, L! g# y7 O* _( {# R& j
/ Z7 o! K# n: y9 J+ `3. 伯努利朴素贝叶斯模型
: d% G# b5 X3 k5 S0 P: R' s( L在伯努利朴素贝叶斯模型中，每个特征的取值是布尔型，或以0和1表示，所以伯努利模型中，每个特征值为0或者1。
! L+ @: @2 ]( ?1 o
在文本分类中，多项分布朴素贝叶斯是以文本的单词为粒度（以特征为粒度）进行计算的，即计算该单词特征在对应的文档中出现的次数，而伯努利朴素贝叶斯是为文本为粒度（以样本为粒度）进行计算的，即计算 存在该单词特征的文档个数。因此伯努利朴素贝叶斯模型在一定程度上忽略了同一个文档中的单词的 词频。
/ e3 w+ x5 n) i0 G+ X+ ]; K8 w; d
& U. z% v5 l  [; ~; L其中 带有平滑的  的计算如下：

.

其中  为 每一个 类别为  的样本的所有特征的特征值的总和。

! \9 u8 [. \% U" ?- i其中  为 在类别为 的样本集中，第 i 个特征的值不为0的样本个数。
" k: l0 R9 I+ I
+ z8 N! O$ w. i4 X9 o同样地： n 为训练集样本特征的维度， 为平滑系数 。
- R) A+ r+ y+ K! a
) ?2 I, c8 G. [4 m( K/ p而  = 类别为的样本的个数 / 训练集所有样本的个数。

8 @/ V* D# ?6 i+ X7 T同样，以 使用带有Laplace平滑的伯努利朴素贝叶斯进行文本分类 为例：
  A  ^& L% l" t- ~% }9 s2 W8 `7 a
如下所示：已知在训练集中有三个文本及其标记，下面进行伯努利朴素贝叶斯模型预测第四个文本的标记。

. |4 q9 K  x; Y. ]9 findex     X                        Y
  1       Yes,Yes,I like it.       1
  2       No,No,I dislike it.      0
  3       No,No.Dislike            0
  e: r- p9 U) t: N( F6 C
  4       No.dislike it.           ?

若使用词袋法，且 以 训练集中的文本 为词汇表，即将训练集中的文本中出现的单词都统计出来作为词典，那么记单词的数目为 n，这代表了文本的n个维度。通过统计，可得 n = 6 。

; y3 v1 n  U# U% a/ e以上四个文本在伯努利贝叶斯模型的这 6 个特征维度上的表示为（特征值为1或0）：
" D# U+ F6 k& |
1 k3 \3 ]' ]5 p8 }& L& J    Yes  No  I  like  dislike  it
9 u' w2 |) z5 f  r6 k# eX1   1   0   1    1      0     1

X2   0   1   1    0      1     1
X3   0   1   0    0      1     0

) U: M5 |9 x5 |& oX4   0   1   0    0      1     1

$ m9 N$ R, F, c: S7 i" j0 O下面需要计算：

伯努利分布朴素贝叶斯是以文本为粒度进行计算的，因此通过训练集的学习可有：

P(Y=1) = 1/3 、 P(Y=0) = 2/3

* T1 Y! D5 l- X- z* D6 r0 ~由Laplace平滑公式计算可得：
. T0 B3 v8 t( S0 N7 @' O  H% A$ l5 d
P(No|Y=1) = (0+1)/(4+6) = 1/10
P(dislike|Y=1) = (0+1)/(4+6) = 1/10
P(it|Y=1) = (1+1)/(4+6) = 1/5
8 x/ ^* O: B& W% h6 ]------------------------------------
- o+ a9 d) k/ D1 IP(No|Y=0) = (2+1)/(6+6) = 1/4
P(dislike|Y=0) = (2+1)/(6+6) = 1/4
P(it|Y=0) = (1+1)/(6+6) = 1/6
3 s1 l1 W( N, x, b4 f# V! t
+ e2 o. x3 a6 x7 y- C所以：可计算出X4的样本标记Y=1和Y=0的相对概率：
0 G  W- ^4 D4 S7 U" s7 j* T2 F* X
) b, ^9 i  X9 t& n2 H. MP(Y=1|X4) = P(Y=1)*P(No|Y=1)*P(dislike|Y=1)*P(it|Y=1)
= 1/3*1/10*1/10*1/5 = 1/1500  
------------------------------------------
3 y0 m0 ~2 x  _+ w" U* z5 U% oP(Y=0|X4) = P(Y=0)*P(No|Y=0)*P(dislike|Y=0)*P(it|Y=0)
    = 2/3*1/4*1/4*1/6 = 1/144

由公式  ：

可知，X4的类别标记应该为 Y=0 。以上即完成了对X4的类别预测。

7 ^. O( Q& @+ x. G3 f1 Q6 e三、词袋法的特征值计算
% |  i2 O, F: o  E0 `词袋法对文本特征抽取：不考虑单词出现的顺序，只将每个出现过的单词视作一列特征，称这些不重复的词汇集合为词汇表。因此每一条训练集的文本都可以映射成一个特征向量。文本的特征值的常见计算方式有以下两种：
; x5 A! d. c. F
1. TF计算
TF计算（词频计算）：对于每一个样本文本，只考虑每个单词在当前这条文本中出现的频率（Term Frequency），为避免词频之间差距较大，使用归一化后的词频，通常用在多项分布朴素贝叶斯模型中。具体的公式如下：
9 L% ~$ I' i& `% u; @
.

其中 为 单词A 在当前文本中出现的次数；N 为 当前文本的单词个数。

例如：

& z& U( v# T: H: d: O% d! B$ r0 _    index     X                        
" X" f' B. I! u5 j2 D9 d2 M) \    1         Yes,Yes,I like it.      
$ O9 S' o- _/ ^7 m    2         No,No,I dislike it.      
7 A: g7 b. \* ^
通过词频的计算和倍数乘积可得到向量：
' n5 r% a( ^4 w1 f3 q
7 N6 s5 A# M" `6 ^5 H9 [    Yes  No  I  like  dislike  it
X1   2   0   1    1      0     1
* J( u/ r1 ]6 f! ^- P/ y. tX2   0   2   1    0      1     1
. P' W4 _  m3 M, `3 m
2. TF-IDF计算
TF-IDF计算：除了计算某一个单词在当前文本中出现的频率（TF），而且还考虑 含有这个单词的文本条数 的倒数（Inverse Document Frequency），通常用在 多项分布朴素贝叶斯模型 中。当一个单词几乎在所有文本中都出现时，这个单词会对我们的分类决策带有干扰。TF-IDF的计算具体如下：
- d/ }6 ]( U' e+ {: T: c4 h

- T8 h1 ~( e6 }
其中  、N 的含义同上；D 表示 所有的文本数目，可理解为所有的样本数；表示 包含单词A的文本数目。
1 n7 g7 b/ P, i, e2 b
1 Z" @# V! I+ I! v- w8 F为了考虑  项为0情况，在分母进行了加 1 操作。
- a' e" r6 s9 ~  `) n
越大，对应的特征值 就越小，对分类决策的影响就越小。

3. 二值计算
$ f& o" A" w2 z1 L: T二值计算：直接使用二值特征来表示：若单词出现在文本中则特征值为1，不出现则特征值为0，通常用在伯努利朴素贝叶斯模型
( k/ m1 t1 Q6 y1 B% K7 Q! e中。

四、代码示例
在去掉停用词的情况下，使用TF计算、TF-IDF计算进行多项式朴素贝叶斯模型的学习，对20类新闻文本数据分类；使用高斯朴素贝叶斯模型进行花的分类。

from sklearn.datasets import fetch_20newsgroups
news = fetch_20newsgroups(subset="all") # 下载近2万条文本
* _7 z4 L. j5 `! u8 A5 N5 N7 Z% X4 }7 q# H' Gprint news.target_names # 文本的20种类别
/ Y1 A( c. s( t( T  ?' g# w
['alt.atheism', 'comp.graphics', 'comp.os.ms-windows.misc', 'comp.sys.ibm.pc.hardware', 'comp.sys.mac.hardware', 'comp.windows.x', 'misc.forsale', 'rec.autos', 'rec.motorcycles', 'rec.sport.baseball', 'rec.sport.hockey', 'sci.crypt', 'sci.electronics', 'sci.med', 'sci.space', 'soc.religion.christian', 'talk.politics.guns', 'talk.politics.mideast', 'talk.politics.misc', 'talk.religion.misc']
9 M- m# N: P! h5 A/ S8 Q
: E7 x& @/ F: c5 ?% N: ?from sklearn.model_selection import train_test_split
/ X! h0 o- |5 b" _X_train,X_test,y_train,y_test = train_test_split(news.data,news.target,test_size=1/4.,random_state=38)
' u/ \& C% T- Sfrom sklearn.feature_extraction.text import CountVectorizer
5 }7 v& P' C9 x9 w& r/ j
4 A# Y& ~; x- G# 过滤停用词
count_vec = CountVectorizer(analyzer='word',stop_words='english')
from sklearn.feature_extraction.text import TfidfVectorizer
1 d9 o/ S: }) t' \% ptfidf_vec = TfidfVectorizer(analyzer='word',stop_words='english')

  f. R  j' @* R9 {$ b# 进行TF特征抽取
X_count_train = count_vec.fit_transform(X_train)
X_count_test = count_vec.transform(X_test)
* X$ f# k) h, D
# 进行TF-IDF特征抽取
- A2 N0 P1 n* e. F2 G8 L3 EX_tfidf_train = tfidf_vec.fit_transform(X_train)
8 v3 u3 C7 ~, |6 e& gX_tfidf_test = tfidf_vec.transform(X_test)
, B4 @8 D* Q7 N, B' \8 }- f
# 使用多项分布朴素贝叶斯
) ~# P" t& ~2 p. h2 d1 efrom sklearn.naive_bayes import MultinomialNB

( p3 h6 x4 M# Q- A# O( w7 Z# alpha :平滑因子
mnb_count_clf = MultinomialNB(alpha=1.)
( o" A; L2 e4 K0 a/ W3 e- D  \mnb_count_clf.fit(X_count_train,y_train)
/ b3 o3 }" j2 }& O# _; G" E" M& nmnb_count_clf.score(X_count_test,y_test)
0.86311544991511036
9 c" w5 Q% _# Z: c( h1 R( r( h0 W( {
8 n% B3 w9 }0 P( }: U. Q; a9 S# alpha :平滑因子
+ U: H% p; A/ q( Z6 \0 }8 G9 zmnb_tfidf_clf = MultinomialNB(alpha=1.)  
8 W' f! ]5 b/ l0 F. Jmnb_tfidf_clf.fit(X_tfidf_train,y_train)
; _& e# }: A9 H5 T% e" mmnb_tfidf_clf.score(X_tfidf_test,y_test)
0.86884550084889645

# 使用伯努利朴素贝叶斯
from sklearn.naive_bayes import BernoulliNB
* r3 n" u% F/ Z  l) qbnb_count_clf = BernoulliNB(alpha=1.)
bnb_count_clf.fit(X_count_train,y_train)
bnb_count_clf.score(X_count_test,y_test)
0.72325976230899836

bnb_tfidf_clf = BernoulliNB(alpha=1.)
bnb_tfidf_clf.fit(X_tfidf_train,y_train)
bnb_tfidf_clf.score(X_tfidf_test,y_test)
- h/ o* N3 m, e, e: p- |! m! B0.72325976230899836
＞　可以发现伯努利朴素贝叶斯模型的学习中特征值只有 非0和0 的区别
& A" T* [. e1 ~$ t7 n" J
. v* o2 p5 d/ O9 M; n3 jfrom sklearn.naive_bayes import GaussianNB
#　使用高斯朴素贝叶斯模型进行花的分类
from sklearn.datasets import load_iris
, s/ N* f5 T' |) }. biris = load_iris()
& y4 T7 J4 T. J, k4 _X_train,X_test,y_train,y_test = train_test_split(iris.data,iris.target,test_size=1/4.,random_state=38)
from sklearn.preprocessing import StandardScaler
' I  d+ j3 F7 h8 u6 bstandard_scaler = StandardScaler();
. e: A* B! W$ u+ ?7 X) CX_train = standard_scaler.fit_transform(X_train)
. F3 m6 B" `' ?% r4 nX_test = standard_scaler.transform(X_test)
  k2 n* K- f# @8 ~- ngnb_clf = GaussianNB()
# y, T) m& K* p) W+ Z& ?gnb_clf.fit(X_train,y_train)
* O( n* F. ]9 I- R6 L1 c6 H2 bgnb_clf.score(X_test,y_test)
————————————————
+ f. P" J$ R" Q! y$ }6 a0 r版权声明：本文为CSDN博主「cxmscb」的原创文章，遵循CC 4.0 BY-SA版权协议，转载请附上原文出处链接及本声明。
: i, r$ X7 a5 N4 j- J* O3 ~原文链接：https://blog.csdn.net/cxmscb/article/details/69267326
* {, R# v9 U' {3 v- M8 N

& R' o9 D3 P  z! b' b: A

---

欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)

Powered by Discuz! X2.5