计算朱莉娅的置信区间

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