R/oddsratioci.mp.R

In shearer/PropCIs: Various Confidence Interval Methods for Proportions

oddsratioci.mp <-
function(b,c,conf.level)
{
  z  <- qchisq(conf.level,1)
  A <- b + c + z
  B <- 2*c + z
  C <- c^2/(b+c)
  l <- (B - sqrt(B^2-4*A*C))/(2*A)
  u <- (B + sqrt(B^2-4*A*C))/(2*A)
  ll <- l/(1-l)
  ul <- u/(1-u)
  cint <- c(ll, ul)
  attr(cint, "conf.level") <- conf.level
  rval <- list(conf.int = cint)
  class(rval) <- "htest"
  return(rval)
}