计算 Julia 的置信区间答案

【问题标题】：Calculate confidence interval in Julia计算 Julia 的置信区间
【发布时间】：2013-12-23 21:32:00
【问题描述】：

基本上，我有一个数组，其中包含来自人群的随机样本。我正在尝试计算真实总体均值的 95% 置信区间。

我研究了两个库，但它们似乎只支持两个样本的 t 检验。（我不擅长统计，所以我可能遗漏了一些东西。）

【问题讨论】：

【解决方案1】：

有一个HypothesisTests.jl 包可以为您完成这项工作。

using Distributions, HypothesisTests
data = rand(Normal(10,10),100)
ci(OneSampleTTest(data))

【讨论】：

我认为 ci 已被弃用，取而代之的是 confint。

【解决方案2】：

这可能包含在某处的函数中，但这是假设正常人口的直接方法：

using Distributions
function t_test(x; conf_level=0.95)
    alpha = (1 - conf_level)
    tstar = quantile(TDist(length(x)-1), 1 - alpha/2)
    SE = std(x)/sqrt(length(x))

    lo, hi = mean(x) + [-1, 1] * tstar * SE
    "($lo, $hi)"
end

【讨论】：

目前我认为这是唯一的方法。
能够定义像nbviewer.ipython.org/gist/StefanKarpinski/b8fe9dbb36c1427b9f22 的单元格56 这样的中缀± 运算符会很好，但我认为笔记本是在julia 的私有分支上完成的。