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中的列表推导浅析
Apr 26 Python
python字符串连接方法分析
Apr 12 Python
Python探索之SocketServer详解
Oct 28 Python
Python实现端口检测的方法
Jul 24 Python
python 定义n个变量方法 (变量声明自动化)
Nov 10 Python
Python3.7 dataclass使用指南小结
Feb 22 Python
Python for循环与range函数的使用详解
Mar 23 Python
Python中文分词库jieba,pkusegwg性能准确度比较
Feb 11 Python
python简单的三元一次方程求解实例
Apr 02 Python
python Matplotlib数据可视化(1):简单入门
Sep 30 Python
python 实现一个简单的线性回归案例
Dec 17 Python
JAVA SpringMVC实现自定义拦截器
Mar 16 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
crontab无法执行php的解决方法
2016/01/25 PHP
ThinkPHP框架实现导出excel数据的方法示例【基于PHPExcel】
2018/05/12 PHP
一页面多XMLHttpRequest对象
2007/01/22 Javascript
javascript实现的鼠标链接提示效果生成器代码
2007/06/28 Javascript
Javascript String对象扩展HTML编码和解码的方法
2009/06/02 Javascript
jQuery 版本的文本输入框检查器Input Check
2009/07/09 Javascript
js 文件引入实现代码
2010/04/23 Javascript
Jquery工作常用实例 使用AJAX使网页进行异步更新
2011/07/26 Javascript
javascript实现2016新年版日历
2016/01/25 Javascript
Javascript实现苹果悬浮虚拟按钮
2016/04/10 Javascript
JQuery核心函数是什么及使用方法介绍
2016/05/03 Javascript
KnockoutJS 3.X API 第四章之表单textInput、hasFocus、checked绑定
2016/10/11 Javascript
JavaScript下拉菜单功能实例代码
2017/03/01 Javascript
详解如何搭建mpvue框架搭配vant组件库的小程序项目
2019/05/16 Javascript
node基于async/await对mysql进行封装
2019/06/20 Javascript
小程序选项卡以及swiper套用(跨页面)
2020/06/19 Javascript
[33:33]完美世界DOTA2联赛PWL S2 FTD.C vs SZ 第二场 11.27
2020/11/30 DOTA
python中的reduce内建函数使用方法指南
2014/08/31 Python
Python函数式编程指南(三):迭代器详解
2015/06/24 Python
python常用知识梳理(必看篇)
2017/03/23 Python
Python实现自动为照片添加日期并分类的方法
2017/09/30 Python
python2.7+selenium2实现淘宝滑块自动认证功能
2018/02/24 Python
python和opencv实现抠图
2018/07/18 Python
解决使用export_graphviz可视化树报错的问题
2019/08/09 Python
Python实现不规则图形填充的思路
2020/02/02 Python
python里反向传播算法详解
2020/11/22 Python
中国旅游网站:同程旅游
2016/09/11 全球购物
英国羊绒服装购物网站:Pure Collection
2018/10/22 全球购物
爱尔兰最大的体育零售商:Life Style Sports
2019/06/12 全球购物
教师简历自我评价
2014/02/03 职场文书
职工趣味运动会方案
2014/02/10 职场文书
党员公开承诺书范文
2014/03/25 职场文书
党建工作先进材料
2014/05/02 职场文书
高一新生军训方案
2014/05/12 职场文书
开学第一天的感想
2015/08/10 职场文书
节约用水广告语60条
2019/11/14 职场文书