Nothing
# Compute Clopper-Pearson interval for a binomial proportion
#
# Copyright (C) 2004 Kaspar Rufibach, 2016 Leonhard Held
clopperPearson <- function(x, n, conf.level = 0.95)
{
stopifnot(is.wholenumber(x), is.wholenumber(n), (x<=n), (n>=1), conf.level<1, conf.level>0)
a <- 1 - conf.level
if (x == 0){
ll <- 0
ul <- 1 - (a / 2) ^ (1 / n)}
else if (x == n){
ll <- (a / 2) ^ (1 / n)
ul <- 1
}
else {
ll <- 1/(1 + (n - x + 1) / (x * qf(a / 2, 2 * x, 2 * (n - x + 1))))
ul <- 1/(1 + (n - x) / ((x + 1) * qf(1 - a / 2, 2 * (x + 1), 2 * (n - x))))
}
res <- c("lower" = ll, "prop" = x / n, "upper" = ul)
return(res)
}
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.