## File Name: mcmc_vcov.R
## File Version: 0.13
#--- variance covariance matrix
mcmc_vcov <- function( mcmcobj, exclude="deviance" )
{
mcmcobj <- mcmc_extract_samples_first_chain(mcmcobj=mcmcobj)
mcmcobj <- mcmcobj[, ! ( colnames(mcmcobj) %in% exclude ) ]
res <- stats::var(mcmcobj)
colnames(mcmcobj) -> colnames(res) -> rownames(res)
return(res)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.