statistics - 计算 Julia 的置信区间

Question

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

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

Answer 1

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

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

Answer 2

有一个HypothesisTests.jl包可以为你做这件事。

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