#' @importFrom stats rnorm
rmvnorm <- function (n = 1, mean, var) {
p <- length(mean)
## if (!all(dim(var)==c(p,p)))
## pStop_("length of ",sQuote("mean")," must equal the dimension of the square matrix ",sQuote("var"),".")
cf <- t(chol(var))
matrix(mean,p,n)+cf%*%matrix(rnorm(p*n),p,n)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.