From 01ebe3b50761991fc73bbb7b23b611cfb38d633c Mon Sep 17 00:00:00 2001
From: corentinchoisy <corentin.choisy@proton.me>
Date: Tue, 9 Jan 2024 14:06:51 +0100
Subject: [PATCH] Added scenarios 15 abd 16

---
 .../DIF/scenarios_DIF_baseline_100.do         | 2985 +++++++++++++++--
 1 file changed, 2733 insertions(+), 252 deletions(-)

diff --git a/Scripts/Scenarios/DIF/scenarios_DIF_baseline_100.do b/Scripts/Scenarios/DIF/scenarios_DIF_baseline_100.do
index 0e69f03..ca27769 100644
--- a/Scripts/Scenarios/DIF/scenarios_DIF_baseline_100.do
+++ b/Scripts/Scenarios/DIF/scenarios_DIF_baseline_100.do
@@ -5100,149 +5100,149 @@ forvalues replication = 1/1000 {
     if `difi'==1 {
       qui local difj = runiformint(2,7)
       if `difj'==2 {
-        mat D= (-0.85 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+        mat D= (-1.85,-0.85,0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
       }
       if `difj'==3 {
-        mat D= (-0.85 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.15)
+        mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
       }
       if `difj'==4 {
-        mat D= (-0.85 \ -0.67 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+        mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
       }
       if `difj'==5 {
-        mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.15)
+        mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
       }
       if `difj'==6 {
-        mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.15)
+        mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
       }
       if `difj'==7 {
-        mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.45)
+        mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
       }
     }
     if `difi'==2 {
       qui local difj = runiformint(2,7)
       if `difj'==2 {
-        mat D= (-0.85 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+        mat D= (-1.85,-0.85,0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
       }
       if `difj'==3 {
-        mat D= (-1.15 \ -0.37 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.15)
+        mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
       }
       if `difj'==4 {
-        mat D= (-1.15 \ -0.37 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+        mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
       }
       if `difj'==5 {
-        mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.15)
+        mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
       }
       if `difj'==6 {
-        mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.15)
+        mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
       }
       if `difj'==7 {
-        mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.45)
+        mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
       }
     }
 
     if `difi'==3 {
       qui local difj = runiformint(2,7)
       if `difj'==2 {
-        mat D= (-0.85 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.15)
+        mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
       }
       if `difj'==3 {
-        mat D= (-1.15 \ -0.37 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.15)
+        mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
       }
       if `difj'==4 {
-        mat D= (-1.15 \ -0.67 \ -0.02 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
       }
       if `difj'==5 {
-        mat D= (-1.15 \ -0.67 \-0.02 \ 0 \ 0.62 \ 0.67 \ 1.15)
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
       }
       if `difj'==6 {
-        mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.97 \ 1.15)
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
       }
       if `difj'==7 {
-        mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.45)
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
       }
     }
     if `difi'==4 {
       qui local difj = runiformint(2,7)
       if `difj'==2 {
-        mat D= (-0.85 \ -0.67 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+        mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
       }
       if `difj'==3 {
-        mat D= (-1.15 \ -0.37 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+        mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
       }
       if `difj'==4 {
-        mat D= (-1.15 \ -0.67 \ -0.02 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
       }
       if `difj'==5 {
-        mat D= (-1.15 \ -0.67 \-0.32 \ 0.3 \ 0.62 \ 0.67 \ 1.15)
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
       }
       if `difj'==6 {
-        mat D= (-1.15 \ -0.67 \ -0.32 \ 0.3 \ 0.32 \ 0.97 \ 1.15)
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
       }
       if `difj'==7 {
-        mat D= (-1.15 \ -0.67 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.45)
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
       }
     }
     if `difi'==5 {
       qui local difj = runiformint(2,7)
       if `difj'==2 {
-        mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.15)
+        mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
       }
       if `difj'==3 {
-        mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.15)
+        mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
       }
       if `difj'==4 {
-        mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.62 \ 0.67 \ 1.15)
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
       }
       if `difj'==5 {
-        mat D= (-1.15 \ -0.67 \-0.32 \ 0.3 \ 0.62 \ 0.67 \ 1.15)
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
       }
       if `difj'==6 {
-        mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.97 \ 1.15)
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
       }
       if `difj'==7 {
-        mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.45)
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
       }
     }
     if `difi'==6 {
       qui local difj = runiformint(2,7)
       if `difj'==2 {
-        mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.15)
+        mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
       }
       if `difj'==3 {
-        mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.15)
+        mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
       }
       if `difj'==4 {
-        mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.97 \ 1.15)
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
       }
       if `difj'==5 {
-        mat D= (-1.15 \ -0.67 \-0.32 \ 0.3 \ 0.32 \ 0.97 \ 1.15)
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
       }
       if `difj'==6 {
-        mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.97 \ 1.15)
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
       }
       if `difj'==7 {
-        mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.45)
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.45,1.45,2.45)
       }
     }
     if `difi'==7 {
       qui local difj = runiformint(2,7)
       if `difj'==2 {
-        mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.45)
+        mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
       }
       if `difj'==3 {
-        mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.45)
+        mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
       }
       if `difj'==4 {
-        mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.45)
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
       }
       if `difj'==5 {
-        mat D= (-1.15 \ -0.67 \-0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.45)
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
       }
       if `difj'==6 {
-        mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.45)
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
       }
       if `difj'==7 {
-        mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.45)
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.45,1.45,2.45)
       }
     }
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
@@ -5290,149 +5290,149 @@ forvalues replication = 1/1000 {
       if `difi'==1 {
         qui local difj = runiformint(2,7)
         if `difj'==2 {
-          mat D= (-0.85 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+          mat D= (-1.85,-0.85,0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
         }
         if `difj'==3 {
-          mat D= (-0.85 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.15)
+          mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
         }
         if `difj'==4 {
-          mat D= (-0.85 \ -0.67 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+          mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
         }
         if `difj'==5 {
-          mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.15)
+          mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
         }
         if `difj'==6 {
-          mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.15)
+          mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
         }
         if `difj'==7 {
-          mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.45)
+          mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
         }
       }
       if `difi'==2 {
         qui local difj = runiformint(2,7)
         if `difj'==2 {
-          mat D= (-0.85 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+          mat D= (-1.85,-0.85,0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
         }
         if `difj'==3 {
-          mat D= (-1.15 \ -0.37 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.15)
+          mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
         }
         if `difj'==4 {
-          mat D= (-1.15 \ -0.37 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+          mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
         }
         if `difj'==5 {
-          mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.15)
+          mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
         }
         if `difj'==6 {
-          mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.15)
+          mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
         }
         if `difj'==7 {
-          mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.45)
+          mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
         }
       }
 
       if `difi'==3 {
         qui local difj = runiformint(2,7)
         if `difj'==2 {
-          mat D= (-0.85 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.15)
+          mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
         }
         if `difj'==3 {
-          mat D= (-1.15 \ -0.37 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.15)
+          mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
         }
         if `difj'==4 {
-          mat D= (-1.15 \ -0.67 \ -0.02 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
         }
         if `difj'==5 {
-          mat D= (-1.15 \ -0.67 \-0.02 \ 0 \ 0.62 \ 0.67 \ 1.15)
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
         }
         if `difj'==6 {
-          mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.97 \ 1.15)
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
         }
         if `difj'==7 {
-          mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.45)
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
         }
       }
       if `difi'==4 {
         qui local difj = runiformint(2,7)
         if `difj'==2 {
-          mat D= (-0.85 \ -0.67 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+          mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
         }
         if `difj'==3 {
-          mat D= (-1.15 \ -0.37 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+          mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
         }
         if `difj'==4 {
-          mat D= (-1.15 \ -0.67 \ -0.02 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
         }
         if `difj'==5 {
-          mat D= (-1.15 \ -0.67 \-0.32 \ 0.3 \ 0.62 \ 0.67 \ 1.15)
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
         }
         if `difj'==6 {
-          mat D= (-1.15 \ -0.67 \ -0.32 \ 0.3 \ 0.32 \ 0.97 \ 1.15)
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
         }
         if `difj'==7 {
-          mat D= (-1.15 \ -0.67 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.45)
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
         }
       }
       if `difi'==5 {
         qui local difj = runiformint(2,7)
         if `difj'==2 {
-          mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.15)
+          mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
         }
         if `difj'==3 {
-          mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.15)
+          mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
         }
         if `difj'==4 {
-          mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.62 \ 0.67 \ 1.15)
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
         }
         if `difj'==5 {
-          mat D= (-1.15 \ -0.67 \-0.32 \ 0.3 \ 0.62 \ 0.67 \ 1.15)
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
         }
         if `difj'==6 {
-          mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.97 \ 1.15)
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
         }
         if `difj'==7 {
-          mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.45)
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
         }
       }
       if `difi'==6 {
         qui local difj = runiformint(2,7)
         if `difj'==2 {
-          mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.15)
+          mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
         }
         if `difj'==3 {
-          mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.15)
+          mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
         }
         if `difj'==4 {
-          mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.97 \ 1.15)
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
         }
         if `difj'==5 {
-          mat D= (-1.15 \ -0.67 \-0.32 \ 0.3 \ 0.32 \ 0.97 \ 1.15)
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
         }
         if `difj'==6 {
-          mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.97 \ 1.15)
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
         }
         if `difj'==7 {
-          mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.45)
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.45,1.45,2.45)
         }
       }
       if `difi'==7 {
         qui local difj = runiformint(2,7)
         if `difj'==2 {
-          mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.45)
+          mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
         }
         if `difj'==3 {
-          mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.45)
+          mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
         }
         if `difj'==4 {
-          mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.45)
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
         }
         if `difj'==5 {
-          mat D= (-1.15 \ -0.67 \-0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.45)
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
         }
         if `difj'==6 {
-          mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.45)
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
         }
         if `difj'==7 {
-          mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.45)
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.45,1.45,2.45)
         }
       }
   		qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear
