[Python] KDE图[作密度图(Kernel Density Estimate,核密度估计)]介绍和使用场景(案例)

news2024/9/29 21:23:30

KDE图是什么?

核密度估计(Kernel Density Estimate,KDE)是一种非参数统计方法,用于估计未知随机变量的概率分布。它通过在每个数据点附近放置一个核函数,并将这些核函数加总起来,得到对概率分布的估计。

KDE的主要思想是通过在每个数据点附近放置一个核函数来估计概率分布。核函数可以是各种形式,常用的有高斯核、均匀核等。核函数在数据点附近产生一个非负的函数值,表示该点附近的概率密度。然后,将所有核函数的值加总起来,得到整个数据集的概率密度函数。

KDE图可以理解为是对直方图的加窗平滑。KDE图解决了一个基本的数据平滑问题,即根据有限的数据样本对总体进行推断。通过KDE分布图,可以查看并对训练数据集和测试数据集中特征变量的分布情况。

KDE图定义

设 (x1,x2,...,xn) 是从某个单变量分布中抽取的独立同分布样本,在任何给定点 x 具有未知密度 f 。我们有兴趣估计这个函数 f 的形状。它的核密度估计量是:

其中 K 是内核——一个非负函数——h > 0 是一个称为带宽的平滑参数。下标为 h 的内核称为缩放内核,定义为 Kh(x) = 1/h K(x/h)。直觉上,人们希望选择数据允许的尽可能小的 h;然而,估计量的偏差与其方差之间总是存在权衡。

通常使用一系列核函数:uniform、triangular、biweight、triweight、Epanechnikov、normal 等。Epanechnikov 内核在均方误差意义上是最优的,尽管前面列出的内核的效率损失很小。 由于其方便的数学性质,常使用正态核,即 K(x) = ϕ(x),其中 ϕ 是标准正态密度函数。核密度估计的构建在密度估计之外的领域找到了解释。 例如,在热力学中,这相当于将热核(热方程的基本解)放置在每个数据点位置 xi 时产生的热量。类似的方法用于在点云上构建离散拉普拉斯算子以进行流形学习(例如扩散图)。

例子

核密度估计与直方图密切相关,但可以通过使用合适的核赋予平滑度或连续性等属性。下图基于这 6 个数据点说明了这种关系:

样品123456
-2.1-1.3-0.41.95.16.2

对于直方图,首先,水平轴被划分为覆盖数据范围的子区间或区间:在这种情况下,每个区间有 6 个宽度为 2 的区间。每当一个数据点落在该区间内时,一个高度为 1 的区间/12 放在那里。如果不止一个数据点落在同一个箱子里,这些箱子就会堆叠在一起。对于核密度估计,将方差为 2.25 的正态核(由红色虚线表示)放置在每个数据点 xi 上。对内核求和以进行内核密度估计(蓝色实线)。核密度估计的平滑性(与直方图的离散性相比)说明了核密度估计如何更快地收敛到连续随机变量的真实基础密度。 

使用相同数据构建的直方图(左)和核密度估计(右)的比较。 六个单独的内核是红色虚线,内核密度估计蓝色曲线。 数据点是横轴上的地毯图。

如上图的右图所示,我理解的红色曲线就是基函数,类似于傅里叶函数的基函数叠加形成上面的蓝色曲线。

KDE图使用场景

KDE可以用于多种场景,其中一些案例包括:

  1. 数据可视化:KDE可以用来绘制密度图,可以清晰地展示数据分布的形态。例如,在金融领域,可以使用KDE绘制股票收益率的密度图,以帮助理解收益率的分布情况。

  2. 数据挖掘:KDE可以用于异常检测。通过估计数据的概率密度分布,可以确定哪些数据点的概率较低,从而识别出异常点。例如,在网络安全领域,可以使用KDE来检测网络流量中的异常活动。

  3. 信号处理:KDE可以用于信号处理中的噪声估计。通过估计信号的概率密度分布,可以确定噪声的分布情况,进而进行噪声滤波等操作。

KDE图使用案例

在这个例子中,我们导入了tips数据集,并使用sns.kdeplot()函数绘制了total_bill变量的KDE图和直方图。可以使用title()和xlabel()函数设置图表标题和轴标签。最后,使用plt.show()函数显示图表。

