diff --git a/Figures_manuscript/scripts/Fig_SUPnote4_simulation_test_boundaries.R b/Figures_manuscript/scripts/Fig_SUPnote4_simulation_test_boundaries.R
index 99d6e4c9fe612cf4a4acfee464c241839ceb6eaf..e46b8368a07558cb0362960cc5be4e43bc2fa112 100644
--- a/Figures_manuscript/scripts/Fig_SUPnote4_simulation_test_boundaries.R
+++ b/Figures_manuscript/scripts/Fig_SUPnote4_simulation_test_boundaries.R
@@ -17,7 +17,7 @@ epsh = matrix(c(h1, 0, 0, h2), ncol=2)/M
 epshlow = matrix(c(0.1, 0, 0, 0.1), ncol=2)/M
 
 epsg1 = matrix(c(h1, pgh, pgh, h2), ncol=2)/ M
-epsg2 = matrix(c(h1, -0.2, -0.2, h2), ncol=2)/ M
+epsg2 = matrix(c(h1, -pgh, -pgh, h2), ncol=2)/ M
 
 # draw genetic effect according to several scenarios
 BG_mat <-list()