Dies ist das von prop.test geschätzte Konfidenzintervall
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
Ich habe versucht, es zu reproduzieren, indem ich den Code unter prop.test gelesen habe. Hier ist ein vereinfachter Weg, um diese beiden Grenzen zu erreichen
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
Können Sie mir erklären, warum die zugrunde liegende Erfolgswahrscheinlichkeit (p) in Zeile 5 verwendet wird? oder vielleicht könnten Sie vorschlagen, wo ich weitere Informationen zu dieser YATES-Korrektur finden kann, die sich auf die SCHÄTZUNG auswirkt.
Vielen Dank
quelle
binom
Paket von R enthält auch den Agresti-Coull CI.