Grouplasso

uses glmnet(..., family = "mgaussian") to fit a multi-response regression, suitable for VAR models.  glmnet uses Group Lasso penalty instead of plain lasso penalty for "mgaussian" families