Il s'agit de l'intervalle de confiance estimé par test prop.
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
J'ai essayé de le reproduire en lisant le code sous prop.test. Voici un moyen simplifié d'obtenir ces deux limites
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
Pouvez-vous m'expliquer pourquoi la probabilité de réussite sous-jacente (p) est utilisée à la ligne 5? ou peut-être pourriez-vous suggérer où puis-je trouver plus d'informations sur cette correction YATES qui affecte l'ESTIMATION.
Je vous remercie
la source
binom
package de R a également le CI Agresti-Coull.