Modeling Confidence Interval (CI) in R: How?

I was wondering how I can verify, through simulations in R, that the 95% confidence interval obtained from a binomial test with 5 successes in 15 trials when TRUE p = .5 has 95% coverage "Probability" in the long run?

Here is the 95% CI for such a test using R (how can the following CI be shown to have 95% coverage in the long run if TRUE p = .5):

as.numeric(binom.test(x = 5, n = 15, p = .5)[[4]])
# > [1] 0.1182411 0.6161963 (in the long-run 95% of the time, ".5" is contained within these
#                            two numbers, how to show this in R?)

      

0


source to share


1 answer


Something like that?

fun <- function(n = 15, p = 0.5){
    x <- rbinom(1, size = n, prob = p)
    res <- binom.test(x, n, p)[[4]]
    c(Lower = res[1], Upper = res[2])
}

set.seed(3183)
R <- 10000
sim <- t(replicate(R, fun()))

      



Note that binom.test

when called with 5 hits, 15 probes and p = 0.5 always return the same value, hence the call rbinom

. The number of successes will vary. We can calculate the proportion of cases when p

is between Lower

and Upper

.

cov <- mean(sim[,1] <= .5 & .5 <= sim[,2])

      

+1


source







All Articles