@@ -5472,149 +5472,149 @@ forvalues replication = 1/1000 {
         if `difi'==1 {
           qui local difj = runiformint(2,7)
           if `difj'==2 {
-            mat D= (-0.85 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+            mat D= (-1.85,-0.85,0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
           }
           if `difj'==3 {
-            mat D= (-0.85 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.15)
+            mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
           }
           if `difj'==4 {
-            mat D= (-0.85 \ -0.67 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+            mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
           }
           if `difj'==5 {
-            mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.15)
+            mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
           }
           if `difj'==6 {
-            mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.15)
+            mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
           }
           if `difj'==7 {
-            mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.45)
+            mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
           }
         }
         if `difi'==2 {
           qui local difj = runiformint(2,7)
           if `difj'==2 {
-            mat D= (-0.85 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+            mat D= (-1.85,-0.85,0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
           }
           if `difj'==3 {
-            mat D= (-1.15 \ -0.37 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.15)
+            mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
           }
           if `difj'==4 {
-            mat D= (-1.15 \ -0.37 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+            mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
           }
           if `difj'==5 {
-            mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.15)
+            mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
           }
           if `difj'==6 {
-            mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.15)
+            mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
           }
           if `difj'==7 {
-            mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.45)
+            mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
           }
         }
 
         if `difi'==3 {
           qui local difj = runiformint(2,7)
           if `difj'==2 {
-            mat D= (-0.85 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.15)
+            mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
           }
           if `difj'==3 {
-            mat D= (-1.15 \ -0.37 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.15)
+            mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
           }
           if `difj'==4 {
-            mat D= (-1.15 \ -0.67 \ -0.02 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
           }
           if `difj'==5 {
-            mat D= (-1.15 \ -0.67 \-0.02 \ 0 \ 0.62 \ 0.67 \ 1.15)
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
           }
           if `difj'==6 {
-            mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.97 \ 1.15)
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
           }
           if `difj'==7 {
-            mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.45)
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
           }
         }
         if `difi'==4 {
           qui local difj = runiformint(2,7)
           if `difj'==2 {
-            mat D= (-0.85 \ -0.67 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+            mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
           }
           if `difj'==3 {
-            mat D= (-1.15 \ -0.37 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+            mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
           }
           if `difj'==4 {
-            mat D= (-1.15 \ -0.67 \ -0.02 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
           }
           if `difj'==5 {
-            mat D= (-1.15 \ -0.67 \-0.32 \ 0.3 \ 0.62 \ 0.67 \ 1.15)
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
           }
           if `difj'==6 {
-            mat D= (-1.15 \ -0.67 \ -0.32 \ 0.3 \ 0.32 \ 0.97 \ 1.15)
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
           }
           if `difj'==7 {
-            mat D= (-1.15 \ -0.67 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.45)
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
           }
         }
         if `difi'==5 {
           qui local difj = runiformint(2,7)
           if `difj'==2 {
-            mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.15)
+            mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
           }
           if `difj'==3 {
-            mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.15)
+            mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
           }
           if `difj'==4 {
-            mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.62 \ 0.67 \ 1.15)
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
           }
           if `difj'==5 {
-            mat D= (-1.15 \ -0.67 \-0.32 \ 0.3 \ 0.62 \ 0.67 \ 1.15)
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
           }
           if `difj'==6 {
-            mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.97 \ 1.15)
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
           }
           if `difj'==7 {
-            mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.45)
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
           }
         }
         if `difi'==6 {
           qui local difj = runiformint(2,7)
           if `difj'==2 {
-            mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.15)
+            mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
           }
           if `difj'==3 {
-            mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.15)
+            mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
           }
           if `difj'==4 {
-            mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.97 \ 1.15)
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
           }
           if `difj'==5 {
-            mat D= (-1.15 \ -0.67 \-0.32 \ 0.3 \ 0.32 \ 0.97 \ 1.15)
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
           }
           if `difj'==6 {
-            mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.97 \ 1.15)
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
           }
           if `difj'==7 {
-            mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.45)
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.45,1.45,2.45)
           }
         }
         if `difi'==7 {
           qui local difj = runiformint(2,7)
           if `difj'==2 {
-            mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.45)
+            mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
           }
           if `difj'==3 {
-            mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.45)
+            mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
           }
           if `difj'==4 {
-            mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.45)
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
           }
           if `difj'==5 {
-            mat D= (-1.15 \ -0.67 \-0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.45)
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
           }
           if `difj'==6 {
-            mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.45)
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
           }
           if `difj'==7 {
-            mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.45)
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.45,1.45,2.45)
           }
         }
     		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
@@ -5669,149 +5669,149 @@ forvalues replication = 1/1000 {
             if `difi'==1 {
               qui local difj = runiformint(2,7)
               if `difj'==2 {
-                mat D= (-0.85 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+                mat D= (-1.85,-0.85,0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
               }
               if `difj'==3 {
-                mat D= (-0.85 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.15)
+                mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
               }
               if `difj'==4 {
-                mat D= (-0.85 \ -0.67 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+                mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
               }
               if `difj'==5 {
-                mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.15)
+                mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
               }
               if `difj'==6 {
-                mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.15)
+                mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
               }
               if `difj'==7 {
-                mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.45)
+                mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
               }
             }
             if `difi'==2 {
               qui local difj = runiformint(2,7)
               if `difj'==2 {
-                mat D= (-0.85 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+                mat D= (-1.85,-0.85,0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
               }
               if `difj'==3 {
-                mat D= (-1.15 \ -0.37 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.15)
+                mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
               }
               if `difj'==4 {
-                mat D= (-1.15 \ -0.37 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+                mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
               }
               if `difj'==5 {
-                mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.15)
+                mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
               }
               if `difj'==6 {
-                mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.15)
+                mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
               }
               if `difj'==7 {
-                mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.45)
+                mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
               }
             }
 
             if `difi'==3 {
               qui local difj = runiformint(2,7)
               if `difj'==2 {
-                mat D= (-0.85 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.15)
+                mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
               }
               if `difj'==3 {
-                mat D= (-1.15 \ -0.37 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.15)
+                mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
               }
               if `difj'==4 {
-                mat D= (-1.15 \ -0.67 \ -0.02 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
               }
               if `difj'==5 {
-                mat D= (-1.15 \ -0.67 \-0.02 \ 0 \ 0.62 \ 0.67 \ 1.15)
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
               }
               if `difj'==6 {
-                mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.97 \ 1.15)
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
               }
               if `difj'==7 {
-                mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.45)
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
               }
             }
             if `difi'==4 {
               qui local difj = runiformint(2,7)
               if `difj'==2 {
-                mat D= (-0.85 \ -0.67 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+                mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
               }
               if `difj'==3 {
-                mat D= (-1.15 \ -0.37 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+                mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
               }
               if `difj'==4 {
-                mat D= (-1.15 \ -0.67 \ -0.02 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
               }
               if `difj'==5 {
-                mat D= (-1.15 \ -0.67 \-0.32 \ 0.3 \ 0.62 \ 0.67 \ 1.15)
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
               }
               if `difj'==6 {
-                mat D= (-1.15 \ -0.67 \ -0.32 \ 0.3 \ 0.32 \ 0.97 \ 1.15)
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
               }
               if `difj'==7 {
-                mat D= (-1.15 \ -0.67 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.45)
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
               }
             }
             if `difi'==5 {
               qui local difj = runiformint(2,7)
               if `difj'==2 {
-                mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.15)
+                mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
               }
               if `difj'==3 {
-                mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.15)
+                mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
               }
               if `difj'==4 {
-                mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.62 \ 0.67 \ 1.15)
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
               }
               if `difj'==5 {
-                mat D= (-1.15 \ -0.67 \-0.32 \ 0.3 \ 0.62 \ 0.67 \ 1.15)
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
               }
               if `difj'==6 {
-                mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.97 \ 1.15)
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
               }
               if `difj'==7 {
-                mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.45)
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
               }
             }
             if `difi'==6 {
               qui local difj = runiformint(2,7)
               if `difj'==2 {
-                mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.15)
+                mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
               }
               if `difj'==3 {
-                mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.15)
+                mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
               }
               if `difj'==4 {
-                mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.97 \ 1.15)
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
               }
               if `difj'==5 {
-                mat D= (-1.15 \ -0.67 \-0.32 \ 0.3 \ 0.32 \ 0.97 \ 1.15)
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
               }
               if `difj'==6 {
-                mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.97 \ 1.15)
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
               }
               if `difj'==7 {
-                mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.45)
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.45,1.45,2.45)
               }
             }
             if `difi'==7 {
               qui local difj = runiformint(2,7)
               if `difj'==2 {
-                mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.45)
+                mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
               }
               if `difj'==3 {
-                mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.45)
+                mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
               }
               if `difj'==4 {
-                mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.45)
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
               }
               if `difj'==5 {
-                mat D= (-1.15 \ -0.67 \-0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.45)
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
               }
               if `difj'==6 {
-                mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.45)
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
               }
               if `difj'==7 {
-                mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.45)
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.45,1.45,2.45)
               }
             }
         		qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear
@@ -5851,149 +5851,149 @@ forvalues replication = 1/1000 {
               if `difi'==1 {
                 qui local difj = runiformint(2,7)
                 if `difj'==2 {
-                  mat D= (-0.85 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+                  mat D= (-1.85,-0.85,0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
                 }
                 if `difj'==3 {
-                  mat D= (-0.85 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.15)
+                  mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
                 }
                 if `difj'==4 {
-                  mat D= (-0.85 \ -0.67 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+                  mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
                 }
                 if `difj'==5 {
-                  mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.15)
+                  mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
                 }
                 if `difj'==6 {
-                  mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.15)
+                  mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
                 }
                 if `difj'==7 {
-                  mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.45)
+                  mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
                 }
               }
               if `difi'==2 {
                 qui local difj = runiformint(2,7)
                 if `difj'==2 {
-                  mat D= (-0.85 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+                  mat D= (-1.85,-0.85,0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
                 }
                 if `difj'==3 {
-                  mat D= (-1.15 \ -0.37 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.15)
+                  mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
                 }
                 if `difj'==4 {
-                  mat D= (-1.15 \ -0.37 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+                  mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
                 }
                 if `difj'==5 {
-                  mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.15)
+                  mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
                 }
                 if `difj'==6 {
-                  mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.15)
+                  mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
                 }
                 if `difj'==7 {
-                  mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.45)
+                  mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
                 }
               }
 
               if `difi'==3 {
                 qui local difj = runiformint(2,7)
                 if `difj'==2 {
