diff --git a/Scripts/Scenarios/NoDIF/scenarios_noDIF_baseline.do b/Scripts/Scenarios/NoDIF/scenarios_noDIF_baseline.do
index 18b8780..d3ab757 100644
--- a/Scripts/Scenarios/NoDIF/scenarios_noDIF_baseline.do
+++ b/Scripts/Scenarios/NoDIF/scenarios_noDIF_baseline.do
@@ -40,13 +40,13 @@ forvalues replication = 1/1000 {
     }
 	  di "replication = `replication'"
 		mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+		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= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
 		qui gen TT = 1
 		tempfile grp1
 		qui save `grp1',replace
@@ -79,13 +79,13 @@ forvalues replication = 1/1000 {
 	  	di "replication = `replication'"
 	    }
   		mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-  		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+  		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= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-  		qui simirt, nbobs(Nn) mu(0.2) cov(1) dim(4) pcm(D) clear
+  		qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear
   		qui gen TT = 1
   		tempfile grp1
   		qui save `grp1',replace
@@ -117,13 +117,13 @@ forvalues replication = 1/1000 {
   	di "replication = `replication'"
     }
   		mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-  		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+  		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= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-  		qui simirt, nbobs(Nn) mu(0.4) cov(1) dim(4) pcm(D) clear
+  		qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear
   		qui gen TT = 1
   		tempfile grp1
   		qui save `grp1',replace
@@ -160,13 +160,13 @@ forvalues replication = 1/1000 {
 		  	di "replication = `replication'"
 		    }
     		mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-    		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+    		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.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-    		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+    		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
     		qui gen TT = 1
     		tempfile grp1
     		qui save `grp1',replace
@@ -198,13 +198,13 @@ forvalues replication = 1/1000 {
 			  	di "replication = `replication'"
 			    }
           mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-      		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+      		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.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-      		qui simirt, nbobs(Nn) mu(0.2) cov(1) dim(4) pcm(D) clear
+      		qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear
       		qui gen TT = 1
       		tempfile grp1
       		qui save `grp1',replace
@@ -235,13 +235,13 @@ forvalues replication = 1/1000 {
 		  	di "replication = `replication'"
 		    }
           mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-      		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+      		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.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-      		qui simirt, nbobs(Nn) mu(0.4) cov(1) dim(4) pcm(D) clear
+      		qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear
       		qui gen TT = 1
       		tempfile grp1
       		qui save `grp1',replace
@@ -278,13 +278,13 @@ forvalues replication = 1/1000 {
   	di "replication = `replication'"
     }
 		mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
-		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) 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) cov(1) dim(7) pcm(D) clear
+		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
 		qui gen TT = 1
 		tempfile grp1
 		qui save `grp1',replace
@@ -316,13 +316,13 @@ forvalues replication = 1/1000 {
 	  	di "replication = `replication'"
 	    }
       mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
-  		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+  		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) 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(7) pcm(D) clear
+  		qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(7) pcm(D) clear
   		qui gen TT = 1
   		tempfile grp1
   		qui save `grp1',replace
@@ -353,13 +353,13 @@ forvalues replication = 1/1000 {
   	di "replication = `replication'"
     }
       mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
-  		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+  		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) 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(7) pcm(D) clear
+  		qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(7) pcm(D) clear
   		qui gen TT = 1
   		tempfile grp1
   		qui save `grp1',replace
@@ -396,13 +396,13 @@ forvalues replication = 1/1000 {
 		  	di "replication = `replication'"
 		    }
     		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.15,1.15,2.15)
-    		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+    		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
     		qui gen TT = 0
     		tempfile grp0
     		qui save `grp0',replace
 
         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.15,1.15,2.15)
-    		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+    		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
     		qui gen TT = 1
     		tempfile grp1
     		qui save `grp1',replace
@@ -434,13 +434,13 @@ forvalues replication = 1/1000 {
 		  	di "replication = `replication'"
 		    }
           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.15,1.15,2.15)
