mean - scipy实现的置信区间与数学公式不符
mean - scipy实现的置信区间与数学公式不符
提问于 2019-06-19 05:20:37
均值的置信区间具有以下解析解：
假设我的数据集是正态分布的，并且我不知道总体标准差，我可以使用t-score来计算均值的CI。所以我就这么做了：
from scipy import stats
import numpy as np
arr = np.array([4, 4, 1, 6, 6, 8, 1, 2, 3, 2, 2, 3, 4, 7, 6, 8, 0, 2, 8, 6, 5])
alpha = 0.05 # significance level = 5%
df = len(arr) - 1 # degress of freedom = 20
t = stats.t.ppf(1 - alpha/2, df) # 95% confidence t-score = 2.086
s = np.std(arr, ddof=1) # sample standard deviation = 2.502
n = len(arr)
lower = np.mean(arr) - (t * s / np.sqrt(n))
upper = np.mean(arr) + (t * s / np.sqrt(n))
print((lower, upper))
>>> (3.0514065531195387, 5.329545827832843)
print(stats.t.interval(1 - alpha/2, df, loc=np.mean(arr), scale=s / np.sqrt(n)))
>>> (2.8672993716475763, 5.513653009304806)并且使用该方程手动计算的间隔I与CI的scipy实现不一致。这个错误来自哪里？
回答 1
Stack Overflow用户
回答已采纳
发布于 2019-06-19 06:25:15
你的显著性水平是0.05，所以置信水平是0.95。将该值传递给stats.t.interval。不要除以2；函数会帮你除以2：
In [62]: print(stats.t.interval(1 - alpha, df, loc=np.mean(arr), scale=s / np.sqrt(n)))
(3.0514065531195387, 5.329545827832843)页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/56660487
复制
![[开发技巧]·Python实现信号滤波（基于scipy）](https://ask.qcloudimg.com/http-save/yehe-2144603/wyoyax4pnv.png)
相似问题
scipy.signal.spectrogram输出与预期不符
Scipy威布尔参数置信区间
如何实现数学公式
scipy.weave.inline与数学库的工作不符合预期
混淆数学公式的实现
添加站长 进交流群
领取专属 10元无门槛券
AI混元助手 在线答疑
关注 腾讯云开发者公众号
洞察 腾讯核心技术
剖析业界实践案例
社区富文本编辑器全新改版！诚邀体验～
全新交互，全新视觉，新增快捷键、悬浮工具栏、高亮块等功能并同时优化现有功能，全面提升创作效率和体验
腾讯云开发者
