## File Name: confint.mlnormal.R
## File Version: 0.08
###########################################
# confidence interval
confint.mlnormal <- function( object, parm, level=.95, ... )
{
c1 <- coef.mlnormal(object)
v1 <- vcov.mlnormal(object)
if ( ! missing(parm) ){
c1 <- c1[parm]
v1 <- v1[ parm, parm]
}
q1 <- ( 1 - level ) / 2
q2 <- 1 - ( 1 - level ) / 2
quant <- stats::qnorm(q2)
se <- sqrt( diag(v1) )
res <- data.frame( "a"=c1 - quant*se, "b"=c1 + quant*se )
colnames(res)[1] <- paste0( round( 100*q1,1 ), " %")
colnames(res)[2] <- paste0( round( 100*q2,1 ), " %")
rownames(res) <- names(c1)
return(res)
}
###############################################
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.