glmregNB function

fit a negative binomial model with lasso (or elastic net), snet and mnet regularization