如何使用方程式中的95％置信区间在r中使用95％置信区间？

Question

在程序r中，如何使用95％置信区间=（人口平均值）+/- 1.96 x（总体平均值的标准误差）方程式来求解置信区间？我的生物识别学教授希望我们使用该方程，但我们不知道如何将其输入到程序R中。

Answer 1

由于您没有提供示例数据，因此以下是使用iris示例数据集的示例。具体来说，让我们计算出萼片宽度平均值的95％CI。

1. 让我们绘制分布Sepal.Width

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

1. 95％的置信区间是由位于均值周围的标准偏差（假设基本正态分布）在1.96单位内的那些观测值定义的。

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

2. 我们确认确实有95％的观察结果位于CI

   之内

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