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 相关文章推荐
Python3搜索及替换文件中文本的方法
May 22 Python
Python的条件语句与运算符优先级详解
Oct 13 Python
使用Django Form解决表单数据无法动态刷新的两种方法
Jul 14 Python
python实现二叉查找树实例代码
Feb 08 Python
浅谈Pandas中map, applymap and apply的区别
Apr 10 Python
Python中.join()和os.path.join()两个函数的用法详解
Jun 11 Python
python中文编码与json中文输出问题详解
Aug 24 Python
Python之lambda匿名函数及map和filter的用法
Mar 05 Python
Python 余弦相似度与皮尔逊相关系数 计算实例
Dec 23 Python
基于spring boot 日志(logback)报错的解决方式
Feb 20 Python
Tensorflow tf.tile()的用法实例分析
May 22 Python
Python 整行读取文本方法并去掉readlines换行\n操作
Sep 03 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变量引用的面试题
2010/08/08 PHP
php使用正则表达式提取字符串中尖括号、小括号、中括号、大括号中的字符串
2020/04/05 PHP
PHP生成图片缩略图类示例
2017/01/12 PHP
php实现的支付宝网页支付功能示例【基于TP5框架】
2019/09/16 PHP
script标签的 charset 属性使用说明
2010/12/04 Javascript
基于jquery的从一个页面跳转到另一个页面的指定位置的实现代码(带平滑移动的效果)
2011/05/24 Javascript
jQuery 联动日历实现代码
2012/05/31 Javascript
javascript模块化是什么及其优缺点介绍
2013/09/02 Javascript
关闭页面时window.location事件未执行的原因分析及解决方案
2014/09/01 Javascript
js电话号码验证方法
2015/09/28 Javascript
js实现贪吃蛇小游戏(容易理解)
2017/01/22 Javascript
微信小程序 弹幕功能简单实例
2017/02/14 Javascript
javaScript+turn.js实现图书翻页效果实例代码
2017/02/16 Javascript
在node中如何使用 ES6
2017/04/22 Javascript
javascript获取指定区间范围随机数的方法
2017/09/08 Javascript
angular指令笔记ng-options的使用方法
2017/09/18 Javascript
nodejs项目windows下开机自启动的方法
2017/11/22 NodeJs
axios进阶实践之利用最优雅的方式写ajax请求
2017/12/20 Javascript
JS实现十分钟倒计时代码实例
2018/10/18 Javascript
利用hasOwnProperty给数组去重的面试题分享
2018/11/05 Javascript
JavaScript静态作用域和动态作用域实例详解
2019/06/17 Javascript
利用python程序帮大家清理windows垃圾
2017/01/15 Python
Tensorflow卷积神经网络实例进阶
2018/05/24 Python
Python适配器模式代码实现解析
2019/08/02 Python
Python django框架输入汉字,数字,字符生成二维码实现详解
2019/09/24 Python
CSS3中颜色线性渐变实战
2015/07/18 HTML / CSS
写自荐信有哪些不宜?
2013/10/17 职场文书
企业消防安全制度
2014/02/02 职场文书
总账会计岗位职责
2014/03/13 职场文书
班主任工作经验交流材料
2014/05/13 职场文书
项目合作意向书模板
2014/07/29 职场文书
亲属关系公证书样本
2015/01/23 职场文书
乡镇一岗双责责任书
2015/01/29 职场文书
机修车间主任岗位职责
2015/04/08 职场文书
植树节新闻稿
2015/07/17 职场文书
windows安装 redis 6.2.6最新步骤详解
2022/04/26 Redis