基本上,我有一个包含来自总体的随机样本的数组。我正在尝试计算真实总体均值的95%置信区间。
我已经查看了两个库,但他们似乎只支持双样本t检验。 (我不擅长统计,所以我可能会遗漏一些东西。)
答案 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))