-      		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+      		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
       		qui gen TT = 0
       		tempfile grp0
       		qui save `grp0',replace
 
           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.15,1.15,2.15)
-      		qui simirt, nbobs(Nn) mu(0.2) cov(1) dim(7) pcm(D) clear
+      		qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(7) pcm(D) clear
       		qui gen TT = 1
       		tempfile grp1
       		qui save `grp1',replace
@@ -471,13 +471,13 @@ forvalues replication = 1/1000 {
 		  	di "replication = `replication'"
 		    }
           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.15,1.15,2.15)
-      		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+      		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
       		qui gen TT = 0
       		tempfile grp0
       		qui save `grp0',replace
 
           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.15,1.15,2.15)
-      		qui simirt, nbobs(Nn) mu(0.4) cov(1) dim(7) pcm(D) clear
+      		qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(7) pcm(D) clear
       		qui gen TT = 1
       		tempfile grp1
       		qui save `grp1',replace
@@ -545,13 +545,13 @@ local Nn = 200
 		  	di "replication = `replication'"
 		    }
         		mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-        		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+        		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= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-        		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+        		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
         		qui gen TT = 1
         		tempfile grp1
         		qui save `grp1',replace
@@ -583,13 +583,13 @@ local Nn = 200
 		  	di "replication = `replication'"
 		    }
           		mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-          		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+          		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= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-          		qui simirt, nbobs(Nn) mu(0.2) cov(1) dim(4) pcm(D) clear
+          		qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear
           		qui gen TT = 1
           		tempfile grp1
           		qui save `grp1',replace
@@ -620,13 +620,13 @@ local Nn = 200
 		  	di "replication = `replication'"
 		    }
           		mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-          		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+          		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= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-          		qui simirt, nbobs(Nn) mu(0.4) cov(1) dim(4) pcm(D) clear
+          		qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear
           		qui gen TT = 1
           		tempfile grp1
           		qui save `grp1',replace
@@ -663,13 +663,13 @@ local Nn = 200
 		  	di "replication = `replication'"
 		    }
             		mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-            		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+            		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.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-            		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+            		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
             		qui gen TT = 1
             		tempfile grp1
             		qui save `grp1',replace
@@ -701,13 +701,13 @@ local Nn = 200
 		  	di "replication = `replication'"
 		    }
                   mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-              		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+              		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.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-              		qui simirt, nbobs(Nn) mu(0.2) cov(1) dim(4) pcm(D) clear
+              		qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear
               		qui gen TT = 1
               		tempfile grp1
               		qui save `grp1',replace
@@ -738,13 +738,13 @@ local Nn = 200
 		  	di "replication = `replication'"
 		    }
                   mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-              		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+              		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.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-              		qui simirt, nbobs(Nn) mu(0.4) cov(1) dim(4) pcm(D) clear
+              		qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear
               		qui gen TT = 1
               		tempfile grp1
               		qui save `grp1',replace
@@ -781,13 +781,13 @@ local Nn = 200
 		  	di "replication = `replication'"
 		    }
         		mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
-        		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+        		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) 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) cov(1) dim(7) pcm(D) clear
+        		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
         		qui gen TT = 1
         		tempfile grp1
         		qui save `grp1',replace
@@ -819,13 +819,13 @@ local Nn = 200
 		  	di "replication = `replication'"
 		    }
               mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
-          		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+          		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) 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(7) pcm(D) clear
+          		qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(7) pcm(D) clear
           		qui gen TT = 1
           		tempfile grp1
           		qui save `grp1',replace
@@ -856,13 +856,13 @@ local Nn = 200
 		  	di "replication = `replication'"
 		    }
               mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
-          		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+          		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) 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(7) pcm(D) clear
+          		qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(7) pcm(D) clear
           		qui gen TT = 1
           		tempfile grp1
           		qui save `grp1',replace
@@ -899,13 +899,13 @@ local Nn = 200
 		  	di "replication = `replication'"
 		    }
             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.15,1.15,2.15)
-            qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+            qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
             qui gen TT = 0
             tempfile grp0
             qui save `grp0',replace
 
             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.15,1.15,2.15)
-            qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+            qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
             qui gen TT = 1
             tempfile grp1
             qui save `grp1',replace
