diff --git a/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_100.do b/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_100.do
index ec9127e..8385e51 100644
--- a/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_100.do
+++ b/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_100.do
@@ -61,7 +61,7 @@ forvalues replication = 1/1000 {
     if `difi'==4 {
     mat E= (0 \ 0 \ 0 \ 0.3)
     }
-    mat D=B+E
+    mat D=B-E
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -466,7 +466,7 @@ forvalues replication = 1/1000 {
     if `difi'==4 {
       mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3)
     }
-    mat D=B+E
+    mat D=B-E
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -871,7 +871,7 @@ forvalues replication = 1/1000 {
     if `difi'==4 {
     mat E= (0 \ 0 \ 0 \ 0.5)
     }
-    mat D=B+E
+    mat D=B-E
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -1277,7 +1277,7 @@ forvalues replication = 1/1000 {
     if `difi'==4 {
       mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5)
     }
-    mat D=B+E
+    mat D=B-E
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -1705,7 +1705,7 @@ forvalues replication = 1/1000 {
     if `difj'==4 {
       mat F=(0 \ 0 \ 0 \ 0.3)
     }
-    mat D=B+E+F
+    mat D=B-E-F
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -2246,7 +2246,7 @@ forvalues replication = 1/1000 {
     if `difj'==4 {
       mat F=(0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3)
     }
-    mat D=B+E+F
+    mat D=B-E-F
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -2787,7 +2787,7 @@ forvalues replication = 1/1000 {
     if `difj'==4 {
       mat F=(0 \ 0 \ 0 \ 0.5)
     }
-    mat D=B+E+F
+    mat D=B-E-F
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -3326,7 +3326,7 @@ forvalues replication = 1/1000 {
     if `difj'==4 {
       mat F=(0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5)
     }
-    mat D=B+E+F
+    mat D=B-E-F
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
diff --git a/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_200.do b/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_200.do
index 78bc3dc..789cf71 100644
--- a/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_200.do
+++ b/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_200.do
@@ -61,7 +61,7 @@ forvalues replication = 1/1000 {
     if `difi'==4 {
     mat E= (0 \ 0 \ 0 \ 0.3)
     }
-    mat D=B+E
+    mat D=B-E
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -466,7 +466,7 @@ forvalues replication = 1/1000 {
     if `difi'==4 {
       mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3)
     }
-    mat D=B+E
+    mat D=B-E
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -871,7 +871,7 @@ forvalues replication = 1/1000 {
     if `difi'==4 {
     mat E= (0 \ 0 \ 0 \ 0.5)
     }
-    mat D=B+E
+    mat D=B-E
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -1277,7 +1277,7 @@ forvalues replication = 1/1000 {
     if `difi'==4 {
       mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5)
     }
-    mat D=B+E
+    mat D=B-E
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -1705,7 +1705,7 @@ forvalues replication = 1/1000 {
     if `difj'==4 {
       mat F=(0 \ 0 \ 0 \ 0.3)
     }
-    mat D=B+E+F
+    mat D=B-E-F
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -2246,7 +2246,7 @@ forvalues replication = 1/1000 {
     if `difj'==4 {
       mat F=(0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3)
     }
-    mat D=B+E+F
+    mat D=B-E-F
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -2787,7 +2787,7 @@ forvalues replication = 1/1000 {
     if `difj'==4 {
       mat F=(0 \ 0 \ 0 \ 0.5)
     }
-    mat D=B+E+F
+    mat D=B-E-F
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -3326,7 +3326,7 @@ forvalues replication = 1/1000 {
     if `difj'==4 {
       mat F=(0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5)
     }
-    mat D=B+E+F
+    mat D=B-E-F
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
diff --git a/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_300.do b/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_300.do
index ee998b1..fefb4fe 100644
--- a/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_300.do
+++ b/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_300.do
@@ -61,7 +61,7 @@ forvalues replication = 1/1000 {
     if `difi'==4 {
     mat E= (0 \ 0 \ 0 \ 0.3)
     }
-    mat D=B+E
+    mat D=B-E
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -466,7 +466,7 @@ forvalues replication = 1/1000 {
     if `difi'==4 {
       mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3)
     }
-    mat D=B+E
+    mat D=B-E
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -871,7 +871,7 @@ forvalues replication = 1/1000 {
     if `difi'==4 {
     mat E= (0 \ 0 \ 0 \ 0.5)
     }
-    mat D=B+E
+    mat D=B-E
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -1277,7 +1277,7 @@ forvalues replication = 1/1000 {
     if `difi'==4 {
       mat E= (0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5)
     }
