diff --git a/Figures_manuscript/scripts/Fig_SUPnote4_simulation_test_boundaries.R b/Figures_manuscript/scripts/Fig_SUPnote4_simulation_test_boundaries.R
index e46b8368a07558cb0362960cc5be4e43bc2fa112..be198b18076f6acb1c4c921a8b93f29e546fb571 100644
--- a/Figures_manuscript/scripts/Fig_SUPnote4_simulation_test_boundaries.R
+++ b/Figures_manuscript/scripts/Fig_SUPnote4_simulation_test_boundaries.R
@@ -73,6 +73,7 @@ for(hypothesis in names(BG_mat)){
         D_pval[simu,"N_samp"] = N
         D_pval[simu,"env_cor"] = ce
         # Draw SNPs :
+	set.seed(1111)
         MAF = runif(M,0.05, 0.95)
         i= 1