-                  mat D= (-0.85 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.15)
+                  mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
                 }
                 if `difj'==3 {
-                  mat D= (-1.15 \ -0.37 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.15)
+                  mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
                 }
                 if `difj'==4 {
-                  mat D= (-1.15 \ -0.67 \ -0.02 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
                 }
                 if `difj'==5 {
-                  mat D= (-1.15 \ -0.67 \-0.02 \ 0 \ 0.62 \ 0.67 \ 1.15)
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
                 }
                 if `difj'==6 {
-                  mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.97 \ 1.15)
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
                 }
                 if `difj'==7 {
-                  mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.45)
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
                 }
               }
               if `difi'==4 {
                 qui local difj = runiformint(2,7)
                 if `difj'==2 {
-                  mat D= (-0.85 \ -0.67 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+                  mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
                 }
                 if `difj'==3 {
-                  mat D= (-1.15 \ -0.37 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+                  mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
                 }
                 if `difj'==4 {
-                  mat D= (-1.15 \ -0.67 \ -0.02 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
                 }
                 if `difj'==5 {
-                  mat D= (-1.15 \ -0.67 \-0.32 \ 0.3 \ 0.62 \ 0.67 \ 1.15)
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
                 }
                 if `difj'==6 {
-                  mat D= (-1.15 \ -0.67 \ -0.32 \ 0.3 \ 0.32 \ 0.97 \ 1.15)
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
                 }
                 if `difj'==7 {
-                  mat D= (-1.15 \ -0.67 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.45)
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
                 }
               }
               if `difi'==5 {
                 qui local difj = runiformint(2,7)
                 if `difj'==2 {
-                  mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.15)
+                  mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
                 }
                 if `difj'==3 {
-                  mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.15)
+                  mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
                 }
                 if `difj'==4 {
-                  mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.62 \ 0.67 \ 1.15)
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
                 }
                 if `difj'==5 {
-                  mat D= (-1.15 \ -0.67 \-0.32 \ 0.3 \ 0.62 \ 0.67 \ 1.15)
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
                 }
                 if `difj'==6 {
-                  mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.97 \ 1.15)
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
                 }
                 if `difj'==7 {
-                  mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.45)
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
                 }
               }
               if `difi'==6 {
                 qui local difj = runiformint(2,7)
                 if `difj'==2 {
-                  mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.15)
+                  mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
                 }
                 if `difj'==3 {
-                  mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.15)
+                  mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
                 }
                 if `difj'==4 {
-                  mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.97 \ 1.15)
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
                 }
                 if `difj'==5 {
-                  mat D= (-1.15 \ -0.67 \-0.32 \ 0.3 \ 0.32 \ 0.97 \ 1.15)
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
                 }
                 if `difj'==6 {
-                  mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.97 \ 1.15)
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
                 }
                 if `difj'==7 {
-                  mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.45)
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.45,1.45,2.45)
                 }
               }
               if `difi'==7 {
                 qui local difj = runiformint(2,7)
                 if `difj'==2 {
-                  mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.45)
+                  mat D= (-1.85,-0.85,0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
                 }
                 if `difj'==3 {
-                  mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.45)
+                  mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
                 }
                 if `difj'==4 {
-                  mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.45)
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
                 }
                 if `difj'==5 {
-                  mat D= (-1.15 \ -0.67 \-0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.45)
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
                 }
                 if `difj'==6 {
-                  mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.45)
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
                 }
                 if `difj'==7 {
-                  mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.45)
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.45,1.45,2.45)
                 }
               }
           		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
@@ -6044,25 +6044,25 @@ forvalues replication = 1/1000 {
     }
 
     if `difi'==1 {
-      mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+      mat D= (-1.85,-0.85, 0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
     }
     if `difi'==2 {
-      mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+      mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
     }
     if `difi'==3 {
-      mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.15)
+      mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
     }
     if `difi'==4 {
-      mat D= (-1.15 \ -0.67 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+      mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
     }
     if `difi'==5 {
-      mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.15)
+      mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
     }
     if `difi'==6 {
-      mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.15)
+      mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
     }
     if `difi'==7 {
-      mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.45)
+      mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
     }
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
     qui gen TT = 0
@@ -6071,25 +6071,25 @@ forvalues replication = 1/1000 {
 
 
     if `difj'==1 {
-      mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+      mat D= (-1.85,-0.85, 0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
     }
     if `difj'==2 {
-      mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+      mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
     }
     if `difj'==3 {
-      mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.15)
+      mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
     }
     if `difj'==4 {
-      mat D= (-1.15 \ -0.67 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+      mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
     }
     if `difj'==5 {
-      mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.15)
+      mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
     }
     if `difj'==6 {
-      mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.15)
+      mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
     }
     if `difj'==7 {
-      mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.45)
+      mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
     }
 		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
     qui gen TT = 1
@@ -6132,25 +6132,25 @@ forvalues replication = 1/1000 {
       }
 
       if `difi'==1 {
-        mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+        mat D= (-1.85,-0.85, 0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
       }
       if `difi'==2 {
-        mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+        mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
       }
       if `difi'==3 {
-        mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.15)
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
       }
       if `difi'==4 {
-        mat D= (-1.15 \ -0.67 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
       }
       if `difi'==5 {
-        mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.15)
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
       }
       if `difi'==6 {
-        mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.15)
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
       }
       if `difi'==7 {
-        mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.45)
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
       }
   		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
       qui gen TT = 0
@@ -6159,25 +6159,25 @@ forvalues replication = 1/1000 {
 
 
       if `difj'==1 {
-        mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+        mat D= (-1.85,-0.85, 0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
       }
       if `difj'==2 {
-        mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+        mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
       }
       if `difj'==3 {
-        mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.15)
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
       }
       if `difj'==4 {
-        mat D= (-1.15 \ -0.67 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
       }
       if `difj'==5 {
-        mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.15)
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
       }
       if `difj'==6 {
-        mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.15)
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
       }
       if `difj'==7 {
-        mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.45)
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
       }
   		qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear
       qui gen TT = 1
@@ -6219,25 +6219,25 @@ forvalues replication = 1/1000 {
           }
 
           if `difi'==1 {
-            mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+            mat D= (-1.85,-0.85, 0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
           }
           if `difi'==2 {
-            mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+            mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
           }
           if `difi'==3 {
-            mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.15)
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
           }
           if `difi'==4 {
-            mat D= (-1.15 \ -0.67 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
           }
           if `difi'==5 {
-            mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.15)
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
           }
           if `difi'==6 {
-            mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.15)
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
           }
           if `difi'==7 {
-            mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.45)
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
           }
       		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
           qui gen TT = 0
@@ -6246,25 +6246,25 @@ forvalues replication = 1/1000 {
 
 
           if `difj'==1 {
-            mat D= (-0.85 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+            mat D= (-1.85,-0.85, 0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
           }
           if `difj'==2 {
-            mat D= (-1.15 \ -0.37 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+            mat D= (-2.15,-1.15,-0.15 \ -1.37,-0.37,0.63 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
           }
           if `difj'==3 {
-            mat D= (-1.15 \ -0.67 \ -0.02 \ 0 \ 0.32 \ 0.67 \ 1.15)
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.02,-0.02,0.98 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
           }
           if `difj'==4 {
-            mat D= (-1.15 \ -0.67 \ -0.32 \ 0.3 \ 0.32 \ 0.67 \ 1.15)
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.7,0.3,1.3 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
           }
           if `difj'==5 {
-            mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.62 \ 0.67 \ 1.15)
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.38,0.62,1.62 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
           }
           if `difj'==6 {
-            mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.97 \ 1.15)
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.03,0.97,1.97 \ 0.15,1.15,2.15)
           }
           if `difj'==7 {
-            mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.45)
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.45,1.45,2.45)
           }
       		qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear
           qui gen TT = 1
@@ -6289,3 +6289,2484 @@ forvalues replication = 1/1000 {
       		}
       	}
       	export delimited using "`path'/scenario_14H_100.csv", replace
