Internal function to fit a GLM with lasso (or elastic net), snet and mnet regularization
Useful links