powered by
symmetrize(a)
jj <- matrix(rnorm(50),10,5) X <- crossprod(jj,jj) # X has a Wishart distribution (and in # particular is positive definite) chol(X) symmetrize(chol(X))
Run the code above in your browser using DataLab