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)