使用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之wxPython菜单使用详解
Sep 28 Python
Python编码类型转换方法详解
Jul 01 Python
解读! Python在人工智能中的作用
Nov 14 Python
Windows下python3.6.4安装教程
Jul 31 Python
python计算两个矩形框重合百分比的实例
Nov 07 Python
使用TensorFlow实现二分类的方法示例
Feb 05 Python
详解python3中用HTMLTestRunner.py报ImportError: No module named 'StringIO'如何解决
Aug 27 Python
Django media static外部访问Django中的图片设置教程
Apr 07 Python
在ipython notebook中使用argparse方式
Apr 20 Python
python中的测试框架
Nov 13 Python
python用分数表示矩阵的方法实例
Jan 11 Python
Python 中如何使用 virtualenv 管理虚拟环境
Jan 21 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与Java进行通信的实现方法
2013/10/21 PHP
Yii2创建控制器(createController)方法详解
2016/07/23 PHP
PHP简单实现合并2个数字键数组值的方法
2017/05/30 PHP
javascript监听鼠标滚轮事件浅析
2014/06/05 Javascript
js创建表单元素并使用submit进行提交
2014/08/14 Javascript
jQuery实现倒计时按钮功能代码分享
2014/09/03 Javascript
AngularJS 防止页面闪烁的方法
2017/03/09 Javascript
jQuery实现可兼容IE6的滚动监听功能
2017/09/20 jQuery
JavaScript 日期时间选择器一些小结
2018/04/02 Javascript
echarts实现地图定时切换散点与多图表级联联动详解
2018/08/07 Javascript
在node中使用jwt签发与验证token的方法
2019/04/03 Javascript
javascript实现扫雷简易版
2020/08/18 Javascript
[05:20]2018DOTA2亚洲邀请赛主赛事第三日战况回顾 LGD率先挺进胜者组决赛
2018/04/06 DOTA
Python基础入门之seed()方法的使用
2015/05/15 Python
pymssql数据库操作MSSQL2005实例分析
2015/05/25 Python
python实现判断数组是否包含指定元素的方法
2015/07/15 Python
python实现汉诺塔方法汇总
2016/07/25 Python
python 将大文件切分为多个小文件的实例
2019/01/14 Python
python zip()函数使用方法解析
2019/10/31 Python
对python中assert、isinstance的用法详解
2019/11/27 Python
深入了解如何基于Python读写Kafka
2019/12/31 Python
python GUI编程(Tkinter) 创建子窗口及在窗口上用图片绘图实例
2020/03/04 Python
python中wx模块的具体使用方法
2020/05/15 Python
Python容器类型公共方法总结
2020/08/19 Python
PyCharm 2020.2下配置Anaconda环境的方法步骤
2020/09/23 Python
详解Python中第三方库Faker
2020/09/25 Python
详解pycharm连接远程linux服务器的虚拟环境的方法
2020/11/13 Python
eDreams澳大利亚:预订机票、酒店和度假产品
2017/04/19 全球购物
Pamela Love官网:纽约设计师Pamela Love的精美、时尚和穿孔珠宝
2020/10/19 全球购物
护校行动方案
2014/05/31 职场文书
欢迎词怎么写
2015/01/23 职场文书
奖学金个人总结
2015/03/04 职场文书
2015年清明节扫墓演讲稿
2015/03/18 职场文书
2015年办公室文员工作总结
2015/04/24 职场文书
防汛通知
2015/04/25 职场文书
合同审查法律意见书
2015/06/04 职场文书