使用python绘制cdf的多种实现方法


Posted in Python onFebruary 25, 2020

首先我们先用随机函数编造一个包含1000个数值的一维numpy数组,如下:

// An highlighted block
rng = np.random.RandomState(seed=12345)
samples = stats.norm.rvs(size=1000, random_state=rng)

接下来我们将使用各种方法画出以上数据的累积分布图

1、matplotlib.pyplot.hist()

def hist(self, x, bins=None, range=None, density=None, weights=None,
       cumulative=False, bottom=None, histtype='bar', align='mid',
       orientation='vertical', rwidth=None, log=False,
       color=None, label=None, stacked=False, normed=None,
       **kwargs):

第一种方法,我们使用matplotlib图形库中的hist函数,熟悉该库的人应该知道这是一个直方图绘制函数,以上是从API中找到的hist函数的所有参数,我们给出一维数组或者列表x,使用hist画出该数据的直方图。

直方图有两种形式,分别是概率分布直方图和累积分布直方图(可能说的不准确- -!),可以通过参数cucumulative来调节,默认为False,画出的是PDF,那么True画出的便是CDF直方图。

PDF(figure1)可以观察到整个数据在横轴范围内的分布,CDF(figure2)则可以看出不同的数据分布间的差异性,也可以观察到整个数据的增长趋势和波动情况。

使用python绘制cdf的多种实现方法

上图是概率分布直方图,纵轴代表概率,如果置参数normed=False,纵轴代表频数

使用python绘制cdf的多种实现方法

如果我们要观察两种数据分布的差异,可能使用直方图就不是很直观,各种直方柱会相互重叠,我们只需更改直方图的图像类型,令histtype=‘step',就会画出一条曲线来(Figure3,实际上就是将直方柱并在一起,除边界外颜色透明),类似于累积分布曲线。这时,我们就能很好地观察到不同数据分布曲线间的差异。

使用python绘制cdf的多种实现方法

2、numpy.histogram

def histogram(a, bins=10, range=None, normed=False, weights=None,
     density=None)

第二种方法我们使用numpy中画直方图的函数histogram,该函数不是一个直接的绘图函数(废话- -!过渡句,哈哈),给定一组数据a,它会返回两个数组hist和bin_edges,默认情况下hist是数据在各个区间上的频率,bin_edges是划分的各个区间的边界,说到这我们大概可以想到其实该函数算是上一个函数的底层函数,我们可以依据得到的这两个数组来画直方图,我们也可以用频率数组来直接画分布曲线(Figure4)

使用python绘制cdf的多种实现方法

这里我只给出了一个最原始的图像,直接用hist数组画的,如果想要变成合格的累积分布曲线图,纵轴为概率(频率乘区间长度),横轴为区间(从bin_edges数组中取n-1个)就可以了

3、stats.relfreq

def relfreq(a, numbins=10, defaultreallimits=None, weights=None)
Returns
-------
frequency : ndarray
  Binned values of relative frequency.
lowerlimit : float
  Lower real limit
binsize : float
  Width of each bin.
extrapoints : int
  Extra points.

第三种方法我们使用stats中的relfreq函数,该函数和第二种的方法类似,也并非是直接画图,而是返回关于直方图的一些数据,这里的frequency直接是概率而非频率,可以直接作为CDF图的纵轴,但是横轴需要自己计算,计算公式:

res.lowerlimit + np.linspace(0,res.binsize*res.frequency.size, res.frequency.size)

这个公式应该很好理解,我就不多说了,当然这些返回值都要依赖我们所给出的bins的数目。下面我给出一段代码,便是使用stats.relfreq画出概率分布直方图和累积分布曲线图。

rng = np.random.RandomState(seed=12345)
samples = stats.norm.rvs(size=1000, random_state=rng)
res = stats.relfreq(samples, numbins=25)
x = res.lowerlimit + np.linspace(0, res.binsize*res.frequency.size,res.frequency.size)
fig = plt.figure(figsize=(5, 4))
ax = fig.add_subplot(1, 1, 1)
ax.bar(x, res.frequency, width=res.binsize)
ax.set_title('Relative frequency histogram')
ax.set_xlim([x.min(), x.max()])
plt.show()

使用python绘制cdf的多种实现方法

rng = np.random.RandomState(seed=12345)
samples = stats.norm.rvs(size=1000, random_state=rng)
res = stats.relfreq(samples, numbins=25)
x = res.lowerlimit + np.linspace(0, res.binsize*res.frequency.size,res.frequency.size)
y=np.cumsum(res.frequency)
plt.plot(x,y)
plt.title('Figure6 累积分布直方图')
plt.show()

使用python绘制cdf的多种实现方法

