Yes, Bill almost always has helpful ideas. Just a comment: If indeed the process is gobbling up too much memory, that might indicate a problem with your function or implementation. I defer to real experts on this, however.
Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Tue, Dec 2, 2014 at 1:59 PM, Ramiro Barrantes <ram...@precisionbioassay.com> wrote: > Thanks so much for your reply. I am using try but nlme never returns!! and > I think the process is getting killed by the system as it is taking over all > the memory. However, I do like William Dunlap's idea of using > R.utils::withTimeout to limit the time. > > Thanks again for your help! > ________________________________________ > From: Bert Gunter [gunter.ber...@gene.com] > Sent: Tuesday, December 02, 2014 4:30 PM > To: Ramiro Barrantes > Cc: r-help@r-project.org > Subject: Re: [R] How to have NLME return when convergence not reached > > ?try > Or > ?tryCatch > > Bert > > Sent from my iPhone -- please excuse typos. > >> On Dec 2, 2014, at 12:57 PM, Ramiro Barrantes <ram...@precisionbioassay.com> >> wrote: >> >> Hello, >> >> I am trying to fit many hundreds of simulated datasets using NLME (it's all >> in a big loop in R). Some don't seem to converge. I am working on >> addressing the issues by perhaps adjusting my simulation, or tweaking >> iteration steps in nlme, etc. However, when it doesn't converge, NLME just >> hangs, and my program either stalls for hours/days or takes over the >> computer memory and everything crashes eventually. Is there a way to tell >> nlme to stop when it doesn't seem to be converging somehow? I have been >> looking at the parameters in nlmeControl() but see nothing obvious. >> >> Thanks in advance, >> >> Ramiro >> >> [[alternative HTML version deleted]] >> >> ______________________________________________ >> R-help@r-project.org 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. ______________________________________________ R-help@r-project.org 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.
