corentinchoisy/simul_these
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Updated display in noDIF scenario code + separated code intro N=100, 200 and 300 files

Browse Source
This commit is contained in:
Corentin Choisy
2024-01-05 10:38:50 +01:00
parent c0cdde3a03
commit a3ffdbb46c
5 changed files with 1672 additions and 271 deletions
Download Patch File Download Diff File Expand all files Collapse all files

252
Scripts/Scenarios/NoDIF/scenarios_noDIF_baseline.do
View File

@ -32,8 +32,13 @@ local Nn = 100
* Scenario 1A : H_0 is TRUE
di "SCENARIO 1A - N=100"
forvalues replication = 1/1000 {
di "replication = `replication'"
if mod(`replication',10)==0 {
di "replication = `replication'"
}
di "replication = `replication'"
mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
@ -67,8 +72,12 @@ forvalues replication = 1/1000 {
* Scenario 1B : H_0 is FALSE / Effect size = 0.2
forvalues replication = 1/1000 {
di "replication = `replication'"
di "SCENARIO 1B - N=100"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
@ -101,8 +110,12 @@ forvalues replication = 1/1000 {
* Scenario 1C : H_0 is FALSE / Effect size = 0.4
forvalues replication = 1/1000 {
di "replication = `replication'"
di "SCENARIO 1C - N=100"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
@ -140,9 +153,12 @@ forvalues replication = 1/1000 {
** Scenario 2: J = 4 items / M = 4 modalities / N=100 per group / TT=treatment variable
* Scenario 2A : H_0 is TRUE
di "SCENARIO 2A - N=100"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
@ -175,9 +191,12 @@ forvalues replication = 1/1000 {
* Scenario 2B : H_0 is FALSE / Effect size = 0.2
di "SCENARIO 2B - N=100"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
@ -209,9 +228,12 @@ forvalues replication = 1/1000 {
export delimited using "`path'/scenario_2B_100.csv", replace
* Scenario 2C : H_0 is FALSE / Effect size = 0.4
di "SCENARIO 2C - N=100"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
@ -249,9 +271,12 @@ forvalues replication = 1/1000 {
** Scenario 3: J = 7 items / M = 2 modalities / N=100 per group / TT=treatment variable
* Scenario 3A : H_0 is TRUE
di "SCENARIO 3A - N=100"
forvalues replication = 1/1000 {
di "replication = `replication'"
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
@ -284,9 +309,12 @@ forvalues replication = 1/1000 {
* Scenario 3B : H_0 is FALSE / Effect size = 0.2
di "SCENARIO 3B - N=100"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
@ -318,9 +346,12 @@ forvalues replication = 1/1000 {
export delimited using "`path'/scenario_3B_100.csv", replace
* Scenario 3C : H_0 is FALSE / Effect size = 0.4
di "SCENARIO 3C - N=100"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
@ -358,9 +389,12 @@ forvalues replication = 1/1000 {
** Scenario 4: J = 7 items / M = 4 modalities / N=100 per group / TT=treatment variable
* Scenario 4A : H_0 is TRUE
di "SCENARIO 4A - N=100"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
@ -393,9 +427,12 @@ forvalues replication = 1/1000 {
* Scenario 4B : H_0 is FALSE / Effect size = 0.2
di "SCENARIO 4B - N=100"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
@ -427,9 +464,12 @@ forvalues replication = 1/1000 {
export delimited using "`path'/scenario_4B_100.csv", replace
* Scenario 4C : H_0 is FALSE / Effect size = 0.4
di "SCENARIO 4C - N=100"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
@ -498,9 +538,12 @@ local Nn = 200
** Scenario 1: J = 4 items / M = 2 modalities / N=100 per group / TT=treatment variable
* Scenario 1A : H_0 is TRUE
di "SCENARIO 1A - N=200"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
@ -533,9 +576,12 @@ local Nn = 200
* Scenario 1B : H_0 is FALSE / Effect size = 0.2
di "SCENARIO 1B - N=200"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
@ -567,9 +613,12 @@ local Nn = 200
export delimited using "`path'/scenario_1B_200.csv", replace
* Scenario 1C : H_0 is FALSE / Effect size = 0.4
di "SCENARIO 1C - N=200"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
@ -607,9 +656,12 @@ local Nn = 200
** Scenario 2: J = 4 items / M = 4 modalities / N=100 per group / TT=treatment variable
* Scenario 2A : H_0 is TRUE
di "SCENARIO 2A - N=200"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
@ -642,9 +694,12 @@ local Nn = 200
* Scenario 2B : H_0 is FALSE / Effect size = 0.2
di "SCENARIO 2B - N=200"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
@ -676,9 +731,12 @@ local Nn = 200
export delimited using "`path'/scenario_2B_200.csv", replace
* Scenario 2C : H_0 is FALSE / Effect size = 0.4
di "SCENARIO 2C - N=200"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
@ -716,9 +774,12 @@ local Nn = 200
** Scenario 3: J = 7 items / M = 2 modalities / N=100 per group / TT=treatment variable
* Scenario 3A : H_0 is TRUE
di "SCENARIO 3A - N=200"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
@ -751,9 +812,12 @@ local Nn = 200
* Scenario 3B : H_0 is FALSE / Effect size = 0.2
di "SCENARIO 3B - N=200"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
@ -785,9 +849,12 @@ local Nn = 200
export delimited using "`path'/scenario_3B_200.csv", replace
* Scenario 3C : H_0 is FALSE / Effect size = 0.4
di "SCENARIO 3C - N=200"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
@ -825,9 +892,12 @@ local Nn = 200
** Scenario 4: J = 7 items / M = 4 modalities / N=100 per group / TT=treatment variable
* Scenario 4A : H_0 is TRUE
di "SCENARIO 4A - N=200"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
@ -860,9 +930,12 @@ local Nn = 200
* Scenario 4B : H_0 is FALSE / Effect size = 0.2
di "SCENARIO 4B - N=200"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
@ -894,9 +967,12 @@ local Nn = 200
export delimited using "`path'/scenario_4B_200.csv", replace
* Scenario 4C : H_0 is FALSE / Effect size = 0.4
di "SCENARIO 4C - N=200"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
@ -965,9 +1041,12 @@ local Nn = 300
** Scenario 1: J = 4 items / M = 2 modalities / N=100 per group / TT=treatment variable
* Scenario 1A : H_0 is TRUE
di "SCENARIO 1A - N=300"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
@ -1000,9 +1079,12 @@ local Nn = 300
* Scenario 1B : H_0 is FALSE / Effect size = 0.2
di "SCENARIO 1B - N=300"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
@ -1034,9 +1116,12 @@ local Nn = 300
export delimited using "`path'/scenario_1B_300.csv", replace
* Scenario 1C : H_0 is FALSE / Effect size = 0.4
di "SCENARIO 1C - N=300"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
@ -1074,9 +1159,12 @@ local Nn = 300
** Scenario 2: J = 4 items / M = 4 modalities / N=100 per group / TT=treatment variable
* Scenario 2A : H_0 is TRUE
di "SCENARIO 2A - N=300"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
@ -1109,9 +1197,12 @@ local Nn = 300
* Scenario 2B : H_0 is FALSE / Effect size = 0.2
di "SCENARIO 2B - N=300"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
@ -1143,9 +1234,12 @@ local Nn = 300
export delimited using "`path'/scenario_2B_300.csv", replace
* Scenario 2C : H_0 is FALSE / Effect size = 0.4
di "SCENARIO 2C - N=300"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
@ -1183,9 +1277,12 @@ local Nn = 300
** Scenario 3: J = 7 items / M = 2 modalities / N=100 per group / TT=treatment variable
* Scenario 3A : H_0 is TRUE
di "SCENARIO 3A - N=300"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
@ -1218,9 +1315,12 @@ local Nn = 300
* Scenario 3B : H_0 is FALSE / Effect size = 0.2
di "SCENARIO 3B - N=300"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
@ -1252,9 +1352,12 @@ local Nn = 300
export delimited using "`path'/scenario_3B_300.csv", replace
* Scenario 3C : H_0 is FALSE / Effect size = 0.4
di "SCENARIO 3C - N=300"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
@ -1292,9 +1395,12 @@ local Nn = 300
** Scenario 4: J = 7 items / M = 4 modalities / N=100 per group / TT=treatment variable
* Scenario 4A : H_0 is TRUE
di "SCENARIO 4A - N=300"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
@ -1327,9 +1433,12 @@ local Nn = 300
* Scenario 4B : H_0 is FALSE / Effect size = 0.2
di "SCENARIO 4B - N=300"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
@ -1361,9 +1470,12 @@ local Nn = 300
export delimited using "`path'/scenario_4B_300.csv", replace
* Scenario 4C : H_0 is FALSE / Effect size = 0.4
di "SCENARIO 4C - N=300"
forvalues replication = 1/1000 {
di "replication = `replication'"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0

501
Scripts/Scenarios/NoDIF/scenarios_noDIF_baseline_100.do Normal file
View File

@ -0,0 +1,501 @@
*=================================================================================================================================================
* Date : 2024-01-04
* Stata version : Stata 18 SE
*
* This program creates dataset without DIF for a randomized controlled trial scenario
*
* ado-files needed : - simirt (version 4.3 August 29, 2019, available on OSF)
*
* outputs : scenario_1,scenario_2,scenario_3,scenario_4, for N=100/200/300
*
*
* Warning : To obtain reproduce the data obtained in the .csv files in this repository, use 'simirt_setseed.ado' instead of 'simirt.ado'
*
*
*================================================================================================================================================
* Load simirt.ado
adopath+"/home/corentin/Documents/These/Recherche/Simulations/Modules/"
* Set data output folder path
local path = "/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF"
* Scenarios with : n = 100
*==========================
local path = "/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N100"
local Nn = 100
** Scenario 1: J = 4 items / M = 2 modalities / N=100 per group / TT=treatment variable
* Scenario 1A : H_0 is TRUE
di "SCENARIO 1A - N=100"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
di "replication = `replication'"
mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui 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_1A_100.csv", replace
* Scenario 1B : H_0 is FALSE / Effect size = 0.2
di "SCENARIO 1B - N=100"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
qui simirt, nbobs(Nn) mu(0.2) cov(1) dim(4) pcm(D) clear
qui 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_1B_100.csv", replace
* Scenario 1C : H_0 is FALSE / Effect size = 0.4
di "SCENARIO 1C - N=100"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
qui simirt, nbobs(Nn) mu(0.4) cov(1) dim(4) pcm(D) clear
qui 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_1C_100.csv", replace
**------------------------------------------------------------------------------------**
** Scenario 2: J = 4 items / M = 4 modalities / N=100 per group / TT=treatment variable
* Scenario 2A : H_0 is TRUE
di "SCENARIO 2A - N=100"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui 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_2A_100.csv", replace
* Scenario 2B : H_0 is FALSE / Effect size = 0.2
di "SCENARIO 2B - N=100"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
qui simirt, nbobs(Nn) mu(0.2) cov(1) dim(4) pcm(D) clear
qui 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_2B_100.csv", replace
* Scenario 2C : H_0 is FALSE / Effect size = 0.4
di "SCENARIO 2C - N=100"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
qui simirt, nbobs(Nn) mu(0.4) cov(1) dim(4) pcm(D) clear
qui 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_2C_100.csv", replace
**------------------------------------------------------------------------------------**
** Scenario 3: J = 7 items / M = 2 modalities / N=100 per group / TT=treatment variable
* Scenario 3A : H_0 is TRUE
di "SCENARIO 3A - N=100"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui 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_3A_100.csv", replace
* Scenario 3B : H_0 is FALSE / Effect size = 0.2
di "SCENARIO 3B - N=100"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
qui simirt, nbobs(Nn) mu(0.2) cov(1) dim(7) pcm(D) clear
qui 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_3B_100.csv", replace
* Scenario 3C : H_0 is FALSE / Effect size = 0.4
di "SCENARIO 3C - N=100"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
qui simirt, nbobs(Nn) mu(0.4) cov(1) dim(7) pcm(D) clear
qui 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_3C_100.csv", replace
**------------------------------------------------------------------------------------**
** Scenario 4: J = 7 items / M = 4 modalities / N=100 per group / TT=treatment variable
* Scenario 4A : H_0 is TRUE
di "SCENARIO 4A - N=100"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui 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_4A_100.csv", replace
* Scenario 4B : H_0 is FALSE / Effect size = 0.2
di "SCENARIO 4B - N=100"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
qui simirt, nbobs(Nn) mu(0.2) cov(1) dim(7) pcm(D) clear
qui 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_4B_100.csv", replace
* Scenario 4C : H_0 is FALSE / Effect size = 0.4
di "SCENARIO 4C - N=100"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
qui simirt, nbobs(Nn) mu(0.4) cov(1) dim(7) pcm(D) clear
qui 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_4C_100.csv", replace

494
Scripts/Scenarios/NoDIF/scenarios_noDIF_baseline_200.do Normal file
View File

@ -0,0 +1,494 @@
*=================================================================================================================================================
* Date : 2024-01-04
* Stata version : Stata 18 SE
*
* This program creates dataset without DIF for a randomized controlled trial scenario
*
* ado-files needed : - simirt (version 4.3 August 29, 2019, available on OSF)
*
* outputs : scenario_1,scenario_2,scenario_3,scenario_4, for N=100/200/300
*
*
* Warning : To obtain reproduce the data obtained in the .csv files in this repository, use 'simirt_setseed.ado' instead of 'simirt.ado'
*
*
*================================================================================================================================================
* Load simirt.ado
adopath+"/home/corentin/Documents/These/Recherche/Simulations/Modules/"
* Set data output folder path
local path = "/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF"
*==========================
local path = "/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N200"
local Nn = 200
** Scenario 1: J = 4 items / M = 2 modalities / N=100 per group / TT=treatment variable
* Scenario 1A : H_0 is TRUE
di "SCENARIO 1A - N=200"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui 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_1A_200.csv", replace
* Scenario 1B : H_0 is FALSE / Effect size = 0.2
di "SCENARIO 1B - N=200"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
qui simirt, nbobs(Nn) mu(0.2) cov(1) dim(4) pcm(D) clear
qui 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_1B_200.csv", replace
* Scenario 1C : H_0 is FALSE / Effect size = 0.4
di "SCENARIO 1C - N=200"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
qui simirt, nbobs(Nn) mu(0.4) cov(1) dim(4) pcm(D) clear
qui 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_1C_200.csv", replace
**------------------------------------------------------------------------------------**
** Scenario 2: J = 4 items / M = 4 modalities / N=100 per group / TT=treatment variable
* Scenario 2A : H_0 is TRUE
di "SCENARIO 2A - N=200"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui 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_2A_200.csv", replace
* Scenario 2B : H_0 is FALSE / Effect size = 0.2
di "SCENARIO 2B - N=200"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
qui simirt, nbobs(Nn) mu(0.2) cov(1) dim(4) pcm(D) clear
qui 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_2B_200.csv", replace
* Scenario 2C : H_0 is FALSE / Effect size = 0.4
di "SCENARIO 2C - N=200"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
qui simirt, nbobs(Nn) mu(0.4) cov(1) dim(4) pcm(D) clear
qui 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_2C_200.csv", replace
**------------------------------------------------------------------------------------**
** Scenario 3: J = 7 items / M = 2 modalities / N=100 per group / TT=treatment variable
* Scenario 3A : H_0 is TRUE
di "SCENARIO 3A - N=200"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui 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_3A_200.csv", replace
* Scenario 3B : H_0 is FALSE / Effect size = 0.2
di "SCENARIO 3B - N=200"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
qui simirt, nbobs(Nn) mu(0.2) cov(1) dim(7) pcm(D) clear
qui 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_3B_200.csv", replace
* Scenario 3C : H_0 is FALSE / Effect size = 0.4
di "SCENARIO 3C - N=200"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
qui simirt, nbobs(Nn) mu(0.4) cov(1) dim(7) pcm(D) clear
qui 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_3C_200.csv", replace
**------------------------------------------------------------------------------------**
** Scenario 4: J = 7 items / M = 4 modalities / N=100 per group / TT=treatment variable
* Scenario 4A : H_0 is TRUE
di "SCENARIO 4A - N=200"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui 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_4A_200.csv", replace
* Scenario 4B : H_0 is FALSE / Effect size = 0.2
di "SCENARIO 4B - N=200"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
qui simirt, nbobs(Nn) mu(0.2) cov(1) dim(7) pcm(D) clear
qui 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_4B_200.csv", replace
* Scenario 4C : H_0 is FALSE / Effect size = 0.4
di "SCENARIO 4C - N=200"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
qui simirt, nbobs(Nn) mu(0.4) cov(1) dim(7) pcm(D) clear
qui 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_4C_200.csv", replace

495
Scripts/Scenarios/NoDIF/scenarios_noDIF_baseline_300.do Normal file
View File

@ -0,0 +1,495 @@
*=================================================================================================================================================
* Date : 2024-01-04
* Stata version : Stata 18 SE
*
* This program creates dataset without DIF for a randomized controlled trial scenario
*
* ado-files needed : - simirt (version 4.3 August 29, 2019, available on OSF)
*
* outputs : scenario_1,scenario_2,scenario_3,scenario_4, for N=100/200/300
*
*
* Warning : To obtain reproduce the data obtained in the .csv files in this repository, use 'simirt_setseed.ado' instead of 'simirt.ado'
*
*
*================================================================================================================================================
* Load simirt.ado
adopath+"/home/corentin/Documents/These/Recherche/Simulations/Modules/"
* Set data output folder path
local path = "/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF"
* Scenarios with : n = 300
*==========================
local path = "/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N300"
local Nn = 300
** Scenario 1: J = 4 items / M = 2 modalities / N=100 per group / TT=treatment variable
* Scenario 1A : H_0 is TRUE
di "SCENARIO 1A - N=300"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui 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_1A_300.csv", replace
* Scenario 1B : H_0 is FALSE / Effect size = 0.2
di "SCENARIO 1B - N=300"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
qui simirt, nbobs(Nn) mu(0.2) cov(1) dim(4) pcm(D) clear
qui 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_1B_300.csv", replace
* Scenario 1C : H_0 is FALSE / Effect size = 0.4
di "SCENARIO 1C - N=300"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-0.84 \ -0.25 \ 0.25 \ 0.84)
qui simirt, nbobs(Nn) mu(0.4) cov(1) dim(4) pcm(D) clear
qui 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_1C_300.csv", replace
**------------------------------------------------------------------------------------**
** Scenario 2: J = 4 items / M = 4 modalities / N=100 per group / TT=treatment variable
* Scenario 2A : H_0 is TRUE
di "SCENARIO 2A - N=300"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui 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_2A_300.csv", replace
* Scenario 2B : H_0 is FALSE / Effect size = 0.2
di "SCENARIO 2B - N=300"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
qui simirt, nbobs(Nn) mu(0.2) cov(1) dim(4) pcm(D) clear
qui 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_2B_300.csv", replace
* Scenario 2C : H_0 is FALSE / Effect size = 0.4
di "SCENARIO 2C - N=300"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(4) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
qui simirt, nbobs(Nn) mu(0.4) cov(1) dim(4) pcm(D) clear
qui 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_2C_300.csv", replace
**------------------------------------------------------------------------------------**
** Scenario 3: J = 7 items / M = 2 modalities / N=100 per group / TT=treatment variable
* Scenario 3A : H_0 is TRUE
di "SCENARIO 3A - N=300"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui 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_3A_300.csv", replace
* Scenario 3B : H_0 is FALSE / Effect size = 0.2
di "SCENARIO 3B - N=300"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
qui simirt, nbobs(Nn) mu(0.2) cov(1) dim(7) pcm(D) clear
qui 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_3B_300.csv", replace
* Scenario 3C : H_0 is FALSE / Effect size = 0.4
di "SCENARIO 3C - N=300"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-1.15 \ -0.67 \ -0.32 \ 0 \ 0.32 \ 0.67 \ 1.15)
qui simirt, nbobs(Nn) mu(0.4) cov(1) dim(7) pcm(D) clear
qui 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_3C_300.csv", replace
**------------------------------------------------------------------------------------**
** Scenario 4: J = 7 items / M = 4 modalities / N=100 per group / TT=treatment variable
* Scenario 4A : H_0 is TRUE
di "SCENARIO 4A - N=300"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui 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_4A_300.csv", replace
* Scenario 4B : H_0 is FALSE / Effect size = 0.2
di "SCENARIO 4B - N=300"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
qui simirt, nbobs(Nn) mu(0.2) cov(1) dim(7) pcm(D) clear
qui 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_4B_300.csv", replace
* Scenario 4C : H_0 is FALSE / Effect size = 0.4
di "SCENARIO 4C - N=300"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
qui simirt, nbobs(Nn) mu(0) cov(1) dim(7) pcm(D) clear
qui gen TT = 0
tempfile grp0
qui save `grp0',replace
mat D= (-2.15,-1.15,-0.15 \ -1.67,-0.67,0.33 \ -1.32,-0.32,0.68 \ -1,0,1 \ -0.68,0.32,1.32 \ -0.33,0.67,1.67 \ 0.15,1.15,2.15)
qui simirt, nbobs(Nn) mu(0.4) cov(1) dim(7) pcm(D) clear
qui 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_4C_300.csv", replace