[R] linear constraint optim with bounds/reparametrization

[Thomas Lumley]

On Mon, 9 Aug 2004, Kahra Hannu wrote:

> >1) constrOptim does not work in this case because it only fits inequality
> >constraints, ie A%*%theta > =  c
>                           --- I was struggling with the same problem a
> few weeks ago in the portfolio optimization context. You can impose
> equality constraints by using inequality constraints >= and <=
> simultaneously. See the example bellow.
>

Ick. You do not want to use constrOptim for equality constraints.
constrOptim is a log-barrier interior-point method, meaning that it adds
a multiple of log(A%*%theta-c) to the objective function. This is a really
bad idea as a way of faking equality constraints.

Use Lagrange multipliers and optim.

	-thomas