[R] nlm recursion problem

[Peter Dalgaard BSA]

John Logsdon <j.logsdon at lancaster.ac.uk> writes:

> nlm is called by many functions and packages so perhaps this is quite a
> common problem - what is the slickest way round it?  The solution required
> is just the intercept and slope of a line where there may be censored data
> so lm() won't do and neither will optimise() (it is only for 1 dimension
> and I require 2).

There really is no way around this, except either to rewrite the
maximiser in a language that allows reentrancy ( = C in practice) or
wait for someone to write a Fortran compiler that does. We've
previously done this for the solvers, and nlm appears to be under
scrutiny at the moment anyway, so perhaps...

-- 
   O__  ---- Peter Dalgaard             Blegdamsvej 3  
  c/ /'_ --- Dept. of Biostatistics     2200 Cph. N   
 (*) \(*) -- University of Copenhagen   Denmark      Ph: (+45) 35327918
~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk)             FAX: (+45) 35327907
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help 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-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._