powered by
stein_eig(l, n)
p <- 5 n <- 10 S <- rWishart(1, n, diag(p))[,,1] l <- eigen(S)$val stein_eig(l, n)
Run the code above in your browser using DataLab