R/extract.solution.r

In jae0/snowcrab: Snow crab assessment suite using aegis* and associated projects and data streams.

extract.solution = function( oms, optim.model.core ) {
  solution = NULL
  if (!is.null ( oms ) ){
    good = which( oms[,"convergence"] == 0 )
    if (length(good) > 1 ) {
      oms = oms[ good, ]
      final = which.min( oms[,"error"] )
      solution = optim.model.core( inits= unlist(oms[ final, "par"]) )
    } else if (length(good) == 1) { 
      oms = oms[ good, ]
      solution = optim.model.core( inits=unlist( oms["par"] ) )
    }
  }
  if (is.null (solution) ) print("Convergence problems")
  return (solution)
}