+
+
+
+
+
+**** Scenarios with : DIF size 0.5 ****
+
+** Scenario 15: J = 7 items / M = 2 modalities / DIF SIZE = 0.5
+
+* Scenario 15A : H_0 is TRUE / DIF on treatment x 2
+
+di "Scenario 15A - N=100"
+
+forvalues replication = 1/1000 {
+  	if mod(`replication',10)==0 {
+  	di "replication = `replication'"
+    }
+    qui local difi = runiformint(1,7)
+
+    mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
+    qui gen TT = 0
+		tempfile grp0
+		qui save `grp0',replace
+
+
+    if `difi'==1 {
+      qui local difj = runiformint(2,7)
+      if `difj'==2 {
+        mat D= (-0.65 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+      }
+      if `difj'==3 {
+        mat D= (-0.65 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.15)
+      }
+      if `difj'==4 {
+        mat D= (-0.65 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+      }
+      if `difj'==5 {
+        mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.15)
+      }
+      if `difj'==6 {
+        mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.15)
+      }
+      if `difj'==7 {
+        mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.65)
+      }
+    }
+    if `difi'==2 {
+      qui local difj = runiformint(2,7)
+      if `difj'==2 {
+        mat D= (-0.65 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+      }
+      if `difj'==3 {
+        mat D= (-1.15 \ -0.17 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.15)
+      }
+      if `difj'==4 {
+        mat D= (-1.15 \ -0.17 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+      }
+      if `difj'==5 {
+        mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.15)
+      }
+      if `difj'==6 {
+        mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.15)
+      }
+      if `difj'==7 {
+        mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.65)
+      }
+    }
+
+    if `difi'==3 {
+      qui local difj = runiformint(2,7)
+      if `difj'==2 {
+        mat D= (-0.65 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.15)
+      }
+      if `difj'==3 {
+        mat D= (-1.15 \ -0.17 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.15)
+      }
+      if `difj'==4 {
+        mat D= (-1.15 \ -0.67 \ 0.18 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+      }
+      if `difj'==5 {
+        mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.82 \ 0.67 \ 1.15)
+      }
+      if `difj'==6 {
+        mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 1.17 \ 1.15)
+      }
+      if `difj'==7 {
+        mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.65)
+      }
+    }
+    if `difi'==4 {
+      qui local difj = runiformint(2,7)
+      if `difj'==2 {
+        mat D= (-0.65 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+      }
+      if `difj'==3 {
+        mat D= (-1.15 \ -0.17 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+      }
+      if `difj'==4 {
+        mat D= (-1.15 \ -0.67 \ 0.18 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+      }
+      if `difj'==5 {
+        mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.82 \ 0.67 \ 1.15)
+      }
+      if `difj'==6 {
+        mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 1.17 \ 1.15)
+      }
+      if `difj'==7 {
+        mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.65)
+      }
+    }
+    if `difi'==5 {
+      qui local difj = runiformint(2,7)
+      if `difj'==2 {
+        mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.15)
+      }
+      if `difj'==3 {
+        mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.15)
+      }
+      if `difj'==4 {
+        mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.82 \ 0.67 \ 1.15)
+      }
+      if `difj'==5 {
+        mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.82 \ 0.67 \ 1.15)
+      }
+      if `difj'==6 {
+        mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 1.17 \ 1.15)
+      }
+      if `difj'==7 {
+        mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.65)
+      }
+    }
+    if `difi'==6 {
+      qui local difj = runiformint(2,7)
+      if `difj'==2 {
+        mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.15)
+      }
+      if `difj'==3 {
+        mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.15)
+      }
+      if `difj'==4 {
+        mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 1.17 \ 1.15)
+      }
+      if `difj'==5 {
+        mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 1.17 \ 1.15)
+      }
+      if `difj'==6 {
+        mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 1.17 \ 1.15)
+      }
+      if `difj'==7 {
+        mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.65)
+      }
+    }
+    if `difi'==7 {
+      qui local difj = runiformint(2,7)
+      if `difj'==2 {
+        mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.65)
+      }
+      if `difj'==3 {
+        mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.65)
+      }
+      if `difj'==4 {
+        mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.65)
+      }
+      if `difj'==5 {
+        mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.65)
+      }
+      if `difj'==6 {
+        mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.65)
+      }
+      if `difj'==7 {
+        mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.65)
+      }
+    }
+		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
+    qui gen TT = 1
+		tempfile grp1
+		qui save `grp1',replace
+
+
+		clear
+		use `grp0'
+    qui append using `grp1'
+		drop id
+		qui gen id = _n
+		order(id)
+		qui gen replication = `replication'
+		if `replication'==1{
+			tempfile data
+			qui save `data'
+		}
+		else{
+			qui append using `data'
+			qui save `data',replace
+		}
+	}
+	export delimited using "`path'/scenario_15A_100.csv", replace
+
+
+  * Scenario 15B : H_1 is TRUE / Effect size 0.2 / DIF on treatment x2
+
+  di "Scenario 15B - N=100"
+
+  forvalues replication = 1/1000 {
+    	if mod(`replication',10)==0 {
+    	di "replication = `replication'"
+      }
+      qui local difi = runiformint(1,7)
+
+      mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+  		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
+      qui gen TT = 0
+  		tempfile grp0
+  		qui save `grp0',replace
+
+
+      if `difi'==1 {
+        qui local difj = runiformint(2,7)
+        if `difj'==2 {
+          mat D= (-0.65 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+        }
+        if `difj'==3 {
+          mat D= (-0.65 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.15)
+        }
+        if `difj'==4 {
+          mat D= (-0.65 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+        }
+        if `difj'==5 {
+          mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.15)
+        }
+        if `difj'==6 {
+          mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.15)
+        }
+        if `difj'==7 {
+          mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.65)
+        }
+      }
+      if `difi'==2 {
+        qui local difj = runiformint(2,7)
+        if `difj'==2 {
+          mat D= (-0.65 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+        }
+        if `difj'==3 {
+          mat D= (-1.15 \ -0.17 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.15)
+        }
+        if `difj'==4 {
+          mat D= (-1.15 \ -0.17 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+        }
+        if `difj'==5 {
+          mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.15)
+        }
+        if `difj'==6 {
+          mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.15)
+        }
+        if `difj'==7 {
+          mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.65)
+        }
+      }
+
+      if `difi'==3 {
+        qui local difj = runiformint(2,7)
+        if `difj'==2 {
+          mat D= (-0.65 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.15)
+        }
+        if `difj'==3 {
+          mat D= (-1.15 \ -0.17 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.15)
+        }
+        if `difj'==4 {
+          mat D= (-1.15 \ -0.67 \ 0.18 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+        }
+        if `difj'==5 {
+          mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.82 \ 0.67 \ 1.15)
+        }
+        if `difj'==6 {
+          mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 1.17 \ 1.15)
+        }
+        if `difj'==7 {
+          mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.65)
+        }
+      }
+      if `difi'==4 {
+        qui local difj = runiformint(2,7)
+        if `difj'==2 {
+          mat D= (-0.65 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+        }
+        if `difj'==3 {
+          mat D= (-1.15 \ -0.17 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+        }
+        if `difj'==4 {
+          mat D= (-1.15 \ -0.67 \ 0.18 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+        }
+        if `difj'==5 {
+          mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.82 \ 0.67 \ 1.15)
+        }
+        if `difj'==6 {
+          mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 1.17 \ 1.15)
+        }
+        if `difj'==7 {
+          mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.65)
+        }
+      }
+      if `difi'==5 {
+        qui local difj = runiformint(2,7)
+        if `difj'==2 {
+          mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.15)
+        }
+        if `difj'==3 {
+          mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.15)
+        }
+        if `difj'==4 {
+          mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.82 \ 0.67 \ 1.15)
+        }
+        if `difj'==5 {
+          mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.82 \ 0.67 \ 1.15)
+        }
+        if `difj'==6 {
+          mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 1.17 \ 1.15)
+        }
+        if `difj'==7 {
+          mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.65)
+        }
+      }
+      if `difi'==6 {
+        qui local difj = runiformint(2,7)
+        if `difj'==2 {
+          mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.15)
+        }
+        if `difj'==3 {
+          mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.15)
+        }
+        if `difj'==4 {
+          mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 1.17 \ 1.15)
+        }
+        if `difj'==5 {
+          mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 1.17 \ 1.15)
+        }
+        if `difj'==6 {
+          mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 1.17 \ 1.15)
+        }
+        if `difj'==7 {
+          mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.65)
+        }
+      }
+      if `difi'==7 {
+        qui local difj = runiformint(2,7)
+        if `difj'==2 {
+          mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.65)
+        }
+        if `difj'==3 {
+          mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.65)
+        }
+        if `difj'==4 {
+          mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.65)
+        }
+        if `difj'==5 {
+          mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.65)
+        }
+        if `difj'==6 {
+          mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.65)
+        }
+        if `difj'==7 {
+          mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.65)
+        }
+      }
+  		qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear
+      qui gen TT = 1
+  		tempfile grp1
+  		qui save `grp1',replace
+
+
+  		clear
+  		use `grp0'
+      qui append using `grp1'
+  		drop id
+  		qui gen id = _n
+  		order(id)
+  		qui gen replication = `replication'
+  		if `replication'==1{
+  			tempfile data
+  			qui save `data'
+  		}
+  		else{
+  			qui append using `data'
+  			qui save `data',replace
+  		}
+  	}
+  	export delimited using "`path'/scenario_15B_100.csv", replace
+
+* Scenario 15C : H_1 is TRUE / Effect size 0.2 / DIF on control x2
+
+    di "Scenario 15C - N=100"
+
+    forvalues replication = 1/1000 {
+      	if mod(`replication',10)==0 {
+      	di "replication = `replication'"
+        }
+        qui local difi = runiformint(1,7)
+
+        if `difi'==1 {
+          qui local difj = runiformint(2,7)
+          if `difj'==2 {
+            mat D= (-0.65 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+          }
+          if `difj'==3 {
+            mat D= (-0.65 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.15)
+          }
+          if `difj'==4 {
+            mat D= (-0.65 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+          }
+          if `difj'==5 {
+            mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.15)
+          }
+          if `difj'==6 {
+            mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.15)
+          }
+          if `difj'==7 {
+            mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.65)
+          }
+        }
+        if `difi'==2 {
+          qui local difj = runiformint(2,7)
+          if `difj'==2 {
+            mat D= (-0.65 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+          }
+          if `difj'==3 {
+            mat D= (-1.15 \ -0.17 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.15)
+          }
+          if `difj'==4 {
+            mat D= (-1.15 \ -0.17 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+          }
+          if `difj'==5 {
+            mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.15)
+          }
+          if `difj'==6 {
+            mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.15)
+          }
+          if `difj'==7 {
+            mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.65)
+          }
+        }
+
+        if `difi'==3 {
+          qui local difj = runiformint(2,7)
+          if `difj'==2 {
+            mat D= (-0.65 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.15)
+          }
+          if `difj'==3 {
+            mat D= (-1.15 \ -0.17 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.15)
+          }
+          if `difj'==4 {
+            mat D= (-1.15 \ -0.67 \ 0.18 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+          }
+          if `difj'==5 {
+            mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.82 \ 0.67 \ 1.15)
+          }
+          if `difj'==6 {
+            mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 1.17 \ 1.15)
+          }
+          if `difj'==7 {
+            mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.65)
+          }
+        }
+        if `difi'==4 {
+          qui local difj = runiformint(2,7)
+          if `difj'==2 {
+            mat D= (-0.65 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+          }
+          if `difj'==3 {
+            mat D= (-1.15 \ -0.17 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+          }
+          if `difj'==4 {
+            mat D= (-1.15 \ -0.67 \ 0.18 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+          }
+          if `difj'==5 {
+            mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.82 \ 0.67 \ 1.15)
+          }
+          if `difj'==6 {
+            mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 1.17 \ 1.15)
+          }
+          if `difj'==7 {
+            mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.65)
+          }
+        }
+        if `difi'==5 {
+          qui local difj = runiformint(2,7)
+          if `difj'==2 {
+            mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.15)
+          }
+          if `difj'==3 {
+            mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.15)
+          }
+          if `difj'==4 {
+            mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.82 \ 0.67 \ 1.15)
+          }
+          if `difj'==5 {
+            mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.82 \ 0.67 \ 1.15)
+          }
+          if `difj'==6 {
+            mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 1.17 \ 1.15)
+          }
+          if `difj'==7 {
+            mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.65)
+          }
+        }
+        if `difi'==6 {
+          qui local difj = runiformint(2,7)
+          if `difj'==2 {
+            mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.15)
+          }
+          if `difj'==3 {
+            mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.15)
+          }
+          if `difj'==4 {
+            mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 1.17 \ 1.15)
+          }
+          if `difj'==5 {
+            mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 1.17 \ 1.15)
+          }
+          if `difj'==6 {
+            mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 1.17 \ 1.15)
+          }
+          if `difj'==7 {
+            mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.65)
+          }
+        }
+        if `difi'==7 {
+          qui local difj = runiformint(2,7)
+          if `difj'==2 {
+            mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.65)
+          }
+          if `difj'==3 {
+            mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.65)
+          }
+          if `difj'==4 {
+            mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.65)
+          }
+          if `difj'==5 {
+            mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.65)
+          }
+          if `difj'==6 {
+            mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.65)
+          }
+          if `difj'==7 {
+            mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.65)
+          }
+        }
+    		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
+        qui gen TT = 0
+    		tempfile grp0
+    		qui save `grp0',replace
+
+
+        mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+    		qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear
+        qui gen TT = 1
+    		tempfile grp1
+    		qui save `grp1',replace
+
+
+    		clear
+    		use `grp0'
+        qui append using `grp1'
+    		drop id
+    		qui gen id = _n
+    		order(id)
+    		qui gen replication = `replication'
+    		if `replication'==1{
+    			tempfile data
+    			qui save `data'
+    		}
+    		else{
+    			qui append using `data'
+    			qui save `data',replace
+    		}
+    	}
+    	export delimited using "`path'/scenario_15C_100.csv", replace
+
+
+* Scenario 15D : H_1 is TRUE / Effect size 0.4 / DIF on treatment x2
+
+        di "Scenario 15D - N=100"
+
+        forvalues replication = 1/1000 {
+          	if mod(`replication',10)==0 {
+          	di "replication = `replication'"
+            }
+            qui local difi = runiformint(1,7)
+
+            mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+        		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
+            qui gen TT = 0
+        		tempfile grp0
+        		qui save `grp0',replace
+
+
+            if `difi'==1 {
+              qui local difj = runiformint(2,7)
+              if `difj'==2 {
+                mat D= (-0.65 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+              }
+              if `difj'==3 {
+                mat D= (-0.65 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.15)
+              }
+              if `difj'==4 {
+                mat D= (-0.65 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+              }
+              if `difj'==5 {
+                mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.15)
+              }
+              if `difj'==6 {
+                mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.15)
+              }
+              if `difj'==7 {
+                mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.65)
+              }
+            }
+            if `difi'==2 {
+              qui local difj = runiformint(2,7)
+              if `difj'==2 {
+                mat D= (-0.65 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+              }
+              if `difj'==3 {
+                mat D= (-1.15 \ -0.17 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.15)
+              }
+              if `difj'==4 {
+                mat D= (-1.15 \ -0.17 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+              }
+              if `difj'==5 {
+                mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.15)
+              }
+              if `difj'==6 {
+                mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.15)
+              }
+              if `difj'==7 {
+                mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.65)
+              }
+            }
+
+            if `difi'==3 {
+              qui local difj = runiformint(2,7)
+              if `difj'==2 {
+                mat D= (-0.65 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.15)
+              }
+              if `difj'==3 {
+                mat D= (-1.15 \ -0.17 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.15)
+              }
+              if `difj'==4 {
+                mat D= (-1.15 \ -0.67 \ 0.18 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+              }
+              if `difj'==5 {
+                mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.82 \ 0.67 \ 1.15)
+              }
+              if `difj'==6 {
+                mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 1.17 \ 1.15)
+              }
+              if `difj'==7 {
+                mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.65)
+              }
+            }
+            if `difi'==4 {
+              qui local difj = runiformint(2,7)
+              if `difj'==2 {
+                mat D= (-0.65 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+              }
+              if `difj'==3 {
+                mat D= (-1.15 \ -0.17 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+              }
+              if `difj'==4 {
+                mat D= (-1.15 \ -0.67 \ 0.18 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+              }
+              if `difj'==5 {
+                mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.82 \ 0.67 \ 1.15)
+              }
+              if `difj'==6 {
+                mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 1.17 \ 1.15)
+              }
+              if `difj'==7 {
+                mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.65)
+              }
+            }
+            if `difi'==5 {
+              qui local difj = runiformint(2,7)
+              if `difj'==2 {
+                mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.15)
+              }
+              if `difj'==3 {
+                mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.15)
+              }
+              if `difj'==4 {
+                mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.82 \ 0.67 \ 1.15)
+              }
+              if `difj'==5 {
+                mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.82 \ 0.67 \ 1.15)
+              }
+              if `difj'==6 {
+                mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 1.17 \ 1.15)
+              }
+              if `difj'==7 {
+                mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.65)
+              }
+            }
+            if `difi'==6 {
+              qui local difj = runiformint(2,7)
+              if `difj'==2 {
+                mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.15)
+              }
+              if `difj'==3 {
+                mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.15)
+              }
+              if `difj'==4 {
+                mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 1.17 \ 1.15)
+              }
+              if `difj'==5 {
+                mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 1.17 \ 1.15)
+              }
+              if `difj'==6 {
+                mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 1.17 \ 1.15)
+              }
+              if `difj'==7 {
+                mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.65)
+              }
+            }
+            if `difi'==7 {
+              qui local difj = runiformint(2,7)
+              if `difj'==2 {
+                mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.65)
+              }
+              if `difj'==3 {
+                mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.65)
+              }
+              if `difj'==4 {
+                mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.65)
+              }
+              if `difj'==5 {
+                mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.65)
+              }
+              if `difj'==6 {
+                mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.65)
+              }
+              if `difj'==7 {
+                mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.65)
+              }
+            }
+        		qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear
+            qui gen TT = 1
+        		tempfile grp1
+        		qui save `grp1',replace
+
+
+        		clear
+        		use `grp0'
+            qui append using `grp1'
+        		drop id
+        		qui gen id = _n
+        		order(id)
+        		qui gen replication = `replication'
+        		if `replication'==1{
+        			tempfile data
+        			qui save `data'
+        		}
+        		else{
+        			qui append using `data'
+        			qui save `data',replace
+        		}
+        	}
+        	export delimited using "`path'/scenario_15D_100.csv", replace
+
+* Scenario 15E : H_1 is TRUE / Effect size 0.4 / DIF on control x2
+
+          di "Scenario 15E - N=100"
+
+          forvalues replication = 1/1000 {
+            	if mod(`replication',10)==0 {
+            	di "replication = `replication'"
+              }
+              qui local difi = runiformint(1,7)
+
+              if `difi'==1 {
+                qui local difj = runiformint(2,7)
+                if `difj'==2 {
+                  mat D= (-0.65 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+                }
+                if `difj'==3 {
+                  mat D= (-0.65 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.15)
+                }
+                if `difj'==4 {
+                  mat D= (-0.65 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+                }
+                if `difj'==5 {
+                  mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.15)
+                }
+                if `difj'==6 {
+                  mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.15)
+                }
+                if `difj'==7 {
+                  mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.65)
+                }
+              }
+              if `difi'==2 {
+                qui local difj = runiformint(2,7)
+                if `difj'==2 {
+                  mat D= (-0.65 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+                }
+                if `difj'==3 {
+                  mat D= (-1.15 \ -0.17 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.15)
+                }
+                if `difj'==4 {
+                  mat D= (-1.15 \ -0.17 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+                }
+                if `difj'==5 {
+                  mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.15)
+                }
+                if `difj'==6 {
+                  mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.15)
+                }
+                if `difj'==7 {
+                  mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.65)
+                }
+              }
+
+              if `difi'==3 {
+                qui local difj = runiformint(2,7)
+                if `difj'==2 {
+                  mat D= (-0.65 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.15)
+                }
+                if `difj'==3 {
+                  mat D= (-1.15 \ -0.17 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.15)
+                }
+                if `difj'==4 {
+                  mat D= (-1.15 \ -0.67 \ 0.18 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+                }
+                if `difj'==5 {
+                  mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.82 \ 0.67 \ 1.15)
+                }
+                if `difj'==6 {
+                  mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 1.17 \ 1.15)
+                }
+                if `difj'==7 {
+                  mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.65)
+                }
+              }
+              if `difi'==4 {
+                qui local difj = runiformint(2,7)
+                if `difj'==2 {
+                  mat D= (-0.65 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+                }
+                if `difj'==3 {
+                  mat D= (-1.15 \ -0.17 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+                }
+                if `difj'==4 {
+                  mat D= (-1.15 \ -0.67 \ 0.18 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+                }
+                if `difj'==5 {
+                  mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.82 \ 0.67 \ 1.15)
+                }
+                if `difj'==6 {
+                  mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 1.17 \ 1.15)
+                }
+                if `difj'==7 {
+                  mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.65)
+                }
+              }
+              if `difi'==5 {
+                qui local difj = runiformint(2,7)
+                if `difj'==2 {
+                  mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.15)
+                }
+                if `difj'==3 {
+                  mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.15)
+                }
+                if `difj'==4 {
+                  mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.82 \ 0.67 \ 1.15)
+                }
+                if `difj'==5 {
+                  mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.82 \ 0.67 \ 1.15)
+                }
+                if `difj'==6 {
+                  mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 1.17 \ 1.15)
+                }
+                if `difj'==7 {
+                  mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.65)
+                }
+              }
+              if `difi'==6 {
+                qui local difj = runiformint(2,7)
+                if `difj'==2 {
+                  mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.15)
+                }
+                if `difj'==3 {
+                  mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.15)
+                }
+                if `difj'==4 {
+                  mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 1.17 \ 1.15)
+                }
+                if `difj'==5 {
+                  mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 1.17 \ 1.15)
+                }
+                if `difj'==6 {
+                  mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 1.17 \ 1.15)
+                }
+                if `difj'==7 {
+                  mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.65)
+                }
+              }
+              if `difi'==7 {
+                qui local difj = runiformint(2,7)
+                if `difj'==2 {
+                  mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.65)
+                }
+                if `difj'==3 {
+                  mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.65)
+                }
+                if `difj'==4 {
+                  mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.65)
+                }
+                if `difj'==5 {
+                  mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.65)
+                }
+                if `difj'==6 {
+                  mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.65)
+                }
+                if `difj'==7 {
+                  mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.65)
+                }
+              }
+          		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
+              qui gen TT = 0
+          		tempfile grp0
+          		qui save `grp0',replace
+
+
+              mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+          		qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear
+              qui gen TT = 1
+          		tempfile grp1
+          		qui save `grp1',replace
+
+
+          		clear
+          		use `grp0'
+              qui append using `grp1'
+          		drop id
+          		qui gen id = _n
+          		order(id)
+          		qui gen replication = `replication'
+          		if `replication'==1{
+          			tempfile data
+          			qui save `data'
+          		}
+          		else{
+          			qui append using `data'
+          			qui save `data',replace
+          		}
+          	}
+          	export delimited using "`path'/scenario_15E_100.csv", replace
+
+
+* Scenario 15F : H_0 is TRUE / DIF on treatment + DIF on control (other item)
+
+di "Scenario 15F - N=100"
+
+forvalues replication = 1/1000 {
+  	if mod(`replication',10)==0 {
+  	di "replication = `replication'"
+    }
+    qui local difi = runiformint(1,7)
+    qui local difj = runiformint(1,7)
+    while `difi'==`difj' {
+      qui local difi = runiformint(1,7)
+      qui local difj = runiformint(1,7)
+    }
+
+    if `difi'==1 {
+      mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+    }
+    if `difi'==2 {
+      mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+    }
+    if `difi'==3 {
+      mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.15)
+    }
+    if `difi'==4 {
+      mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+    }
+    if `difi'==5 {
+      mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.15)
+    }
+    if `difi'==6 {
+      mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.15)
+    }
+    if `difi'==7 {
+      mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.65)
+    }
+		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
+    qui gen TT = 0
+		tempfile grp0
+		qui save `grp0',replace
+
+
+    if `difj'==1 {
+      mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+    }
+    if `difj'==2 {
+      mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+    }
+    if `difj'==3 {
+      mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.15)
+    }
+    if `difj'==4 {
+      mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+    }
+    if `difj'==5 {
+      mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.15)
+    }
+    if `difj'==6 {
+      mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.15)
+    }
+    if `difj'==7 {
+      mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.65)
+    }
+		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
+    qui gen TT = 1
+		tempfile grp1
+		qui save `grp1',replace
+
+
+		clear
+		use `grp0'
+    qui append using `grp1'
+		drop id
+		qui gen id = _n
+		order(id)
+		qui gen replication = `replication'
+		if `replication'==1{
+			tempfile data
+			qui save `data'
+		}
+		else{
+			qui append using `data'
+			qui save `data',replace
+		}
+	}
+	export delimited using "`path'/scenario_15F_100.csv", replace
+
+
+* Scenario 15G : H_1 is TRUE / Effect size 0.2 / DIF on treatment + DIF on control (other item)
+
+  di "Scenario 15G - N=100"
+
+  forvalues replication = 1/1000 {
+    	if mod(`replication',10)==0 {
+    	di "replication = `replication'"
+      }
+      qui local difi = runiformint(1,7)
+      qui local difj = runiformint(1,7)
+      while `difi'==`difj' {
+        qui local difi = runiformint(1,7)
+        qui local difj = runiformint(1,7)
+      }
+
+      if `difi'==1 {
+        mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+      }
+      if `difi'==2 {
+        mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+      }
+      if `difi'==3 {
+        mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.15)
+      }
+      if `difi'==4 {
+        mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+      }
+      if `difi'==5 {
+        mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.15)
+      }
+      if `difi'==6 {
+        mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.15)
+      }
+      if `difi'==7 {
+        mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.65)
+      }
+  		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
+      qui gen TT = 0
+  		tempfile grp0
+  		qui save `grp0',replace
+
+
+      if `difj'==1 {
+        mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+      }
+      if `difj'==2 {
+        mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+      }
+      if `difj'==3 {
+        mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.15)
+      }
+      if `difj'==4 {
+        mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+      }
+      if `difj'==5 {
+        mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.15)
+      }
+      if `difj'==6 {
+        mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.15)
+      }
+      if `difj'==7 {
+        mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.65)
+      }
+  		qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear
+      qui gen TT = 1
+  		tempfile grp1
+  		qui save `grp1',replace
+
+
+  		clear
+  		use `grp0'
+      qui append using `grp1'
+  		drop id
+  		qui gen id = _n
+  		order(id)
+  		qui gen replication = `replication'
+  		if `replication'==1{
+  			tempfile data
+  			qui save `data'
+  		}
+  		else{
+  			qui append using `data'
+  			qui save `data',replace
+  		}
+  	}
+  	export delimited using "`path'/scenario_15G_100.csv", replace
+
+* Scenario 15H : H_1 is TRUE / Effect size 0.4 / DIF on treatment + DIF on control (other item)
+
+      di "Scenario 15H - N=100"
+
+      forvalues replication = 1/1000 {
+        	if mod(`replication',10)==0 {
+        	di "replication = `replication'"
+          }
+          qui local difi = runiformint(1,7)
+          qui local difj = runiformint(1,7)
+          while `difi'==`difj' {
+            qui local difi = runiformint(1,7)
+            qui local difj = runiformint(1,7)
+          }
+
+          if `difi'==1 {
+            mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+          }
+          if `difi'==2 {
+            mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+          }
+          if `difi'==3 {
+            mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.15)
+          }
+          if `difi'==4 {
+            mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+          }
+          if `difi'==5 {
+            mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.15)
+          }
+          if `difi'==6 {
+            mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.15)
+          }
+          if `difi'==7 {
+            mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.65)
+          }
+      		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
+          qui gen TT = 0
+      		tempfile grp0
+      		qui save `grp0',replace
+
+
+          if `difj'==1 {
+            mat D= (-0.65 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+          }
+          if `difj'==2 {
+            mat D= (-1.15 \ -0.17 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+          }
+          if `difj'==3 {
+            mat D= (-1.15 \ -0.67 \ 0.18 \ 0 \ 0.32 \ 0.67 \ 1.15)
+          }
+          if `difj'==4 {
+            mat D= (-1.15 \ -0.67 \ -0.32 \ 0.5 \ 0.32 \ 0.67 \ 1.15)
+          }
+          if `difj'==5 {
+            mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.82 \ 0.67 \ 1.15)
+          }
+          if `difj'==6 {
+            mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 1.17 \ 1.15)
+          }
+          if `difj'==7 {
+            mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.65)
+          }
+      		qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear
+          qui gen TT = 1
+      		tempfile grp1
+      		qui save `grp1',replace
+
+
+      		clear
+      		use `grp0'
+          qui append using `grp1'
+      		drop id
+      		qui gen id = _n
+      		order(id)
+      		qui gen replication = `replication'
+      		if `replication'==1{
+      			tempfile data
+      			qui save `data'
+      		}
+      		else{
+      			qui append using `data'
+      			qui save `data',replace
+      		}
+      	}
+      	export delimited using "`path'/scenario_15H_100.csv", replace
+
+
+
+
+
+
+
+
+
+** Scenario 16: J = 7 items / M = 4 modalities / DIF SIZE = 0.5
+
+* Scenario 16A : H_0 is TRUE / DIF on treatment x 2
+
+di "Scenario 16A - N=100"
+
+forvalues replication = 1/1000 {
+  	if mod(`replication',10)==0 {
+  	di "replication = `replication'"
+    }
+    qui local difi = runiformint(1,7)
+
+    mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
+    qui gen TT = 0
+		tempfile grp0
+		qui save `grp0',replace
+
+
+    if `difi'==1 {
+      qui local difj = runiformint(2,7)
+      if `difj'==2 {
+        mat D= (-1.65,-0.65,0.35 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+      }
+      if `difj'==3 {
+        mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+      }
+      if `difj'==4 {
+        mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+      }
+      if `difj'==5 {
+        mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+      }
+      if `difj'==6 {
+        mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+      }
+      if `difj'==7 {
+        mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+      }
+    }
+    if `difi'==2 {
+      qui local difj = runiformint(2,7)
+      if `difj'==2 {
+        mat D= (-1.65,-0.65,0.35 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+      }
+      if `difj'==3 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+      }
+      if `difj'==4 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+      }
+      if `difj'==5 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+      }
+      if `difj'==6 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+      }
+      if `difj'==7 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+      }
+    }
+
+    if `difi'==3 {
+      qui local difj = runiformint(2,7)
+      if `difj'==2 {
+        mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+      }
+      if `difj'==3 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+      }
+      if `difj'==4 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+      }
+      if `difj'==5 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+      }
+      if `difj'==6 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+      }
+      if `difj'==7 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+      }
+    }
+    if `difi'==4 {
+      qui local difj = runiformint(2,7)
+      if `difj'==2 {
+        mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+      }
+      if `difj'==3 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+      }
+      if `difj'==4 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+      }
+      if `difj'==5 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+      }
+      if `difj'==6 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+      }
+      if `difj'==7 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+      }
+    }
+    if `difi'==5 {
+      qui local difj = runiformint(2,7)
+      if `difj'==2 {
+        mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+      }
+      if `difj'==3 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+      }
+      if `difj'==4 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+      }
+      if `difj'==5 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+      }
+      if `difj'==6 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+      }
+      if `difj'==7 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+      }
+    }
+
+    if `difi'==6 {
+      qui local difj = runiformint(2,7)
+      if `difj'==2 {
+        mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+      }
+      if `difj'==3 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+      }
+      if `difj'==4 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+      }
+      if `difj'==5 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+      }
+      if `difj'==6 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+      }
+      if `difj'==7 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.65,1.65,2.65)
+      }
+    }
+    if `difi'==7 {
+      qui local difj = runiformint(2,7)
+      if `difj'==2 {
+        mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+      }
+      if `difj'==3 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+      }
+      if `difj'==4 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+      }
+      if `difj'==5 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+      }
+      if `difj'==6 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+      }
+      if `difj'==7 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.65,1.65,2.65)
+      }
+    }
+		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
+    qui gen TT = 1
+		tempfile grp1
+		qui save `grp1',replace
+
+
+		clear
+		use `grp0'
+    qui append using `grp1'
+		drop id
+		qui gen id = _n
+		order(id)
+		qui gen replication = `replication'
+		if `replication'==1{
+			tempfile data
+			qui save `data'
+		}
+		else{
+			qui append using `data'
+			qui save `data',replace
+		}
+	}
+	export delimited using "`path'/scenario_16A_100.csv", replace
+
+
+  * Scenario 16B : H_1 is TRUE / Effect size 0.2 / DIF on treatment x2
+
+  di "Scenario 16B - N=100"
+
+  forvalues replication = 1/1000 {
+    	if mod(`replication',10)==0 {
+    	di "replication = `replication'"
+      }
+      qui local difi = runiformint(1,7)
+
+      mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+  		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
+      qui gen TT = 0
+  		tempfile grp0
+  		qui save `grp0',replace
+
+
+      if `difi'==1 {
+        qui local difj = runiformint(2,7)
+        if `difj'==2 {
+          mat D= (-1.65,-0.65,0.35 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+        }
+        if `difj'==3 {
+          mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+        }
+        if `difj'==4 {
+          mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+        }
+        if `difj'==5 {
+          mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+        }
+        if `difj'==6 {
+          mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+        }
+        if `difj'==7 {
+          mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+        }
+      }
+      if `difi'==2 {
+        qui local difj = runiformint(2,7)
+        if `difj'==2 {
+          mat D= (-1.65,-0.65,0.35 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+        }
+        if `difj'==3 {
+          mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+        }
+        if `difj'==4 {
+          mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+        }
+        if `difj'==5 {
+          mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+        }
+        if `difj'==6 {
+          mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+        }
+        if `difj'==7 {
+          mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+        }
+      }
+
+      if `difi'==3 {
+        qui local difj = runiformint(2,7)
+        if `difj'==2 {
+          mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+        }
+        if `difj'==3 {
+          mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+        }
+        if `difj'==4 {
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+        }
+        if `difj'==5 {
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+        }
+        if `difj'==6 {
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+        }
+        if `difj'==7 {
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+        }
+      }
+      if `difi'==4 {
+        qui local difj = runiformint(2,7)
+        if `difj'==2 {
+          mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+        }
+        if `difj'==3 {
+          mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+        }
+        if `difj'==4 {
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+        }
+        if `difj'==5 {
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+        }
+        if `difj'==6 {
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+        }
+        if `difj'==7 {
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+        }
+      }
+      if `difi'==5 {
+        qui local difj = runiformint(2,7)
+        if `difj'==2 {
+          mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+        }
+        if `difj'==3 {
+          mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+        }
+        if `difj'==4 {
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+        }
+        if `difj'==5 {
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+        }
+        if `difj'==6 {
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+        }
+        if `difj'==7 {
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+        }
+      }
+
+      if `difi'==6 {
+        qui local difj = runiformint(2,7)
+        if `difj'==2 {
+          mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+        }
+        if `difj'==3 {
+          mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+        }
+        if `difj'==4 {
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+        }
+        if `difj'==5 {
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+        }
+        if `difj'==6 {
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+        }
+        if `difj'==7 {
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.65,1.65,2.65)
+        }
+      }
+      if `difi'==7 {
+        qui local difj = runiformint(2,7)
+        if `difj'==2 {
+          mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+        }
+        if `difj'==3 {
+          mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+        }
+        if `difj'==4 {
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+        }
+        if `difj'==5 {
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+        }
+        if `difj'==6 {
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+        }
+        if `difj'==7 {
+          mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.65,1.65,2.65)
+        }
+      }
+  		qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear
+      qui gen TT = 1
+  		tempfile grp1
+  		qui save `grp1',replace
+
+
+  		clear
+  		use `grp0'
+      qui append using `grp1'
+  		drop id
+  		qui gen id = _n
+  		order(id)
+  		qui gen replication = `replication'
+  		if `replication'==1{
+  			tempfile data
+  			qui save `data'
+  		}
+  		else{
+  			qui append using `data'
+  			qui save `data',replace
+  		}
+  	}
+  	export delimited using "`path'/scenario_16B_100.csv", replace
+
+* Scenario 16C : H_1 is TRUE / Effect size 0.2 / DIF on control x2
+
+    di "Scenario 16C - N=100"
+
+    forvalues replication = 1/1000 {
+      	if mod(`replication',10)==0 {
+      	di "replication = `replication'"
+        }
+        qui local difi = runiformint(1,7)
+
+        if `difi'==1 {
+          qui local difj = runiformint(2,7)
+          if `difj'==2 {
+            mat D= (-1.65,-0.65,0.35 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+          }
+          if `difj'==3 {
+            mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+          }
+          if `difj'==4 {
+            mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+          }
+          if `difj'==5 {
+            mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+          }
+          if `difj'==6 {
+            mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+          }
+          if `difj'==7 {
+            mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+          }
+        }
+        if `difi'==2 {
+          qui local difj = runiformint(2,7)
+          if `difj'==2 {
+            mat D= (-1.65,-0.65,0.35 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+          }
+          if `difj'==3 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+          }
+          if `difj'==4 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+          }
+          if `difj'==5 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+          }
+          if `difj'==6 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+          }
+          if `difj'==7 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+          }
+        }
+
+        if `difi'==3 {
+          qui local difj = runiformint(2,7)
+          if `difj'==2 {
+            mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+          }
+          if `difj'==3 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+          }
+          if `difj'==4 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+          }
+          if `difj'==5 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+          }
+          if `difj'==6 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+          }
+          if `difj'==7 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+          }
+        }
+        if `difi'==4 {
+          qui local difj = runiformint(2,7)
+          if `difj'==2 {
+            mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+          }
+          if `difj'==3 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+          }
+          if `difj'==4 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+          }
+          if `difj'==5 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+          }
+          if `difj'==6 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+          }
+          if `difj'==7 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+          }
+        }
+        if `difi'==5 {
+          qui local difj = runiformint(2,7)
+          if `difj'==2 {
+            mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+          }
+          if `difj'==3 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+          }
+          if `difj'==4 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+          }
+          if `difj'==5 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+          }
+          if `difj'==6 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+          }
+          if `difj'==7 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+          }
+        }
+
+        if `difi'==6 {
+          qui local difj = runiformint(2,7)
+          if `difj'==2 {
+            mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+          }
+          if `difj'==3 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+          }
+          if `difj'==4 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+          }
+          if `difj'==5 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+          }
+          if `difj'==6 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+          }
+          if `difj'==7 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.65,1.65,2.65)
+          }
+        }
+        if `difi'==7 {
+          qui local difj = runiformint(2,7)
+          if `difj'==2 {
+            mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+          }
+          if `difj'==3 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+          }
+          if `difj'==4 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+          }
+          if `difj'==5 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+          }
+          if `difj'==6 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+          }
+          if `difj'==7 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.65,1.65,2.65)
+          }
+        }
+    		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
+        qui gen TT = 0
+    		tempfile grp0
+    		qui save `grp0',replace
+
+
+        mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+    		qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear
+        qui gen TT = 1
+    		tempfile grp1
+    		qui save `grp1',replace
+
+
+    		clear
+    		use `grp0'
+        qui append using `grp1'
+    		drop id
+    		qui gen id = _n
+    		order(id)
+    		qui gen replication = `replication'
+    		if `replication'==1{
+    			tempfile data
+    			qui save `data'
+    		}
+    		else{
+    			qui append using `data'
+    			qui save `data',replace
+    		}
+    	}
+    	export delimited using "`path'/scenario_16C_100.csv", replace
+
+
+* Scenario 16D : H_1 is TRUE / Effect size 0.4 / DIF on treatment x2
+
+        di "Scenario 16D - N=100"
+
+        forvalues replication = 1/1000 {
+          	if mod(`replication',10)==0 {
+          	di "replication = `replication'"
+            }
+            qui local difi = runiformint(1,7)
+
+            mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+        		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
+            qui gen TT = 0
+        		tempfile grp0
+        		qui save `grp0',replace
+
+
+            if `difi'==1 {
+              qui local difj = runiformint(2,7)
+              if `difj'==2 {
+                mat D= (-1.65,-0.65,0.35 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+              }
+              if `difj'==3 {
+                mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+              }
+              if `difj'==4 {
+                mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+              }
+              if `difj'==5 {
+                mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+              }
+              if `difj'==6 {
+                mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+              }
+              if `difj'==7 {
+                mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+              }
+            }
+            if `difi'==2 {
+              qui local difj = runiformint(2,7)
+              if `difj'==2 {
+                mat D= (-1.65,-0.65,0.35 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+              }
+              if `difj'==3 {
+                mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+              }
+              if `difj'==4 {
+                mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+              }
+              if `difj'==5 {
+                mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+              }
+              if `difj'==6 {
+                mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+              }
+              if `difj'==7 {
+                mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+              }
+            }
+
+            if `difi'==3 {
+              qui local difj = runiformint(2,7)
+              if `difj'==2 {
+                mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+              }
+              if `difj'==3 {
+                mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+              }
+              if `difj'==4 {
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+              }
+              if `difj'==5 {
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+              }
+              if `difj'==6 {
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+              }
+              if `difj'==7 {
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+              }
+            }
+            if `difi'==4 {
+              qui local difj = runiformint(2,7)
+              if `difj'==2 {
+                mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+              }
+              if `difj'==3 {
+                mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+              }
+              if `difj'==4 {
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+              }
+              if `difj'==5 {
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+              }
+              if `difj'==6 {
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+              }
+              if `difj'==7 {
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+              }
+            }
+            if `difi'==5 {
+              qui local difj = runiformint(2,7)
+              if `difj'==2 {
+                mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+              }
+              if `difj'==3 {
+                mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+              }
+              if `difj'==4 {
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+              }
+              if `difj'==5 {
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+              }
+              if `difj'==6 {
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+              }
+              if `difj'==7 {
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+              }
+            }
+
+            if `difi'==6 {
+              qui local difj = runiformint(2,7)
+              if `difj'==2 {
+                mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+              }
+              if `difj'==3 {
+                mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+              }
+              if `difj'==4 {
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+              }
+              if `difj'==5 {
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+              }
+              if `difj'==6 {
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+              }
+              if `difj'==7 {
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.65,1.65,2.65)
+              }
+            }
+            if `difi'==7 {
+              qui local difj = runiformint(2,7)
+              if `difj'==2 {
+                mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+              }
+              if `difj'==3 {
+                mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+              }
+              if `difj'==4 {
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+              }
+              if `difj'==5 {
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+              }
+              if `difj'==6 {
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+              }
+              if `difj'==7 {
+                mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.65,1.65,2.65)
+              }
+            }
+        		qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear
+            qui gen TT = 1
+        		tempfile grp1
+        		qui save `grp1',replace
+
+
+        		clear
+        		use `grp0'
+            qui append using `grp1'
+        		drop id
+        		qui gen id = _n
+        		order(id)
+        		qui gen replication = `replication'
+        		if `replication'==1{
+        			tempfile data
+        			qui save `data'
+        		}
+        		else{
+        			qui append using `data'
+        			qui save `data',replace
+        		}
+        	}
+        	export delimited using "`path'/scenario_16D_100.csv", replace
+
+* Scenario 16E : H_1 is TRUE / Effect size 0.4 / DIF on control x2
+
+          di "Scenario 16E - N=100"
+
+          forvalues replication = 1/1000 {
+            	if mod(`replication',10)==0 {
+            	di "replication = `replication'"
+              }
+              qui local difi = runiformint(1,7)
+
+              if `difi'==1 {
+                qui local difj = runiformint(2,7)
+                if `difj'==2 {
+                  mat D= (-1.65,-0.65,0.35 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+                }
+                if `difj'==3 {
+                  mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+                }
+                if `difj'==4 {
+                  mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+                }
+                if `difj'==5 {
+                  mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+                }
+                if `difj'==6 {
+                  mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+                }
+                if `difj'==7 {
+                  mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+                }
+              }
+              if `difi'==2 {
+                qui local difj = runiformint(2,7)
+                if `difj'==2 {
+                  mat D= (-1.65,-0.65,0.35 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+                }
+                if `difj'==3 {
+                  mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+                }
+                if `difj'==4 {
+                  mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+                }
+                if `difj'==5 {
+                  mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+                }
+                if `difj'==6 {
+                  mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+                }
+                if `difj'==7 {
+                  mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+                }
+              }
+
+              if `difi'==3 {
+                qui local difj = runiformint(2,7)
+                if `difj'==2 {
+                  mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+                }
+                if `difj'==3 {
+                  mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+                }
+                if `difj'==4 {
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+                }
+                if `difj'==5 {
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+                }
+                if `difj'==6 {
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+                }
+                if `difj'==7 {
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+                }
+              }
+              if `difi'==4 {
+                qui local difj = runiformint(2,7)
+                if `difj'==2 {
+                  mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+                }
+                if `difj'==3 {
+                  mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+                }
+                if `difj'==4 {
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+                }
+                if `difj'==5 {
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+                }
+                if `difj'==6 {
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+                }
+                if `difj'==7 {
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+                }
+              }
+              if `difi'==5 {
+                qui local difj = runiformint(2,7)
+                if `difj'==2 {
+                  mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+                }
+                if `difj'==3 {
+                  mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+                }
+                if `difj'==4 {
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+                }
+                if `difj'==5 {
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+                }
+                if `difj'==6 {
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+                }
+                if `difj'==7 {
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+                }
+              }
+
+              if `difi'==6 {
+                qui local difj = runiformint(2,7)
+                if `difj'==2 {
+                  mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+                }
+                if `difj'==3 {
+                  mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+                }
+                if `difj'==4 {
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+                }
+                if `difj'==5 {
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+                }
+                if `difj'==6 {
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+                }
+                if `difj'==7 {
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.65,1.65,2.65)
+                }
+              }
+              if `difi'==7 {
+                qui local difj = runiformint(2,7)
+                if `difj'==2 {
+                  mat D= (-1.65,-0.65,0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+                }
+                if `difj'==3 {
+                  mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+                }
+                if `difj'==4 {
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+                }
+                if `difj'==5 {
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+                }
+                if `difj'==6 {
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+                }
+                if `difj'==7 {
+                  mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.65,1.65,2.65)
+                }
+              }
+          		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
+              qui gen TT = 0
+          		tempfile grp0
+          		qui save `grp0',replace
+
+
+              mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
+          		qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear
+              qui gen TT = 1
+          		tempfile grp1
+          		qui save `grp1',replace
+
+
+          		clear
+          		use `grp0'
+              qui append using `grp1'
+          		drop id
+          		qui gen id = _n
+          		order(id)
+          		qui gen replication = `replication'
+          		if `replication'==1{
+          			tempfile data
+          			qui save `data'
+          		}
+          		else{
+          			qui append using `data'
+          			qui save `data',replace
+          		}
+          	}
+          	export delimited using "`path'/scenario_16E_100.csv", replace
+
+
+* Scenario 16F : H_0 is TRUE / DIF on treatment + DIF on control (other item)
+
+di "Scenario 16F - N=100"
+
+forvalues replication = 1/1000 {
+  	if mod(`replication',10)==0 {
+  	di "replication = `replication'"
+    }
+    qui local difi = runiformint(1,7)
+    qui local difj = runiformint(1,7)
+    while `difi'==`difj' {
+      qui local difi = runiformint(1,7)
+      qui local difj = runiformint(1,7)
+    }
+
+    if `difi'==1 {
+      mat D= (-1.65,-0.65, 0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+    }
+    if `difi'==2 {
+      mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+    }
+    if `difi'==3 {
+      mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+    }
+    if `difi'==4 {
+      mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+    }
+    if `difi'==5 {
+      mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+    }
+    if `difi'==6 {
+      mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+    }
+    if `difi'==7 {
+      mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+    }
+		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
+    qui gen TT = 0
+		tempfile grp0
+		qui save `grp0',replace
+
+
+    if `difj'==1 {
+      mat D= (-1.65,-0.65, 0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+    }
+    if `difj'==2 {
+      mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+    }
+    if `difj'==3 {
+      mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+    }
+    if `difj'==4 {
+      mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+    }
+    if `difj'==5 {
+      mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+    }
+    if `difj'==6 {
+      mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+    }
+    if `difj'==7 {
+      mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+    }
+		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
+    qui gen TT = 1
+		tempfile grp1
+		qui save `grp1',replace
+
+
+		clear
+		use `grp0'
+    qui append using `grp1'
+		drop id
+		qui gen id = _n
+		order(id)
+		qui gen replication = `replication'
+		if `replication'==1{
+			tempfile data
+			qui save `data'
+		}
+		else{
+			qui append using `data'
+			qui save `data',replace
+		}
+	}
+	export delimited using "`path'/scenario_16F_100.csv", replace
+
+
+* Scenario 16G : H_1 is TRUE / Effect size 0.2 / DIF on treatment + DIF on control (other item)
+
+  di "Scenario 16G - N=100"
+
+  forvalues replication = 1/1000 {
+    	if mod(`replication',10)==0 {
+    	di "replication = `replication'"
+      }
+      qui local difi = runiformint(1,7)
+      qui local difj = runiformint(1,7)
+      while `difi'==`difj' {
+        qui local difi = runiformint(1,7)
+        qui local difj = runiformint(1,7)
+      }
+
+      if `difi'==1 {
+        mat D= (-1.65,-0.65, 0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+      }
+      if `difi'==2 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+      }
+      if `difi'==3 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+      }
+      if `difi'==4 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+      }
+      if `difi'==5 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+      }
+      if `difi'==6 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+      }
+      if `difi'==7 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+      }
+  		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
+      qui gen TT = 0
+  		tempfile grp0
+  		qui save `grp0',replace
+
+
+      if `difj'==1 {
+        mat D= (-1.65,-0.65, 0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+      }
+      if `difj'==2 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+      }
+      if `difj'==3 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+      }
+      if `difj'==4 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+      }
+      if `difj'==5 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+      }
+      if `difj'==6 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+      }
+      if `difj'==7 {
+        mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+      }
+  		qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear
+      qui gen TT = 1
+  		tempfile grp1
+  		qui save `grp1',replace
+
+
+  		clear
+  		use `grp0'
+      qui append using `grp1'
+  		drop id
+  		qui gen id = _n
+  		order(id)
+  		qui gen replication = `replication'
+  		if `replication'==1{
+  			tempfile data
+  			qui save `data'
+  		}
+  		else{
+  			qui append using `data'
+  			qui save `data',replace
+  		}
+  	}
+  	export delimited using "`path'/scenario_16G_100.csv", replace
+
+* Scenario 16H : H_1 is TRUE / Effect size 0.4 / DIF on treatment + DIF on control (other item)
+
+      di "Scenario 16H - N=100"
+
+      forvalues replication = 1/1000 {
+        	if mod(`replication',10)==0 {
+        	di "replication = `replication'"
+          }
+          qui local difi = runiformint(1,7)
+          qui local difj = runiformint(1,7)
+          while `difi'==`difj' {
+            qui local difi = runiformint(1,7)
+            qui local difj = runiformint(1,7)
+          }
+
+          if `difi'==1 {
+            mat D= (-1.65,-0.65, 0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+          }
+          if `difi'==2 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+          }
+          if `difi'==3 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+          }
+          if `difi'==4 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+          }
+          if `difi'==5 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+          }
+          if `difi'==6 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+          }
+          if `difi'==7 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+          }
+      		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
+          qui gen TT = 0
+      		tempfile grp0
+      		qui save `grp0',replace
+
+
+          if `difj'==1 {
+            mat D= (-1.65,-0.65, 0.35 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+          }
+          if `difj'==2 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.17,-0.17,0.83 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+          }
+          if `difj'==3 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -0.82,0.18,1.18 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+          }
+          if `difj'==4 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -0.5,0.5,1.5 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+          }
+          if `difj'==5 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.18,0.82,1.82 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
+          }
+          if `difj'==6 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ 0.17,1.17,2.17 \ 0.15,1.15,2.15)
+          }
+          if `difj'==7 {
+            mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.65,1.65,2.65)
+          }
+      		qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear
+          qui gen TT = 1
+      		tempfile grp1
+      		qui save `grp1',replace
+
+
+      		clear
+      		use `grp0'
+          qui append using `grp1'
+      		drop id
+      		qui gen id = _n
+      		order(id)
+      		qui gen replication = `replication'
+      		if `replication'==1{
+      			tempfile data
+      			qui save `data'
+      		}
+      		else{
+      			qui append using `data'
+      			qui save `data',replace
+      		}
+      	}
+      	export delimited using "`path'/scenario_16H_100.csv", replace
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+set obs 7
+gen idif=_n
+sample 3,count
+levelsof idif
+local difi = r(levels)