Hi Ivo,

The only problem is that uniroot() actually uses Brent's algorithm, and is based on the C code from netlib (there is also Fortran code available - zeroin.f).  Brent's algorithm is more sophisticated than a simple bisection approach that you have vectorized.  It combines bisection and secant methods along with some quadratic interpolation.  The idea behind this hybrid approach (which by the way is a fundamental theme in all of numerical analysis) is to get faster convergence while not sacrificing the global convergence of bisection.

So, the existing uniroot() will not be deprecated unless you can vectorize Brent's hybrid root-finding approach!

Best,
Ravi

Ravi Varadhan, Ph.D.
Assistant Professor
The Center on Aging and Health
Division of Geriatric Medicine & Gerontology
Johns Hopkins University
rvaradhan@jhmi.edu<mailto:rvaradhan@jhmi.edu>
410-502-2619


	[[alternative HTML version deleted]]


