powered by
rwishart(N, df, Sigma)
rmultnorm
rgamma
x <- matrix(rnorm(100), 50, 2) XX <- crossprod(x) tmp <- rwishart(1, 50, XX)
Run the code above in your browser using DataLab