【机器学习】朴素贝叶斯算法:从贝叶斯定理到手写数字识别实战

朴素贝叶斯算法:从贝叶斯定理到手写数字识别实战

简介

所谓的贝叶斯定理源于他生前为解决一个"逆概"问题写的一篇文章,而这篇文章是在他死后才由他的一位朋友发表出来的。在贝叶斯写这篇文章之前,人们已经能够计算"正向概率",如"假设袋子里面有 N 个白球,M 个黑球,你伸手进去摸一把,摸出黑球的概率是多大"。而一个自然而然的问题是反过来:"如果我们事先并不知道袋子里面黑白球的比例,而是闭着眼睛摸出一个(或好几个)球,观察这些取出来的球的颜色之后,那么我们可以就此对袋子里面的黑白球的比例作出什么样的推测"。这个问题,就是所谓的逆向概率问题。

贝叶斯定理的思想出现在18世纪,但真正大规模派上用途还得等到计算机的出现。因为这个定理需要大规模的数据计算推理才能凸显效果,它在很多计算机应用领域中都大有作为,如自然语言处理,机器学习,推荐系统,图像识别,博弈论等等。我们今天说的也就是最简单的贝叶斯在机器学习上的使用------朴素贝叶斯。

一、贝叶斯算法原理详解

贝叶斯定理

正向概率

假设袋子里有 10 个白球,90 个黑球,然后从袋子里面拿出一个球,拿出的球是白球的概率是多少?

P(摸出白球) = 10/100 = 1/10

逆向概率

如果我们事先并不知道袋子里白球、黑球的比例,然后通过多次试验,根据拿出来球的颜色推测袋子里白球、黑球的比例。

经典题目

在一个学校里面,男生占 60%,女生占 40%;男生:100% 穿长裤,女生:50% 穿长裤,50% 穿裙子;此时迎面走来一个穿长裤的学生,刚好你是高度近视,那么要判断一下这个穿长裤的学生是女生的概率该怎么计算?

实质

穿长裤的是女生的概率 = 女生中穿长裤的人数 / 穿长裤的总人数

假设

全校人数(U):1000 人

男生中穿长裤的人数:

U × P(男生) × P(穿长裤 | 男生) = 1000 × 60% × 100% = 600(人)

女生中穿长裤的人数:

U × P(女生) × P(穿长裤 | 女生) = 1000 × 40% × 50% = 200(人)

穿长裤且是女生的概率:

200 / (200 + 600) = 25%

贝叶斯算法

实质

穿长裤的是女生的概率 = 女生中穿长裤的人数 / 穿长裤的总人数

定义

A --- 男生,B --- 女生,C --- 穿长裤

那么:

P(B|C) = U × P(B) × P(C\|B) / U × P(A) × P(C\|A) + U × P(B) × P(C\|B)

= P(B) × P(C\|B) / P©

= (40% × 50%) / (60% + 40% × 50%)

= 25%

P© = P(B) × P(C\|B) / P©

贝叶斯如何将样本数据进行计算实现预测

示例:一个包含 3 个特征(A、B 和 C)和一个类别标签(D)的简单数据集

特征 A 特征 B 特征 C 类别 D
1 0 1 0
0 1 1 1
0 1 0 0
1 0 0 1
1 0 0 0
0 1 1 1

目标:使用朴素贝叶斯分类器来预测一个新的样本(特征 A=1,特征 B=1,特征 C=0)的类别。特征 A=0,特征 B=0,特征 C=1 属于?

核心:把输入当作条件,结果为 0 的概率是多大?结果为 1 的概率是多大?

公式

首先,我们计算类别的先验概率:

P(D=0) = 3/6 = 0.5

P(D=1) = 3/6 = 0.5

然后,我们计算每个特征在每个类别下的条件概率:

P(A=1|D=0) = 2/3

P(A=1|D=1) = 1/3

P(B=1|D=0) = 1/3

P(B=1|D=1) = 2/3

P(C=0|D=0) = 1/3

P(C=0|D=1) = 1/3

现在,我们使用朴素贝叶斯分类器来计算给定特征值下每个类别的后验概率:

P(D=0|A=1,B=1,C=0) ∝ P(D=0) × P(A=1|D=0) × P(B=1|D=0) × P(C=0|D=0)

= 0.5 × (2/3) × (1/3) × (1/3) = 0.5 × 2/27 = 1/27

P(D=1|A=1,B=1,C=0) ∝ P(D=1) × P(A=1|D=1) × P(B=1|D=1) × P(C=0|D=1)

= 0.5 × (1/3) × (2/3) × (1/3) = 0.5 × 2/27 = 1/27

所以特征(A=1,特征 B=1,特征 C=0)的类别 D=0,由此可见贝叶斯算法是用于分类的算法。