-    mat D=B+E
+    mat D=B-E
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -1705,7 +1705,7 @@ forvalues replication = 1/1000 {
     if `difj'==4 {
       mat F=(0 \ 0 \ 0 \ 0.3)
     }
-    mat D=B+E+F
+    mat D=B-E-F
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -2246,7 +2246,7 @@ forvalues replication = 1/1000 {
     if `difj'==4 {
       mat F=(0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3)
     }
-    mat D=B+E+F
+    mat D=B-E-F
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -2787,7 +2787,7 @@ forvalues replication = 1/1000 {
     if `difj'==4 {
       mat F=(0 \ 0 \ 0 \ 0.5)
     }
-    mat D=B+E+F
+    mat D=B-E-F
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -3326,7 +3326,7 @@ forvalues replication = 1/1000 {
     if `difj'==4 {
       mat F=(0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5)
     }
-    mat D=B+E+F
+    mat D=B-E-F
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
diff --git a/Scripts/Scenarios/DIF/scenarios_DIF_J7_baseline_100.do b/Scripts/Scenarios/DIF/scenarios_DIF_J7_baseline_100.do
index 601dee4..789dc33 100644
--- a/Scripts/Scenarios/DIF/scenarios_DIF_J7_baseline_100.do
+++ b/Scripts/Scenarios/DIF/scenarios_DIF_J7_baseline_100.do
@@ -103,7 +103,7 @@ forvalues replication = 1/1000 {
     if `difj'==7 {
       mat F=(0 \ 0 \ 0 \ 0 \ 0 \ 0 \ 0.3)
     }
-    mat D=B+E+F
+    mat D=B-E-F
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -774,7 +774,7 @@ forvalues replication = 1/1000 {
     if `difj'==7 {
       mat F=(0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3)
     }
-    mat D=B+E+F
+    mat D=B-E-F
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -1440,7 +1440,7 @@ forvalues replication = 1/1000 {
     if `difj'==7 {
       mat F=(0 \ 0 \ 0 \ 0 \ 0 \ 0 \ 0.5)
     }
-    mat D=B+E+F
+    mat D=B-E-F
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -2109,7 +2109,7 @@ forvalues replication = 1/1000 {
     if `difj'==7 {
       mat F=(0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5)
     }
-    mat D=B+E+F
+    mat D=B-E-F
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -2801,7 +2801,7 @@ forvalues replication = 1/1000 {
     if `difk'==7 {
       mat G=(0 \ 0 \ 0 \ 0 \ 0 \ 0 \ 0.3)
     }
-    mat D=B+E+F+G
+    mat D=B-E-F-G
     di
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
     qui gen TT = 1
@@ -3661,7 +3661,7 @@ forvalues replication = 1/1000 {
     if `difk'==7 {
       mat G=(0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3)
     }
-    mat D=B+E+F+G
+    mat D=B-E-F-G
     di
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
     qui gen TT = 1
@@ -4532,7 +4532,7 @@ forvalues replication = 1/1000 {
     if `difk'==7 {
       mat G=(0 \ 0 \ 0 \ 0 \ 0 \ 0 \ 0.5)
     }
-    mat D=B+E+F+G
+    mat D=B-E-F-G
     di
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
     qui gen TT = 1
@@ -5368,7 +5368,7 @@ forvalues replication = 1/1000 {
     if `difk'==7 {
       mat G=(0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5)
     }
-    mat D=B+E+F+G
+    mat D=B-E-F-G
     di
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
     qui gen TT = 1
diff --git a/Scripts/Scenarios/DIF/scenarios_DIF_J7_baseline_200.do b/Scripts/Scenarios/DIF/scenarios_DIF_J7_baseline_200.do
index 54ca56d..2586c8b 100644
--- a/Scripts/Scenarios/DIF/scenarios_DIF_J7_baseline_200.do
+++ b/Scripts/Scenarios/DIF/scenarios_DIF_J7_baseline_200.do
@@ -103,7 +103,7 @@ forvalues replication = 1/1000 {
     if `difj'==7 {
       mat F=(0 \ 0 \ 0 \ 0 \ 0 \ 0 \ 0.3)
     }
-    mat D=B+E+F
+    mat D=B-E-F
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -774,7 +774,7 @@ forvalues replication = 1/1000 {
     if `difj'==7 {
       mat F=(0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3)
     }