以上就是本人整理出来的关于画cdf直方图和曲线的三种方法,整理这方面东西的初忠是在发现在进行数据分析的时候,概率分布直方图只能观察到数据大概的分布情况,而在不同的数据样本进行比较时却很难直观滴反映其差异性,通过看论文发现cdf可以做到这一点。

本人并不是数学专业出身,想要表达其意义,但有些描述和用词不当,大家借鉴就好。希望大家多多支持三水点靠木!

Python 相关文章推荐
Python 函数基础知识汇总
Mar 09 Python
Django框架首页和登录页分离操作示例
May 28 Python
在python中画正态分布图像的实例
Jul 08 Python
python3.7 的新特性详解
Jul 25 Python
Python解决pip install时出现的Could not fetch URL问题
Aug 01 Python
python 并发编程 阻塞IO模型原理解析
Aug 20 Python
Django实现简单网页弹出警告代码
Nov 15 Python
opencv-python 读取图像并转换颜色空间实例
Dec 09 Python
Python tensorflow实现mnist手写数字识别示例【非卷积与卷积实现】
Dec 19 Python
Django 项目通过加载不同env文件来区分不同环境
Feb 17 Python
Python实现列表拼接和去重的三种方式
Jul 02 Python
Python MNIST手写体识别详解与试练
Nov 07 Python
python GUI库图形界面开发之PyQt5开发环境配置与基础使用
Feb 25 #Python
python GUI库图形界面开发之PyQt5信号与槽基本操作
Feb 25 #Python
python GUI库图形界面开发之PyQt5信号与槽机制、自定义信号基础介绍
Feb 25 #Python
python模拟点击网页按钮实现方法
Feb 25 #Python
python GUI库图形界面开发之PyQt5动态加载QSS样式文件
Feb 25 #Python
python 计算概率密度、累计分布、逆函数的例子
Feb 25 #Python
python GUI库图形界面开发之PyQt5窗口背景与不规则窗口实例
Feb 25 #Python
You might like
php防止伪造的数据从URL提交方法
2014/06/27 PHP
php设置静态内容缓存时间的方法
2014/12/01 PHP
php实现在新浪云中使用imagick生成缩略图并上传的方法
2016/09/26 PHP
用JavaScript隐藏控件的方法
2009/09/21 Javascript
侧栏跟随滚动的简单实现代码
2013/03/18 Javascript
计算新浪Weibo消息长度(还可以输入119字)
2013/07/02 Javascript
js常用自定义公共函数汇总
2014/01/15 Javascript
深入分析Javascript跨域问题
2015/04/17 Javascript
总结JavaScript三种数据存储方式之间的区别
2016/05/03 Javascript
jQuery获取同级元素的简单代码
2016/07/09 Javascript
jQuery bt气泡实现悬停显示及移开隐藏功能的方法
2016/07/12 Javascript
JavaScript将DOM事件处理程序封装为event.js 出现的低级错误问题
2016/08/03 Javascript
javascript的几种写法总结
2016/09/30 Javascript
jQuery插件ajaxFileUpload使用实例解析
2016/10/19 Javascript
JS动态给对象添加属性和值的实现方法
2016/10/21 Javascript
基于BootStrap与jQuery.validate实现表单提交校验功能
2016/12/22 Javascript
对layui中table组件工具栏的使用详解
2019/09/19 Javascript
Angular封装表单控件及思想总结
2019/12/11 Javascript
Vue 构造选项 - 进阶使用说明
2020/08/14 Javascript
微信小程序自定义底部弹出框功能
2020/11/18 Javascript
js前端对于大量数据的展示方式及处理方法
2020/12/02 Javascript
python安装twisted的问题解析
2018/08/21 Python
wxPython的安装与使用教程
2018/08/31 Python
Python JSON格式数据的提取和保存的实现
2019/03/22 Python
python tools实现视频的每一帧提取并保存
2020/03/20 Python
Django 中自定义 Admin 样式与功能的实现方法
2019/07/04 Python
django 类视图的使用方法详解
2019/07/24 Python
使用 tf.nn.dynamic_rnn 展开时间维度方式
2020/01/21 Python
在keras 中获取张量 tensor 的维度大小实例
2020/06/10 Python
Python局部变量与全局变量区别原理解析
2020/07/14 Python
澳大利亚牛仔裤商店:Just Jeans
2016/10/13 全球购物
Sephora丝芙兰印尼官方网站:购买化妆品和护肤品
2018/07/02 全球购物
AJax面试题
2014/11/25 面试题
森林防火工作方案
2014/02/14 职场文书
遗产继承公证书
2014/04/09 职场文书
运动会宣传口号
2014/06/09 职场文书