From 04fe611abaf97ba157aa903d384c10f8f1622cbe Mon Sep 17 00:00:00 2001
From: corentinchoisy <corentin.choisy@proton.me>
Date: Fri, 12 Jan 2024 10:54:33 +0100
Subject: [PATCH] Updated catalogue and scripts

---
 .../DIF/scenarios_DIF_J4_baseline_100.do      | 140 ++++++++++----
 .../DIF/scenarios_DIF_J4_baseline_200.do      | 140 ++++++++++----
 .../DIF/scenarios_DIF_J4_baseline_300.do      | 140 ++++++++++----
 .../DIF/scenarios_DIF_J7_baseline_100.do      | 180 ++++++++++++++----
 .../DIF/scenarios_DIF_J7_baseline_200.do      | 180 ++++++++++++++----
 .../DIF/scenarios_DIF_J7_baseline_300.do      | 180 ++++++++++++++----
 catalogue.md                                  |  16 +-
 7 files changed, 728 insertions(+), 248 deletions(-)

diff --git a/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_100.do b/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_100.do
index d42adbd..163209a 100644
--- a/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_100.do
+++ b/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_100.do
@@ -37,7 +37,7 @@ local Nn = 100
 di "Scenario 5A - N=100"
 
 forvalues replication = 1/1000 {
-  	if mod(`replication',10)==0 {
+  	if mod(`replication',100)==0 {
   	di "replication = `replication'"
     }
     qui local difi = runiformint(1,4)
@@ -72,6 +72,7 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -92,7 +93,7 @@ forvalues replication = 1/1000 {
   di "Scenario 5B - N=100"
 
   forvalues replication = 1/1000 {
-    	if mod(`replication',10)==0 {
+    	if mod(`replication',100)==0 {
     	di "replication = `replication'"
       }
       qui local difi = runiformint(1,4)
@@ -128,6 +129,7 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -147,7 +149,7 @@ forvalues replication = 1/1000 {
     di "Scenario 5C - N=100"
 
     forvalues replication = 1/1000 {
-      	if mod(`replication',10)==0 {
+      	if mod(`replication',100)==0 {
       	di "replication = `replication'"
         }
         qui local difi = runiformint(1,4)
@@ -183,6 +185,7 @@ forvalues replication = 1/1000 {
     		use `grp0'
         qui append using `grp1'
     		drop id
+        qui gen dif1 = `difi'
     		qui gen id = _n
     		order(id)
     		qui gen replication = `replication'
@@ -203,7 +206,7 @@ forvalues replication = 1/1000 {
         di "Scenario 5D - N=100"
 
         forvalues replication = 1/1000 {
-          	if mod(`replication',10)==0 {
+          	if mod(`replication',100)==0 {
           	di "replication = `replication'"
             }
             qui local difi = runiformint(1,4)
@@ -239,6 +242,7 @@ forvalues replication = 1/1000 {
         		use `grp0'
             qui append using `grp1'
         		drop id
+            qui gen dif1 = `difi'
         		qui gen id = _n
         		order(id)
         		qui gen replication = `replication'
@@ -258,7 +262,7 @@ forvalues replication = 1/1000 {
           di "Scenario 5E - N=100"
 
           forvalues replication = 1/1000 {
-            	if mod(`replication',10)==0 {
+            	if mod(`replication',100)==0 {
             	di "replication = `replication'"
               }
               qui local difi = runiformint(1,4)
@@ -294,6 +298,7 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -326,7 +331,7 @@ forvalues replication = 1/1000 {
 di "Scenario 6A - N=100"
 
 forvalues replication = 1/1000 {
-  	if mod(`replication',10)==0 {
+  	if mod(`replication',100)==0 {
   	di "replication = `replication'"
     }
     qui local difi = runiformint(1,4)
@@ -361,6 +366,7 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -380,7 +386,7 @@ forvalues replication = 1/1000 {
   di "Scenario 6B - N=100"
 
   forvalues replication = 1/1000 {
-    	if mod(`replication',10)==0 {
+    	if mod(`replication',100)==0 {
     	di "replication = `replication'"
       }
       qui local difi = runiformint(1,4)
@@ -416,6 +422,7 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -435,7 +442,7 @@ forvalues replication = 1/1000 {
       di "Scenario 6C - N=100"
 
       forvalues replication = 1/1000 {
-        	if mod(`replication',10)==0 {
+        	if mod(`replication',100)==0 {
         	di "replication = `replication'"
           }
           qui local difi = runiformint(1,4)
@@ -471,6 +478,7 @@ forvalues replication = 1/1000 {
       		use `grp0'
           qui append using `grp1'
       		drop id
+          qui gen dif1 = `difi'
       		qui gen id = _n
       		order(id)
       		qui gen replication = `replication'
@@ -490,7 +498,7 @@ forvalues replication = 1/1000 {
           di "Scenario 6D - N=100"
 
           forvalues replication = 1/1000 {
-            	if mod(`replication',10)==0 {
+            	if mod(`replication',100)==0 {
             	di "replication = `replication'"
               }
               qui local difi = runiformint(1,4)
@@ -526,6 +534,7 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -545,7 +554,7 @@ forvalues replication = 1/1000 {
               di "Scenario 6E - N=100"
 
               forvalues replication = 1/1000 {
-                	if mod(`replication',10)==0 {
+                	if mod(`replication',100)==0 {
                 	di "replication = `replication'"
                   }
                   qui local difi = runiformint(1,4)
@@ -581,6 +590,7 @@ forvalues replication = 1/1000 {
               		use `grp0'
                   qui append using `grp1'
               		drop id
+                  qui gen dif1 = `difi'
               		qui gen id = _n
               		order(id)
               		qui gen replication = `replication'
@@ -612,7 +622,7 @@ forvalues replication = 1/1000 {
 di "Scenario 7A - N=100"
 
 forvalues replication = 1/1000 {
-  	if mod(`replication',10)==0 {
+  	if mod(`replication',100)==0 {
   	di "replication = `replication'"
     }
     qui local difi = runiformint(1,4)
@@ -648,6 +658,7 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -668,7 +679,7 @@ forvalues replication = 1/1000 {
   di "Scenario 7B - N=100"
 
   forvalues replication = 1/1000 {
-    	if mod(`replication',10)==0 {
+    	if mod(`replication',100)==0 {
     	di "replication = `replication'"
       }
       qui local difi = runiformint(1,4)
@@ -704,6 +715,7 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -723,7 +735,7 @@ forvalues replication = 1/1000 {
     di "Scenario 7C - N=100"
 
     forvalues replication = 1/1000 {
-      	if mod(`replication',10)==0 {
+      	if mod(`replication',100)==0 {
       	di "replication = `replication'"
         }
         qui local difi = runiformint(1,4)
@@ -759,6 +771,7 @@ forvalues replication = 1/1000 {
     		use `grp0'
         qui append using `grp1'
     		drop id
+        qui gen dif1 = `difi'
     		qui gen id = _n
     		order(id)
     		qui gen replication = `replication'
@@ -779,7 +792,7 @@ forvalues replication = 1/1000 {
         di "Scenario 7D - N=100"
 
         forvalues replication = 1/1000 {
-          	if mod(`replication',10)==0 {
+          	if mod(`replication',100)==0 {
           	di "replication = `replication'"
             }
             qui local difi = runiformint(1,4)
@@ -815,6 +828,7 @@ forvalues replication = 1/1000 {
         		use `grp0'
             qui append using `grp1'
         		drop id
+            qui gen dif1 = `difi'
         		qui gen id = _n
         		order(id)
         		qui gen replication = `replication'
@@ -834,7 +848,7 @@ forvalues replication = 1/1000 {
           di "Scenario 7E - N=100"
 
           forvalues replication = 1/1000 {
-            	if mod(`replication',10)==0 {
+            	if mod(`replication',100)==0 {
             	di "replication = `replication'"
               }
               qui local difi = runiformint(1,4)
@@ -870,6 +884,7 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -902,7 +917,7 @@ forvalues replication = 1/1000 {
 di "Scenario 8A - N=100"
 
 forvalues replication = 1/1000 {
-  	if mod(`replication',10)==0 {
+  	if mod(`replication',100)==0 {
   	di "replication = `replication'"
     }
     qui local difi = runiformint(1,4)
@@ -938,6 +953,7 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -957,7 +973,7 @@ forvalues replication = 1/1000 {
   di "Scenario 8B - N=100"
 
   forvalues replication = 1/1000 {
-    	if mod(`replication',10)==0 {
+    	if mod(`replication',100)==0 {
     	di "replication = `replication'"
       }
       qui local difi = runiformint(1,4)
@@ -993,6 +1009,7 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -1012,7 +1029,7 @@ forvalues replication = 1/1000 {
       di "Scenario 8C - N=100"
 
       forvalues replication = 1/1000 {
-        	if mod(`replication',10)==0 {
+        	if mod(`replication',100)==0 {
         	di "replication = `replication'"
           }
           qui local difi = runiformint(1,4)
@@ -1048,6 +1065,7 @@ forvalues replication = 1/1000 {
       		use `grp0'
           qui append using `grp1'
       		drop id
+          qui gen dif1 = `difi'
       		qui gen id = _n
       		order(id)
       		qui gen replication = `replication'
@@ -1067,7 +1085,7 @@ forvalues replication = 1/1000 {
           di "Scenario 8D - N=100"
 
           forvalues replication = 1/1000 {
-            	if mod(`replication',10)==0 {
+            	if mod(`replication',100)==0 {
             	di "replication = `replication'"
               }
               qui local difi = runiformint(1,4)
@@ -1103,6 +1121,7 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -1122,7 +1141,7 @@ forvalues replication = 1/1000 {
               di "Scenario 8E - N=100"
 
               forvalues replication = 1/1000 {
-                	if mod(`replication',10)==0 {
+                	if mod(`replication',100)==0 {
                 	di "replication = `replication'"
                   }
                   qui local difi = runiformint(1,4)
@@ -1158,6 +1177,7 @@ forvalues replication = 1/1000 {
               		use `grp0'
                   qui append using `grp1'
               		drop id
+                  qui gen dif1 = `difi'
               		qui gen id = _n
               		order(id)
               		qui gen replication = `replication'
@@ -1197,7 +1217,7 @@ di "Scenario 9A - N=100"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 4
@@ -1251,6 +1271,8 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -1272,7 +1294,7 @@ forvalues replication = 1/1000 {
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 4
@@ -1326,6 +1348,8 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -1346,7 +1370,7 @@ forvalues replication = 1/1000 {
 
     forvalues replication = 1/1000 {
       clear
-      if mod(`replication',10)==0 {
+      if mod(`replication',100)==0 {
       di "replication = `replication'"
       }
       qui set obs 4
@@ -1400,6 +1424,8 @@ forvalues replication = 1/1000 {
     		use `grp0'
         qui append using `grp1'
     		drop id
+        qui gen dif1 = `difi'
+        qui gen dif2 = `difj'
     		qui gen id = _n
     		order(id)
     		qui gen replication = `replication'
@@ -1421,7 +1447,7 @@ forvalues replication = 1/1000 {
 
         forvalues replication = 1/1000 {
           clear
-          if mod(`replication',10)==0 {
+          if mod(`replication',100)==0 {
           di "replication = `replication'"
           }
           qui set obs 4
@@ -1475,6 +1501,8 @@ forvalues replication = 1/1000 {
         		use `grp0'
             qui append using `grp1'
         		drop id
+            qui gen dif1 = `difi'
+            qui gen dif2 = `difj'
         		qui gen id = _n
         		order(id)
         		qui gen replication = `replication'
@@ -1495,7 +1523,7 @@ forvalues replication = 1/1000 {
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 4
@@ -1549,6 +1577,8 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -1576,7 +1606,7 @@ di "Scenario 10A - N=100"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 4
@@ -1630,6 +1660,8 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -1651,7 +1683,7 @@ forvalues replication = 1/1000 {
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 4
@@ -1705,6 +1737,8 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -1725,7 +1759,7 @@ forvalues replication = 1/1000 {
 
     forvalues replication = 1/1000 {
       clear
-      if mod(`replication',10)==0 {
+      if mod(`replication',100)==0 {
       di "replication = `replication'"
       }
       qui set obs 4
@@ -1779,6 +1813,8 @@ forvalues replication = 1/1000 {
     		use `grp0'
         qui append using `grp1'
     		drop id
+        qui gen dif1 = `difi'
+        qui gen dif2 = `difj'
     		qui gen id = _n
     		order(id)
     		qui gen replication = `replication'
@@ -1800,7 +1836,7 @@ forvalues replication = 1/1000 {
 
         forvalues replication = 1/1000 {
           clear
-          if mod(`replication',10)==0 {
+          if mod(`replication',100)==0 {
           di "replication = `replication'"
           }
           qui set obs 4
@@ -1854,6 +1890,8 @@ forvalues replication = 1/1000 {
         		use `grp0'
             qui append using `grp1'
         		drop id
+            qui gen dif1 = `difi'
+            qui gen dif2 = `difj'
         		qui gen id = _n
         		order(id)
         		qui gen replication = `replication'
@@ -1874,7 +1912,7 @@ forvalues replication = 1/1000 {
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 4
@@ -1928,6 +1966,8 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -1952,7 +1992,7 @@ di "Scenario 11A - N=100"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 4
@@ -2006,6 +2046,8 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -2027,7 +2069,7 @@ forvalues replication = 1/1000 {
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 4
@@ -2081,6 +2123,8 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -2101,7 +2145,7 @@ forvalues replication = 1/1000 {
 
     forvalues replication = 1/1000 {
       clear
-      if mod(`replication',10)==0 {
+      if mod(`replication',100)==0 {
       di "replication = `replication'"
       }
       qui set obs 4
@@ -2155,6 +2199,8 @@ forvalues replication = 1/1000 {
     		use `grp0'
         qui append using `grp1'
     		drop id
+        qui gen dif1 = `difi'
+        qui gen dif2 = `difj'
     		qui gen id = _n
     		order(id)
     		qui gen replication = `replication'
@@ -2176,7 +2222,7 @@ forvalues replication = 1/1000 {
 
         forvalues replication = 1/1000 {
           clear
-          if mod(`replication',10)==0 {
+          if mod(`replication',100)==0 {
           di "replication = `replication'"
           }
           qui set obs 4
@@ -2230,6 +2276,8 @@ forvalues replication = 1/1000 {
         		use `grp0'
             qui append using `grp1'
         		drop id
+            qui gen dif1 = `difi'
+            qui gen dif2 = `difj'
         		qui gen id = _n
         		order(id)
         		qui gen replication = `replication'
@@ -2250,7 +2298,7 @@ forvalues replication = 1/1000 {
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 4
@@ -2304,6 +2352,8 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -2331,7 +2381,7 @@ di "Scenario 12A - N=100"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 4
@@ -2385,6 +2435,8 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -2406,7 +2458,7 @@ forvalues replication = 1/1000 {
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 4
@@ -2460,6 +2512,8 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -2480,7 +2534,7 @@ forvalues replication = 1/1000 {
 
     forvalues replication = 1/1000 {
       clear
-      if mod(`replication',10)==0 {
+      if mod(`replication',100)==0 {
       di "replication = `replication'"
       }
       qui set obs 4
@@ -2534,6 +2588,8 @@ forvalues replication = 1/1000 {
     		use `grp0'
         qui append using `grp1'
     		drop id
+        qui gen dif1 = `difi'
+        qui gen dif2 = `difj'
     		qui gen id = _n
     		order(id)
     		qui gen replication = `replication'
@@ -2555,7 +2611,7 @@ forvalues replication = 1/1000 {
 
         forvalues replication = 1/1000 {
           clear
-          if mod(`replication',10)==0 {
+          if mod(`replication',100)==0 {
           di "replication = `replication'"
           }
           qui set obs 4
@@ -2609,6 +2665,8 @@ forvalues replication = 1/1000 {
         		use `grp0'
             qui append using `grp1'
         		drop id
+            qui gen dif1 = `difi'
+            qui gen dif2 = `difj'
         		qui gen id = _n
         		order(id)
         		qui gen replication = `replication'
@@ -2629,7 +2687,7 @@ forvalues replication = 1/1000 {
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 4
@@ -2683,6 +2741,8 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
diff --git a/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_200.do b/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_200.do
index d72c96d..345129c 100644
--- a/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_200.do
+++ b/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_200.do
@@ -37,7 +37,7 @@ local Nn = 200
 di "Scenario 5A - N=200"
 
 forvalues replication = 1/1000 {
-  	if mod(`replication',10)==0 {
+  	if mod(`replication',100)==0 {
   	di "replication = `replication'"
     }
     qui local difi = runiformint(1,4)
@@ -72,6 +72,7 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -92,7 +93,7 @@ forvalues replication = 1/1000 {
   di "Scenario 5B - N=200"
 
   forvalues replication = 1/1000 {
-    	if mod(`replication',10)==0 {
+    	if mod(`replication',100)==0 {
     	di "replication = `replication'"
       }
       qui local difi = runiformint(1,4)
@@ -128,6 +129,7 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -147,7 +149,7 @@ forvalues replication = 1/1000 {
     di "Scenario 5C - N=200"
 
     forvalues replication = 1/1000 {
-      	if mod(`replication',10)==0 {
+      	if mod(`replication',100)==0 {
       	di "replication = `replication'"
         }
         qui local difi = runiformint(1,4)
@@ -183,6 +185,7 @@ forvalues replication = 1/1000 {
     		use `grp0'
         qui append using `grp1'
     		drop id
+        qui gen dif1 = `difi'
     		qui gen id = _n
     		order(id)
     		qui gen replication = `replication'
@@ -203,7 +206,7 @@ forvalues replication = 1/1000 {
         di "Scenario 5D - N=200"
 
         forvalues replication = 1/1000 {
-          	if mod(`replication',10)==0 {
+          	if mod(`replication',100)==0 {
           	di "replication = `replication'"
             }
             qui local difi = runiformint(1,4)
@@ -239,6 +242,7 @@ forvalues replication = 1/1000 {
         		use `grp0'
             qui append using `grp1'
         		drop id
+            qui gen dif1 = `difi'
         		qui gen id = _n
         		order(id)
         		qui gen replication = `replication'
@@ -258,7 +262,7 @@ forvalues replication = 1/1000 {
           di "Scenario 5E - N=200"
 
           forvalues replication = 1/1000 {
-            	if mod(`replication',10)==0 {
+            	if mod(`replication',100)==0 {
             	di "replication = `replication'"
               }
               qui local difi = runiformint(1,4)
@@ -294,6 +298,7 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -326,7 +331,7 @@ forvalues replication = 1/1000 {
 di "Scenario 6A - N=200"
 
 forvalues replication = 1/1000 {
-  	if mod(`replication',10)==0 {
+  	if mod(`replication',100)==0 {
   	di "replication = `replication'"
     }
     qui local difi = runiformint(1,4)
@@ -361,6 +366,7 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -380,7 +386,7 @@ forvalues replication = 1/1000 {
   di "Scenario 6B - N=200"
 
   forvalues replication = 1/1000 {
-    	if mod(`replication',10)==0 {
+    	if mod(`replication',100)==0 {
     	di "replication = `replication'"
       }
       qui local difi = runiformint(1,4)
@@ -416,6 +422,7 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -435,7 +442,7 @@ forvalues replication = 1/1000 {
       di "Scenario 6C - N=200"
 
       forvalues replication = 1/1000 {
-        	if mod(`replication',10)==0 {
+        	if mod(`replication',100)==0 {
         	di "replication = `replication'"
           }
           qui local difi = runiformint(1,4)
@@ -471,6 +478,7 @@ forvalues replication = 1/1000 {
       		use `grp0'
           qui append using `grp1'
       		drop id
+          qui gen dif1 = `difi'
       		qui gen id = _n
       		order(id)
       		qui gen replication = `replication'
@@ -490,7 +498,7 @@ forvalues replication = 1/1000 {
           di "Scenario 6D - N=200"
 
           forvalues replication = 1/1000 {
-            	if mod(`replication',10)==0 {
+            	if mod(`replication',100)==0 {
             	di "replication = `replication'"
               }
               qui local difi = runiformint(1,4)
@@ -526,6 +534,7 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -545,7 +554,7 @@ forvalues replication = 1/1000 {
               di "Scenario 6E - N=200"
 
               forvalues replication = 1/1000 {
-                	if mod(`replication',10)==0 {
+                	if mod(`replication',100)==0 {
                 	di "replication = `replication'"
                   }
                   qui local difi = runiformint(1,4)
@@ -581,6 +590,7 @@ forvalues replication = 1/1000 {
               		use `grp0'
                   qui append using `grp1'
               		drop id
+                  qui gen dif1 = `difi'
               		qui gen id = _n
               		order(id)
               		qui gen replication = `replication'
@@ -612,7 +622,7 @@ forvalues replication = 1/1000 {
 di "Scenario 7A - N=200"
 
 forvalues replication = 1/1000 {
-  	if mod(`replication',10)==0 {
+  	if mod(`replication',100)==0 {
   	di "replication = `replication'"
     }
     qui local difi = runiformint(1,4)
@@ -648,6 +658,7 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -668,7 +679,7 @@ forvalues replication = 1/1000 {
   di "Scenario 7B - N=200"
 
   forvalues replication = 1/1000 {
-    	if mod(`replication',10)==0 {
+    	if mod(`replication',100)==0 {
     	di "replication = `replication'"
       }
       qui local difi = runiformint(1,4)
@@ -704,6 +715,7 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -723,7 +735,7 @@ forvalues replication = 1/1000 {
     di "Scenario 7C - N=200"
 
     forvalues replication = 1/1000 {
-      	if mod(`replication',10)==0 {
+      	if mod(`replication',100)==0 {
       	di "replication = `replication'"
         }
         qui local difi = runiformint(1,4)
@@ -759,6 +771,7 @@ forvalues replication = 1/1000 {
     		use `grp0'
         qui append using `grp1'
     		drop id
+        qui gen dif1 = `difi'
     		qui gen id = _n
     		order(id)
     		qui gen replication = `replication'
@@ -779,7 +792,7 @@ forvalues replication = 1/1000 {
         di "Scenario 7D - N=200"
 
         forvalues replication = 1/1000 {
-          	if mod(`replication',10)==0 {
+          	if mod(`replication',100)==0 {
           	di "replication = `replication'"
             }
             qui local difi = runiformint(1,4)
@@ -815,6 +828,7 @@ forvalues replication = 1/1000 {
         		use `grp0'
             qui append using `grp1'
         		drop id
+            qui gen dif1 = `difi'
         		qui gen id = _n
         		order(id)
         		qui gen replication = `replication'
@@ -834,7 +848,7 @@ forvalues replication = 1/1000 {
           di "Scenario 7E - N=200"
 
           forvalues replication = 1/1000 {
-            	if mod(`replication',10)==0 {
+            	if mod(`replication',100)==0 {
             	di "replication = `replication'"
               }
               qui local difi = runiformint(1,4)
@@ -870,6 +884,7 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -902,7 +917,7 @@ forvalues replication = 1/1000 {
 di "Scenario 8A - N=200"
 
 forvalues replication = 1/1000 {
-  	if mod(`replication',10)==0 {
+  	if mod(`replication',100)==0 {
   	di "replication = `replication'"
     }
     qui local difi = runiformint(1,4)
@@ -938,6 +953,7 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -957,7 +973,7 @@ forvalues replication = 1/1000 {
   di "Scenario 8B - N=200"
 
   forvalues replication = 1/1000 {
-    	if mod(`replication',10)==0 {
+    	if mod(`replication',100)==0 {
     	di "replication = `replication'"
       }
       qui local difi = runiformint(1,4)
@@ -993,6 +1009,7 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -1012,7 +1029,7 @@ forvalues replication = 1/1000 {
       di "Scenario 8C - N=200"
 
       forvalues replication = 1/1000 {
-        	if mod(`replication',10)==0 {
+        	if mod(`replication',100)==0 {
         	di "replication = `replication'"
           }
           qui local difi = runiformint(1,4)
@@ -1048,6 +1065,7 @@ forvalues replication = 1/1000 {
       		use `grp0'
           qui append using `grp1'
       		drop id
+          qui gen dif1 = `difi'
       		qui gen id = _n
       		order(id)
       		qui gen replication = `replication'
@@ -1067,7 +1085,7 @@ forvalues replication = 1/1000 {
           di "Scenario 8D - N=200"
 
           forvalues replication = 1/1000 {
-            	if mod(`replication',10)==0 {
+            	if mod(`replication',100)==0 {
             	di "replication = `replication'"
               }
               qui local difi = runiformint(1,4)
@@ -1103,6 +1121,7 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -1122,7 +1141,7 @@ forvalues replication = 1/1000 {
               di "Scenario 8E - N=200"
 
               forvalues replication = 1/1000 {
-                	if mod(`replication',10)==0 {
+                	if mod(`replication',100)==0 {
                 	di "replication = `replication'"
                   }
                   qui local difi = runiformint(1,4)
@@ -1158,6 +1177,7 @@ forvalues replication = 1/1000 {
               		use `grp0'
                   qui append using `grp1'
               		drop id
+                  qui gen dif1 = `difi'
               		qui gen id = _n
               		order(id)
               		qui gen replication = `replication'
@@ -1197,7 +1217,7 @@ di "Scenario 9A - N=200"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 4
@@ -1251,6 +1271,8 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -1272,7 +1294,7 @@ forvalues replication = 1/1000 {
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 4
@@ -1326,6 +1348,8 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -1346,7 +1370,7 @@ forvalues replication = 1/1000 {
 
     forvalues replication = 1/1000 {
       clear
-      if mod(`replication',10)==0 {
+      if mod(`replication',100)==0 {
       di "replication = `replication'"
       }
       qui set obs 4
@@ -1400,6 +1424,8 @@ forvalues replication = 1/1000 {
     		use `grp0'
         qui append using `grp1'
     		drop id
+        qui gen dif1 = `difi'
+        qui gen dif2 = `difj'
     		qui gen id = _n
     		order(id)
     		qui gen replication = `replication'
@@ -1421,7 +1447,7 @@ forvalues replication = 1/1000 {
 
         forvalues replication = 1/1000 {
           clear
-          if mod(`replication',10)==0 {
+          if mod(`replication',100)==0 {
           di "replication = `replication'"
           }
           qui set obs 4
@@ -1475,6 +1501,8 @@ forvalues replication = 1/1000 {
         		use `grp0'
             qui append using `grp1'
         		drop id
+            qui gen dif1 = `difi'
+            qui gen dif2 = `difj'
         		qui gen id = _n
         		order(id)
         		qui gen replication = `replication'
@@ -1495,7 +1523,7 @@ forvalues replication = 1/1000 {
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 4
@@ -1549,6 +1577,8 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -1576,7 +1606,7 @@ di "Scenario 10A - N=200"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 4
@@ -1630,6 +1660,8 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -1651,7 +1683,7 @@ forvalues replication = 1/1000 {
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 4
@@ -1705,6 +1737,8 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -1725,7 +1759,7 @@ forvalues replication = 1/1000 {
 
     forvalues replication = 1/1000 {
       clear
-      if mod(`replication',10)==0 {
+      if mod(`replication',100)==0 {
       di "replication = `replication'"
       }
       qui set obs 4
@@ -1779,6 +1813,8 @@ forvalues replication = 1/1000 {
     		use `grp0'
         qui append using `grp1'
     		drop id
+        qui gen dif1 = `difi'
+        qui gen dif2 = `difj'
     		qui gen id = _n
     		order(id)
     		qui gen replication = `replication'
@@ -1800,7 +1836,7 @@ forvalues replication = 1/1000 {
 
         forvalues replication = 1/1000 {
           clear
-          if mod(`replication',10)==0 {
+          if mod(`replication',100)==0 {
           di "replication = `replication'"
           }
           qui set obs 4
@@ -1854,6 +1890,8 @@ forvalues replication = 1/1000 {
         		use `grp0'
             qui append using `grp1'
         		drop id
+            qui gen dif1 = `difi'
+            qui gen dif2 = `difj'
         		qui gen id = _n
         		order(id)
         		qui gen replication = `replication'
@@ -1874,7 +1912,7 @@ forvalues replication = 1/1000 {
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 4
@@ -1928,6 +1966,8 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -1952,7 +1992,7 @@ di "Scenario 11A - N=200"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 4
@@ -2006,6 +2046,8 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -2027,7 +2069,7 @@ forvalues replication = 1/1000 {
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 4
@@ -2081,6 +2123,8 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -2101,7 +2145,7 @@ forvalues replication = 1/1000 {
 
     forvalues replication = 1/1000 {
       clear
-      if mod(`replication',10)==0 {
+      if mod(`replication',100)==0 {
       di "replication = `replication'"
       }
       qui set obs 4
@@ -2155,6 +2199,8 @@ forvalues replication = 1/1000 {
     		use `grp0'
         qui append using `grp1'
     		drop id
+        qui gen dif1 = `difi'
+        qui gen dif2 = `difj'
     		qui gen id = _n
     		order(id)
     		qui gen replication = `replication'
@@ -2176,7 +2222,7 @@ forvalues replication = 1/1000 {
 
         forvalues replication = 1/1000 {
           clear
-          if mod(`replication',10)==0 {
+          if mod(`replication',100)==0 {
           di "replication = `replication'"
           }
           qui set obs 4
@@ -2230,6 +2276,8 @@ forvalues replication = 1/1000 {
         		use `grp0'
             qui append using `grp1'
         		drop id
+            qui gen dif1 = `difi'
+            qui gen dif2 = `difj'
         		qui gen id = _n
         		order(id)
         		qui gen replication = `replication'
@@ -2250,7 +2298,7 @@ forvalues replication = 1/1000 {
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 4
@@ -2304,6 +2352,8 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -2331,7 +2381,7 @@ di "Scenario 12A - N=200"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 4
@@ -2385,6 +2435,8 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -2406,7 +2458,7 @@ forvalues replication = 1/1000 {
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 4
@@ -2460,6 +2512,8 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -2480,7 +2534,7 @@ forvalues replication = 1/1000 {
 
     forvalues replication = 1/1000 {
       clear
-      if mod(`replication',10)==0 {
+      if mod(`replication',100)==0 {
       di "replication = `replication'"
       }
       qui set obs 4
@@ -2534,6 +2588,8 @@ forvalues replication = 1/1000 {
     		use `grp0'
         qui append using `grp1'
     		drop id
+        qui gen dif1 = `difi'
+        qui gen dif2 = `difj'
     		qui gen id = _n
     		order(id)
     		qui gen replication = `replication'
@@ -2555,7 +2611,7 @@ forvalues replication = 1/1000 {
 
         forvalues replication = 1/1000 {
           clear
-          if mod(`replication',10)==0 {
+          if mod(`replication',100)==0 {
           di "replication = `replication'"
           }
           qui set obs 4
@@ -2609,6 +2665,8 @@ forvalues replication = 1/1000 {
         		use `grp0'
             qui append using `grp1'
         		drop id
+            qui gen dif1 = `difi'
+            qui gen dif2 = `difj'
         		qui gen id = _n
         		order(id)
         		qui gen replication = `replication'
@@ -2629,7 +2687,7 @@ forvalues replication = 1/1000 {
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 4
@@ -2683,6 +2741,8 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
diff --git a/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_300.do b/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_300.do
index 00b4e19..41a11e9 100644
--- a/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_300.do
+++ b/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_300.do
@@ -37,7 +37,7 @@ local Nn = 300
 di "Scenario 5A - N=300"
 
 forvalues replication = 1/1000 {
-  	if mod(`replication',10)==0 {
+  	if mod(`replication',100)==0 {
   	di "replication = `replication'"
     }
     qui local difi = runiformint(1,4)
@@ -72,6 +72,7 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -92,7 +93,7 @@ forvalues replication = 1/1000 {
   di "Scenario 5B - N=300"
 
   forvalues replication = 1/1000 {
-    	if mod(`replication',10)==0 {
+    	if mod(`replication',100)==0 {
     	di "replication = `replication'"
       }
       qui local difi = runiformint(1,4)
@@ -128,6 +129,7 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -147,7 +149,7 @@ forvalues replication = 1/1000 {
     di "Scenario 5C - N=300"
 
     forvalues replication = 1/1000 {
-      	if mod(`replication',10)==0 {
+      	if mod(`replication',100)==0 {
       	di "replication = `replication'"
         }
         qui local difi = runiformint(1,4)
@@ -183,6 +185,7 @@ forvalues replication = 1/1000 {
     		use `grp0'
         qui append using `grp1'
     		drop id
+        qui gen dif1 = `difi'
     		qui gen id = _n
     		order(id)
     		qui gen replication = `replication'
@@ -203,7 +206,7 @@ forvalues replication = 1/1000 {
         di "Scenario 5D - N=300"
 
         forvalues replication = 1/1000 {
-          	if mod(`replication',10)==0 {
+          	if mod(`replication',100)==0 {
           	di "replication = `replication'"
             }
             qui local difi = runiformint(1,4)
@@ -239,6 +242,7 @@ forvalues replication = 1/1000 {
         		use `grp0'
             qui append using `grp1'
         		drop id
+            qui gen dif1 = `difi'
         		qui gen id = _n
         		order(id)
         		qui gen replication = `replication'
@@ -258,7 +262,7 @@ forvalues replication = 1/1000 {
           di "Scenario 5E - N=300"
 
           forvalues replication = 1/1000 {
-            	if mod(`replication',10)==0 {
+            	if mod(`replication',100)==0 {
             	di "replication = `replication'"
               }
               qui local difi = runiformint(1,4)
@@ -294,6 +298,7 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -326,7 +331,7 @@ forvalues replication = 1/1000 {
 di "Scenario 6A - N=300"
 
 forvalues replication = 1/1000 {
-  	if mod(`replication',10)==0 {
+  	if mod(`replication',100)==0 {
   	di "replication = `replication'"
     }
     qui local difi = runiformint(1,4)
@@ -361,6 +366,7 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -380,7 +386,7 @@ forvalues replication = 1/1000 {
   di "Scenario 6B - N=300"
 
   forvalues replication = 1/1000 {
-    	if mod(`replication',10)==0 {
+    	if mod(`replication',100)==0 {
     	di "replication = `replication'"
       }
       qui local difi = runiformint(1,4)
@@ -416,6 +422,7 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -435,7 +442,7 @@ forvalues replication = 1/1000 {
       di "Scenario 6C - N=300"
 
       forvalues replication = 1/1000 {
-        	if mod(`replication',10)==0 {
+        	if mod(`replication',100)==0 {
         	di "replication = `replication'"
           }
           qui local difi = runiformint(1,4)
@@ -471,6 +478,7 @@ forvalues replication = 1/1000 {
       		use `grp0'
           qui append using `grp1'
       		drop id
+          qui gen dif1 = `difi'
       		qui gen id = _n
       		order(id)
       		qui gen replication = `replication'
@@ -490,7 +498,7 @@ forvalues replication = 1/1000 {
           di "Scenario 6D - N=300"
 
           forvalues replication = 1/1000 {
-            	if mod(`replication',10)==0 {
+            	if mod(`replication',100)==0 {
             	di "replication = `replication'"
               }
               qui local difi = runiformint(1,4)
@@ -526,6 +534,7 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -545,7 +554,7 @@ forvalues replication = 1/1000 {
               di "Scenario 6E - N=300"
 
               forvalues replication = 1/1000 {
-                	if mod(`replication',10)==0 {
+                	if mod(`replication',100)==0 {
                 	di "replication = `replication'"
                   }
                   qui local difi = runiformint(1,4)
@@ -581,6 +590,7 @@ forvalues replication = 1/1000 {
               		use `grp0'
                   qui append using `grp1'
               		drop id
+                  qui gen dif1 = `difi'
               		qui gen id = _n
               		order(id)
               		qui gen replication = `replication'
@@ -612,7 +622,7 @@ forvalues replication = 1/1000 {
 di "Scenario 7A - N=300"
 
 forvalues replication = 1/1000 {
-  	if mod(`replication',10)==0 {
+  	if mod(`replication',100)==0 {
   	di "replication = `replication'"
     }
     qui local difi = runiformint(1,4)
@@ -648,6 +658,7 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -668,7 +679,7 @@ forvalues replication = 1/1000 {
   di "Scenario 7B - N=300"
 
   forvalues replication = 1/1000 {
-    	if mod(`replication',10)==0 {
+    	if mod(`replication',100)==0 {
     	di "replication = `replication'"
       }
       qui local difi = runiformint(1,4)
@@ -704,6 +715,7 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -723,7 +735,7 @@ forvalues replication = 1/1000 {
     di "Scenario 7C - N=300"
 
     forvalues replication = 1/1000 {
-      	if mod(`replication',10)==0 {
+      	if mod(`replication',100)==0 {
       	di "replication = `replication'"
         }
         qui local difi = runiformint(1,4)
@@ -759,6 +771,7 @@ forvalues replication = 1/1000 {
     		use `grp0'
         qui append using `grp1'
     		drop id
+        qui gen dif1 = `difi'
     		qui gen id = _n
     		order(id)
     		qui gen replication = `replication'
@@ -779,7 +792,7 @@ forvalues replication = 1/1000 {
         di "Scenario 7D - N=300"
 
         forvalues replication = 1/1000 {
-          	if mod(`replication',10)==0 {
+          	if mod(`replication',100)==0 {
           	di "replication = `replication'"
             }
             qui local difi = runiformint(1,4)
@@ -815,6 +828,7 @@ forvalues replication = 1/1000 {
         		use `grp0'
             qui append using `grp1'
         		drop id
+            qui gen dif1 = `difi'
         		qui gen id = _n
         		order(id)
         		qui gen replication = `replication'
@@ -834,7 +848,7 @@ forvalues replication = 1/1000 {
           di "Scenario 7E - N=300"
 
           forvalues replication = 1/1000 {
-            	if mod(`replication',10)==0 {
+            	if mod(`replication',100)==0 {
             	di "replication = `replication'"
               }
               qui local difi = runiformint(1,4)
@@ -870,6 +884,7 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -902,7 +917,7 @@ forvalues replication = 1/1000 {
 di "Scenario 8A - N=300"
 
 forvalues replication = 1/1000 {
-  	if mod(`replication',10)==0 {
+  	if mod(`replication',100)==0 {
   	di "replication = `replication'"
     }
     qui local difi = runiformint(1,4)
@@ -938,6 +953,7 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -957,7 +973,7 @@ forvalues replication = 1/1000 {
   di "Scenario 8B - N=300"
 
   forvalues replication = 1/1000 {
-    	if mod(`replication',10)==0 {
+    	if mod(`replication',100)==0 {
     	di "replication = `replication'"
       }
       qui local difi = runiformint(1,4)
@@ -993,6 +1009,7 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -1012,7 +1029,7 @@ forvalues replication = 1/1000 {
       di "Scenario 8C - N=300"
 
       forvalues replication = 1/1000 {
-        	if mod(`replication',10)==0 {
+        	if mod(`replication',100)==0 {
         	di "replication = `replication'"
           }
           qui local difi = runiformint(1,4)
@@ -1048,6 +1065,7 @@ forvalues replication = 1/1000 {
       		use `grp0'
           qui append using `grp1'
       		drop id
+          qui gen dif1 = `difi'
       		qui gen id = _n
       		order(id)
       		qui gen replication = `replication'
@@ -1067,7 +1085,7 @@ forvalues replication = 1/1000 {
           di "Scenario 8D - N=300"
 
           forvalues replication = 1/1000 {
-            	if mod(`replication',10)==0 {
+            	if mod(`replication',100)==0 {
             	di "replication = `replication'"
               }
               qui local difi = runiformint(1,4)
@@ -1103,6 +1121,7 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -1122,7 +1141,7 @@ forvalues replication = 1/1000 {
               di "Scenario 8E - N=300"
 
               forvalues replication = 1/1000 {
-                	if mod(`replication',10)==0 {
+                	if mod(`replication',100)==0 {
                 	di "replication = `replication'"
                   }
                   qui local difi = runiformint(1,4)
@@ -1158,6 +1177,7 @@ forvalues replication = 1/1000 {
               		use `grp0'
                   qui append using `grp1'
               		drop id
+                  qui gen dif1 = `difi'
               		qui gen id = _n
               		order(id)
               		qui gen replication = `replication'
@@ -1197,7 +1217,7 @@ di "Scenario 9A - N=300"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 4
@@ -1251,6 +1271,8 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -1272,7 +1294,7 @@ forvalues replication = 1/1000 {
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 4
@@ -1326,6 +1348,8 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -1346,7 +1370,7 @@ forvalues replication = 1/1000 {
 
     forvalues replication = 1/1000 {
       clear
-      if mod(`replication',10)==0 {
+      if mod(`replication',100)==0 {
       di "replication = `replication'"
       }
       qui set obs 4
@@ -1400,6 +1424,8 @@ forvalues replication = 1/1000 {
     		use `grp0'
         qui append using `grp1'
     		drop id
+        qui gen dif1 = `difi'
+        qui gen dif2 = `difj'
     		qui gen id = _n
     		order(id)
     		qui gen replication = `replication'
@@ -1421,7 +1447,7 @@ forvalues replication = 1/1000 {
 
         forvalues replication = 1/1000 {
           clear
-          if mod(`replication',10)==0 {
+          if mod(`replication',100)==0 {
           di "replication = `replication'"
           }
           qui set obs 4
@@ -1475,6 +1501,8 @@ forvalues replication = 1/1000 {
         		use `grp0'
             qui append using `grp1'
         		drop id
+            qui gen dif1 = `difi'
+            qui gen dif2 = `difj'
         		qui gen id = _n
         		order(id)
         		qui gen replication = `replication'
@@ -1495,7 +1523,7 @@ forvalues replication = 1/1000 {
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 4
@@ -1549,6 +1577,8 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -1576,7 +1606,7 @@ di "Scenario 10A - N=300"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 4
@@ -1630,6 +1660,8 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -1651,7 +1683,7 @@ forvalues replication = 1/1000 {
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 4
@@ -1705,6 +1737,8 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -1725,7 +1759,7 @@ forvalues replication = 1/1000 {
 
     forvalues replication = 1/1000 {
       clear
-      if mod(`replication',10)==0 {
+      if mod(`replication',100)==0 {
       di "replication = `replication'"
       }
       qui set obs 4
@@ -1779,6 +1813,8 @@ forvalues replication = 1/1000 {
     		use `grp0'
         qui append using `grp1'
     		drop id
+        qui gen dif1 = `difi'
+        qui gen dif2 = `difj'
     		qui gen id = _n
     		order(id)
     		qui gen replication = `replication'
@@ -1800,7 +1836,7 @@ forvalues replication = 1/1000 {
 
         forvalues replication = 1/1000 {
           clear
-          if mod(`replication',10)==0 {
+          if mod(`replication',100)==0 {
           di "replication = `replication'"
           }
           qui set obs 4
@@ -1854,6 +1890,8 @@ forvalues replication = 1/1000 {
         		use `grp0'
             qui append using `grp1'
         		drop id
+            qui gen dif1 = `difi'
+            qui gen dif2 = `difj'
         		qui gen id = _n
         		order(id)
         		qui gen replication = `replication'
@@ -1874,7 +1912,7 @@ forvalues replication = 1/1000 {
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 4
@@ -1928,6 +1966,8 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -1952,7 +1992,7 @@ di "Scenario 11A - N=300"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 4
@@ -2006,6 +2046,8 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -2027,7 +2069,7 @@ forvalues replication = 1/1000 {
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 4
@@ -2081,6 +2123,8 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -2101,7 +2145,7 @@ forvalues replication = 1/1000 {
 
     forvalues replication = 1/1000 {
       clear
-      if mod(`replication',10)==0 {
+      if mod(`replication',100)==0 {
       di "replication = `replication'"
       }
       qui set obs 4
@@ -2155,6 +2199,8 @@ forvalues replication = 1/1000 {
     		use `grp0'
         qui append using `grp1'
     		drop id
+        qui gen dif1 = `difi'
+        qui gen dif2 = `difj'
     		qui gen id = _n
     		order(id)
     		qui gen replication = `replication'
@@ -2176,7 +2222,7 @@ forvalues replication = 1/1000 {
 
         forvalues replication = 1/1000 {
           clear
-          if mod(`replication',10)==0 {
+          if mod(`replication',100)==0 {
           di "replication = `replication'"
           }
           qui set obs 4
@@ -2230,6 +2276,8 @@ forvalues replication = 1/1000 {
         		use `grp0'
             qui append using `grp1'
         		drop id
+            qui gen dif1 = `difi'
+            qui gen dif2 = `difj'
         		qui gen id = _n
         		order(id)
         		qui gen replication = `replication'
@@ -2250,7 +2298,7 @@ forvalues replication = 1/1000 {
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 4
@@ -2304,6 +2352,8 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -2331,7 +2381,7 @@ di "Scenario 12A - N=300"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 4
@@ -2385,6 +2435,8 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -2406,7 +2458,7 @@ forvalues replication = 1/1000 {
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 4
@@ -2460,6 +2512,8 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -2480,7 +2534,7 @@ forvalues replication = 1/1000 {
 
     forvalues replication = 1/1000 {
       clear
-      if mod(`replication',10)==0 {
+      if mod(`replication',100)==0 {
       di "replication = `replication'"
       }
       qui set obs 4
@@ -2534,6 +2588,8 @@ forvalues replication = 1/1000 {
     		use `grp0'
         qui append using `grp1'
     		drop id
+        qui gen dif1 = `difi'
+        qui gen dif2 = `difj'
     		qui gen id = _n
     		order(id)
     		qui gen replication = `replication'
@@ -2555,7 +2611,7 @@ forvalues replication = 1/1000 {
 
         forvalues replication = 1/1000 {
           clear
-          if mod(`replication',10)==0 {
+          if mod(`replication',100)==0 {
           di "replication = `replication'"
           }
           qui set obs 4
@@ -2609,6 +2665,8 @@ forvalues replication = 1/1000 {
         		use `grp0'
             qui append using `grp1'
         		drop id
+            qui gen dif1 = `difi'
+            qui gen dif2 = `difj'
         		qui gen id = _n
         		order(id)
         		qui gen replication = `replication'
@@ -2629,7 +2687,7 @@ forvalues replication = 1/1000 {
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 4
@@ -2683,6 +2741,8 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
diff --git a/Scripts/Scenarios/DIF/scenarios_DIF_J7_baseline_100.do b/Scripts/Scenarios/DIF/scenarios_DIF_J7_baseline_100.do
index 7f5bd3a..d8c8c94 100644
--- a/Scripts/Scenarios/DIF/scenarios_DIF_J7_baseline_100.do
+++ b/Scripts/Scenarios/DIF/scenarios_DIF_J7_baseline_100.do
@@ -39,7 +39,7 @@ di "Scenario 13A - N=100"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 7
@@ -114,6 +114,8 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -135,7 +137,7 @@ forvalues replication = 1/1000 {
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 7
@@ -208,6 +210,8 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -228,7 +232,7 @@ forvalues replication = 1/1000 {
 
     forvalues replication = 1/1000 {
       clear
-      if mod(`replication',10)==0 {
+      if mod(`replication',100)==0 {
       di "replication = `replication'"
       }
       qui set obs 7
@@ -301,6 +305,8 @@ forvalues replication = 1/1000 {
     		use `grp0'
         qui append using `grp1'
     		drop id
+        qui gen dif1 = `difi'
+        qui gen dif2 = `difj'
     		qui gen id = _n
     		order(id)
     		qui gen replication = `replication'
@@ -322,7 +328,7 @@ forvalues replication = 1/1000 {
 
         forvalues replication = 1/1000 {
           clear
-          if mod(`replication',10)==0 {
+          if mod(`replication',100)==0 {
           di "replication = `replication'"
           }
           qui set obs 7
@@ -394,6 +400,8 @@ forvalues replication = 1/1000 {
         		use `grp0'
             qui append using `grp1'
         		drop id
+            qui gen dif1 = `difi'
+            qui gen dif2 = `difj'
         		qui gen id = _n
         		order(id)
         		qui gen replication = `replication'
@@ -414,7 +422,7 @@ forvalues replication = 1/1000 {
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 7
@@ -486,6 +494,8 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -511,7 +521,7 @@ di "Scenario 14A - N=100"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 7
@@ -583,6 +593,8 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -604,7 +616,7 @@ forvalues replication = 1/1000 {
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 7
@@ -676,6 +688,8 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -696,7 +710,7 @@ forvalues replication = 1/1000 {
 
     forvalues replication = 1/1000 {
       clear
-      if mod(`replication',10)==0 {
+      if mod(`replication',100)==0 {
       di "replication = `replication'"
       }
       qui set obs 7
@@ -768,6 +782,8 @@ forvalues replication = 1/1000 {
     		use `grp0'
         qui append using `grp1'
     		drop id
+        qui gen dif1 = `difi'
+        qui gen dif2 = `difj'
     		qui gen id = _n
     		order(id)
     		qui gen replication = `replication'
@@ -789,7 +805,7 @@ forvalues replication = 1/1000 {
 
         forvalues replication = 1/1000 {
           clear
-          if mod(`replication',10)==0 {
+          if mod(`replication',100)==0 {
           di "replication = `replication'"
           }
           qui set obs 7
@@ -861,6 +877,8 @@ forvalues replication = 1/1000 {
         		use `grp0'
             qui append using `grp1'
         		drop id
+            qui gen dif1 = `difi'
+            qui gen dif2 = `difj'
         		qui gen id = _n
         		order(id)
         		qui gen replication = `replication'
@@ -881,7 +899,7 @@ forvalues replication = 1/1000 {
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 7
@@ -953,6 +971,8 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -978,7 +998,7 @@ di "Scenario 15A - N=100"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 7
@@ -1049,6 +1069,8 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -1070,7 +1092,7 @@ forvalues replication = 1/1000 {
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 7
@@ -1142,6 +1164,8 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -1162,7 +1186,7 @@ forvalues replication = 1/1000 {
 
     forvalues replication = 1/1000 {
       clear
-      if mod(`replication',10)==0 {
+      if mod(`replication',100)==0 {
       di "replication = `replication'"
       }
       qui set obs 7
@@ -1234,6 +1258,8 @@ forvalues replication = 1/1000 {
     		use `grp0'
         qui append using `grp1'
     		drop id
+        qui gen dif1 = `difi'
+        qui gen dif2 = `difj'
     		qui gen id = _n
     		order(id)
     		qui gen replication = `replication'
@@ -1255,7 +1281,7 @@ forvalues replication = 1/1000 {
 
         forvalues replication = 1/1000 {
           clear
-          if mod(`replication',10)==0 {
+          if mod(`replication',100)==0 {
           di "replication = `replication'"
           }
           qui set obs 7
@@ -1327,6 +1353,8 @@ forvalues replication = 1/1000 {
         		use `grp0'
             qui append using `grp1'
         		drop id
+            qui gen dif1 = `difi'
+            qui gen dif2 = `difj'
         		qui gen id = _n
         		order(id)
         		qui gen replication = `replication'
@@ -1347,7 +1375,7 @@ forvalues replication = 1/1000 {
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 7
@@ -1419,6 +1447,8 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -1444,7 +1474,7 @@ di "Scenario 16A - N=100"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 7
@@ -1516,6 +1546,8 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -1537,7 +1569,7 @@ forvalues replication = 1/1000 {
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 7
@@ -1609,6 +1641,8 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -1629,7 +1663,7 @@ forvalues replication = 1/1000 {
 
     forvalues replication = 1/1000 {
       clear
-      if mod(`replication',10)==0 {
+      if mod(`replication',100)==0 {
       di "replication = `replication'"
       }
       qui set obs 7
@@ -1701,6 +1735,8 @@ forvalues replication = 1/1000 {
     		use `grp0'
         qui append using `grp1'
     		drop id
+        qui gen dif1 = `difi'
+        qui gen dif2 = `difj'
     		qui gen id = _n
     		order(id)
     		qui gen replication = `replication'
@@ -1722,7 +1758,7 @@ forvalues replication = 1/1000 {
 
         forvalues replication = 1/1000 {
           clear
-          if mod(`replication',10)==0 {
+          if mod(`replication',100)==0 {
           di "replication = `replication'"
           }
           qui set obs 7
@@ -1794,6 +1830,8 @@ forvalues replication = 1/1000 {
         		use `grp0'
             qui append using `grp1'
         		drop id
+            qui gen dif1 = `difi'
+            qui gen dif2 = `difj'
         		qui gen id = _n
         		order(id)
         		qui gen replication = `replication'
@@ -1814,7 +1852,7 @@ forvalues replication = 1/1000 {
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 7
@@ -1886,6 +1924,8 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -1915,7 +1955,7 @@ di "Scenario 17A - N=100"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 7
@@ -2011,6 +2051,9 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
+    qui gen dif3 = `difk'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -2032,7 +2075,7 @@ di "Scenario 17B - N=100"
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 7
@@ -2128,6 +2171,9 @@ di "Scenario 17B - N=100"
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
+      qui gen dif3 = `difk'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -2158,7 +2204,7 @@ di "Scenario 17B - N=100"
 
       forvalues replication = 1/1000 {
         clear
-        if mod(`replication',10)==0 {
+        if mod(`replication',100)==0 {
         di "replication = `replication'"
         }
         qui set obs 7
@@ -2253,6 +2299,9 @@ di "Scenario 17B - N=100"
       		use `grp0'
           qui append using `grp1'
       		drop id
+          qui gen dif1 = `difi'
+          qui gen dif2 = `difj'
+          qui gen dif3 = `difk'
       		qui gen id = _n
       		order(id)
       		qui gen replication = `replication'
@@ -2279,7 +2328,7 @@ di "Scenario 17B - N=100"
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 7
@@ -2375,6 +2424,9 @@ di "Scenario 17B - N=100"
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
+              qui gen dif3 = `difk'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -2405,7 +2457,7 @@ di "Scenario 17B - N=100"
 
               forvalues replication = 1/1000 {
                   clear
-                	if mod(`replication',10)==0 {
+                	if mod(`replication',100)==0 {
                 	di "replication = `replication'"
                   }
                   qui set obs 7
@@ -2498,6 +2550,9 @@ di "Scenario 17B - N=100"
               		use `grp0'
                   qui append using `grp1'
               		drop id
+                  qui gen dif1 = `difi'
+                  qui gen dif2 = `difj'
+                  qui gen dif3 = `difk'
               		qui gen id = _n
               		order(id)
               		qui gen replication = `replication'
@@ -2525,7 +2580,7 @@ di "Scenario 18A - N=100"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 7
@@ -2621,6 +2676,9 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
+    qui gen dif3 = `difk'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -2642,7 +2700,7 @@ di "Scenario 18B - N=100"
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 7
@@ -2738,6 +2796,9 @@ di "Scenario 18B - N=100"
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
+      qui gen dif3 = `difk'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -2768,7 +2829,7 @@ di "Scenario 18B - N=100"
 
       forvalues replication = 1/1000 {
         clear
-        if mod(`replication',10)==0 {
+        if mod(`replication',100)==0 {
         di "replication = `replication'"
         }
         qui set obs 7
@@ -2863,6 +2924,9 @@ di "Scenario 18B - N=100"
       		use `grp0'
           qui append using `grp1'
       		drop id
+          qui gen dif1 = `difi'
+          qui gen dif2 = `difj'
+          qui gen dif3 = `difk'
       		qui gen id = _n
       		order(id)
       		qui gen replication = `replication'
@@ -2889,7 +2953,7 @@ di "Scenario 18B - N=100"
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 7
@@ -2985,6 +3049,9 @@ di "Scenario 18B - N=100"
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
+              qui gen dif3 = `difk'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -3015,7 +3082,7 @@ di "Scenario 18B - N=100"
 
               forvalues replication = 1/1000 {
                   clear
-                	if mod(`replication',10)==0 {
+                	if mod(`replication',100)==0 {
                 	di "replication = `replication'"
                   }
                   qui set obs 7
@@ -3108,6 +3175,9 @@ di "Scenario 18B - N=100"
               		use `grp0'
                   qui append using `grp1'
               		drop id
+                  qui gen dif1 = `difi'
+                  qui gen dif2 = `difj'
+                  qui gen dif3 = `difk'
               		qui gen id = _n
               		order(id)
               		qui gen replication = `replication'
@@ -3143,7 +3213,7 @@ di "Scenario 19A - N=100"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 7
@@ -3239,6 +3309,9 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
+    qui gen dif3 = `difk'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -3260,7 +3333,7 @@ di "Scenario 19B - N=100"
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 7
@@ -3356,6 +3429,9 @@ di "Scenario 19B - N=100"
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
+      qui gen dif3 = `difk'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -3378,7 +3454,7 @@ di "Scenario 19B - N=100"
 
       forvalues replication = 1/1000 {
         clear
-        if mod(`replication',10)==0 {
+        if mod(`replication',100)==0 {
         di "replication = `replication'"
         }
         qui set obs 7
@@ -3473,6 +3549,9 @@ di "Scenario 19B - N=100"
       		use `grp0'
           qui append using `grp1'
       		drop id
+          qui gen dif1 = `difi'
+          qui gen dif2 = `difj'
+          qui gen dif3 = `difk'
       		qui gen id = _n
       		order(id)
       		qui gen replication = `replication'
@@ -3495,7 +3574,7 @@ di "Scenario 19B - N=100"
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 7
@@ -3591,6 +3670,9 @@ di "Scenario 19B - N=100"
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
+              qui gen dif3 = `difk'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -3613,7 +3695,7 @@ di "Scenario 19B - N=100"
 
               forvalues replication = 1/1000 {
                   clear
-                	if mod(`replication',10)==0 {
+                	if mod(`replication',100)==0 {
                 	di "replication = `replication'"
                   }
                   qui set obs 7
@@ -3706,6 +3788,9 @@ di "Scenario 19B - N=100"
               		use `grp0'
                   qui append using `grp1'
               		drop id
+                  qui gen dif1 = `difi'
+                  qui gen dif2 = `difj'
+                  qui gen dif3 = `difk'
               		qui gen id = _n
               		order(id)
               		qui gen replication = `replication'
@@ -3732,7 +3817,7 @@ di "Scenario 20A - N=100"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 7
@@ -3828,6 +3913,9 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
+    qui gen dif3 = `difk'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -3849,7 +3937,7 @@ di "Scenario 20B - N=100"
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 7
@@ -3945,6 +4033,9 @@ di "Scenario 20B - N=100"
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
+      qui gen dif3 = `difk'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -3966,7 +4057,7 @@ di "Scenario 20B - N=100"
 
       forvalues replication = 1/1000 {
         clear
-        if mod(`replication',10)==0 {
+        if mod(`replication',100)==0 {
         di "replication = `replication'"
         }
         qui set obs 7
@@ -4061,6 +4152,9 @@ di "Scenario 20B - N=100"
       		use `grp0'
           qui append using `grp1'
       		drop id
+          qui gen dif1 = `difi'
+          qui gen dif2 = `difj'
+          qui gen dif3 = `difk'
       		qui gen id = _n
       		order(id)
       		qui gen replication = `replication'
@@ -4083,7 +4177,7 @@ di "Scenario 20B - N=100"
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 7
@@ -4179,6 +4273,9 @@ di "Scenario 20B - N=100"
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
+              qui gen dif3 = `difk'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -4201,7 +4298,7 @@ di "Scenario 20B - N=100"
 
               forvalues replication = 1/1000 {
                   clear
-                	if mod(`replication',10)==0 {
+                	if mod(`replication',100)==0 {
                 	di "replication = `replication'"
                   }
                   qui set obs 7
@@ -4294,6 +4391,9 @@ di "Scenario 20B - N=100"
               		use `grp0'
                   qui append using `grp1'
               		drop id
+                  qui gen dif1 = `difi'
+                  qui gen dif2 = `difj'
+                  qui gen dif3 = `difk'
               		qui gen id = _n
               		order(id)
               		qui gen replication = `replication'
diff --git a/Scripts/Scenarios/DIF/scenarios_DIF_J7_baseline_200.do b/Scripts/Scenarios/DIF/scenarios_DIF_J7_baseline_200.do
index d1d5fe9..c38cacc 100644
--- a/Scripts/Scenarios/DIF/scenarios_DIF_J7_baseline_200.do
+++ b/Scripts/Scenarios/DIF/scenarios_DIF_J7_baseline_200.do
@@ -39,7 +39,7 @@ di "Scenario 13A - N=200"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 7
@@ -114,6 +114,8 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -135,7 +137,7 @@ forvalues replication = 1/1000 {
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 7
@@ -208,6 +210,8 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -228,7 +232,7 @@ forvalues replication = 1/1000 {
 
     forvalues replication = 1/1000 {
       clear
-      if mod(`replication',10)==0 {
+      if mod(`replication',100)==0 {
       di "replication = `replication'"
       }
       qui set obs 7
@@ -301,6 +305,8 @@ forvalues replication = 1/1000 {
     		use `grp0'
         qui append using `grp1'
     		drop id
+        qui gen dif1 = `difi'
+        qui gen dif2 = `difj'
     		qui gen id = _n
     		order(id)
     		qui gen replication = `replication'
@@ -322,7 +328,7 @@ forvalues replication = 1/1000 {
 
         forvalues replication = 1/1000 {
           clear
-          if mod(`replication',10)==0 {
+          if mod(`replication',100)==0 {
           di "replication = `replication'"
           }
           qui set obs 7
@@ -394,6 +400,8 @@ forvalues replication = 1/1000 {
         		use `grp0'
             qui append using `grp1'
         		drop id
+            qui gen dif1 = `difi'
+            qui gen dif2 = `difj'
         		qui gen id = _n
         		order(id)
         		qui gen replication = `replication'
@@ -414,7 +422,7 @@ forvalues replication = 1/1000 {
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 7
@@ -486,6 +494,8 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -511,7 +521,7 @@ di "Scenario 14A - N=200"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 7
@@ -583,6 +593,8 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -604,7 +616,7 @@ forvalues replication = 1/1000 {
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 7
@@ -676,6 +688,8 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -696,7 +710,7 @@ forvalues replication = 1/1000 {
 
     forvalues replication = 1/1000 {
       clear
-      if mod(`replication',10)==0 {
+      if mod(`replication',100)==0 {
       di "replication = `replication'"
       }
       qui set obs 7
@@ -768,6 +782,8 @@ forvalues replication = 1/1000 {
     		use `grp0'
         qui append using `grp1'
     		drop id
+        qui gen dif1 = `difi'
+        qui gen dif2 = `difj'
     		qui gen id = _n
     		order(id)
     		qui gen replication = `replication'
@@ -789,7 +805,7 @@ forvalues replication = 1/1000 {
 
         forvalues replication = 1/1000 {
           clear
-          if mod(`replication',10)==0 {
+          if mod(`replication',100)==0 {
           di "replication = `replication'"
           }
           qui set obs 7
@@ -861,6 +877,8 @@ forvalues replication = 1/1000 {
         		use `grp0'
             qui append using `grp1'
         		drop id
+            qui gen dif1 = `difi'
+            qui gen dif2 = `difj'
         		qui gen id = _n
         		order(id)
         		qui gen replication = `replication'
@@ -881,7 +899,7 @@ forvalues replication = 1/1000 {
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 7
@@ -953,6 +971,8 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -978,7 +998,7 @@ di "Scenario 15A - N=200"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 7
@@ -1049,6 +1069,8 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -1070,7 +1092,7 @@ forvalues replication = 1/1000 {
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 7
@@ -1142,6 +1164,8 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -1162,7 +1186,7 @@ forvalues replication = 1/1000 {
 
     forvalues replication = 1/1000 {
       clear
-      if mod(`replication',10)==0 {
+      if mod(`replication',100)==0 {
       di "replication = `replication'"
       }
       qui set obs 7
@@ -1234,6 +1258,8 @@ forvalues replication = 1/1000 {
     		use `grp0'
         qui append using `grp1'
     		drop id
+        qui gen dif1 = `difi'
+        qui gen dif2 = `difj'
     		qui gen id = _n
     		order(id)
     		qui gen replication = `replication'
@@ -1255,7 +1281,7 @@ forvalues replication = 1/1000 {
 
         forvalues replication = 1/1000 {
           clear
-          if mod(`replication',10)==0 {
+          if mod(`replication',100)==0 {
           di "replication = `replication'"
           }
           qui set obs 7
@@ -1327,6 +1353,8 @@ forvalues replication = 1/1000 {
         		use `grp0'
             qui append using `grp1'
         		drop id
+            qui gen dif1 = `difi'
+            qui gen dif2 = `difj'
         		qui gen id = _n
         		order(id)
         		qui gen replication = `replication'
@@ -1347,7 +1375,7 @@ forvalues replication = 1/1000 {
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 7
@@ -1419,6 +1447,8 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -1444,7 +1474,7 @@ di "Scenario 16A - N=200"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 7
@@ -1516,6 +1546,8 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -1537,7 +1569,7 @@ forvalues replication = 1/1000 {
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 7
@@ -1609,6 +1641,8 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -1629,7 +1663,7 @@ forvalues replication = 1/1000 {
 
     forvalues replication = 1/1000 {
       clear
-      if mod(`replication',10)==0 {
+      if mod(`replication',100)==0 {
       di "replication = `replication'"
       }
       qui set obs 7
@@ -1701,6 +1735,8 @@ forvalues replication = 1/1000 {
     		use `grp0'
         qui append using `grp1'
     		drop id
+        qui gen dif1 = `difi'
+        qui gen dif2 = `difj'
     		qui gen id = _n
     		order(id)
     		qui gen replication = `replication'
@@ -1722,7 +1758,7 @@ forvalues replication = 1/1000 {
 
         forvalues replication = 1/1000 {
           clear
-          if mod(`replication',10)==0 {
+          if mod(`replication',100)==0 {
           di "replication = `replication'"
           }
           qui set obs 7
@@ -1794,6 +1830,8 @@ forvalues replication = 1/1000 {
         		use `grp0'
             qui append using `grp1'
         		drop id
+            qui gen dif1 = `difi'
+            qui gen dif2 = `difj'
         		qui gen id = _n
         		order(id)
         		qui gen replication = `replication'
@@ -1814,7 +1852,7 @@ forvalues replication = 1/1000 {
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 7
@@ -1886,6 +1924,8 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -1915,7 +1955,7 @@ di "Scenario 17A - N=200"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 7
@@ -2011,6 +2051,9 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
+    qui gen dif3 = `difk'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -2032,7 +2075,7 @@ di "Scenario 17B - N=200"
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 7
@@ -2128,6 +2171,9 @@ di "Scenario 17B - N=200"
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
+      qui gen dif3 = `difk'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -2158,7 +2204,7 @@ di "Scenario 17B - N=200"
 
       forvalues replication = 1/1000 {
         clear
-        if mod(`replication',10)==0 {
+        if mod(`replication',100)==0 {
         di "replication = `replication'"
         }
         qui set obs 7
@@ -2253,6 +2299,9 @@ di "Scenario 17B - N=200"
       		use `grp0'
           qui append using `grp1'
       		drop id
+          qui gen dif1 = `difi'
+          qui gen dif2 = `difj'
+          qui gen dif3 = `difk'
       		qui gen id = _n
       		order(id)
       		qui gen replication = `replication'
@@ -2279,7 +2328,7 @@ di "Scenario 17B - N=200"
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 7
@@ -2375,6 +2424,9 @@ di "Scenario 17B - N=200"
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
+              qui gen dif3 = `difk'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -2405,7 +2457,7 @@ di "Scenario 17B - N=200"
 
               forvalues replication = 1/1000 {
                   clear
-                	if mod(`replication',10)==0 {
+                	if mod(`replication',100)==0 {
                 	di "replication = `replication'"
                   }
                   qui set obs 7
@@ -2498,6 +2550,9 @@ di "Scenario 17B - N=200"
               		use `grp0'
                   qui append using `grp1'
               		drop id
+                  qui gen dif1 = `difi'
+                  qui gen dif2 = `difj'
+                  qui gen dif3 = `difk'
               		qui gen id = _n
               		order(id)
               		qui gen replication = `replication'
@@ -2525,7 +2580,7 @@ di "Scenario 18A - N=200"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 7
@@ -2621,6 +2676,9 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
+    qui gen dif3 = `difk'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -2642,7 +2700,7 @@ di "Scenario 18B - N=200"
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 7
@@ -2738,6 +2796,9 @@ di "Scenario 18B - N=200"
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
+      qui gen dif3 = `difk'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -2768,7 +2829,7 @@ di "Scenario 18B - N=200"
 
       forvalues replication = 1/1000 {
         clear
-        if mod(`replication',10)==0 {
+        if mod(`replication',100)==0 {
         di "replication = `replication'"
         }
         qui set obs 7
@@ -2863,6 +2924,9 @@ di "Scenario 18B - N=200"
       		use `grp0'
           qui append using `grp1'
       		drop id
+          qui gen dif1 = `difi'
+          qui gen dif2 = `difj'
+          qui gen dif3 = `difk'
       		qui gen id = _n
       		order(id)
       		qui gen replication = `replication'
@@ -2889,7 +2953,7 @@ di "Scenario 18B - N=200"
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 7
@@ -2985,6 +3049,9 @@ di "Scenario 18B - N=200"
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
+              qui gen dif3 = `difk'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -3015,7 +3082,7 @@ di "Scenario 18B - N=200"
 
               forvalues replication = 1/1000 {
                   clear
-                	if mod(`replication',10)==0 {
+                	if mod(`replication',100)==0 {
                 	di "replication = `replication'"
                   }
                   qui set obs 7
@@ -3108,6 +3175,9 @@ di "Scenario 18B - N=200"
               		use `grp0'
                   qui append using `grp1'
               		drop id
+                  qui gen dif1 = `difi'
+                  qui gen dif2 = `difj'
+                  qui gen dif3 = `difk'
               		qui gen id = _n
               		order(id)
               		qui gen replication = `replication'
@@ -3143,7 +3213,7 @@ di "Scenario 19A - N=200"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 7
@@ -3239,6 +3309,9 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
+    qui gen dif3 = `difk'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -3260,7 +3333,7 @@ di "Scenario 19B - N=200"
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 7
@@ -3356,6 +3429,9 @@ di "Scenario 19B - N=200"
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
+      qui gen dif3 = `difk'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -3378,7 +3454,7 @@ di "Scenario 19B - N=200"
 
       forvalues replication = 1/1000 {
         clear
-        if mod(`replication',10)==0 {
+        if mod(`replication',100)==0 {
         di "replication = `replication'"
         }
         qui set obs 7
@@ -3473,6 +3549,9 @@ di "Scenario 19B - N=200"
       		use `grp0'
           qui append using `grp1'
       		drop id
+          qui gen dif1 = `difi'
+          qui gen dif2 = `difj'
+          qui gen dif3 = `difk'
       		qui gen id = _n
       		order(id)
       		qui gen replication = `replication'
@@ -3495,7 +3574,7 @@ di "Scenario 19B - N=200"
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 7
@@ -3591,6 +3670,9 @@ di "Scenario 19B - N=200"
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
+              qui gen dif3 = `difk'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -3613,7 +3695,7 @@ di "Scenario 19B - N=200"
 
               forvalues replication = 1/1000 {
                   clear
-                	if mod(`replication',10)==0 {
+                	if mod(`replication',100)==0 {
                 	di "replication = `replication'"
                   }
                   qui set obs 7
@@ -3706,6 +3788,9 @@ di "Scenario 19B - N=200"
               		use `grp0'
                   qui append using `grp1'
               		drop id
+                  qui gen dif1 = `difi'
+                  qui gen dif2 = `difj'
+                  qui gen dif3 = `difk'
               		qui gen id = _n
               		order(id)
               		qui gen replication = `replication'
@@ -3732,7 +3817,7 @@ di "Scenario 20A - N=200"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 7
@@ -3828,6 +3913,9 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
+    qui gen dif3 = `difk'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -3849,7 +3937,7 @@ di "Scenario 20B - N=200"
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 7
@@ -3945,6 +4033,9 @@ di "Scenario 20B - N=200"
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
+      qui gen dif3 = `difk'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -3966,7 +4057,7 @@ di "Scenario 20B - N=200"
 
       forvalues replication = 1/1000 {
         clear
-        if mod(`replication',10)==0 {
+        if mod(`replication',100)==0 {
         di "replication = `replication'"
         }
         qui set obs 7
@@ -4061,6 +4152,9 @@ di "Scenario 20B - N=200"
       		use `grp0'
           qui append using `grp1'
       		drop id
+          qui gen dif1 = `difi'
+          qui gen dif2 = `difj'
+          qui gen dif3 = `difk'
       		qui gen id = _n
       		order(id)
       		qui gen replication = `replication'
@@ -4083,7 +4177,7 @@ di "Scenario 20B - N=200"
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 7
@@ -4179,6 +4273,9 @@ di "Scenario 20B - N=200"
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
+              qui gen dif3 = `difk'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -4201,7 +4298,7 @@ di "Scenario 20B - N=200"
 
               forvalues replication = 1/1000 {
                   clear
-                	if mod(`replication',10)==0 {
+                	if mod(`replication',100)==0 {
                 	di "replication = `replication'"
                   }
                   qui set obs 7
@@ -4294,6 +4391,9 @@ di "Scenario 20B - N=200"
               		use `grp0'
                   qui append using `grp1'
               		drop id
+                  qui gen dif1 = `difi'
+                  qui gen dif2 = `difj'
+                  qui gen dif3 = `difk'
               		qui gen id = _n
               		order(id)
               		qui gen replication = `replication'
diff --git a/Scripts/Scenarios/DIF/scenarios_DIF_J7_baseline_300.do b/Scripts/Scenarios/DIF/scenarios_DIF_J7_baseline_300.do
index 093e9bb..a17d90c 100644
--- a/Scripts/Scenarios/DIF/scenarios_DIF_J7_baseline_300.do
+++ b/Scripts/Scenarios/DIF/scenarios_DIF_J7_baseline_300.do
@@ -39,7 +39,7 @@ di "Scenario 13A - N=300"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 7
@@ -114,6 +114,8 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -135,7 +137,7 @@ forvalues replication = 1/1000 {
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 7
@@ -208,6 +210,8 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -228,7 +232,7 @@ forvalues replication = 1/1000 {
 
     forvalues replication = 1/1000 {
       clear
-      if mod(`replication',10)==0 {
+      if mod(`replication',100)==0 {
       di "replication = `replication'"
       }
       qui set obs 7
@@ -301,6 +305,8 @@ forvalues replication = 1/1000 {
     		use `grp0'
         qui append using `grp1'
     		drop id
+        qui gen dif1 = `difi'
+        qui gen dif2 = `difj'
     		qui gen id = _n
     		order(id)
     		qui gen replication = `replication'
@@ -322,7 +328,7 @@ forvalues replication = 1/1000 {
 
         forvalues replication = 1/1000 {
           clear
-          if mod(`replication',10)==0 {
+          if mod(`replication',100)==0 {
           di "replication = `replication'"
           }
           qui set obs 7
@@ -394,6 +400,8 @@ forvalues replication = 1/1000 {
         		use `grp0'
             qui append using `grp1'
         		drop id
+            qui gen dif1 = `difi'
+            qui gen dif2 = `difj'
         		qui gen id = _n
         		order(id)
         		qui gen replication = `replication'
@@ -414,7 +422,7 @@ forvalues replication = 1/1000 {
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 7
@@ -486,6 +494,8 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -511,7 +521,7 @@ di "Scenario 14A - N=300"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 7
@@ -583,6 +593,8 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -604,7 +616,7 @@ forvalues replication = 1/1000 {
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 7
@@ -676,6 +688,8 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -696,7 +710,7 @@ forvalues replication = 1/1000 {
 
     forvalues replication = 1/1000 {
       clear
-      if mod(`replication',10)==0 {
+      if mod(`replication',100)==0 {
       di "replication = `replication'"
       }
       qui set obs 7
@@ -768,6 +782,8 @@ forvalues replication = 1/1000 {
     		use `grp0'
         qui append using `grp1'
     		drop id
+        qui gen dif1 = `difi'
+        qui gen dif2 = `difj'
     		qui gen id = _n
     		order(id)
     		qui gen replication = `replication'
@@ -789,7 +805,7 @@ forvalues replication = 1/1000 {
 
         forvalues replication = 1/1000 {
           clear
-          if mod(`replication',10)==0 {
+          if mod(`replication',100)==0 {
           di "replication = `replication'"
           }
           qui set obs 7
@@ -861,6 +877,8 @@ forvalues replication = 1/1000 {
         		use `grp0'
             qui append using `grp1'
         		drop id
+            qui gen dif1 = `difi'
+            qui gen dif2 = `difj'
         		qui gen id = _n
         		order(id)
         		qui gen replication = `replication'
@@ -881,7 +899,7 @@ forvalues replication = 1/1000 {
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 7
@@ -953,6 +971,8 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -978,7 +998,7 @@ di "Scenario 15A - N=300"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 7
@@ -1049,6 +1069,8 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -1070,7 +1092,7 @@ forvalues replication = 1/1000 {
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 7
@@ -1142,6 +1164,8 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -1162,7 +1186,7 @@ forvalues replication = 1/1000 {
 
     forvalues replication = 1/1000 {
       clear
-      if mod(`replication',10)==0 {
+      if mod(`replication',100)==0 {
       di "replication = `replication'"
       }
       qui set obs 7
@@ -1234,6 +1258,8 @@ forvalues replication = 1/1000 {
     		use `grp0'
         qui append using `grp1'
     		drop id
+        qui gen dif1 = `difi'
+        qui gen dif2 = `difj'
     		qui gen id = _n
     		order(id)
     		qui gen replication = `replication'
@@ -1255,7 +1281,7 @@ forvalues replication = 1/1000 {
 
         forvalues replication = 1/1000 {
           clear
-          if mod(`replication',10)==0 {
+          if mod(`replication',100)==0 {
           di "replication = `replication'"
           }
           qui set obs 7
@@ -1327,6 +1353,8 @@ forvalues replication = 1/1000 {
         		use `grp0'
             qui append using `grp1'
         		drop id
+            qui gen dif1 = `difi'
+            qui gen dif2 = `difj'
         		qui gen id = _n
         		order(id)
         		qui gen replication = `replication'
@@ -1347,7 +1375,7 @@ forvalues replication = 1/1000 {
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 7
@@ -1419,6 +1447,8 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -1444,7 +1474,7 @@ di "Scenario 16A - N=300"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 7
@@ -1516,6 +1546,8 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -1537,7 +1569,7 @@ forvalues replication = 1/1000 {
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 7
@@ -1609,6 +1641,8 @@ forvalues replication = 1/1000 {
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -1629,7 +1663,7 @@ forvalues replication = 1/1000 {
 
     forvalues replication = 1/1000 {
       clear
-      if mod(`replication',10)==0 {
+      if mod(`replication',100)==0 {
       di "replication = `replication'"
       }
       qui set obs 7
@@ -1701,6 +1735,8 @@ forvalues replication = 1/1000 {
     		use `grp0'
         qui append using `grp1'
     		drop id
+        qui gen dif1 = `difi'
+        qui gen dif2 = `difj'
     		qui gen id = _n
     		order(id)
     		qui gen replication = `replication'
@@ -1722,7 +1758,7 @@ forvalues replication = 1/1000 {
 
         forvalues replication = 1/1000 {
           clear
-          if mod(`replication',10)==0 {
+          if mod(`replication',100)==0 {
           di "replication = `replication'"
           }
           qui set obs 7
@@ -1794,6 +1830,8 @@ forvalues replication = 1/1000 {
         		use `grp0'
             qui append using `grp1'
         		drop id
+            qui gen dif1 = `difi'
+            qui gen dif2 = `difj'
         		qui gen id = _n
         		order(id)
         		qui gen replication = `replication'
@@ -1814,7 +1852,7 @@ forvalues replication = 1/1000 {
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 7
@@ -1886,6 +1924,8 @@ forvalues replication = 1/1000 {
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -1915,7 +1955,7 @@ di "Scenario 17A - N=300"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 7
@@ -2011,6 +2051,9 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
+    qui gen dif3 = `difk'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -2032,7 +2075,7 @@ di "Scenario 17B - N=300"
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 7
@@ -2128,6 +2171,9 @@ di "Scenario 17B - N=300"
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
+      qui gen dif3 = `difk'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -2158,7 +2204,7 @@ di "Scenario 17B - N=300"
 
       forvalues replication = 1/1000 {
         clear
-        if mod(`replication',10)==0 {
+        if mod(`replication',100)==0 {
         di "replication = `replication'"
         }
         qui set obs 7
@@ -2253,6 +2299,9 @@ di "Scenario 17B - N=300"
       		use `grp0'
           qui append using `grp1'
       		drop id
+          qui gen dif1 = `difi'
+          qui gen dif2 = `difj'
+          qui gen dif3 = `difk'
       		qui gen id = _n
       		order(id)
       		qui gen replication = `replication'
@@ -2279,7 +2328,7 @@ di "Scenario 17B - N=300"
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 7
@@ -2375,6 +2424,9 @@ di "Scenario 17B - N=300"
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
+              qui gen dif3 = `difk'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -2405,7 +2457,7 @@ di "Scenario 17B - N=300"
 
               forvalues replication = 1/1000 {
                   clear
-                	if mod(`replication',10)==0 {
+                	if mod(`replication',100)==0 {
                 	di "replication = `replication'"
                   }
                   qui set obs 7
@@ -2498,6 +2550,9 @@ di "Scenario 17B - N=300"
               		use `grp0'
                   qui append using `grp1'
               		drop id
+                  qui gen dif1 = `difi'
+                  qui gen dif2 = `difj'
+                  qui gen dif3 = `difk'
               		qui gen id = _n
               		order(id)
               		qui gen replication = `replication'
@@ -2525,7 +2580,7 @@ di "Scenario 18A - N=300"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 7
@@ -2621,6 +2676,9 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
+    qui gen dif3 = `difk'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -2642,7 +2700,7 @@ di "Scenario 18B - N=300"
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 7
@@ -2738,6 +2796,9 @@ di "Scenario 18B - N=300"
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
+      qui gen dif3 = `difk'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -2768,7 +2829,7 @@ di "Scenario 18B - N=300"
 
       forvalues replication = 1/1000 {
         clear
-        if mod(`replication',10)==0 {
+        if mod(`replication',100)==0 {
         di "replication = `replication'"
         }
         qui set obs 7
@@ -2863,6 +2924,9 @@ di "Scenario 18B - N=300"
       		use `grp0'
           qui append using `grp1'
       		drop id
+          qui gen dif1 = `difi'
+          qui gen dif2 = `difj'
+          qui gen dif3 = `difk'
       		qui gen id = _n
       		order(id)
       		qui gen replication = `replication'
@@ -2889,7 +2953,7 @@ di "Scenario 18B - N=300"
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 7
@@ -2985,6 +3049,9 @@ di "Scenario 18B - N=300"
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
+              qui gen dif3 = `difk'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -3015,7 +3082,7 @@ di "Scenario 18B - N=300"
 
               forvalues replication = 1/1000 {
                   clear
-                	if mod(`replication',10)==0 {
+                	if mod(`replication',100)==0 {
                 	di "replication = `replication'"
                   }
                   qui set obs 7
@@ -3108,6 +3175,9 @@ di "Scenario 18B - N=300"
               		use `grp0'
                   qui append using `grp1'
               		drop id
+                  qui gen dif1 = `difi'
+                  qui gen dif2 = `difj'
+                  qui gen dif3 = `difk'
               		qui gen id = _n
               		order(id)
               		qui gen replication = `replication'
@@ -3143,7 +3213,7 @@ di "Scenario 19A - N=300"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 7
@@ -3239,6 +3309,9 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
+    qui gen dif3 = `difk'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -3260,7 +3333,7 @@ di "Scenario 19B - N=300"
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 7
@@ -3356,6 +3429,9 @@ di "Scenario 19B - N=300"
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
+      qui gen dif3 = `difk'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -3378,7 +3454,7 @@ di "Scenario 19B - N=300"
 
       forvalues replication = 1/1000 {
         clear
-        if mod(`replication',10)==0 {
+        if mod(`replication',100)==0 {
         di "replication = `replication'"
         }
         qui set obs 7
@@ -3473,6 +3549,9 @@ di "Scenario 19B - N=300"
       		use `grp0'
           qui append using `grp1'
       		drop id
+          qui gen dif1 = `difi'
+          qui gen dif2 = `difj'
+          qui gen dif3 = `difk'
       		qui gen id = _n
       		order(id)
       		qui gen replication = `replication'
@@ -3495,7 +3574,7 @@ di "Scenario 19B - N=300"
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 7
@@ -3591,6 +3670,9 @@ di "Scenario 19B - N=300"
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
+              qui gen dif3 = `difk'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -3613,7 +3695,7 @@ di "Scenario 19B - N=300"
 
               forvalues replication = 1/1000 {
                   clear
-                	if mod(`replication',10)==0 {
+                	if mod(`replication',100)==0 {
                 	di "replication = `replication'"
                   }
                   qui set obs 7
@@ -3706,6 +3788,9 @@ di "Scenario 19B - N=300"
               		use `grp0'
                   qui append using `grp1'
               		drop id
+                  qui gen dif1 = `difi'
+                  qui gen dif2 = `difj'
+                  qui gen dif3 = `difk'
               		qui gen id = _n
               		order(id)
               		qui gen replication = `replication'
@@ -3732,7 +3817,7 @@ di "Scenario 20A - N=300"
 
 forvalues replication = 1/1000 {
   clear
-  if mod(`replication',10)==0 {
+  if mod(`replication',100)==0 {
   di "replication = `replication'"
   }
   qui set obs 7
@@ -3828,6 +3913,9 @@ forvalues replication = 1/1000 {
 		use `grp0'
     qui append using `grp1'
 		drop id
+    qui gen dif1 = `difi'
+    qui gen dif2 = `difj'
+    qui gen dif3 = `difk'
 		qui gen id = _n
 		order(id)
 		qui gen replication = `replication'
@@ -3849,7 +3937,7 @@ di "Scenario 20B - N=300"
 
   forvalues replication = 1/1000 {
     clear
-    if mod(`replication',10)==0 {
+    if mod(`replication',100)==0 {
     di "replication = `replication'"
     }
     qui set obs 7
@@ -3945,6 +4033,9 @@ di "Scenario 20B - N=300"
   		use `grp0'
       qui append using `grp1'
   		drop id
+      qui gen dif1 = `difi'
+      qui gen dif2 = `difj'
+      qui gen dif3 = `difk'
   		qui gen id = _n
   		order(id)
   		qui gen replication = `replication'
@@ -3966,7 +4057,7 @@ di "Scenario 20B - N=300"
 
       forvalues replication = 1/1000 {
         clear
-        if mod(`replication',10)==0 {
+        if mod(`replication',100)==0 {
         di "replication = `replication'"
         }
         qui set obs 7
@@ -4061,6 +4152,9 @@ di "Scenario 20B - N=300"
       		use `grp0'
           qui append using `grp1'
       		drop id
+          qui gen dif1 = `difi'
+          qui gen dif2 = `difj'
+          qui gen dif3 = `difk'
       		qui gen id = _n
       		order(id)
       		qui gen replication = `replication'
@@ -4083,7 +4177,7 @@ di "Scenario 20B - N=300"
 
           forvalues replication = 1/1000 {
             clear
-            if mod(`replication',10)==0 {
+            if mod(`replication',100)==0 {
             di "replication = `replication'"
             }
             qui set obs 7
@@ -4179,6 +4273,9 @@ di "Scenario 20B - N=300"
           		use `grp0'
               qui append using `grp1'
           		drop id
+              qui gen dif1 = `difi'
+              qui gen dif2 = `difj'
+              qui gen dif3 = `difk'
           		qui gen id = _n
           		order(id)
           		qui gen replication = `replication'
@@ -4201,7 +4298,7 @@ di "Scenario 20B - N=300"
 
               forvalues replication = 1/1000 {
                   clear
-                	if mod(`replication',10)==0 {
+                	if mod(`replication',100)==0 {
                 	di "replication = `replication'"
                   }
                   qui set obs 7
@@ -4294,6 +4391,9 @@ di "Scenario 20B - N=300"
               		use `grp0'
                   qui append using `grp1'
               		drop id
+                  qui gen dif1 = `difi'
+                  qui gen dif2 = `difj'
+                  qui gen dif3 = `difk'
               		qui gen id = _n
               		order(id)
               		qui gen replication = `replication'
diff --git a/catalogue.md b/catalogue.md
index 22e7cd4..8916761 100644
--- a/catalogue.md
+++ b/catalogue.md
@@ -26,14 +26,14 @@ Chaque scénario est séparé en 3 sous-scénarios:
   | **Scenario 10** | X |   |   | X |  | X | | X | |
   | **Scenario 11** | X |  | X |   | | X | |  | X |
   | **Scenario 12** | X |  |   | X | | X | |  | X |
-  | **Scenario 13** | X |   | X |   |
-  | **Scenario 14** | X |   |   | X |
-  | **Scenario 15** |   | X | X |   |
-  | **Scenario 16** |   | X |   | X |
-  | **Scenario 17** | X |   | X |   |
-  | **Scenario 18** | X |   |   | X |
-  | **Scenario 19** |   | X | X |   |
-  | **Scenario 20** |   | X |   | X |
+  | **Scenario 13** |  | X | X |   |  | X | | X | |
+  | **Scenario 14** |  | X |   | X |  | X | | X | |
+  | **Scenario 15** |   | X | X |   | | X | |  | X |
+  | **Scenario 16** |   | X |   | X | | X | |  | X |
+  | **Scenario 17** |   | X | X |   |  |  | X | X | |
+  | **Scenario 18** |   | X |   | X |  |  | X | X | |
+  | **Scenario 19** |   | X | X |   |  |  | X |  | X |
+  | **Scenario 20** |   | X |   | X |  |  | X | | X |
 
   Chaque scénario est séparé en 3 sous-scénarios:
     * A - H<sub>0</sub> vraie