@@ -937,13 +937,13 @@ local Nn = 200
 		  	di "replication = `replication'"
 		    }
             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.15,1.15,2.15)
-            qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+            qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
             qui gen TT = 0
             tempfile grp0
             qui save `grp0',replace
 
             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.15,1.15,2.15)
-            qui simirt, nbobs(Nn) mu(0.2) cov(1) dim(7) pcm(D) clear
+            qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(7) pcm(D) clear
             qui gen TT = 1
             tempfile grp1
             qui save `grp1',replace
@@ -974,13 +974,13 @@ local Nn = 200
 		  	di "replication = `replication'"
 		    }
             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.15,1.15,2.15)
-            qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+            qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
             qui gen TT = 0
             tempfile grp0
             qui save `grp0',replace
 
             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.15,1.15,2.15)
-            qui simirt, nbobs(Nn) mu(0.4) cov(1) dim(7) pcm(D) clear
+            qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(7) pcm(D) clear
             qui gen TT = 1
             tempfile grp1
             qui save `grp1',replace
@@ -1048,13 +1048,13 @@ local Nn = 300
 		  	di "replication = `replication'"
 		    }
             mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-            qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+            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= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-            qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+            qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
             qui gen TT = 1
             tempfile grp1
             qui save `grp1',replace
@@ -1086,13 +1086,13 @@ local Nn = 300
 		  	di "replication = `replication'"
 		    }
             mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-            qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+            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= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-            qui simirt, nbobs(Nn) mu(0.2) cov(1) dim(4) pcm(D) clear
+            qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear
             qui gen TT = 1
             tempfile grp1
             qui save `grp1',replace
@@ -1123,13 +1123,13 @@ local Nn = 300
 		  	di "replication = `replication'"
 		    }
             mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-            qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+            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= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-            qui simirt, nbobs(Nn) mu(0.4) cov(1) dim(4) pcm(D) clear
+            qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear
             qui gen TT = 1
             tempfile grp1
             qui save `grp1',replace
@@ -1166,13 +1166,13 @@ local Nn = 300
 		  	di "replication = `replication'"
 		    }
                             		mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-                            		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+                            		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.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-                            		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+                            		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
                             		qui gen TT = 1
                             		tempfile grp1
                             		qui save `grp1',replace
@@ -1204,13 +1204,13 @@ local Nn = 300
 															  	di "replication = `replication'"
 															    }
                                   mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-                              		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+                              		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.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-                              		qui simirt, nbobs(Nn) mu(0.2) cov(1) dim(4) pcm(D) clear
+                              		qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear
                               		qui gen TT = 1
                               		tempfile grp1
                               		qui save `grp1',replace
@@ -1241,13 +1241,13 @@ local Nn = 300
 		  	di "replication = `replication'"
 		    }
                                   mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-                              		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+                              		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.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-                              		qui simirt, nbobs(Nn) mu(0.4) cov(1) dim(4) pcm(D) clear
+                              		qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear
                               		qui gen TT = 1
                               		tempfile grp1
                               		qui save `grp1',replace
@@ -1284,13 +1284,13 @@ local Nn = 300
 		  	di "replication = `replication'"
 		    }
                         		mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
-                        		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+                        		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) 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) cov(1) dim(7) pcm(D) clear
+                        		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
                         		qui gen TT = 1
                         		tempfile grp1
                         		qui save `grp1',replace
@@ -1322,13 +1322,13 @@ local Nn = 300
 		  	di "replication = `replication'"
 		    }
                               mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
-                          		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+                          		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) 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(7) pcm(D) clear
+                          		qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(7) pcm(D) clear
                           		qui gen TT = 1
                           		tempfile grp1
                           		qui save `grp1',replace
@@ -1359,13 +1359,13 @@ local Nn = 300
 		  	di "replication = `replication'"
 		    }
                               mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
-                          		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+                          		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) 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(7) pcm(D) clear
+                          		qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(7) pcm(D) clear
                           		qui gen TT = 1
                           		tempfile grp1
                           		qui save `grp1',replace
@@ -1402,13 +1402,13 @@ local Nn = 300
 		  	di "replication = `replication'"
 		    }
                             		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.15,1.15,2.15)
