run_extra_optimization function

Run extra optimization on an already fitted object

Description

run_extra_optimization(object, nlminb_loops = 0, newton_loops = 1)

Arguments

• object: An object from sdmTMB().

• nlminb_loops: How many extra times to run stats::nlminb()

  optimization. Sometimes restarting the optimizer at the previous best values aids convergence.

• newton_loops: How many extra Newton optimization loops to try with stats::optimHess(). Sometimes aids convergence.

Returns

An updated model fit of class sdmTMB.

Examples

# Run extra optimization steps to help convergence:
# (Not typically needed)
fit <- sdmTMB(density ~ 0 + poly(depth, 2) + as.factor(year),
  data = pcod_2011, mesh = pcod_mesh_2011, family = tweedie())
fit_1 <- run_extra_optimization(fit, newton_loops = 1)
max(fit$gradients)
max(fit_1$gradients)