File: zclopperpearson.R

package info (click to toggle)
r-cran-epir 2.0.80%2Bdfsg-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 2,332 kB
  • sloc: makefile: 5
file content (12 lines) | stat: -rw-r--r-- 350 bytes parent folder | download | duplicates (2) 
1
2
3
4
5
6
7
8
9
10
11
12
 
zclopperpearson <- function(dat, conf.level){
  # From RSurveillance function binom.cp:
  a <- dat[,1]
  n <- dat[,2]
  p <- a / n
  
  tails <- 2
  low <- stats::qbeta((1 - conf.level) / tails, a, n - a + 1)
  upp <- stats::qbeta(1 - (1 - conf.level) / tails, a + 1, n - a)
  rval <- data.frame(est = p, lower = low, upper = upp)
  rval
}