Sample from the multivariate normal distribution
rmvnorm(mu, Sig2)
mu
Sig2
a p x 1 vector
Peter Hoff
rmvnorm( c(0,0,0),diag(rep(3,1)) )