朴素贝叶斯

目录

1.模型概述
2.后验最大化
3.参数估计
  3.1 最大似然估计
  3.2 贝叶斯估计
4.代码实现

模型概述

朴素贝叶斯(naive bayes) 是一种基于贝叶斯定理以及假设特征之间独立的分类算法。属于监督学习的一种。对于此模型,重点是基于两个理论,贝叶斯定理后验最大概率估计
贝叶斯定理:P(Y|X) = \frac{P(X|Y)P(Y)}{P(X)}
最大后验估计:y = {\underset {c_k}{\operatorname {arg\max} }}\ P(c_k|X = x)
给定x条件下,y的条件分布当中找一个使得y出现的概率最大的类别,然后用该类别作为预测结果,在这个预测过程中,分别计算出y等于每一个类别的概率,在求这些概率的时候,分母对于不同的类别取值是不变的,最大化条件概率就是最大化分子的值。P(Y|X) \rightarrow P(X|Y)P(Y)贝叶斯公式里面,分子这一项有一个条件独立性的假设,来拆成N项的乘积(假设有N个特征)P(X|Y) = \prod_{i=1}^{N}P(X^{(1)} = x^{(1)}|Y=c_k)最终通过后验概率最大来进行分类,在这个模型里面,需要估计的参数是这两部分P(X|Y),P(Y),方法有极大似然估计,和贝叶斯估计
极大似然估计里面可能出现一个比较尴尬的问题就是,某个条件的情况为0,训练集上数据比较少,分类的类别又比较多,其中一种类别,训练集没有相对应的实例。
贝叶斯估计为了使得不为0,在所有分子加上一项\lambda。对于贝叶斯估计,对每一个参数给一个先验分布,Dirichlet分布 ,是beta分布在多维上的一个推广。

后验最大化

朴素贝叶斯中利用到了最大后验估计的原理,这里证明y = {\underset {c_k}{\operatorname {arg\max} }}\ P(c_k|X = x) 的等号成立的过程。
假设模型使用的是0-1损失函数L(Y,f(X)) = \begin{cases} 1, & Y \neq f(X)\\ 0, & Y {=} f(X) \end{cases}
分类器的任务是使得期望风险函数最小min E \ L(Y,f(X))

推导过程

参数估计

在模型里面,需要估计的是两部分P(Y = c_k)P(X = x|Y=c_k)

极大似然估计

P(Y = c_i) = \frac{\# \{y_j = c_i\}}{N},i = 1,2,...k
P(X^{(i)} = x^{(i)}|y = c_i) = \frac{\# \{y_j = c_i,X^{(i)} = x^{(i)} \}}{\# \{y_j = c_i \}}

贝叶斯估计

极大似然估计的化可能会出现某些参数的值为0。使用贝叶斯估计可以避免这一情况,但是相应的也会复杂一些。
Y \in \{ c_1,c_2,...,c_k\}
对于每一个类别的概率的参数记作\theta ,\sum_{i = 1}^k \theta = 1


当时,也叫拉普拉斯修正。

代码实现

class NaiveBayes():
    def __init__(self):
        self.lam = lam
    
    def fit(self,X,y):
        #prior 
        self.prior = {}
        counts_y = list(Counter(y).values())
        values = counts_y / np.sum(counts_y)
        keys = list(Counter(y).keys())
        for i in range(len(values)):
            self.prior[keys[i]] = values[i]
                
        #likelihood
        self.likelihood = {}
        for i in range(np.shape(X)[1]):
            x = X[:,i]
            counts_x = list(Counter(x).values())
            values = counts_x / (np.sum(counts_x) 
            keys = list(Counter(x).keys())
            for i in range(len(values)):
                self.likelihood[keys[i]] = values[i]
        
    def predict(self,X):
        posterioris = []
        
        for j in range(len(self.prior)):
            posteriori = 1
            for i in range(len(X)):     
                posteriori *= self.likelihood[X[i]]
            posteriori *= list(self.prior.values())[j]
            posterioris.append(posteriori)
            
        index = posterioris.index(max(posterioris))
        return list(prior.keys())[index]
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 218,607评论 6 507
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 93,239评论 3 395
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 164,960评论 0 355
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 58,750评论 1 294
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 67,764评论 6 392
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,604评论 1 305
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,347评论 3 418
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 39,253评论 0 276
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,702评论 1 315
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,893评论 3 336
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 40,015评论 1 348
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,734评论 5 346
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,352评论 3 330
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,934评论 0 22
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 33,052评论 1 270
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 48,216评论 3 371
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,969评论 2 355

推荐阅读更多精彩内容