Python求解正态分布置信区间教程


Posted in Python onNovember 20, 2019

正态分布和置信区间

正态分布(Normal Distribution)又叫高斯分布,是一种非常重要的概率分布。其概率密度函数的数学表达如下:

Python求解正态分布置信区间教程

置信区间是对该区间能包含未知参数的可置信的程度的描述。

使用SciPy求解置信区间

import numpy as np
import matplotlib.pyplot as plt
from scipy import stats

N = 10000
x = np.random.normal(0, 1, N)
# ddof取值为1是因为在统计学中样本的标准偏差除的是(N-1)而不是N,统计学中的标准偏差除的是N
# SciPy中的std计算默认是采用统计学中标准差的计算方式
mean, std = x.mean(), x.std(ddof=1)
print(mean, std)
# 计算置信区间
# 这里的0.9是置信水平
conf_intveral = stats.norm.interval(0.9, loc=mean, scale=std)
print(conf_intveral)

输出如下:

0.0033541207210673997 0.9986647964318905
(-1.639303291798682, 1.6460115332408163)

这里的-1.639303291798682是置信上界,1.6460115332408163是置信下界,两个数值构成的区间就是置信区间

使用Matplotlib绘制正态分布密度曲线

# 绘制概率密度分布图
x = np.arange(-5, 5, 0.001)
# PDF是概率密度函数
y = stats.norm.pdf(x, loc=mean, scale=std)
plt.plot(x, y)
plt.show()

这里的pdf()函数是Probability density function,就是本文最开始的那个公式

最后的输出图像如下,可以看到结果跟理论上的正太分布还是比较像的:

Python求解正态分布置信区间教程

正态分布置信区间规律

函数曲线下68.268949%的面积在平均数左右的一个标准差范围内

函数曲线下95.449974%的面积在平均数左右两个标准差的范围内

函数曲线下99.730020%的面积在平均数左右三个标准差的范围内

函数曲线下99.993666%的面积在平均数左右四个标准差的范围内

以上这篇Python求解正态分布置信区间教程就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持三水点靠木。

Python 相关文章推荐
Python实现国外赌场热门游戏Craps(双骰子)
Mar 31 Python
Python算法应用实战之栈详解
Feb 04 Python
解决python "No module named pip" 的问题
Oct 13 Python
Python3实现腾讯云OCR识别
Nov 27 Python
浅谈python 读excel数值为浮点型的问题
Dec 25 Python
python面向对象实现名片管理系统文件版
Apr 26 Python
Python编写带选项的命令行程序方法
Aug 13 Python
详解Ubuntu环境下部署Django+uwsgi+nginx总结
Apr 02 Python
Python CSS选择器爬取京东网商品信息过程解析
Jun 01 Python
解决python 执行sql语句时所传参数含有单引号的问题
Jun 06 Python
浅析PyCharm 的初始设置(知道)
Oct 12 Python
Python与C++中梯度方向直方图的实现
Mar 17 Python
Python pandas RFM模型应用实例详解
Nov 20 #Python
使用Python实现正态分布、正态分布采样
Nov 20 #Python
Python pandas自定义函数的使用方法示例
Nov 20 #Python
Python求正态分布曲线下面积实例
Nov 20 #Python
复化梯形求积分实例——用Python进行数值计算
Nov 20 #Python
Python实现数值积分方式
Nov 20 #Python
基于Python批量生成指定尺寸缩略图代码实例
Nov 20 #Python
You might like
php 伪造本地文件包含漏洞的代码
2011/11/03 PHP
PHP ? EasyUI DataGrid 资料取的方式介绍
2012/11/07 PHP
php使用json_encode对变量json编码
2014/04/07 PHP
PHP SPL标准库之数据结构栈(SplStack)介绍
2015/05/12 PHP
基础的WordPress插件制作教程
2015/11/24 PHP
使用php从身份证号中获取一系列线索(星座、生肖、生日等)
2016/05/11 PHP
thinkphp5.0自定义验证规则使用方法
2017/11/16 PHP
thinkPHP5框架实现基于ajax的分页功能示例
2018/06/12 PHP
laravel-admin 管理平台获取当前登陆用户信息的例子
2019/10/08 PHP
从JavaScript的函数重名看其初始化方式
2007/03/08 Javascript
JQuery 将元素显示在屏幕的中央的代码
2010/02/27 Javascript
JQuery的Ajax跨域请求原理概述及实例
2013/04/26 Javascript
js中array的sort()方法使用介绍
2014/02/20 Javascript
jQuery实现点击该行即可删除HTML表格行
2014/10/17 Javascript
详解Wondows下Node.js使用MongoDB的环境配置
2016/03/01 Javascript
jQuery soColorPacker 网页拾色器
2016/06/22 Javascript
Vue.js双向绑定操作技巧(初级入门)
2016/12/27 Javascript
jQuery的三种bind/One/Live/On事件绑定使用方法
2017/02/23 Javascript
Express URL跳转(重定向)的实现方法
2017/04/07 Javascript
解决Angular.js中使用Swiper插件不能滑动的问题
2018/02/26 Javascript
基于axios 解决跨域cookie丢失的问题
2018/09/26 Javascript
微信小程序整合使用富文本编辑器的方法详解
2019/04/25 Javascript
Vue.js自定义指令学习使用详解
2019/10/19 Javascript
JS实现点餐自动选择框(案例分析)
2019/12/10 Javascript
微信小程序实现下滑到底部自动翻页功能
2020/03/07 Javascript
一文读懂vue动态属性数据绑定(v-bind指令)
2020/07/20 Javascript
JavaScript实现矩形块大小任意缩放
2020/08/25 Javascript
Python列表(list)、字典(dict)、字符串(string)基本操作小结
2014/11/28 Python
CentOS 6.5下安装Python 3.5.2(与Python2并存)
2017/06/05 Python
python中的插值 scipy-interp的实现代码
2018/07/23 Python
Django中的ajax请求
2018/10/19 Python
Python3模拟登录操作实例分析
2019/03/12 Python
解决Django加载静态资源失败的问题
2019/07/28 Python
python将dict中的unicode打印成中文实例
2020/05/11 Python
惠普香港官方商店:HP香港
2019/04/30 全球购物
小学评语大全
2014/04/22 职场文书