import matplotlib.pyplot as plt
import seaborn as sns

# 导入数据集
tips = sns.load_dataset('tips')
tips

 

# 绘制KDE图
sns.kdeplot(data=tips, x='total_bill')

# 设置图表标题和轴标签
plt.title('KDE Plot of Total Bill')
plt.xlabel('Total Bill')

# 显示图表
plt.show()

# 绘制直方图
sns.histplot(data=tips, x='total_bill')

# 设置图表标题和轴标签
plt.title('Hist Plot of Total Bill')
plt.xlabel('Total Bill')

# 显示图表
plt.show()

结论

总之,KDE是一种非参数统计方法,可以用于估计未知随机变量的概率分布。它在数据可视化、数据挖掘和信号处理等领域都有广泛的应用。

参考资料

KDE: Kernel density estimation 核密度估计

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/1410175.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

用可视化案例讲Rust编程4. 用泛型和特性实现自适配shapefile的读取

本节已经涉及Rust学习曲线上的一个大坑:泛型和特性了,属于语言的深水区,如果初学者,建议看一眼知道有这个功能即可。 如果我们立足于功能实现,那么做到像上一节那样就可以了,从原理上来说,每个…

javascript跳转微信: 扫码,分享链接,下载二维码

文章目录 微信后台配置扫码前端配置扫描二维码分享链接下载二维码 总结 微信后台配置扫码 参考: 微信二维码内容获取 在开发设置中配置扫码路由以及跳转路径 完成该配置后,可以通过扫描对应地址的二维码进入小程序的指定页面 前端配置 扫描二维码 …

【github】使用github action 拉取国外docker镜像

使用github action 拉取国外docker镜像 k8s部署经常用到国外镜像,如果本地无法拉取可以考虑使用github action环境 github action的ci服务器在国外,不受中国防火墙影响github action 自带docker命令运行时直接将你仓库代码拉取下来 步骤 你的国内dock…

测试人年终总结:入行三年,下一步怎么走,思想碰撞

