# 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))
k <- 5
v_up <- est_up(E, k = k)
# }
Run the code above in your browser using DataLab