位置: IT常识 - 正文

机器学习--使用朴素贝叶斯进行垃圾邮件分类

编辑:rootadmin
机器学习--使用朴素贝叶斯进行垃圾邮件分类 一、学习背景垃圾邮件的问题一直困扰着人们,传统的垃圾邮件分类的方法主要有"关键词法"和"校验码法"等,然而这两种方法效果并不理想。其中,如果使用的是“关键词”法,垃圾邮件中如果这个关键词被拆开则可能识别不了,比如,“中奖”如果被拆成“中 --- 奖”可能会识别不了。后来,直到提出了使用“贝叶斯”的方法才使得垃圾邮件的分类达到一个较好的效果,而且随着邮件数目越来越多,贝叶斯分类的效果会更加好。我们想采用的分类方法是通过多个词来判断是否为垃圾邮件,但这个概率难以估计,通过贝叶斯公式,可以转化为求垃圾邮件中这些词出现的概率。二、贝叶斯公式贝叶斯定理 由英国数学家贝叶斯 ( Thomas Bayes 1702-1761 ) 发展,用来描述两个条件 概率 之间的关系

已知两个独立事件A和B,事件B发生的前提下,事件A发生的概率可以表示为P(A|B),即上图中橙色部分占红色部分的比例,即:

推荐整理分享机器学习--使用朴素贝叶斯进行垃圾邮件分类,希望有所帮助,仅作参考,欢迎阅读内容。

文章相关热门搜索词:,内容如对您有帮助,希望把文章链接给更多的朋友!

P(A) 称为”先验概率”,即在B事件发生之前,我们对A事件概率的一个判断。如:正常收到一封邮件,该邮件为垃圾邮件的概率就是“先验概率”P(A|B)称为”后验概率”, 即在B事件发生之后,我们对A事件概率的重新评估。如:邮件中含有“中奖”这个词,该邮件为垃圾邮件的概率就是“后验概率”P(B|A)/P(B)是可能性函数,这是一个调整因子,使得预估概率更接近真实概率。条件概率就是:后验概率 = 先验概率 x 调整因子因为要计算两次概率,关于它们的分母,是这个样本的属性在全部样本中的概率。而这两次计算,它们的分母是不变的,所以我们只要计算分子就行。于是有了下面的结论:即:

        朴素贝叶斯分类器(Naïve Bayes Classifier)采用了“属性条件独立性假设”,即每个属性独立地对分类结果发生影响。为方便公式标记,不妨记P(C=c|X=x)为P(c|x)。在假设每个属性都独立的情况下,贝叶斯公式可以修改为:分母是相同的,于是去掉分母,得:

最终我们利用这个公式,在代码中实现概率的计算来对样本进行分类。三、使用朴素贝叶斯进行垃圾邮件分类1.算法思路分类标准:当 P(垃圾邮件|文字内容)> P(正常邮件|文字内容)时,我们认为该邮件为垃圾邮件,但是单凭单个词而做出判断误差肯定相当大,因此我们可以将所有的词一起进行联合判断。假设我们进行判断的词语有“获奖”、“贷款”、“无利息”,则需要判断P(垃圾邮件|获奖,贷款,无利息)与P(正常|获奖,贷款,无利息),使用贝叶斯公式,P(垃圾邮件|获奖,贷款,无利息)可以变为:

假设所有词语独立同分布,可以得到机器学习--使用朴素贝叶斯进行垃圾邮件分类

同理可得P(正常|获奖,贷款,无利息)因此,对P(垃圾邮件|获奖,贷款,无利息)与P(正常|获奖,贷款,无利息)的比较,只需要对分子进行对比。但是如果对多个词的P(内容|正常/垃圾)进行相乘时,可能会因为某个词的概率很小从而导致最后的结果为0(超出计算机的精度),因此可以对P(内容|正常/垃圾)取自然对数,即ln P(内容|正常/垃圾)。因此可以变为

 

2.数据集准备数据来源于国外的一些垃圾邮件和正常邮件

   

        其中ham是正常邮件,spam是垃圾邮件。 3.数据导入import osimport reimport stringimport mathDATA_DIR = 'enron'target_names = ['ham', 'spam']def get_data(DATA_DIR): subfolders = ['enron%d' % i for i in range(1,7)] #获得enron下面的文件夹 data = [] target = [] for subfolder in subfolders: #垃圾邮件 spam spam_files = os.listdir(os.path.join(DATA_DIR, subfolder, 'spam')) #将文件夹路径进行组合 for spam_file in spam_files: #遍历所有垃圾文件 with open(os.path.join(DATA_DIR, subfolder, 'spam', spam_file), encoding="latin-1") as f: data.append(f.read()) target.append(1) #正常邮件 pam ham_files = os.listdir(os.path.join(DATA_DIR, subfolder, 'ham')) for ham_file in ham_files: with open(os.path.join(DATA_DIR, subfolder, 'ham', ham_file), encoding="latin-1") as f: data.append(f.read()) target.append(0) return data, targetX, y = get_data(DATA_DIR)#print(X,y)邮件内容存储在data中,标签存储在target当中,“1”表示为垃圾邮件,“0”表示为正常邮件。下图为截取部分结果。

