Python求解正态分布置信区间教程
正态分布和置信区间
正态分布(Normal Distribution)又叫高斯分布,是一种非常重要的概率分布。其概率密度函数的数学表达如下:
置信区间是对该区间能包含未知参数的可置信的程度的描述。
使用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,就是本文最开始的那个公式
最后的输出图像如下,可以看到结果跟理论上的正太分布还是比较像的:
正态分布置信区间规律
函数曲线下68.268949%的面积在平均数左右的一个标准差范围内
函数曲线下95.449974%的面积在平均数左右两个标准差的范围内
函数曲线下99.730020%的面积在平均数左右三个标准差的范围内
函数曲线下99.993666%的面积在平均数左右四个标准差的范围内
以上这篇Python求解正态分布置信区间教程就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持脚本之家。
相关文章
详解Python3网络爬虫(二):利用urllib.urlopen向有道翻译发送数据获得翻译结果
这篇文章主要介绍了Python3网络爬虫(二):利用urllib.urlopen向有道翻译发送数据获得翻译结果,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧2019-05-05python tkiner实现 一个小小的图片翻页功能的示例代码
这篇文章主要介绍了python tkiner实现 一个小小的图片翻页功能,需要的朋友可以参考下2020-06-06python命令行执行脚本找不到模块ModuleNotFoundError问题
这篇文章主要介绍了python命令行执行脚本找不到模块ModuleNotFoundError问题,具有很好的参考价值,希望对大家有所帮助。如有错误或未考虑完全的地方,望不吝赐教2023-06-06Python复制Excel带有条件格式的单元格sheet实现步骤
这篇文章主要为大家介绍了Python复制Excel带有条件格式的单元格sheet实现步骤,有需要的朋友可以借鉴参考下,希望能够有所帮助,祝大家多多进步,早日升职加薪2023-07-07
最新评论