Function rand returns a matrix containing pseudorandom values
drawn from the standard uniform distribution (using runif).
Function randn returns a matrix containing pseudorandom values
drawn from the standard normal distribution (using rnorm).
References
P. Carbonetto and M. Stephens (2012). Scalable variational
inference for Bayesian variable selection in regression, and its
accuracy in genetic association studies. Bayesian Analysis7,
73--108.