R/iconfianza.percentil.kc.R

Defines functions iconfianza.percentil.kc

Documented in iconfianza.percentil.kc

#' confidence interval por the median using KC Method
#'
#' @keywords internal
iconfianza.percentil.kc <- function(datos, q = 0.50, nivel = 0.95, ic = T, colas = 2, use.t = F) {
  x <- datos[!is.na(datos)]
  n <- length(x)
  if (n != 0) {
    # pnor<-qnorm((1-nivel)/colas,lower.tail=FALSE)
    if (use.t) pnor <- qt(1 - (1 - nivel) / colas, n - 1) else pnor <- qnorm(1 - (1 - nivel) / colas)
    j <- max(1, floor(n * q - pnor * sqrt(n * q * (1 - q))))
    k <- min(n, ceiling(n * q + pnor * sqrt(n * q * (1 - q))))
    if (floor(n * q) == n * q) {
      p <- (sort(x)[n * q] + sort(x)[n * q + 1]) / 2
    } else {
      p <- sort(x)[1 + floor(n * q)]
    }
    med <- p
    l.i <- sort(x)[j]
    l.s <- sort(x)[k]
  } else {
    med <- NA
    l.i <- NA
    l.s <- NA
  }
  if (ic) iconfres <- c(l.i, med, l.s) else iconfres <- rep(med, 3)
  iconfres
}
lozalojo/mem documentation built on July 7, 2023, 10 a.m.