[R] Constrained least squares

Ben Bolker bolker at ufl.edu
Mon Mar 13 14:23:09 CET 2006

Domenico Vistocco <vistocco <at> unicas.it> writes:

> Is there a function in R for constrained linear least squares?
> I used the matlab function LSQLIN: my aim is to obtain
> non-negative regression coefficients which sum 1.
> Thanks in advance,
> domenico vistocco

  I haven't tried it, but it looks like constrOptim()
[plus your own objective function that calculates sums of squares,
something like

f <- function(param1,param2) {
   y.exp <- ... [function of param1 and param2]

(or substitute -sum(dnorm(y,mean=y.exp,sd=param3,log=TRUE)) for
the last line)

More information about the R-help mailing list