计算朱莉娅的置信区间

时间:2013-12-06 13:45:14

标签: statistics julia

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

我已经查看了两个库,但他们似乎只支持双样本t检验。 (我不擅长统计,所以我可能会遗漏一些东西。)

2 个答案:

答案 0 :(得分:5)

可能这在一个函数中被包含在某个地方,但是假设一个正常的人口,这是一个直接的方法:

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

答案 1 :(得分:3)

有一个HypothesisTests.jl包可以帮到你。

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