confint.hci <-
function(object, parm = NULL, level = 0.95, ...) {
if (!any(class(object) == 'hci')) stop("Object is not of class hci")
parm <- NULL
a <- (1 - level) / 2
a <- c(a, 1 - a)
fac <- qnorm(a)
ci <- object$concentration_index + fac * sqrt(object$variance)
return(ci)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.