Nothing
## binom_power <- function(n, p_alt, alpha = 0.05) {
## p_null <- 0.50
## # reject if X <= critical_lo or X >= critical_hi
## critical_lo <- qbinom( alpha/2, size = n, prob = p_null) - 1
## critical_hi <- qbinom(1 - alpha/2, size = n, prob = p_null)
## pbinom(critical_lo, n, p_alt) + 1 - pbinom(critical_hi - 1, n, p_alt)
## }
##
## PowerData <-
## expand.grid(n = seq(5, 10000, by = 5), p_alt = c(0.52, 0.55, 0.60)) %>%
## mutate(
## power = binom_power(n, p_alt),
## plab = paste("alt prob =", as.character(p_alt))
## )
##
## gf_line(power ~ n | plab, data = PowerData, size = 1) %>%
## gf_labs(y = "power", x = "number of coin tosses") %>%
## gf_lims(y = c(0, 1.1))
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.