Sample random values from the multivariate normal distribution

Usage

sample_mvnorm(mu, sigma)

Arguments

mu

mean vector

sigma

covariance matrix

Samples multivariate normal variables by multiplying univariate random normal variables by the cholesky decomposition of the covariance matrix.

If mu is length 1 then just uses rnorm instead.