Доверительный интервал для нормального распределения - R

Question

Я хочу вычислить доверительный интервал для вектора нормально распределенных значений в R. Выполнение этого вручную не будет слишком большой проблемой, но, безусловно, для этого есть встроенная функция.Кажется, я не могу найти ни одного.Есть идеи?

Answer 1

Нет встроенной функции именно для этого и только для этой цели.Я думаю, причина в том, что это действительно очень легко вычислить вручную.Пусть, как и в ответе @ bobbel,

set.seed(1)
x <- rnorm(20)

Тогда один из способов извлечь интервал будет

t.test(x)$conf.int
# [1] -0.2368920  0.6179398
# attr(,"conf.level")
# [1] 0.95

, что не является плохой идеей, так как часто вы захотите сделатьпроверить в любом случае.Что касается одного только интервала, вам нужно определить свою собственную функцию, например,

normConfInt <- function(x, alpha = 0.05)
  mean(x) + qt(1 - alpha / 2, length(x) - 1) * sd(x) / sqrt(length(x)) * c(-1, 1)
normConfInt(x)
# [1] -0.2368920  0.6179398

Answer 2

Я не уверен, есть ли функция для этого, но простой способ сделать это - установить модель только с перехватом и извлечь из нее доверительный интервал с помощью функции confint.

set.seed(1)
x <- rnorm(20)
m <- lm(x~1)

confint(m)
                2.5 %    97.5 %
(Intercept) -0.236892 0.6179398