如何使用方程式中的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() 

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