#' confidence interval for a point (using arithmetic mean)
#'
#' @keywords internal
iconfianza.x <- function(datos, nivel = 0.95, ic = T, colas = 2, use.t = F) {
datos[datos == -Inf] <- NA
datos[datos == Inf] <- NA
n <- sum(!is.na(datos))
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)
med <- mean(datos, na.rm = T)
std <- sqrt(var(datos, na.rm = T))
l.i <- med - pnor * std
l.s <- med + pnor * std
} else {
med <- NA
l.i <- NA
l.s <- NA
}
if (ic) iconfres <- c(l.i, med, l.s) else iconfres <- rep(med, 3)
iconfres
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.