set.seed(123)
n <- 100
p <- 5
L <- 5
Y <- matrix(rnorm(n), n, 1)
W <- matrix(rnorm(n * p * (L - 1)), n, p * (L - 1))
# Call qrglasso with default parameters
result <- qrglasso(Y = Y, W = W, p = p)
# Visualize the BIC results
plot(result)
Run the code above in your browser using DataLab