原贴地址:入行三年,下一步怎么走,思想碰撞 TesterHome 熟悉环境,进步缓慢;停止思考,举步不前(为什么会有这篇文章why 初心变质:计算机系毕业,毕业时的打算是从测试进&a…

如何实现人口、GDP等的空间化效果?

现在假设你遇到了某个情况,手头有某市十几个区县的人口总数但是你想要更加精细化的空间化数据产品,以方便后续进行核密低分析或者进行其它的空间、数据分析。 那我们应该如何通过十几个点的人口数据进行空间化呢,详细步骤如下: 第…

PHP“引用”漏洞

今日例题&#xff1a; <?php highlight_file(__FILE__); error_reporting(0); include("flag.php"); class just4fun { var $enter; var $secret; } if (isset($_GET[pass])) { $pass $_GET[pass]; $passstr_replace(*,\*,$pass); } $o unser…

Flume介绍

一、介绍 Apache Flume 是一种分布式、可靠且可用的系统&#xff0c;用于有效地收集、汇总大量日志数据&#xff0c;并将其从多个不同来源转移到集中式数据存储区。 Apache Flume 的使用不仅限于日志数据聚合。由于数据源是可定制的&#xff0c;Flume 可用于传输大量事件数据&a…

OpenCV书签 #互信息的原理与相似图片搜索实验

1. 介绍 互信息&#xff08;Mutual Information&#xff09; 是信息论中的一个概念&#xff0c;用于衡量两个随机变量之间的关联程度。在图像处理和计算机视觉中&#xff0c;互信息常被用来度量两幅图像之间的相似性。 互信息可以看成是一个随机变量中包含的关于另一个随机变…

腾讯云tsf平台-部署微服务项目

腾讯云tsf平台-部署微服务项目 一、腾讯云tsf平台简介二、部署准备0&#xff08;数据库、中间件等部署&#xff09;三、部署准备1&#xff08;创建集群和命名空间&#xff09;1、准备部署资源--集群2、使用容器部署微服务步骤 1&#xff1a;创建容器集群步骤 2&#xff1a;创建…

数据结构:非完全二叉树(递归实现)

非完全二叉树是指在二叉树中&#xff0c;除了叶子节点&#xff08;无子节点&#xff09;外&#xff0c;其他节点的子节点个数可以不同&#xff0c;即不一定是每个节点都有两个子节点&#xff0c;有右孩子时也不一定有左孩子。 tree.h /* * 文件名称&#xff1a;tree.h * …

SpringBootAdmin邮件通知

在上一篇中我们学习到了 Admin Service Clustering 分布式缓存配置 &#xff0c;这一篇我们来学习&#xff0c;客户端离线&#xff0c;出现故障的时候&#xff0c;我们这么能及时知道呢&#xff0c;发现邮件通知也许是最简单的方式了&#xff01; 邮件通知 邮件通知将作为使用…

代码随想录 Leetcode111. 二叉树的最小深度

题目&#xff1a; 代码(首刷自解 2024年1月24日&#xff09;&#xff1a; class Solution { public:int minDepth(TreeNode* root) {if(root nullptr) return 0;queue<TreeNode*> que;TreeNode* cur root;que.push(cur);int size 0;int depth 0;while (!que.empty()…

【方法】如何打开7Z格式压缩文件?

7Z是一种常见的压缩文件格式&#xff0c;那遇到7Z压缩文件&#xff0c;要如何打开呢&#xff1f; 想要打开7Z压缩包&#xff0c;需要用到解压缩软件&#xff0c;比如WinRAR、7-Zip等软件。安装软件后&#xff0c;只需要用鼠标双击7Z压缩包&#xff0c;就可以打开压缩包&#x…

高数总结(3

目录 1.总结&#xff1a;小结&#xff1a; 1.总结&#xff1a; 小结&#xff1a; 关注我给大家分享更多有趣的知识&#xff0c;以下是个人公众号&#xff0c;提供 ||代码兼职|| ||代码问题求解|| 由于本号流量还不足以发表推广&#xff0c;搜我的公众号即可&#xff1a;

ECharts实现简单饼图和柱状图

1.JAVA 1.饼图 前端使用vue&#xff0c;后端使用SpringBoot <template><div><div class"card" style"padding: 15px">数据可视化分析&#xff1a;图书类型销量</div><div style"display: flex; margin: 10px 0"&g…

autosar学习笔记 之SecOC

SecOC 接下来SecOC标准就更复杂一点,它不单单是做了通讯校验。 SecOC是基于对称密钥加密的一套机制,需要对ECU间的通讯作身份认证处理,来更好的防止伪装攻击,谈起对称或非对称加密,就会涉及到密钥的存储和Mac值的计算。 因此SECOC机制对于密钥的硬件存储,也有一定的要求…

从零学习Hession RPC

为什么学习Hessian RPC&#xff1f; 存粹的RPC&#xff0c;只解决PRC的四个核心问题&#xff08;1.网络通信2.协议 3.序列化 4.代理&#xff09;Java写的HessianRPC落伍了&#xff0c;但是它的序列化方式还保存着&#xff0c;被Dubbo(Hessian Lite)使用。 被落伍&#xff0c;只…

K8S图像化工具rancher

Rancher是一个开源的企业级多集群的k8s管理平台 Rancher和k8s的区别 都是为了容器的调度和编排系统&#xff0c;但是rancher不仅能够调度&#xff0c;还能挂历k8s集群&#xff0c;自带监控&#xff08;普罗米修斯&#xff09;&#xff0c;你哪怕不知带k8s是什么&#xff0c;一样…

如何在Kali系统配置启动SSH并结合内网穿透实现远程连接

文章目录 1. 启动kali ssh 服务2. kali 安装cpolar 内网穿透3. 配置kali ssh公网地址4. 远程连接5. 固定连接SSH公网地址6. SSH固定地址连接测试 简单几步通过[cpolar 内网穿透](cpolar官网-安全的内网穿透工具 | 无需公网ip | 远程访问 | 搭建网站)软件实现ssh 远程连接kali! …

【Kafka】Kafka安装:Linux本地和Docker

目录 Linux本地安装kafkajava环境配置Zookeeper的安装配置Kafka的安装与配置生产与消费 Docker安装kafkaZookeeper安装Kafka安装 Linux本地安装kafka java环境配置 1、上传jdk-8u261-linux-x64.rpm到服务器并安装&#xff1a; rpm -ivh jdk-8u261-linux-x64.rpm2、配置环境变…