R/WilsonCI.R

Defines functions WilsonCI

Documented in WilsonCI

WilsonCI <-
function(x,n,alpha)  {
  phat <- x/n
  nz2 <- n + dnorm(alpha/2)^2
  firstterm <- phat*n/nz2
  secondterm <- 0.5*dnorm(alpha/2)/nz2
  commonterm <- phat*(1-phat)/n
  commonterm <- commonterm * (n^2) * (dnorm(alpha/2)^2) / (nz2^2)
  commonterm <- commonterm + (0.25 * (dnorm(alpha/2)^4) )/ (nz2^2)
  commonterm <- sqrt(commonterm)
  return(c(firstterm+secondterm-commonterm,firstterm+secondterm+commonterm))
}

Try the ACSWR package in your browser

Any scripts or data that you put into this service are public.

ACSWR documentation built on May 29, 2017, 8:41 p.m.