4.分析数据并进行算法训练定义一个类对数据进行预处理class SpamDetector_1(object): #清除空格 def clean(self, s): translator = str.maketrans("", "", string.punctuation) return s.translate(translator) #分开每个单词 def tokenize(self, text): text = self.clean(text).lower() return re.split("\W+", text) #计算某个单词出现的次数 def get_word_counts(self, words): word_counts = {} for word in words: word_counts[word] = word_counts.get(word, 0.0) + 1.0 return word_counts计算P(垃圾邮件)和P(正常邮件);词汇表(即正常邮件和垃圾邮件中出现的所有单词,方便进行拉普拉斯平滑);垃圾邮件和非垃圾邮件的词频,即给定词在垃圾邮件和非垃圾邮件中出现的次数。选取了第100封之后的邮件作为训练集,前面一百封邮件作为测试集。class SpamDetector_2(SpamDetector_1): # X:data,Y:target标签(垃圾邮件或正常邮件) def fit(self, X, Y): self.num_messages = {} self.log_class_priors = {} self.word_counts = {} # 建立一个集合存储所有出现的单词 self.vocab = set() # 统计spam和ham邮件的个数 self.num_messages['spam'] = sum(1 for label in Y if label == 1) self.num_messages['ham'] = sum(1 for label in Y if label == 0) # 计算先验概率,即所有的邮件中,垃圾邮件和正常邮件所占的比例 self.log_class_priors['spam'] = math.log( self.num_messages['spam'] / (self.num_messages['spam'] + self.num_messages['ham'])) self.log_class_priors['ham'] = math.log( self.num_messages['ham'] / (self.num_messages['spam'] + self.num_messages['ham'])) self.word_counts['spam'] = {} self.word_counts['ham'] = {} for x, y in zip(X, Y): c = 'spam' if y == 1 else 'ham' # 构建一个字典存储单封邮件中的单词以及其个数 counts = self.get_word_counts(self.tokenize(x)) for word, count in counts.items(): if word not in self.vocab: self.vocab.add(word)#确保self.vocab中含有所有邮件中的单词 # 下面语句是为了计算垃圾邮件和非垃圾邮件的词频,即给定词在垃圾邮件和非垃圾邮件中出现的次数。 # c是0或1,垃圾邮件的标签 if word not in self.word_counts[c]: self.word_counts[c][word] = 0.0 self.word_counts[c][word] += countMNB = SpamDetector_2()MNB.fit(X[100:], y[100:])5.测试算法 对测试集进行测试,判断是垃圾邮件还是正常邮件,并计算出准确度class SpamDetector(SpamDetector_2): def predict(self, X): result = [] flag_1 = 0 # 遍历所有的测试集 for x in X: counts = self.get_word_counts(self.tokenize(x)) # 生成可以记录单词以及该单词出现的次数的字典 spam_score = 0 ham_score = 0 flag_2 = 0 for word, _ in counts.items(): if word not in self.vocab: continue #下面计算P(内容|垃圾邮件)和P(内容|正常邮件),所有的单词都要进行拉普拉斯平滑 else: # 该单词存在于正常邮件的训练集和垃圾邮件的训练集当中 if word in self.word_counts['spam'].keys() and word in self.word_counts['ham'].keys(): log_w_given_spam = math.log( (self.word_counts['spam'][word] + 1) / (sum(self.word_counts['spam'].values()) + len(self.vocab))) log_w_given_ham = math.log( (self.word_counts['ham'][word] + 1) / (sum(self.word_counts['ham'].values()) + len( self.vocab))) # 该单词存在于垃圾邮件的训练集当中,但不存在于正常邮件的训练集当中 if word in self.word_counts['spam'].keys() and word not in self.word_counts['ham'].keys(): log_w_given_spam = math.log( (self.word_counts['spam'][word] + 1) / (sum(self.word_counts['spam'].values()) + len(self.vocab))) log_w_given_ham = math.log( 1 / (sum(self.word_counts['ham'].values()) + len( self.vocab))) # 该单词存在于正常邮件的训练集当中,但不存在于垃圾邮件的训练集当中 if word not in self.word_counts['spam'].keys() and word in self.word_counts['ham'].keys(): log_w_given_spam = math.log( 1 / (sum(self.word_counts['spam'].values()) + len(self.vocab))) log_w_given_ham = math.log( (self.word_counts['ham'][word] + 1) / (sum(self.word_counts['ham'].values()) + len( self.vocab))) # 把计算到的P(内容|垃圾邮件)和P(内容|正常邮件)加起来 spam_score += log_w_given_spam ham_score += log_w_given_ham flag_2 += 1 # 最后,还要把先验加上去,即P(垃圾邮件)和P(正常邮件) spam_score += self.log_class_priors['spam'] ham_score += self.log_class_priors['ham'] # 最后进行预测,如果spam_score > ham_score则标志为1,即垃圾邮件 if spam_score > ham_score: result.append(1) else: result.append(0) flag_1 += 1 return resultMNB = SpamDetector()MNB.fit(X[100:], y[100:])pred = MNB.predict(X[:100])true = y[:100]accuracy = 0for i in range(100): if pred[i] == true[i]: accuracy += 1print(accuracy) 最后得到准确率为99%

