Как использовать 95% доверительный интервал в r, используя уравнение для 95% доверительного интервала?

Question

Как бы я использовал уравнение 95% доверительного интервала = (среднее значение по совокупности) +/- 1,96 x (стандартная ошибка среднего по совокупности) в программе r для определения доверительного интервала? Мой профессор биометрии хочет, чтобы мы использовали уравнение, но мы не можем понять, как ввести его в программу R.

Answer 1

Поскольку вы не предоставляете пример данных, вот пример использования набора данных iris sample. В частности, давайте вычислим 95% ДИ среднего значения ширины чашелистика.

1. Построим распределение Sepal.Width

   library(ggplot)
   ggplot(iris, aes(Sepal.Width)) + geom_density()
   

2. 95% доверительный интервал определяется теми наблюдениями, которые лежат в пределах 1,96 единиц стандартного отклонения от среднего значения (при условии нормального нормального распределения).

   CI <- with(iris, mean(Sepal.Width) + c(-1, 1) * 1.96 * sd(Sepal.Width))
   CI
   #[1] 2.203035 3.911631
   

3. Мы подтверждаем, что действительно 95% наблюдений лежат в пределах CI

   with(iris, sum(Sepal.Width >= CI[1] & Sepal.Width <= CI[2]) / length(Sepal.Width))
   #[1] 0.9466667