if (requireNamespace("huge", quietly = TRUE)) {
library(huge)
set.seed(10010)
p <- 40 ; n <- 100
dat <- huge.generator(n, p, "hub", verbose=FALSE, v=.1, u=.3)
## Theoretical lambda max is the maximum abs value of the empirical covariance matrix
maxCov <- getMaxCov(dat$data)
lams <- getLamPath(maxCov, 5e-2*maxCov, len=40)
}
Run the code above in your browser using DataLab