adding jitter to mapping

flr · Apr 16, 2024 · 6df0e2b · 6df0e2b
1 parent cae8693
commit 6df0e2b
Showing 1 changed file with 6 additions and 6 deletions.
diff --git a/R/monteCarloStock.R b/R/monteCarloStock.R
@@ -45,12 +45,12 @@ monteCarloStock <- function(stck,tun,sam,realisations,return.sam=FALSE,saveParsD
       parS[[set.pars.internal$SAM[i]]][set.pars.internal$number[i]+1] <- set.pars.internal$value[i]
       map[[set.pars.internal$SAM[[i]]]][set.pars.internal$number[i]+1] <- NA
     }
-    map=list(logSdLogN=as.factor(map$logSdLogN),
-             logSdLogP=as.factor(map$logSdLogP),
-             logFpar=as.factor(map$logFpar),
-             logQpow=as.factor(map$logQpow),
-             logSdLogFsta=as.factor(map$logSdLogFsta),
-             logSdLogObs=as.factor(map$logSdLogObs))
+    map=list(logSdLogN=as.factor(jitter(map$logSdLogN,0.1)),
+             logSdLogP=as.factor(jitter(map$logSdLogP,0.1)),
+             logFpar=as.factor(jitter(map$logFpar,0.1)),
+             logQpow=as.factor(jitter(map$logQpow,0.1)),
+             logSdLogFsta=as.factor(jitter(map$logSdLogFsta,0.1)),
+             logSdLogObs=as.factor(jitter(map$logSdLogObs,0.1)))
     map <- map[which(names(map) %in% set.pars$SAM)]
 
     object            <- FLSAM(stcks,tun,sam@control,set.pars=set.pars,return.fit=T)