powered by
Simulates a random Wishart-distributed matrix
rwish(S0, nu = dim(S0)[1] + 2)
a positive definite matrix
a positive integer
Peter Hoff
## The expectation is S0*nu S0<-rwish(diag(3)) SS<-matrix(0,3,3) for(s in 1:1000) { SS<-SS+rwish(S0,5) } SS/s S0*5
Run the code above in your browser using DataLab