ensure correct parameter is changed for relaxed rate model

tece-lab · Oct 20, 2023 · 441ceef · 441ceef
1 parent ddf3861
commit 441ceef
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/R/run_daisie_ml.R b/R/run_daisie_ml.R
@@ -109,7 +109,10 @@ run_daisie_ml <- function(daisie_data,
     # prevent initial parameters being greater than upper parameter limit
     pick <- which(c("cladogenesis", "extinction", "carrying_capacity",
                     "immigration", "anagenesis") == cs_version$relaxed_par)
-    initparsopt[initparsopt[pick] > par_upper_bound] <- par_upper_bound / 2
+    if (initparsopt[pick] > par_upper_bound) {
+      initparsopt[pick] <- par_upper_bound / 2
+    }
+
   }
 
   ##### ML Optimization ####