I got several code 6 when running simulations. Some of those cases have a near singular estimated hessian, though the calculated hessian is fine. Is this the source of non-convergence? Would it work to re- run the data with the parameter estimates as initial values? Is it possible to have the optimizer use the calculated hessian all the time?