Nothing
Sims <-
expand.grid(
n = c(35, 100),
pi = c(0.2, 0.3, 0.4, 0.5),
method = c("Wald", "Wilson", "score"),
rep = 1:2000) %>%
group_by(n, pi, method, rep) %>%
do(
confint(
binom.test(rbinom(1, .$n, .$pi), n = .$n, ci.method = .$method)))
Sims %>%
group_by(n, pi, method) %>%
summarise(cover = prop(lower <= pi & pi <= upper))
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.