In R, we can use the following code to find the coverage probability for two sided binomial test.
library(binom) binom.coverage(0.1, 22, conf.level = 0.95, method="bayes")$coverage
So how can we get a coverage probability for one sided binomial test in R? (let's say we only care about the upper limit.
Thanks for your thought!