これは、prop.testによって推定された信頼区間です
n <- 600; x <- 276; p <- 0.40
prop.test(x, n, p, alternative="two.sided", conf.level=0.95, correct=T)
95 percent confidence interval:
0.4196787 0.5008409
prop.testの下のコードを読んで、それを再現しようとしました。これらの2つの制限を取得する簡単な方法を次に示します
ESTIMATE <- x/n
YATES <- 0.5
conf.level <- 0.95
z <- qnorm((1 + conf.level)/2)
YATES <- min(YATES, abs(x - n * p))
z22n <- z^2/(2 * n)
p.c <- ESTIMATE + YATES/n
(p.c + z22n + z * sqrt(p.c * (1 - p.c)/n + z22n/(2 * n)))/(1 + 2 * z22n)
[1] 0.5008409
p.c <- ESTIMATE - YATES/n
(p.c + z22n - z * sqrt(p.c * (1 - p.c)/n + z22n/(2 * n)))/(1 + 2 * z22n)
[1] 0.4196787
成功の潜在的な確率(p)が5行目で使用されている理由を教えていただけますか?または、推定に影響を与えるこのYATES修正についての詳細情報をどこで見つけることができるかを提案できます。
ありがとうございました
binom
パッケージにはAgresti-Coull CIも含まれています。