Skip to content

Arima Models - Error and jump error

1 message · John C Nash

#
The error is what it says: Your program has asked for an approximation to a derivative
(you, or whoever wrote the package you are using,  didn't provide an analytic gradient, so
it's using an approximation), and the result was returned Inf.

This is VERY common -- the BFGS option of optim() is a variant of Roger Fletcher's 1970
variable metric code that I worked out with him in Jan 1976. It's a bit long in the tooth,
and does not include parameter bounds -- we only had about 8K maximum for program and data
at the time. That's K, not M or G.

Some options:
- if you can change to L-BFGS-B variant of optim() you can use box constraints, still with
numerical gradients. nlminb() also has them, as do several routines in optimx. Note that a
bound of 0 may still get you into trouble with numerical gradient approximations -- the
step can go over the "cliff" because the bounds constraints don't always get checked for
the numerical  (f(x+delta) - f(x)) / delta. In fact, there is a lot of work still to be
done to make sure all the holes are plugged. I'm working on putting many such checks in
optimx, but it will be a few weeks yet before they are even on the R-forge version.

- you may be able to change your objective function with penalty or barrier functions, or
using log(parameter) to avoid zero divides etc.

Unfortunately, providing "automatic" minimization tools often means and automatic
generation of lots of error messages as these tools follow downhill directions until they
are in very weird parts of parameter spaces.

JN
On 10/22/2011 06:00 AM, r-help-request at r-project.org wrote: