|
|
|
@ -6,7 +6,7 @@
|
|
|
|
|
*
|
|
|
|
|
* ado-files needed : - simirt (version 4.3 August 29, 2019, available on OSF)
|
|
|
|
|
*
|
|
|
|
|
* outputs : scenario_1,scenario_2,scenario_3,scenario_4, for N=200/200/300
|
|
|
|
|
* outputs : scenario_1,scenario_2,scenario_3,scenario_5, for N=200/200/300
|
|
|
|
|
*
|
|
|
|
|
*
|
|
|
|
|
* Warning : To obtain reproduce the data obtained in the .csv files in this repository, use 'simirt_setseed.ado' instead of 'simirt.ado'
|
|
|
|
@ -19,7 +19,7 @@ adopath+"/home/corentin/Documents/These/Recherche/Simulations/Modules/"
|
|
|
|
|
|
|
|
|
|
* Set data output folder path
|
|
|
|
|
|
|
|
|
|
local path = "/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N200"
|
|
|
|
|
local path = "/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N100"
|
|
|
|
|
local Nn = 200
|
|
|
|
|
|
|
|
|
|
|
|
|
|
@ -33,11 +33,11 @@ local Nn = 200
|
|
|
|
|
|
|
|
|
|
**** Scenarios with : DIF size 0.3 ****
|
|
|
|
|
|
|
|
|
|
** Scenario 4: J = 4 items / M = 2 modalities / DIF SIZE = 0.3
|
|
|
|
|
** Scenario 5: J = 4 items / M = 2 modalities / DIF SIZE = 0.3
|
|
|
|
|
|
|
|
|
|
* Scenario 4A : H_0 is TRUE / DIF on treatment
|
|
|
|
|
* Scenario 5A : H_0 is TRUE / DIF on treatment
|
|
|
|
|
|
|
|
|
|
di "SCENARIO 4A - N=200"
|
|
|
|
|
di "Scenario 5A - N=200"
|
|
|
|
|
|
|
|
|
|
forvalues replication = 1/1000 {
|
|
|
|
|
if mod(`replication',10)==0 {
|
|
|
|
@ -86,12 +86,12 @@ forvalues replication = 1/1000 {
|
|
|
|
|
qui save `data',replace
|
|
|
|
|
}
|
|
|
|
|
}
|
|
|
|
|
export delimited using "`path'/scenario_4A_200.csv", replace
|
|
|
|
|
export delimited using "`path'/scenario_5A_100.csv", replace
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
* Scenario 4B : H_1 is TRUE / Effect size 0.2 / DIF on treatment
|
|
|
|
|
* Scenario 5B : H_1 is TRUE / Effect size 0.2 / DIF on treatment
|
|
|
|
|
|
|
|
|
|
di "SCENARIO 4B - N=200"
|
|
|
|
|
di "Scenario 5B - N=200"
|
|
|
|
|
|
|
|
|
|
forvalues replication = 1/1000 {
|
|
|
|
|
if mod(`replication',10)==0 {
|
|
|
|
@ -140,11 +140,11 @@ forvalues replication = 1/1000 {
|
|
|
|
|
qui save `data',replace
|
|
|
|
|
}
|
|
|
|
|
}
|
|
|
|
|
export delimited using "`path'/scenario_4B_200.csv", replace
|
|
|
|
|
export delimited using "`path'/scenario_5B_100.csv", replace
|
|
|
|
|
|
|
|
|
|
* Scenario 4C : H_1 is TRUE / Effect size 0.2 / DIF on control
|
|
|
|
|
* Scenario 5C : H_1 is TRUE / Effect size 0.2 / DIF on control
|
|
|
|
|
|
|
|
|
|
di "SCENARIO 4C - N=200"
|
|
|
|
|
di "Scenario 5C - N=200"
|
|
|
|
|
|
|
|
|
|
forvalues replication = 1/1000 {
|
|
|
|
|
if mod(`replication',10)==0 {
|
|
|
|
@ -193,12 +193,12 @@ forvalues replication = 1/1000 {
|
|
|
|
|
qui save `data',replace
|
|
|
|
|
}
|
|
|
|
|
}
|
|
|
|
|
export delimited using "`path'/scenario_4C_200.csv", replace
|
|
|
|
|
export delimited using "`path'/scenario_5C_100.csv", replace
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
* Scenario 4D : H_1 is TRUE / Effect size 0.4 / DIF on treatment
|
|
|
|
|
* Scenario 5D : H_1 is TRUE / Effect size 0.4 / DIF on treatment
|
|
|
|
|
|
|
|
|
|
di "SCENARIO 4D - N=200"
|
|
|
|
|
di "Scenario 5D - N=200"
|
|
|
|
|
|
|
|
|
|
forvalues replication = 1/1000 {
|
|
|
|
|
if mod(`replication',10)==0 {
|
|
|
|
@ -247,11 +247,11 @@ forvalues replication = 1/1000 {
|
|
|
|
|
qui save `data',replace
|
|
|
|
|
}
|
|
|
|
|
}
|
|
|
|
|
export delimited using "`path'/scenario_4D_200.csv", replace
|
|
|
|
|
export delimited using "`path'/scenario_5D_100.csv", replace
|
|
|
|
|
|
|
|
|
|
* Scenario 4E : H_1 is TRUE / Effect size 0.4 / DIF on control
|
|
|
|
|
* Scenario 5E : H_1 is TRUE / Effect size 0.4 / DIF on control
|
|
|
|
|
|
|
|
|
|
di "SCENARIO 4E - N=200"
|
|
|
|
|
di "Scenario 5E - N=200"
|
|
|
|
|
|
|
|
|
|
forvalues replication = 1/1000 {
|
|
|
|
|
if mod(`replication',10)==0 {
|
|
|
|
@ -300,7 +300,7 @@ forvalues replication = 1/1000 {
|
|
|
|
|
qui save `data',replace
|
|
|
|
|
}
|
|
|
|
|
}
|
|
|
|
|
export delimited using "`path'/scenario_4E_200.csv", replace
|
|
|
|
|
export delimited using "`path'/scenario_5E_100.csv", replace
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
@ -313,11 +313,11 @@ forvalues replication = 1/1000 {
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
** Scenario 5: J = 4 items / M = 4 modalities / DIF SIZE = 0.3
|
|
|
|
|
** Scenario 6: J = 4 items / M = 4 modalities / DIF SIZE = 0.3
|
|
|
|
|
|
|
|
|
|
* Scenario 5A : H_0 is TRUE / DIF on treatment
|
|
|
|
|
* Scenario 6A : H_0 is TRUE / DIF on treatment
|
|
|
|
|
|
|
|
|
|
di "SCENARIO 5A - N=200"
|
|
|
|
|
di "Scenario 6A - N=200"
|
|
|
|
|
|
|
|
|
|
forvalues replication = 1/1000 {
|
|
|
|
|
if mod(`replication',10)==0 {
|
|
|
|
@ -366,11 +366,11 @@ forvalues replication = 1/1000 {
|
|
|
|
|
qui save `data',replace
|
|
|
|
|
}
|
|
|
|
|
}
|
|
|
|
|
export delimited using "`path'/scenario_5A_200.csv", replace
|
|
|
|
|
export delimited using "`path'/scenario_6A_100.csv", replace
|
|
|
|
|
|
|
|
|
|
* Scenario 5B : H_1 is TRUE / Effect size 0.2 / DIF on treatment
|
|
|
|
|
* Scenario 6B : H_1 is TRUE / Effect size 0.2 / DIF on treatment
|
|
|
|
|
|
|
|
|
|
di "SCENARIO 5B - N=200"
|
|
|
|
|
di "Scenario 6B - N=200"
|
|
|
|
|
|
|
|
|
|
forvalues replication = 1/1000 {
|
|
|
|
|
if mod(`replication',10)==0 {
|
|
|
|
@ -419,11 +419,11 @@ forvalues replication = 1/1000 {
|
|
|
|
|
qui save `data',replace
|
|
|
|
|
}
|
|
|
|
|
}
|
|
|
|
|
export delimited using "`path'/scenario_5B_200.csv", replace
|
|
|
|
|
export delimited using "`path'/scenario_6B_100.csv", replace
|
|
|
|
|
|
|
|
|
|
* Scenario 5C : H_1 is TRUE / Effect size 0.2 / DIF on control
|
|
|
|
|
* Scenario 6C : H_1 is TRUE / Effect size 0.2 / DIF on control
|
|
|
|
|
|
|
|
|
|
di "SCENARIO 5C - N=200"
|
|
|
|
|
di "Scenario 6C - N=200"
|
|
|
|
|
|
|
|
|
|
forvalues replication = 1/1000 {
|
|
|
|
|
if mod(`replication',10)==0 {
|
|
|
|
@ -472,11 +472,11 @@ forvalues replication = 1/1000 {
|
|
|
|
|
qui save `data',replace
|
|
|
|
|
}
|
|
|
|
|
}
|
|
|
|
|
export delimited using "`path'/scenario_5C_200.csv", replace
|
|
|
|
|
export delimited using "`path'/scenario_6C_100.csv", replace
|
|
|
|
|
|
|
|
|
|
* Scenario 5D : H_1 is TRUE / Effect size 0.4 / DIF on treatment
|
|
|
|
|
* Scenario 6D : H_1 is TRUE / Effect size 0.4 / DIF on treatment
|
|
|
|
|
|
|
|
|
|
di "SCENARIO 5D - N=200"
|
|
|
|
|
di "Scenario 6D - N=200"
|
|
|
|
|
|
|
|
|
|
forvalues replication = 1/1000 {
|
|
|
|
|
if mod(`replication',10)==0 {
|
|
|
|
@ -525,11 +525,11 @@ forvalues replication = 1/1000 {
|
|
|
|
|
qui save `data',replace
|
|
|
|
|
}
|
|
|
|
|
}
|
|
|
|
|
export delimited using "`path'/scenario_5D_200.csv", replace
|
|
|
|
|
export delimited using "`path'/scenario_6D_100.csv", replace
|
|
|
|
|
|
|
|
|
|
* Scenario 5E : H_1 is TRUE / Effect size 0.4 / DIF on control
|
|
|
|
|
* Scenario 6E : H_1 is TRUE / Effect size 0.4 / DIF on control
|
|
|
|
|
|
|
|
|
|
di "SCENARIO 5E - N=200"
|
|
|
|
|
di "Scenario 6E - N=200"
|
|
|
|
|
|
|
|
|
|
forvalues replication = 1/1000 {
|
|
|
|
|
if mod(`replication',10)==0 {
|
|
|
|
@ -578,7 +578,7 @@ forvalues replication = 1/1000 {
|
|
|
|
|
qui save `data',replace
|
|
|
|
|
}
|
|
|
|
|
}
|
|
|
|
|
export delimited using "`path'/scenario_5E_200.csv", replace
|
|
|
|
|
export delimited using "`path'/scenario_6E_100.csv", replace
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
@ -590,11 +590,11 @@ forvalues replication = 1/1000 {
|
|
|
|
|
|
|
|
|
|
**** Scenarios with : DIF size 0.5 ****
|
|
|
|
|
|
|
|
|
|
** Scenario 6: J = 4 items / M = 2 modalities / DIF SIZE = 0.5
|
|
|
|
|
** Scenario 7: J = 4 items / M = 2 modalities / DIF SIZE = 0.5
|
|
|
|
|
|
|
|
|
|
* Scenario 6A : H_0 is TRUE / DIF on treatment
|
|
|
|
|
* Scenario 7A : H_0 is TRUE / DIF on treatment
|
|
|
|
|
|
|
|
|
|
di "SCENARIO 6A - N=200"
|
|
|
|
|
di "Scenario 7A - N=200"
|
|
|
|
|
|
|
|
|
|
forvalues replication = 1/1000 {
|
|
|
|
|
if mod(`replication',10)==0 {
|
|
|
|
@ -643,12 +643,12 @@ forvalues replication = 1/1000 {
|
|
|
|
|
qui save `data',replace
|
|
|
|
|
}
|
|
|
|
|
}
|
|
|
|
|
export delimited using "`path'/scenario_6A_200.csv", replace
|
|
|
|
|
export delimited using "`path'/scenario_7A_100.csv", replace
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
* Scenario 6B : H_1 is TRUE / Effect size 0.2 / DIF on treatment
|
|
|
|
|
* Scenario 7B : H_1 is TRUE / Effect size 0.2 / DIF on treatment
|
|
|
|
|
|
|
|
|
|
di "SCENARIO 6B - N=200"
|
|
|
|
|
di "Scenario 7B - N=200"
|
|
|
|
|
|
|
|
|
|
forvalues replication = 1/1000 {
|
|
|
|
|
if mod(`replication',10)==0 {
|
|
|
|
@ -697,11 +697,11 @@ forvalues replication = 1/1000 {
|
|
|
|
|
qui save `data',replace
|
|
|
|
|
}
|
|
|
|
|
}
|
|
|
|
|
export delimited using "`path'/scenario_6B_200.csv", replace
|
|
|
|
|
export delimited using "`path'/scenario_7B_100.csv", replace
|
|
|
|
|
|
|
|
|
|
* Scenario 6C : H_1 is TRUE / Effect size 0.2 / DIF on control
|
|
|
|
|
* Scenario 7C : H_1 is TRUE / Effect size 0.2 / DIF on control
|
|
|
|
|
|
|
|
|
|
di "SCENARIO 6C - N=200"
|
|
|
|
|
di "Scenario 7C - N=200"
|
|
|
|
|
|
|
|
|
|
forvalues replication = 1/1000 {
|
|
|
|
|
if mod(`replication',10)==0 {
|
|
|
|
@ -750,12 +750,12 @@ forvalues replication = 1/1000 {
|
|
|
|
|
qui save `data',replace
|
|
|
|
|
}
|
|
|
|
|
}
|
|
|
|
|
export delimited using "`path'/scenario_6C_200.csv", replace
|
|
|
|
|
export delimited using "`path'/scenario_7C_100.csv", replace
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
* Scenario 6D : H_1 is TRUE / Effect size 0.4 / DIF on treatment
|
|
|
|
|
* Scenario 7D : H_1 is TRUE / Effect size 0.4 / DIF on treatment
|
|
|
|
|
|
|
|
|
|
di "SCENARIO 6D - N=200"
|
|
|
|
|
di "Scenario 7D - N=200"
|
|
|
|
|
|
|
|
|
|
forvalues replication = 1/1000 {
|
|
|
|
|
if mod(`replication',10)==0 {
|
|
|
|
@ -804,11 +804,11 @@ forvalues replication = 1/1000 {
|
|
|
|
|
qui save `data',replace
|
|
|
|
|
}
|
|
|
|
|
}
|
|
|
|
|
export delimited using "`path'/scenario_6D_200.csv", replace
|
|
|
|
|
export delimited using "`path'/scenario_7D_100.csv", replace
|
|
|
|
|
|
|
|
|
|
* Scenario 6E : H_1 is TRUE / Effect size 0.4 / DIF on control
|
|
|
|
|
* Scenario 7E : H_1 is TRUE / Effect size 0.4 / DIF on control
|
|
|
|
|
|
|
|
|
|
di "SCENARIO 6E - N=200"
|
|
|
|
|
di "Scenario 7E - N=200"
|
|
|
|
|
|
|
|
|
|
forvalues replication = 1/1000 {
|
|
|
|
|
if mod(`replication',10)==0 {
|
|
|
|
@ -857,7 +857,7 @@ forvalues replication = 1/1000 {
|
|
|
|
|
qui save `data',replace
|
|
|
|
|
}
|
|
|
|
|
}
|
|
|
|
|
export delimited using "`path'/scenario_6E_200.csv", replace
|
|
|
|
|
export delimited using "`path'/scenario_7E_100.csv", replace
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
@ -870,11 +870,11 @@ forvalues replication = 1/1000 {
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
** Scenario 7: J = 4 items / M = 4 modalities / DIF SIZE = 0.5
|
|
|
|
|
** Scenario 8: J = 4 items / M = 4 modalities / DIF SIZE = 0.5
|
|
|
|
|
|
|
|
|
|
* Scenario 7A : H_0 is TRUE / DIF on treatment
|
|
|
|
|
* Scenario 8A : H_0 is TRUE / DIF on treatment
|
|
|
|
|
|
|
|
|
|
di "SCENARIO 7A - N=200"
|
|
|
|
|
di "Scenario 8A - N=200"
|
|
|
|
|
|
|
|
|
|
forvalues replication = 1/1000 {
|
|
|
|
|
if mod(`replication',10)==0 {
|
|
|
|
@ -923,11 +923,11 @@ forvalues replication = 1/1000 {
|
|
|
|
|
qui save `data',replace
|
|
|
|
|
}
|
|
|
|
|
}
|
|
|
|
|
export delimited using "`path'/scenario_7A_200.csv", replace
|
|
|
|
|
export delimited using "`path'/scenario_8A_100.csv", replace
|
|
|
|
|
|
|
|
|
|
* Scenario 7B : H_1 is TRUE / Effect size 0.2 / DIF on treatment
|
|
|
|
|
* Scenario 8B : H_1 is TRUE / Effect size 0.2 / DIF on treatment
|
|
|
|
|
|
|
|
|
|
di "SCENARIO 7B - N=200"
|
|
|
|
|
di "Scenario 8B - N=200"
|
|
|
|
|
|
|
|
|
|
forvalues replication = 1/1000 {
|
|
|
|
|
if mod(`replication',10)==0 {
|
|
|
|
@ -976,11 +976,11 @@ forvalues replication = 1/1000 {
|
|
|
|
|
qui save `data',replace
|
|
|
|
|
}
|
|
|
|
|
}
|
|
|
|
|
export delimited using "`path'/scenario_7B_200.csv", replace
|
|
|
|
|
export delimited using "`path'/scenario_8B_100.csv", replace
|
|
|
|
|
|
|
|
|
|
* Scenario 7C : H_1 is TRUE / Effect size 0.2 / DIF on control
|
|
|
|
|
* Scenario 8C : H_1 is TRUE / Effect size 0.2 / DIF on control
|
|
|
|
|
|
|
|
|
|
di "SCENARIO 7C - N=200"
|
|
|
|
|
di "Scenario 8C - N=200"
|
|
|
|
|
|
|
|
|
|
forvalues replication = 1/1000 {
|
|
|
|
|
if mod(`replication',10)==0 {
|
|
|
|
@ -1029,11 +1029,11 @@ forvalues replication = 1/1000 {
|
|
|
|
|
qui save `data',replace
|
|
|
|
|
}
|
|
|
|
|
}
|
|
|
|
|
export delimited using "`path'/scenario_7C_200.csv", replace
|
|
|
|
|
export delimited using "`path'/scenario_8C_100.csv", replace
|
|
|
|
|
|
|
|
|
|
* Scenario 7D : H_1 is TRUE / Effect size 0.4 / DIF on treatment
|
|
|
|
|
* Scenario 8D : H_1 is TRUE / Effect size 0.4 / DIF on treatment
|
|
|
|
|
|
|
|
|
|
di "SCENARIO 7D - N=200"
|
|
|
|
|
di "Scenario 8D - N=200"
|
|
|
|
|
|
|
|
|
|
forvalues replication = 1/1000 {
|
|
|
|
|
if mod(`replication',10)==0 {
|
|
|
|
@ -1082,11 +1082,11 @@ forvalues replication = 1/1000 {
|
|
|
|
|
qui save `data',replace
|
|
|
|
|
}
|
|
|
|
|
}
|
|
|
|
|
export delimited using "`path'/scenario_7D_200.csv", replace
|
|
|
|
|
export delimited using "`path'/scenario_8D_100.csv", replace
|
|
|
|
|
|
|
|
|
|
* Scenario 7E : H_1 is TRUE / Effect size 0.4 / DIF on control
|
|
|
|
|
* Scenario 8E : H_1 is TRUE / Effect size 0.4 / DIF on control
|
|
|
|
|
|
|
|
|
|
di "SCENARIO 7E - N=200"
|
|
|
|
|
di "Scenario 8E - N=200"
|
|
|
|
|
|
|
|
|
|
forvalues replication = 1/1000 {
|
|
|
|
|
if mod(`replication',10)==0 {
|
|
|
|
@ -1135,4 +1135,4 @@ forvalues replication = 1/1000 {
|
|
|
|
|
qui save `data',replace
|
|
|
|
|
}
|
|
|
|
|
}
|
|
|
|
|
export delimited using "`path'/scenario_7E_200.csv", replace
|
|
|
|
|
export delimited using "`path'/scenario_8E_100.csv", replace
|
|
|
|
|