R 均值的二项式置信区间

Question

我有 4 个不同的 data.frames，它们的观察结果遵循二项分布，我需要为每个数据帧计算与第二列的平均值相关的置信区间 (Flow)。

在Success 和total number of trials = 85 列中报告成功次数。

如何计算置信区间？ 我怎样才能用 R 做到这一点？

这是我的 data.frames 的示例：

df1 <- read.table(text = 'Flow Success
725.661   4
25.54     4
318.481   4
230.556   4
2.823     3
12.6      3
9.891     3
11.553    1', header = TRUE)

> mean(df1$Flow)
[1] 167.1381


df2 <- read.table(text = 'Flow Success
725.661    3
25.54      3
318.481    3
230.556    2
2.823      2
12.6       1', header = TRUE)

> mean(df2$Flow)
[1] 219.2768

df3 <- read.table(text = 'Flow Success
725.661     2
25.54       2
318.481     1', header = TRUE)

> mean(df3$Flow)
[1] 356.5607

df4 <- read.table(text = 'Flow Success
725.661    2
25.54      2', header = TRUE)

> mean(df4$Flow)
[1] 375.6005

我需要计算上述均值的置信区间。

如果需要，我可以为您提供有关数据的更多信息。

感谢任何愿意帮助我的人。

Answer 1

包binom 提供了计算二项式置信区间的方法。可以选择使用所有可用的方法，或指定一种方法。

x 给出成功次数，n 给出伯努利试验次数。

library(binom)

binom.confint(x = 5, n = 10)
          method x  n mean     lower     upper
1  agresti-coull 5 10  0.5 0.2365931 0.7634069
2     asymptotic 5 10  0.5 0.1901025 0.8098975
3          bayes 5 10  0.5 0.2235287 0.7764713
4        cloglog 5 10  0.5 0.1836056 0.7531741
5          exact 5 10  0.5 0.1870860 0.8129140
6          logit 5 10  0.5 0.2245073 0.7754927
7         probit 5 10  0.5 0.2186390 0.7813610
8        profile 5 10  0.5 0.2176597 0.7823403
9            lrt 5 10  0.5 0.2176212 0.7823788
10     prop.test 5 10  0.5 0.2365931 0.7634069
11        wilson 5 10  0.5 0.2365931 0.7634069

binom.confint(x = 5, n = 10, method = "exact")
  method x  n mean    lower    upper
1  exact 5 10  0.5 0.187086 0.812914