大师兄的数据分析学习笔记(十六):分类模型(二)

大师兄的数据分析学习笔记(十五):分类模型(一)
大师兄的数据分析学习笔记(十七):分类模型(三)

二、朴素贝叶斯

1. 回顾概率
  • 概率P(A)是可能性判别的大小,概率值越大事件A越可能发生,反之则越不可能发生。
  • 条件概率P(A|B),在一定条件(事件B发生的情况)下,事件A发生的概率。
  • 联合概率P(A,B),事件A和事件B共同发生的概率。
  • 条件概率联合概率有如下关系:P(A,B) = P(A|B)P(B) = P(B|A)P(A)
  • 全概率公式P(B) = \sum{P(A_i)P(B|A_i)}
  • 综合以上情况,可以得到贝叶斯公式P(B_i|A) = \frac{P(B_i)P(A|B_i)}{\sum{P(B_j)P(A|B_j)}} -> P(B|A) = \frac{P(A|B)P(B)}{P(A)}
  • 贝叶斯公式用来描述两个条件概率之间的关系。
  • 从上面的公式可以看出,贝叶斯公式的分子可以看作联合概率,而分母可以看作全概率
2. 朴素贝叶斯的基本思想
  • 首先朴素贝叶斯Naive Bayes朴素Naive的含义是,数据中的特征相互独立。
  • 假设1000个微信好友中,真实账号比例C_0 = 0.89,虚假账号比例C_1=0.11,且已知特征:
  • F_1:日志数量/注册天数,离散化:F_1<=0.05, 0.05<F_1<0.2, F_1>=0.2
  • F_2:好友数量/注册天数,离散化:F_2<=0.1, 0.1<F_2<0.8, F_2>0.8
  • F_3:是否使用真实头像,离散化: True = 1/False = 0
  • 如果某账号状态:F_1 落入区间1, F_2落入区间1, F_3落入区间0
  • 已知条件:

真实账号落入F_1的比例为0.5:P(F1=1|C_0)=0.5 / 虚假账号落入F_1的比例为0.1:P(F_1|C_1)=0.1
真实账号落入F_2的比例为0.7:P(F1=1|C_0)=0.7 / 虚假账号落入F_1的比例为0.2:P(F_1|C_1)=0.2
真实账号落入F_3的比例为0.2:P(F1=1|C_0)=0.2 / 虚假账号落入F_1的比例为0.2:P(F_1|C_1)=0.9

  • 根据贝叶斯公式:P(C|F_1=1,F_2=1,F_3=0)=\frac{P(F_1=1,F_2=1,F_3=0|C)P(C)}{P(F_1=1,F_2=1,F_3=0)}
  • 由于朴素贝叶斯特征是相互独立的,内部转化后:P(C|F_1=1,F_2=1,F_3=0)=\frac{P(F_1=1|C)P(F_2=1|C)P(F_3=0|C)P(C)}{P(F_1=1,F_2=1,F_3=0|C)P(C)}
  • 将已知条件带入后,可以分别获得真实账号和虚假账号的值:

P(F_1=1|C_0)P(F_2=1|C_0)P(F_3=0|C_0)P(C_0) = 0.5 \times 0.7 \times 0.2 \times 0.89 = 0.0623
P(F_1=1|C_1)P(F_2=1|C_1)P(F_3=0|C_1)P(C_1) = 0.1 \times 0.2 \times 0.9 \times 0.11 = 0.00198

  • 由于真实账号的值大于虚假账号的值,所以更倾向认为账号是真实账号。
3. 拉普拉斯平滑
  • 回到朴素贝叶斯公式:P(C|F_1=1,F_2=1,F_3=0)=\frac{P(F_1=1|C)P(F_2=1|C)P(F_3=0|C)P(C)}{P(F_1=1,F_2=1,F_3=0|C)P(C)}
  • 假设条件概率F_3=0|C_0:0/F_3=0|C_1:0导致P(F_3=0|C_0)=0/P(F_3=0|C_1)=0 并造成整个公式为0。
  • 为了避免这种情况,需要将全部条件概率加1。
4. 代码实现
>>>import os
>>>import pandas as pd
>>>import numpy as np
>>>from sklearn.model_selection import train_test_split
>>>from sklearn.naive_bayes import GaussianNB,BernoulliNB
>>>from sklearn.metrics import  accuracy_score,recall_score,f1_score

>>>models = []
>>>models.append(("GaussianNB",GaussianNB()))
>>>models.append(("BernoulliNB",BernoulliNB()))

>>>df = pd.read_csv(os.path.join(".", "data", "WA_Fn-UseC_-HR-Employee-Attrition.csv"))
>>>X_tt,X_validation,Y_tt,Y_validation = train_test_split(df.JobLevel,df.JobSatisfaction,test_size=0.2)
>>>X_train,X_test,Y_train,Y_test = train_test_split(X_tt,Y_tt,test_size=0.25)

>>>data = df[["JobSatisfaction","JobLevel"]]

>>>for clf_name,clf in models:
>>>    clf.fit(np.array(X_train).reshape(-1,1),np.array(Y_train).reshape(-1,1))
>>>    xy_lst = [(X_train,Y_train),(X_validation,Y_validation),(X_test,Y_test)]
>>>    for i in range(len(xy_lst)):
>>>        X_part = xy_lst[i][0]
>>>        Y_part = xy_lst[i][1]
>>>        Y_pred = clf.predict(np.array(X_part).reshape(-1,1))
>>>        print(i)
>>>        print(clf_name,"-ACC",accuracy_score(Y_part,Y_pred))
>>>        print(clf_name,"-REC",recall_score(Y_part,Y_pred,average='macro'))
>>>        print(clf_name,"-F1",f1_score(Y_part,Y_pred,average='macro'))
>>>       print("="*40)
0
GaussianNB -ACC 0.3253968253968254
GaussianNB -REC 0.25
GaussianNB -F1 0.12275449101796408
========================================
1
GaussianNB -ACC 0.2755102040816326
GaussianNB -REC 0.25
GaussianNB -F1 0.10799999999999998
========================================
2
GaussianNB -ACC 0.30952380952380953
GaussianNB -REC 0.25
GaussianNB -F1 0.11818181818181818
========================================
0
BernoulliNB -ACC 0.3253968253968254
BernoulliNB -REC 0.25
BernoulliNB -F1 0.12275449101796408
========================================
1
BernoulliNB -ACC 0.2755102040816326
BernoulliNB -REC 0.25
BernoulliNB -F1 0.10799999999999998
========================================
2
BernoulliNB -ACC 0.30952380952380953
BernoulliNB -REC 0.25
BernoulliNB -F1 0.11818181818181818
========================================
5. 生成模型与判别模型
  • 生成模型:通过求输入与输出的联合概率分布,再求解类别归类的概率,比如朴素贝叶斯模型。
  • 判别模型:不通过联合概率分布,直接可以获得输出对应最大分类的概率,比如KNN。
  • 生成模型相对判别模型对数据的要求更高,速度也更快。
  • 判别模型相对生成模型对数据的容忍程度更高,使用范围更广。
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 213,752评论 6 493
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 91,100评论 3 387
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 159,244评论 0 349
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 57,099评论 1 286
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 66,210评论 6 385
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,307评论 1 292
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,346评论 3 412
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,133评论 0 269
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,546评论 1 306
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 36,849评论 2 328
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,019评论 1 341
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,702评论 4 337
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,331评论 3 319
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,030评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,260评论 1 267
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 46,871评论 2 365
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 43,898评论 2 351

推荐阅读更多精彩内容