lsie_ln

solve linear least square problem (min ||A%%x-b||)
with inequality constraints u%%x&gt;=co and equality constraints e%*%x=ce
Method:
reduce the pb to lsi_ln on the null-space of e

We solve non linear least squares problems with optional
equality and/or inequality constraints. Non linear iterations are
globalized with back-tracking method. Linear problems are solved by
dense QR decomposition from 'LAPACK' which can limit the size of
treated problems. On the other side, we avoid condition number
degradation which happens in classical quadratic programming approach.
Inequality constraints treatment on each non
linear iteration is based on 'NNLS' method (by Lawson and Hanson).
We provide an original function 'lsi_ln' for solving linear least squares
problem with inequality constraints in least norm sens. Thus if Jacobian of
the problem is rank deficient a solution still can be provided.
However, truncation errors are probable in this case.
Equality constraints are treated by using a basis of Null-space.
User defined function calculating residuals must return a list having
residual vector (not their squared sum) and Jacobian. If Jacobian is
not in the returned list, package 'numDeriv' is used to calculated
finite difference version of Jacobian. The 'NLSIC' method was fist
published in Sokol et al. (2012) <doi:10.1093/bioinformatics/btr716>.

Serguei Sokol

nlsic

Non Linear Least Squares with Inequality Constraints

lsie_ln function

<dl><dt>a</dt>
<dd>dense matrix A or its QR decomposition</dd>
<dt>b</dt>
<dd>right hand side vector</dd>
<dt>u</dt>
<dd>dense matrix of inequality constraints</dd>
<dt>co</dt>
<dd>right hand side vector of inequality constraints</dd>
<dt>e</dt>
<dd>dense matrix of equality constraints</dd>
<dt>ce</dt>
<dd>right hand side vector of equality constraints</dd>
<dt>rcond</dt>
<dd>maximal condition number for determining rank deficient matrix</dd>
<dt>mnorm</dt>
<dd>matrix to be used to minimize the norm ||mnorm%*%x|| in
case of under-dermined system. If mnorm is NULL (default), it is</dd>
<dt>x0</dt>
<dd>vector from which minimal norm should be searched when required</dd></dl>

Arguments

Linear Least Squares problem with inequality and equality constraints, least norm solution — lsie_ln

<dl>

<dt>a</dt>
<dd>dense matrix A or its QR decomposition</dd>


<dt>b</dt>
<dd>right hand side vector</dd>


<dt>u</dt>
<dd>dense matrix of inequality constraints</dd>


<dt>co</dt>
<dd>right hand side vector of inequality constraints</dd>


<dt>e</dt>
<dd>dense matrix of equality constraints</dd>


<dt>ce</dt>
<dd>right hand side vector of equality constraints</dd>


<dt>rcond</dt>
<dd>maximal condition number for determining rank deficient matrix</dd>


<dt>mnorm</dt>
<dd>matrix to be used to minimize the norm ||mnorm%*%x|| in
case of under-dermined system. If mnorm is NULL (default), it is</dd>


<dt>x0</dt>
<dd>vector from which minimal norm should be searched when required</dd>

</dl>

lsie_ln: Linear Least Squares problem with inequality and equality constraints, least norm solution

Description

Usage

Value

Arguments

See Also