本文链接地址:https://www.jiuchutong.com/zhishi/298743.html 转载请保留说明!

上一篇:普通二本,去过阿里外包,到现在年薪40W+的高级测试工程师,我的两年转行心酸经历...(考到二本)

下一篇:vue系列(三)——手把手教你搭建一个vue3管理后台基础模板(vue系列教程)

  • 家居企业做微博营销的第一步要做的是什么?(家居企业做微博怎么写)

    家居企业做微博营销的第一步要做的是什么?(家居企业做微博怎么写)

  • iqoo8怎么换锁屏壁纸(vivoiqoo锁屏怎么换)

    iqoo8怎么换锁屏壁纸(vivoiqoo锁屏怎么换)

  • qq邮箱怎么绑定手机号呢(微信和qq邮箱怎么绑定)

    qq邮箱怎么绑定手机号呢(微信和qq邮箱怎么绑定)

  • 安卓手机怎么看苹果耳机的电量(安卓手机怎么看AirPods pro电量)

    安卓手机怎么看苹果耳机的电量(安卓手机怎么看AirPods pro电量)

  • 怎么发视频到qq邮箱里(怎么发视频到QQ空间)

    怎么发视频到qq邮箱里(怎么发视频到QQ空间)

  • pause键什么意思(pause键是什么意思啊)

    pause键什么意思(pause键是什么意思啊)

  • 订单编号可以查出什么(订单编号可以查到地址吗)

    订单编号可以查出什么(订单编号可以查到地址吗)

  • 英特尔10400f用什么主板(英特尔10400t)

    英特尔10400f用什么主板(英特尔10400t)

  • i34160相当于i5几代(i34160属于什么级别的)

    i34160相当于i5几代(i34160属于什么级别的)

  • qq音乐付费包和会员有什么区别(qq音乐付费包和绿钻)

    qq音乐付费包和会员有什么区别(qq音乐付费包和绿钻)

  • 未签名的描述文件怎么删除(未签名的描述文件怎么信任)

    未签名的描述文件怎么删除(未签名的描述文件怎么信任)

  • 打印机开不了机(打印机开不了机怎么办)

    打印机开不了机(打印机开不了机怎么办)

  • 手机屏幕出现一条绿色的线(手机屏幕出现一根绿线能修好吗)

    手机屏幕出现一条绿色的线(手机屏幕出现一根绿线能修好吗)

  • oppoby是什么型号(oppodesigned by oppo是什么型号)

    oppoby是什么型号(oppodesigned by oppo是什么型号)

  • 手机照片怎样压缩打包(手机照片怎样压缩)

    手机照片怎样压缩打包(手机照片怎样压缩)

  • 网易云注销后能恢复吗(网易云注销后能再注册吗)

    网易云注销后能恢复吗(网易云注销后能再注册吗)

  • oppoa5的呼吸灯怎么设置(oppoa59s呼吸灯)

    oppoa5的呼吸灯怎么设置(oppoa59s呼吸灯)

  • 小米8dc调光怎么开(小米8青春版调光)

    小米8dc调光怎么开(小米8青春版调光)

  • 手机卡消磁了怎么修复(手机卡消磁了怎么办妙招)

    手机卡消磁了怎么修复(手机卡消磁了怎么办妙招)

  • 火山小视频怎么发小火苗(火山小视频怎么找回原来的账号)

    火山小视频怎么发小火苗(火山小视频怎么找回原来的账号)

  • 荣耀v20防水吗(荣耀v20手机防水级别)

    荣耀v20防水吗(荣耀v20手机防水级别)

  • vivoz5x反向充电需要线吗(vivoz5x反向充电怎么用)

    vivoz5x反向充电需要线吗(vivoz5x反向充电怎么用)

  • b站永久封禁怎么解除(b站永久封禁怎么申诉)

    b站永久封禁怎么解除(b站永久封禁怎么申诉)

  • wps循环引用怎么处理(wps循环引用是什么意思)

    wps循环引用怎么处理(wps循环引用是什么意思)

  • 鸿蒙系统如何自定义图标样式?鸿蒙系统自定义图标样式教程(鸿蒙系统如何自定义图标)

    鸿蒙系统如何自定义图标样式?鸿蒙系统自定义图标样式教程(鸿蒙系统如何自定义图标)

  • PCTVoice.exe是什么进程 PCTVoice进程查询(pctversion)

    PCTVoice.exe是什么进程 PCTVoice进程查询(pctversion)

  • 印花税征税范围是什么
  • 未按期纳税申报怎么处理
  • 人头税是什么税
  • 出口货物退免税管理办法
  • 本年计算扣除限额的基数怎么算
  • 个体工商户生产经营所得税税率表2021
  • 研发支出月末结转转去哪里?
  • etc发票开票中
  • 本年利润每个月都转吗
  • 五联发票一般用于哪种单位
  • 上月未记账本月可以记账吗
  • 房地产企业开发间接费用包括哪些
  • 软件企业收到先征后返的增值税税款是否要缴纳企业所得税
  • 手撕发票怎么区分地区开具
  • 一次性医保补助多少钱
  • 代理进出口公司结售汇
  • 可抵扣进项税怎么抵扣
  • 公园门票可以开发票吗
  • 支付工资资产减少了吗
  • 离职后绩效奖金应该按照整月发吗
  • 上市公司现金分红怎么派发
  • 企业合并中或有负债确认后期调整的项目
  • 公司股权转让的印花税谁来承担
  • 系统安装工具有什么用
  • 怎么修改以前年度的账
  • 应收账款与主营业务收入的比率
  • 以前年度应付账款转收入
  • php 字符串 数组
  • 公司整体收购如何交税
  • 没有实际出资的股东怎么退出
  • 应交增值税专栏下的增值税项目有哪些?
  • 货物品种不太多而数量又相对较大
  • coco2017数据集介绍
  • 污水处理增值税税率
  • linux搭建web服务器
  • 十天学会css教程
  • 竣工结算审计费用在线计算器
  • wordpress标签tag文章
  • 工业总产值和营业收入关系
  • 卖家运费险是不是退货免费
  • 现金盘盈盘亏怎么处理
  • 公交卡充值发票报销单怎么填
  • 母子公司资产划转的涉税问题
  • 个人工程款收入算是工资流水吗
  • 政府补助增值税即征即退
  • 纸质承兑上的印花是什么
  • 一般纳税人跨月发票怎么作废
  • sqlserver2005数据库备份
  • 零申报是怎么回事
  • 假如出口收汇的钱怎么办
  • 企业培训费怎么入账
  • 销售返利如何做会计分录
  • 社保局退回生育津贴分录摘要
  • 会计费用类科目明细
  • 租赁存在的原因有哪些
  • 加油卡充值如何开发票
  • 社保缴费基数调整后对个人的影响
  • 收取加盟费的账户是什么
  • 专票不报销公司查得到吗
  • 公司筹建期间允许投标吗
  • 企业为职工购买互助保障项目经费可以有什么开支
  • 应收的货款
  • windows隐藏功能
  • ubuntu鼠标不能移动
  • vi编辑器是干什么的
  • win10怎么取消禁用
  • w7提高开机速度
  • win7如何获取管理员密码
  • linux内核文件夹
  • opengl画点
  • 让动画不再僵硬的软件
  • androidstudiosdk
  • win7怎么设定关机时间
  • unity怎么新建项目
  • python黑魔法手册
  • 全面解析俄乌武器对比
  • javascript面向对象编程指南 pdf
  • 陕西省税务系统
  • 上海市国家税务局地址
  • 税务检查工作方法有哪些
  • 免责声明:网站部分图片文字素材来源于网络,如有侵权,请及时告知,我们会第一时间删除,谢谢! 邮箱:opceo@qq.com

    鄂ICP备2023003026号

    网站地图: 企业信息 工商信息 财税知识 网络常识 编程技术

    友情链接: 武汉网站建设