We might want to consider the use of `mem.save=TRUE` option in `glmnetPlus::glmnet()` for gaussian family. https://github.com/junyangq/glmnetPlus/pull/5/files