-                            		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+                            		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
                             		qui gen TT = 0
                             		tempfile grp0
                             		qui save `grp0',replace
 
                                 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.15,1.15,2.15)
-                            		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+                            		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
                             		qui gen TT = 1
                             		tempfile grp1
                             		qui save `grp1',replace
@@ -1440,13 +1440,13 @@ local Nn = 300
 		  	di "replication = `replication'"
 		    }
                                   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.15,1.15,2.15)
-                              		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+                              		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
                               		qui gen TT = 0
                               		tempfile grp0
                               		qui save `grp0',replace
 
                                   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.15,1.15,2.15)
-                              		qui simirt, nbobs(Nn) mu(0.2) cov(1) dim(7) pcm(D) clear
+                              		qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(7) pcm(D) clear
                               		qui gen TT = 1
                               		tempfile grp1
                               		qui save `grp1',replace
@@ -1477,13 +1477,13 @@ local Nn = 300
 		  	di "replication = `replication'"
 		    }
                                   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.15,1.15,2.15)
-                              		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+                              		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
                               		qui gen TT = 0
                               		tempfile grp0
                               		qui save `grp0',replace
 
                                   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.15,1.15,2.15)
-                              		qui simirt, nbobs(Nn) mu(0.4) cov(1) dim(7) pcm(D) clear
+                              		qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(7) pcm(D) clear
                               		qui gen TT = 1
                               		tempfile grp1
                               		qui save `grp1',replace
diff --git a/Scripts/Scenarios/NoDIF/scenarios_noDIF_baseline_100.do b/Scripts/Scenarios/NoDIF/scenarios_noDIF_baseline_100.do
index 6ccb475..39da9b3 100644
--- a/Scripts/Scenarios/NoDIF/scenarios_noDIF_baseline_100.do
+++ b/Scripts/Scenarios/NoDIF/scenarios_noDIF_baseline_100.do
@@ -40,13 +40,13 @@ forvalues replication = 1/1000 {
     }
 	  di "replication = `replication'"
 		mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+		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= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
 		qui gen TT = 1
 		tempfile grp1
 		qui save `grp1',replace
@@ -79,13 +79,13 @@ forvalues replication = 1/1000 {
 	  	di "replication = `replication'"
 	    }
   		mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-  		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+  		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= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-  		qui simirt, nbobs(Nn) mu(0.2) cov(1) dim(4) pcm(D) clear
+  		qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear
   		qui gen TT = 1
   		tempfile grp1
   		qui save `grp1',replace
@@ -117,13 +117,13 @@ forvalues replication = 1/1000 {
   	di "replication = `replication'"
     }
   		mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-  		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+  		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= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-  		qui simirt, nbobs(Nn) mu(0.4) cov(1) dim(4) pcm(D) clear
+  		qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear
   		qui gen TT = 1
   		tempfile grp1
   		qui save `grp1',replace
@@ -160,13 +160,13 @@ forvalues replication = 1/1000 {
 		  	di "replication = `replication'"
 		    }
     		mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-    		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+    		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.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-    		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+    		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
     		qui gen TT = 1
     		tempfile grp1
     		qui save `grp1',replace
@@ -198,13 +198,13 @@ forvalues replication = 1/1000 {
 			  	di "replication = `replication'"
 			    }
           mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-      		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+      		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.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-      		qui simirt, nbobs(Nn) mu(0.2) cov(1) dim(4) pcm(D) clear
+      		qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear
       		qui gen TT = 1
       		tempfile grp1
       		qui save `grp1',replace
@@ -235,13 +235,13 @@ forvalues replication = 1/1000 {
 		  	di "replication = `replication'"
 		    }
           mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-      		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+      		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.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-      		qui simirt, nbobs(Nn) mu(0.4) cov(1) dim(4) pcm(D) clear
+      		qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear
       		qui gen TT = 1
       		tempfile grp1
       		qui save `grp1',replace
@@ -278,13 +278,13 @@ forvalues replication = 1/1000 {
   	di "replication = `replication'"
     }
 		mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
