Lasso fit with non-negative coefficients

It seems like GLMNet.jl supports the constraints keyword that you can use.