Nothing
a0_gibbs_rate <-
function(thetas,covmat,a0){
cholcov <- chol(covmat)
invcholcov <- MASS::ginv(cholcov)
tmp <- (1/2)*(colSums((crossprod(invcholcov,thetas))^2))
return(sum(tmp)/a0)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.