[S] uniroot -- doesn't work recursively
Adelchi Azzalini <adelchi.azzalini@jrc.it> writes:
Martin Maechler wrote: ) However, note that R currently has the same problem as S-plus: ) Uniroot() doesn't work reliably, recursively. ) When you found it to be better, you were just lucky. Right. I realised later on that R seemed to behave better but it could not work properly since it calls a Fortran (66? 77?) routine. It just happen to produce less obviously wrong answers than Splus. In fact more tests confirm wrong answers, but not as non-sensical as those from Splus.
The development version now has a reentrant uniroot(). Those of you outside R-core who are interested can get it in the next development snapshot. (and perhaps remind me what is the proper typecast to avoid the following error message...) optimize.c: In function `do_zeroin': optimize.c:205: warning: passing arg 3 of `zeroin' from incompatible pointer type Binaries are a bit further into the future, I'm afraid.
O__ ---- Peter Dalgaard Blegdamsvej 3 c/ /'_ --- Dept. of Biostatistics 2200 Cph. N (*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918 ~~~~~~~~~~ - (p.dalgaard@biostat.ku.dk) FAX: (+45) 35327907 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-devel-request@stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._