data(SimData)
x <- sim.data$x
classes <- sim.data$classes
lambda <- msgl.lambda.seq(x, classes, alpha = .5, d = 50, lambda.min = 0.05)
fit <- msgl(x, classes, alpha = .5, lambda = lambda)
# the nonzero parameters of model 1, 10 and 25
parameters(fit)[c(1,10,25)]
# count the number of nonzero parameters in each model
sapply(parameters(fit), sum)
Run the code above in your browser using DataLab