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递归删除指定目录及其所有内容的方法
Jan 13 Python
python itchat实现微信自动回复的示例代码
Aug 14 Python
Python实现的随机森林算法与简单总结
Jan 30 Python
Django 路由系统URLconf的使用
Oct 11 Python
python运行时强制刷新缓冲区的方法
Jan 14 Python
python中字符串数组逆序排列方法总结
Jun 23 Python
详解pandas数据合并与重塑(pd.concat篇)
Jul 09 Python
python增加图像对比度的方法
Jul 12 Python
安装PyInstaller失败问题解决
Dec 14 Python
Django 解决上传文件时,request.FILES为空的问题
May 20 Python
什么是Python中的匿名函数
Jun 02 Python
python 爬取小说并下载的示例
Dec 07 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 定界符 使用技巧
2009/06/14 PHP
解决php接收shell返回的结果中文乱码问题
2014/01/23 PHP
php通过前序遍历树实现无需递归的无限极分类
2015/07/10 PHP
php验证码实现代码(3种)
2015/09/07 PHP
php封装的图片(缩略图)处理类完整实例
2016/10/19 PHP
ThinkPHP中调用PHPExcel的实现代码
2017/04/08 PHP
PHP 图片处理
2020/09/16 PHP
document.onreadystatechange事件的用法分析
2009/10/17 Javascript
JavaScript高级程序设计 阅读笔记(十七) js事件
2012/08/14 Javascript
jquery ajax请求实例深入解析
2012/11/26 Javascript
JS中prototype关键字的功能介绍及使用示例
2013/07/21 Javascript
javascript定义变量时有var和没有var的区别探讨
2014/07/21 Javascript
Javascript window对象详解
2014/11/12 Javascript
Nodejs如何搭建Web服务器
2016/03/28 NodeJs
javascript简单判断输入内容是否合法的方法
2016/05/11 Javascript
移动端基础事件总结与应用
2017/01/12 Javascript
Javascript DOM事件操作小结(监听鼠标点击、释放,悬停、离开等)
2017/01/20 Javascript
将input框中输入内容显示在相应的div中【三种方法可选】
2017/05/08 Javascript
NodeJs使用Mysql模块实现事务处理实例
2017/05/31 NodeJs
Windows下快速搭建NodeJS本地服务器的步骤
2017/08/09 NodeJs
jQuery.Sumoselect插件实现下拉复选框效果
2017/11/09 jQuery
vue cli3.0结合echarts3.0与地图的使用方法示例
2019/03/26 Javascript
jquery+php后台实现省市区联动功能示例
2019/05/23 jQuery
vue中使用GraphQL的实例代码
2019/11/04 Javascript
vue+vant-UI框架实现购物车的复选框全选和反选功能
2019/11/05 Javascript
vue 解决路由只变化参数页面组件不更新问题
2019/11/05 Javascript
[02:33]2018DOTA2亚洲邀请赛赛前采访——LGD
2018/04/04 DOTA
python二分法实现实例
2013/11/21 Python
Python中if __name__ == "__main__"详细解释
2014/10/21 Python
小小聊天室Python代码实现
2016/08/17 Python
当当网软件测试笔试题
2015/11/24 面试题
环卫工作个人总结
2015/03/04 职场文书
民事答辩状格式范文
2015/05/21 职场文书
2015年小学语文教师工作总结
2015/10/23 职场文书
总结高并发下Nginx性能如何优化
2021/11/01 Servers
教你如何让spark sql写mysql的时候支持update操作
2022/02/15 MySQL