平均二項分布の信頼区間R

Question

私は二項分布に従う観測値を持つ4つの異なるデータフレームを得ており、それぞれについて第2列の平均値（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は、2項の信頼区間を計算する方法を提供します。使用可能なすべてのメソッドを使用するか、または単一のメソッドを指定するかを選択できます。

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