-		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) 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) cov(1) dim(7) pcm(D) clear
+		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
 		qui gen TT = 1
 		tempfile grp1
 		qui save `grp1',replace
@@ -316,13 +316,13 @@ forvalues replication = 1/1000 {
 	  	di "replication = `replication'"
 	    }
       mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
-  		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+  		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) 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(7) pcm(D) clear
+  		qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(7) pcm(D) clear
   		qui gen TT = 1
   		tempfile grp1
   		qui save `grp1',replace
@@ -353,13 +353,13 @@ forvalues replication = 1/1000 {
   	di "replication = `replication'"
     }
       mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
-  		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+  		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) 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(7) pcm(D) clear
+  		qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(7) pcm(D) clear
   		qui gen TT = 1
   		tempfile grp1
   		qui save `grp1',replace
@@ -396,13 +396,13 @@ forvalues replication = 1/1000 {
 		  	di "replication = `replication'"
 		    }
     		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.15,1.15,2.15)
-    		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+    		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
     		qui gen TT = 0
     		tempfile grp0
     		qui save `grp0',replace
 
         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.15,1.15,2.15)
-    		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+    		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
     		qui gen TT = 1
     		tempfile grp1
     		qui save `grp1',replace
@@ -434,13 +434,13 @@ forvalues replication = 1/1000 {
 		  	di "replication = `replication'"
 		    }
           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.15,1.15,2.15)
-      		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+      		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
       		qui gen TT = 0
       		tempfile grp0
       		qui save `grp0',replace
 
           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.15,1.15,2.15)
-      		qui simirt, nbobs(Nn) mu(0.2) cov(1) dim(7) pcm(D) clear
+      		qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(7) pcm(D) clear
       		qui gen TT = 1
       		tempfile grp1
       		qui save `grp1',replace
@@ -471,13 +471,13 @@ forvalues replication = 1/1000 {
 		  	di "replication = `replication'"
 		    }
           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.15,1.15,2.15)
-      		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+      		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
       		qui gen TT = 0
       		tempfile grp0
       		qui save `grp0',replace
 
           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.15,1.15,2.15)
-      		qui simirt, nbobs(Nn) mu(0.4) cov(1) dim(7) pcm(D) clear
+      		qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(7) pcm(D) clear
       		qui gen TT = 1
       		tempfile grp1
       		qui save `grp1',replace
diff --git a/Scripts/Scenarios/NoDIF/scenarios_noDIF_baseline_200.do b/Scripts/Scenarios/NoDIF/scenarios_noDIF_baseline_200.do
index 00f0f78..ac78a2b 100644
--- a/Scripts/Scenarios/NoDIF/scenarios_noDIF_baseline_200.do
+++ b/Scripts/Scenarios/NoDIF/scenarios_noDIF_baseline_200.do
@@ -35,13 +35,13 @@ local Nn = 200
 		  	di "replication = `replication'"
 		    }
         		mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-        		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+        		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= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-        		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+        		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
         		qui gen TT = 1
         		tempfile grp1
         		qui save `grp1',replace
@@ -73,13 +73,13 @@ local Nn = 200
 		  	di "replication = `replication'"
 		    }
           		mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-          		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+          		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= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-          		qui simirt, nbobs(Nn) mu(0.2) cov(1) dim(4) pcm(D) clear
+          		qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear
           		qui gen TT = 1
           		tempfile grp1
           		qui save `grp1',replace
@@ -110,13 +110,13 @@ local Nn = 200
 		  	di "replication = `replication'"
 		    }
           		mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-          		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+          		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= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-          		qui simirt, nbobs(Nn) mu(0.4) cov(1) dim(4) pcm(D) clear
+          		qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear
           		qui gen TT = 1
           		tempfile grp1
           		qui save `grp1',replace
@@ -153,13 +153,13 @@ local Nn = 200
 		  	di "replication = `replication'"
 		    }
             		mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-            		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+            		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.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-            		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+            		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
             		qui gen TT = 1
             		tempfile grp1
             		qui save `grp1',replace
@@ -191,13 +191,13 @@ local Nn = 200
 		  	di "replication = `replication'"
 		    }
                   mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-              		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+              		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.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-              		qui simirt, nbobs(Nn) mu(0.2) cov(1) dim(4) pcm(D) clear
