diff --git a/R/estimate_nhmm.R b/R/estimate_nhmm.R
index 473e7ad..a60e066 100644
--- a/R/estimate_nhmm.R
+++ b/R/estimate_nhmm.R
@@ -23,7 +23,7 @@
 #'  
 #' By default, the convergence is claimed when either the absolute or relative 
 #' change of the objective function is less than `1e-8`, or the absolute or 
-#' relative change of the parameters is less than `1e-4`. These can be changed
+#' relative change of the parameters is less than `1e-6`. These can be changed
 #' by passing arguments `ftol_abs`, `ftol_rel`, `xtol_abs`, and `xtol_rel` via
 #' `...`.  These, as well as argument `maxeval` (maximum number of iterations, 
 #' 1e4 by default), and `print_level` (default is `0`, no console output of 
diff --git a/R/fit_mnhmm.R b/R/fit_mnhmm.R
index 6c9dc8a..e40f5cb 100644
--- a/R/fit_mnhmm.R
+++ b/R/fit_mnhmm.R
@@ -18,8 +18,8 @@ fit_mnhmm <- function(model, inits, init_sd, restarts, lambda, method,
     list(
       ftol_abs = 1e-8,
       ftol_rel = 1e-8,
-      xtol_abs = 1e-4,
-      xtol_rel = 1e-4,
+      xtol_abs = 1e-6,
+      xtol_rel = 1e-6,
       maxeval = 1e4,
       print_level = 0,
       algorithm = "NLOPT_LD_LBFGS",
diff --git a/R/fit_nhmm.R b/R/fit_nhmm.R
index 375bb05..5ce19e8 100644
--- a/R/fit_nhmm.R
+++ b/R/fit_nhmm.R
@@ -18,8 +18,8 @@ fit_nhmm <- function(model, inits, init_sd, restarts, lambda, method, pseudocoun
     list(
       ftol_abs = 1e-8,
       ftol_rel = 1e-8,
-      xtol_abs = 1e-4,
-      xtol_rel = 1e-4,
+      xtol_abs = 1e-6,
+      xtol_rel = 1e-6,
       maxeval = 1e4,
       print_level = 0,
       algorithm = "NLOPT_LD_LBFGS",