Have you plotted f(x)? Plot(s) might help reveal why "the solution is very unstable." If f is a function of a univariate x, this is trivial. If f is a function of a bivariate x, use something like contour or persp. If higher dimensions, I might use something like

          fit <- optim(... hessian=TRUE)
          eigen(fit$hessian)

          hope this helps.
          spencer graves

Martyn Plummer wrote:

On Wed, 2005-06-08 at 17:08 +0200, Fredrik Thuring wrote:

Hi!

Im need a function that solves the equation f(x) = 0 (i.e. the root of
the function) when f is a nonlinear function. Is there any? Ive tried nlm
and optim on the square of the function but the solution is very unstable.


Your colleague Martin Englund asked the same question.  If x is scalar,
you can use uniroot(). If it's a polynomial then polyroot() finds all
the zeros. I hope that helps.

Martyn

______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Reply via email to