+              		qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear
               		qui gen TT = 1
               		tempfile grp1
               		qui save `grp1',replace
@@ -228,13 +228,13 @@ local Nn = 200
 		  	di "replication = `replication'"
 		    }
                   mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-              		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+              		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.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-              		qui simirt, nbobs(Nn) mu(0.4) cov(1) dim(4) pcm(D) clear
+              		qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear
               		qui gen TT = 1
               		tempfile grp1
               		qui save `grp1',replace
@@ -271,13 +271,13 @@ local Nn = 200
 		  	di "replication = `replication'"
 		    }
         		mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
-        		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+        		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) 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) cov(1) dim(7) pcm(D) clear
+        		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
         		qui gen TT = 1
         		tempfile grp1
         		qui save `grp1',replace
@@ -309,13 +309,13 @@ local Nn = 200
 		  	di "replication = `replication'"
 		    }
               mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
-          		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+          		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) 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(7) pcm(D) clear
+          		qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(7) pcm(D) clear
           		qui gen TT = 1
           		tempfile grp1
           		qui save `grp1',replace
@@ -346,13 +346,13 @@ local Nn = 200
 		  	di "replication = `replication'"
 		    }
               mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
-          		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+          		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) 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(7) pcm(D) clear
+          		qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(7) pcm(D) clear
           		qui gen TT = 1
           		tempfile grp1
           		qui save `grp1',replace
@@ -389,13 +389,13 @@ local Nn = 200
 		  	di "replication = `replication'"
 		    }
             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.15,1.15,2.15)
-            qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+            qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
             qui gen TT = 0
             tempfile grp0
             qui save `grp0',replace
 
             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.15,1.15,2.15)
-            qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+            qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
             qui gen TT = 1
             tempfile grp1
             qui save `grp1',replace
@@ -427,13 +427,13 @@ local Nn = 200
 		  	di "replication = `replication'"
 		    }
             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.15,1.15,2.15)
-            qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+            qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
             qui gen TT = 0
             tempfile grp0
             qui save `grp0',replace
 
             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.15,1.15,2.15)
-            qui simirt, nbobs(Nn) mu(0.2) cov(1) dim(7) pcm(D) clear
+            qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(7) pcm(D) clear
             qui gen TT = 1
             tempfile grp1
             qui save `grp1',replace
@@ -464,13 +464,13 @@ local Nn = 200
 		  	di "replication = `replication'"
 		    }
             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.15,1.15,2.15)
-            qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+            qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
             qui gen TT = 0
             tempfile grp0
             qui save `grp0',replace
 
             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.15,1.15,2.15)
-            qui simirt, nbobs(Nn) mu(0.4) cov(1) dim(7) pcm(D) clear
+            qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(7) pcm(D) clear
             qui gen TT = 1
             tempfile grp1
             qui save `grp1',replace
diff --git a/Scripts/Scenarios/NoDIF/scenarios_noDIF_baseline_300.do b/Scripts/Scenarios/NoDIF/scenarios_noDIF_baseline_300.do
index 3cab9dd..612f506 100644
--- a/Scripts/Scenarios/NoDIF/scenarios_noDIF_baseline_300.do
+++ b/Scripts/Scenarios/NoDIF/scenarios_noDIF_baseline_300.do
@@ -36,13 +36,13 @@ local Nn = 300
 		  	di "replication = `replication'"
 		    }
             mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-            qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+            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= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-            qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+            qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
             qui gen TT = 1
             tempfile grp1
             qui save `grp1',replace
@@ -74,13 +74,13 @@ local Nn = 300
 		  	di "replication = `replication'"
 		    }
             mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-            qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+            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= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-            qui simirt, nbobs(Nn) mu(0.2) cov(1) dim(4) pcm(D) clear
+            qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear
             qui gen TT = 1
             tempfile grp1
             qui save `grp1',replace
@@ -111,13 +111,13 @@ local Nn = 300
 		  	di "replication = `replication'"
 		    }
             mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-            qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+            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= (-0.84 \ -0.25 \ 0.25 \ 0.84)
