# NOT RUN {
n <- 30
q <- 100
Sigma <- Simu_Sigma(q = q, diag = FALSE, equal = TRUE)
Matrix::image(Sigma)
E <- matrix(rnorm(n * q), ncol = q) %*% chol(as.matrix(Sigma))
random_eigen <- PA(E, times = 10)
# }
Run the code above in your browser using DataLab