Merge pull request #446 from RobFryer/VDS_cl_seed

Imposex assessment - ensure reproducibility
osparcomm · Mar 5, 2024 · 20fb051 · 20fb051
2 parents 2e7b88d + a851edc
commit 20fb051
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/R/proportional_odds_functions.R b/R/proportional_odds_functions.R
@@ -175,6 +175,8 @@ ctsm.VDS.cl <- function(fit, nsim = 1000) {
 
   indexID <- setdiff(names(fit$par), cutsID)
 
+  set.seed(fit$seed)
+
   data <- MASS::mvrnorm(nsim, fit$par, fit$vcov)
 
   data.cuts <- data[, cutsID, drop = FALSE]