Produce samples from a multivariate normal density using the Cholesky decomposition

rmvtn(n, mu, sigma)

Arguments

n

number of samples

mu

mean vector

sigma

covariance matrix

(source: http://gallery.rcpp.org/articles/simulate-multivariate-normal/)