python中如何使用樸素貝葉斯算法

這里再重復(fù)一下標(biāo)題為什么是"使用"而不是"實(shí)現(xiàn)"：
首先，專(zhuān)業(yè)人士提供的算法比我們自己寫(xiě)的算法無(wú)論是效率還是正確率上都要高。
其次，對(duì)于數(shù)學(xué)不好的人來(lái)說(shuō)，為了實(shí)現(xiàn)算法而去研究一堆公式是很痛苦的事情。
再次，除非他人提供的算法滿(mǎn)足不了自己的需求，否則沒(méi)必要"重復(fù)造輪子"。
下面言歸正傳，不了解貝葉斯算法的可以去查一下相關(guān)資料，這里只是簡(jiǎn)單介紹一下：
1.貝葉斯公式：
P(A|B)=P(AB)/P(B)
2.貝葉斯推斷：
P(A|B)=P(A)×P(B|A)/P(B)
用文字表述：
后驗(yàn)概率=先驗(yàn)概率×相似度/標(biāo)準(zhǔn)化常量
而貝葉斯算法要解決的問(wèn)題就是如何求出相似度，即：P(B|A)的值
3. 在scikit-learn包中提供了三種常用的樸素貝葉斯算法，下面依次說(shuō)明：
1）高斯樸素貝葉斯：假設(shè)屬性/特征是服從正態(tài)分布的(如下圖)，主要應(yīng)用于數(shù)值型特征。
使用scikit-learn包中自帶的數(shù)據(jù)，代碼及說(shuō)明如下：    
>>>from sklearn import datasets ##導(dǎo)入包中的數(shù)據(jù)
>>> iris=datasets.load_iris() ##加載數(shù)據(jù)
>>> iris.feature_names  ##顯示特征名字
 ['sepal length (cm)', 'sepal width (cm)', 'petal length (cm)', 'petal width (cm)']
>>> iris.data   ##顯示數(shù)據(jù)
 array([[ 5.1, 3.5, 1.4, 0.2],[ 4.9, 3. , 1.4, 0.2],[ 4.7, 3.2, 1.3, 0.2]............
>>> iris.data.size  ##數(shù)據(jù)大小 ---600個(gè)
>>> iris.target_names  ##顯示分類(lèi)的名字
 array(['setosa', 'versicolor', 'virginica'], dtype='<U10')
>>> from sklearn.naive_bayes import GaussianNB ##導(dǎo)入高斯樸素貝葉斯算法
>>> clf = GaussianNB()    ##給算法賦一個(gè)變量，主要是為了方便使用
>>> clf.fit(iris.data, iris.target)  ##開(kāi)始分類(lèi)。對(duì)于量特別大的樣本，可以使用函數(shù)partial_fit分類(lèi)，避免一次加載過(guò)多數(shù)據(jù)到內(nèi)存
 
>>> clf.predict(iris.data[0].reshape(1,-1)) ##驗(yàn)證分類(lèi)。標(biāo)紅部分特別說(shuō)明：因?yàn)閜redict的參數(shù)是數(shù)組，data[0]是列表，所以需要轉(zhuǎn)換一下
array([0])
>>> data=np.array([6,4,6,2])   ##驗(yàn)證分類(lèi)
>>> clf.predict(data.reshape(1,-1))
array([2])
這里涉及到一個(gè)問(wèn)題：如何判斷數(shù)據(jù)符合正態(tài)分布？ R語(yǔ)言里面有相關(guān)函數(shù)判斷，或者直接繪圖也可以看出來(lái)，但是都是P(x,y)這種可以在坐標(biāo)系里面直接
畫(huà)出來(lái)的情況，而例子中的數(shù)據(jù)如何確定，目前還沒(méi)有搞明白，這部分后續(xù)會(huì)補(bǔ)上。
2）多項(xiàng)式分布樸素貝葉斯：常用于文本分類(lèi)，特征是單詞，值是單詞出現(xiàn)的次數(shù)。
  
##示例來(lái)在官方文檔，詳細(xì)說(shuō)明見(jiàn)第一個(gè)例子
>>> import numpy as np
>>> X = np.random.randint(5, size=(6, 100)) ##返回隨機(jī)整數(shù)值：范圍[0,5) 大小6*100 6行100列
>>> y = np.array([1, 2, 3, 4, 5, 6])
>>> from sklearn.naive_bayes import MultinomialNB
>>> clf = MultinomialNB()
>>> clf.fit(X, y)
MultinomialNB(alpha=1.0, class_prior=None, fit_prior=True)
>>> print(clf.predict(X[2]))
[3]

3）伯努力樸素貝葉斯：每個(gè)特征都是是布爾型，得出的結(jié)果是0或1，即出現(xiàn)沒(méi)出現(xiàn)    
##示例來(lái)在官方文檔，詳細(xì)說(shuō)明見(jiàn)第一個(gè)例子
>>> import numpy as np
>>> X = np.random.randint(2, size=(6, 100))
>>> Y = np.array([1, 2, 3, 4, 4, 5])
>>> from sklearn.naive_bayes import BernoulliNB
>>> clf = BernoulliNB()
>>> clf.fit(X, Y)
BernoulliNB(alpha=1.0, binarize=0.0, class_prior=None, fit_prior=True)
>>> print(clf.predict(X[2]))
[3]
補(bǔ)充說(shuō)明：此文還不完善，示例一中也有部分說(shuō)明需要寫(xiě)，最近事情較多，后續(xù)會(huì)逐漸完善。