使用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当中的字符串和编码
Apr 25 Python
介绍Python的@property装饰器的用法
Apr 28 Python
对pandas的dataframe绘图并保存的实现方法
Aug 05 Python
python3.5 email实现发送邮件功能
May 22 Python
安装好Pycharm后如何配置Python解释器简易教程
Jun 28 Python
python3 中的字符串(单引号、双引号、三引号)以及字符串与数字的运算
Jul 18 Python
Django实现auth模块下的登录注册与注销功能
Oct 10 Python
python hashlib加密实现代码
Oct 17 Python
Python中生成一个指定长度的随机字符串实现示例
Nov 06 Python
Pytorch使用MNIST数据集实现基础GAN和DCGAN详解
Jan 10 Python
python 子类调用父类的构造函数实例
Mar 12 Python
python 实现图片批量压缩的示例
Dec 18 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
CPU步进是什么意思?i3-9100F B0步进和U0步进区别知识科普
2020/03/17 数码科技
PHP移动文件指针ftell()、fseek()、rewind()函数总结
2014/11/18 PHP
php实现上传图片文件代码
2015/07/19 PHP
thinkPHP5.0框架应用请求生命周期分析
2017/03/25 PHP
Laravel服务容器绑定的几种方法总结
2020/06/14 PHP
jQuery结合ajax实现动态加载文本内容
2015/05/19 Javascript
对象题目的一个坑 理解Javascript对象
2015/12/22 Javascript
原生js配合cookie制作保存路径的拖拽
2015/12/29 Javascript
ECharts仪表盘实例代码(附源码下载)
2016/02/18 Javascript
ES6新特征数字、数组、字符串
2016/10/01 Javascript
JS DOMReady事件的六种实现方法总结
2016/11/23 Javascript
JavaScript原生编写《飞机大战坦克》游戏完整实例
2017/01/04 Javascript
tab栏切换原理
2017/03/22 Javascript
动手写一个angular版本的Message组件的方法
2017/12/16 Javascript
js正则取值的结果数组调试方法
2018/10/10 Javascript
使用p5.js临摹动态图片
2019/11/04 Javascript
基于JavaScript获取url参数2种方法
2020/04/17 Javascript
[58:54]EG vs RNG 2019国际邀请赛小组赛 BO2 第一场 8.16
2019/08/18 DOTA
Python实现二分查找与bisect模块详解
2017/01/13 Python
不可错过的十本Python好书
2017/07/06 Python
浅谈python和C语言混编的几种方式(推荐)
2017/09/27 Python
Python简单计算文件MD5值的方法示例
2018/04/11 Python
使用Selenium破解新浪微博的四宫格验证码
2018/10/19 Python
Python3简单实现串口通信的方法
2019/06/12 Python
浅谈CSS3 动画卡顿解决方案
2019/01/02 HTML / CSS
canvas 绘图时位置偏离的问题解决
2020/09/16 HTML / CSS
演讲稿的格式及范文
2014/08/22 职场文书
入党政审材料范文
2014/12/24 职场文书
2015年南京大屠杀纪念日活动总结
2015/03/24 职场文书
原告代理词范文
2015/05/25 职场文书
新闻稿件写作技巧
2015/07/18 职场文书
二胎满月酒致辞
2015/07/29 职场文书
详解Laravel服务容器的优势
2021/05/29 PHP
MySQL5.7并行复制原理及实现
2021/06/03 MySQL
SpringBoot连接MySQL获取数据写后端接口的操作方法
2021/11/02 MySQL
httpclient调用远程接口的方法
2022/08/14 Java/Android