@pd: you know that a System of equations with more variables than equations is
always solvable
and if a unique solution is desired one of mimimal norm can be used.
According to "Methods for nonlinear least squares problems" by Madsen, Nielsen
and Tingleff the LM-algorithm
solves Systems of the form
[J(x)'J(x)+\mu*I]x=...
with \mu>0 so that the Matrix on the left is always positive definite,
especially nonsingular.
______________________________________________
[email protected] mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