二、朴素贝叶斯的 API

python 复制代码
class sklearn.naive_bayes.MultinomialNB(alpha=1.0, fit_prior=True, class_prior=None)

参数详解

参数 说明
alpha 平滑参数,默认 1.0。alpha=1.0 为拉普拉斯平滑;0 < alpha < 1 为利德斯通平滑;alpha=0 不进行平滑
fit_prior 布尔值,默认 True。True 时从数据中学习先验概率;False 时使用等先验
class_prior 手动指定先验概率,默认 None。若为 None,则由 fit_prior 控制

三、案例分析:手写数字识别

对于朴素贝叶斯算法,该算法的准确率相对于没有那么高,我们重点学习这个算法的思想,下面我就以手写识别为一个案例简单说明一下朴素贝叶斯的使用。

使用 scikit-learn 库实现了基于多项式朴素贝叶斯(MultinomialNB)的手写数字识别分类任务。

1. 导入必要的库

python 复制代码
from sklearn.datasets import load_digits
from sklearn.naive_bayes import MultinomialNB
from sklearn.model_selection import train_test_split
from sklearn import metrics

2. 加载并准备数据

python 复制代码
digits = load_digits()   # 加载手写数字数据集
x = digits.data          # 获取特征数据(图像的像素特征)
y = digits.target        # 获取标签数据(对应的数字类别,0-9)

数据集包含 8×8 像素的手写数字图像,每个图像被转换为 64 维的特征向量(8×8)。

3. 分割数据集

python 复制代码
x_train, x_test, y_train, y_test = train_test_split(
    x, y, test_size=0.2, random_state=0
)

将数据分为训练集(80%)和测试集(20%),test_size=0.2 表示测试集占比 20%,random_state=0 固定随机种子,确保结果可复现。

4. 创建并训练模型

python 复制代码
estimator = MultinomialNB(alpha=1)   # 创建多项式朴素贝叶斯模型,使用拉普拉斯平滑
estimator.fit(x_train, y_train)      # 用训练数据拟合模型

alpha=1 表示使用拉普拉斯平滑,避免零概率问题。

5. 模型预测与评估

python 复制代码
test_predicted = estimator.predict(x_test)   # 用训练好的模型预测测试集
scores = estimator.score(x_test, y_test)     # 计算模型在测试集上的准确率
print(scores)                                # 输出准确率
print(metrics.classification_report(y_test, test_predicted))   # 输出详细分类报告

输出结果说明:

第一个输出是模型的准确率(accuracy),即正确分类的样本占总样本的比例

第二个输出是详细的分类报告,包含:

每个类别的精确率(precision)、召回率(recall)、F1 分数(F1-score)

每个类别的支持样本数(support)

整体的宏平均(macro avg)和加权平均(weighted avg)指标

总结

朴素贝叶斯算法基于贝叶斯定理,通过计算后验概率进行分类。其核心假设是特征之间相互独立(条件独立性假设),虽然这一假设在实际中往往不成立,但朴素贝叶斯在许多场景下仍能取得不错的分类效果。该算法计算效率高、对缺失数据不敏感,适合文本分类、垃圾邮件过滤等任务。但同时也需要关注数据稀疏带来的零概率问题,可通过拉普拉斯平滑等技术进行缓解。

相关推荐
皓月斯语1 小时前
P5736 【深基7.例2】质数筛
c++·算法·题解
腾讯云大数据1 小时前
从多模态数据处理到模型训练:腾讯云EMR-Ray打通Data+AI全流程
人工智能·云计算·腾讯云·mapreduce·腾讯云大数据
louyu6668881 小时前
国产楼宇自控系统哪家靠谱?
大数据·人工智能
FriendshipT1 小时前
Ubuntu 20.04 下使用 Ollama 本地部署 AI 大模型
linux·人工智能·python·深度学习·ubuntu
大模型丫丫1 小时前
Loop Engineering:从强化学习视角看 Agent 循环的本质
java·人工智能·学习
乱世刀疤1 小时前
AI Weekly 7.27-8.2
人工智能
倒流时光三十年1 小时前
第五阶段 41 · 相关性调优(boosting、function_score、filter vs query 上下文)
机器学习·集成学习·boosting
天天爱吃肉82182 小时前
# 商用车多体动力学实战笔记|第7篇:制动系统与制动热衰退、ABS滞环控制
大数据·人工智能·笔记·python·嵌入式硬件·汽车
牧艺2 小时前
别让 Agent 猜需求:前端用「一页 Spec」把返工砍掉一半
人工智能·agent·vibecoding
时空节拍AI数字人2 小时前
多模态交互数字人:语音+视觉+触控如何融合
人工智能·microsoft·ai·aigc·交互·语音识别