How do I modify uniroot function to return .0001 if error ?

Hans W Borchers · 2011-04-04T03:25:10Z

eric aol.com> writes: > > I am calling the uniroot function from inside another function using these > lines (last two lines of the function) : > > d return(d$root) > > The problem is that on occasion there's a problem with the values I'm > passing to uniroot. In those instances uniroot stops and sends a message > that it can't calculate the root because f.upper * f.lower is greater than > zero. All I'd like to do in those cases is be a

Hans W Borchers

Sun, Apr 3, 2011 8:25 PM

eric <ericstrom <at> aol.com> writes:

Do not modify uniroot(). Use 'try' or 'tryCatch', for example

    	e <- try( d <- uniroot(k, c(.001, 250), tol=.05), silent = TRUE )
	if (class(e) == "try-error") {
		return(0.0001)
	} else {
		return(d$root)	
	}

--Hans Werner

Thread (2 messages)

eric How do I modify uniroot function to return .0001 if error ? Apr 3 Hans W Borchers How do I modify uniroot function to return .0001 if error ? Apr 3