-            qui simirt, nbobs(Nn) mu(0.4) cov(1) dim(4) pcm(D) clear
+            qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear
             qui gen TT = 1
             tempfile grp1
             qui save `grp1',replace
@@ -154,13 +154,13 @@ local Nn = 300
 		  	di "replication = `replication'"
 		    }
                             		mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-                            		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+                            		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.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-                            		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+                            		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(4) pcm(D) clear
                             		qui gen TT = 1
                             		tempfile grp1
                             		qui save `grp1',replace
@@ -192,13 +192,13 @@ local Nn = 300
 															  	di "replication = `replication'"
 															    }
                                   mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-                              		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+                              		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.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-                              		qui simirt, nbobs(Nn) mu(0.2) cov(1) dim(4) pcm(D) clear
+                              		qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(4) pcm(D) clear
                               		qui gen TT = 1
                               		tempfile grp1
                               		qui save `grp1',replace
@@ -229,13 +229,13 @@ local Nn = 300
 		  	di "replication = `replication'"
 		    }
                                   mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-                              		qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
+                              		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.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
-                              		qui simirt, nbobs(Nn) mu(0.4) cov(1) dim(4) pcm(D) clear
+                              		qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(4) pcm(D) clear
                               		qui gen TT = 1
                               		tempfile grp1
                               		qui save `grp1',replace
@@ -272,13 +272,13 @@ local Nn = 300
 		  	di "replication = `replication'"
 		    }
                         		mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
-                        		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+                        		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) 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) cov(1) dim(7) pcm(D) clear
+                        		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
                         		qui gen TT = 1
                         		tempfile grp1
                         		qui save `grp1',replace
@@ -310,13 +310,13 @@ local Nn = 300
 		  	di "replication = `replication'"
 		    }
                               mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
-                          		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+                          		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) 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(7) pcm(D) clear
+                          		qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(7) pcm(D) clear
                           		qui gen TT = 1
                           		tempfile grp1
                           		qui save `grp1',replace
@@ -347,13 +347,13 @@ local Nn = 300
 		  	di "replication = `replication'"
 		    }
                               mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
-                          		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+                          		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) 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(7) pcm(D) clear
+                          		qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(7) pcm(D) clear
                           		qui gen TT = 1
                           		tempfile grp1
                           		qui save `grp1',replace
@@ -390,13 +390,13 @@ local Nn = 300
 		  	di "replication = `replication'"
 		    }
                             		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.15,1.15,2.15)
-                            		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+                            		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
                             		qui gen TT = 0
                             		tempfile grp0
                             		qui save `grp0',replace
 
                                 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.15,1.15,2.15)
-                            		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+                            		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
                             		qui gen TT = 1
                             		tempfile grp1
                             		qui save `grp1',replace
@@ -428,13 +428,13 @@ local Nn = 300
 		  	di "replication = `replication'"
 		    }
                                   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.15,1.15,2.15)
-                              		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+                              		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
                               		qui gen TT = 0
                               		tempfile grp0
                               		qui save `grp0',replace
 
                                   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.15,1.15,2.15)
-                              		qui simirt, nbobs(Nn) mu(0.2) cov(1) dim(7) pcm(D) clear
+                              		qui simirt, nbobs(`Nn') mu(0.2) cov(1) dim(7) pcm(D) clear
                               		qui gen TT = 1
                               		tempfile grp1
                               		qui save `grp1',replace
@@ -465,13 +465,13 @@ local Nn = 300
 		  	di "replication = `replication'"
 		    }
                                   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.15,1.15,2.15)
-                              		qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
+                              		qui simirt, nbobs(`Nn') mu(0) cov(1) dim(7) pcm(D) clear
                               		qui gen TT = 0
                               		tempfile grp0
                               		qui save `grp0',replace
 
                                   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.15,1.15,2.15)
-                              		qui simirt, nbobs(Nn) mu(0.4) cov(1) dim(7) pcm(D) clear
+                              		qui simirt, nbobs(`Nn') mu(0.4) cov(1) dim(7) pcm(D) clear
                               		qui gen TT = 1
                               		tempfile grp1
                               		qui save `grp1',replace