I think pcls will do what you want. You'll need to explicitly set up the
constraint matrix, since pcls is designed to deal with general linear
inequality constraints (i.e. it will do more than just non-negativity).
Also the penalty matrix/matrices for pcls would be t(L)%*%L, I think.
btw, there's
Hi all,
I am trying to do an inversion of electromagnetic data with non-negative
least squares method (Tikhonov regularisation) and have got it
programmed in S-Plus. However I am trying to move all my scripts from
S-Plus to R.
Is there an equivalent to nnls.fit in R?
I think this can be done wi