-    mat D=B+E+F
+    mat D=B-E-F
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -1440,7 +1440,7 @@ forvalues replication = 1/1000 {
     if `difj'==7 {
       mat F=(0 \ 0 \ 0 \ 0 \ 0 \ 0 \ 0.5)
     }
-    mat D=B+E+F
+    mat D=B-E-F
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -2109,7 +2109,7 @@ forvalues replication = 1/1000 {
     if `difj'==7 {
       mat F=(0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5)
     }
-    mat D=B+E+F
+    mat D=B-E-F
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -2801,7 +2801,7 @@ forvalues replication = 1/1000 {
     if `difk'==7 {
       mat G=(0 \ 0 \ 0 \ 0 \ 0 \ 0 \ 0.3)
     }
-    mat D=B+E+F+G
+    mat D=B-E-F-G
     di
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
     qui gen TT = 1
@@ -3661,7 +3661,7 @@ forvalues replication = 1/1000 {
     if `difk'==7 {
       mat G=(0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3)
     }
-    mat D=B+E+F+G
+    mat D=B-E-F-G
     di
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
     qui gen TT = 1
@@ -4532,7 +4532,7 @@ forvalues replication = 1/1000 {
     if `difk'==7 {
       mat G=(0 \ 0 \ 0 \ 0 \ 0 \ 0 \ 0.5)
     }
-    mat D=B+E+F+G
+    mat D=B-E-F-G
     di
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
     qui gen TT = 1
@@ -5368,7 +5368,7 @@ forvalues replication = 1/1000 {
     if `difk'==7 {
       mat G=(0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5)
     }
-    mat D=B+E+F+G
+    mat D=B-E-F-G
     di
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
     qui gen TT = 1
diff --git a/Scripts/Scenarios/DIF/scenarios_DIF_J7_baseline_300.do b/Scripts/Scenarios/DIF/scenarios_DIF_J7_baseline_300.do
index 415bc47..2ce98b4 100644
--- a/Scripts/Scenarios/DIF/scenarios_DIF_J7_baseline_300.do
+++ b/Scripts/Scenarios/DIF/scenarios_DIF_J7_baseline_300.do
@@ -103,7 +103,7 @@ forvalues replication = 1/1000 {
     if `difj'==7 {
       mat F=(0 \ 0 \ 0 \ 0 \ 0 \ 0 \ 0.3)
     }
-    mat D=B+E+F
+    mat D=B-E-F
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -774,7 +774,7 @@ forvalues replication = 1/1000 {
     if `difj'==7 {
       mat F=(0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3)
     }
-    mat D=B+E+F
+    mat D=B-E-F
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -1440,7 +1440,7 @@ forvalues replication = 1/1000 {
     if `difj'==7 {
       mat F=(0 \ 0 \ 0 \ 0 \ 0 \ 0 \ 0.5)
     }
-    mat D=B+E+F
+    mat D=B-E-F
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -2109,7 +2109,7 @@ forvalues replication = 1/1000 {
     if `difj'==7 {
       mat F=(0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5)
     }
-    mat D=B+E+F
+    mat D=B-E-F
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
     qui gen TT = 1
 		tempfile grp1
@@ -2801,7 +2801,7 @@ forvalues replication = 1/1000 {
     if `difk'==7 {
       mat G=(0 \ 0 \ 0 \ 0 \ 0 \ 0 \ 0.3)
     }
-    mat D=B+E+F+G
+    mat D=B-E-F-G
     di
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
     qui gen TT = 1
@@ -3661,7 +3661,7 @@ forvalues replication = 1/1000 {
     if `difk'==7 {
       mat G=(0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0.3,0.3,0.3)
     }
-    mat D=B+E+F+G
+    mat D=B-E-F-G
     di
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
     qui gen TT = 1
@@ -4532,7 +4532,7 @@ forvalues replication = 1/1000 {
     if `difk'==7 {
       mat G=(0 \ 0 \ 0 \ 0 \ 0 \ 0 \ 0.5)
     }
-    mat D=B+E+F+G
+    mat D=B-E-F-G
     di
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
     qui gen TT = 1
@@ -5368,7 +5368,7 @@ forvalues replication = 1/1000 {
     if `difk'==7 {
       mat G=(0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0,0,0 \ 0.5,0.5,0.5)
     }
-    mat D=B+E+F+G
+    mat D=B-E-F-G
     di
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
     qui gen TT = 1