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

Updated all scenarios

Browse Source
main
Corentin Choisy 11 months ago
parent 832f836e55
commit c2b7c57293
9 changed files with 8672 additions and 2622 deletions
Show Stats Download Patch File Download Diff File

24
Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_100.do
Unescape Escape View File

@ -1123,7 +1123,7 @@ forvalues replication = 1/1000 {
}
export delimited using "`path'/scenario_7E_100.csv", replace
* Scenario 7F : H_1 is TRUE / Effect size 0.4 / DIF negative
* Scenario 7F : H_1 is TRUE / Effect size -0.2 / DIF negative
di "Scenario 7F - N=100"
@ -1179,7 +1179,7 @@ di "Scenario 7F - N=100"
}
export delimited using "`path'/scenario_7F_100.csv", replace
* Scenario 7G : H_1 is TRUE / Effect size 0.4 / DIF negative
* Scenario 7G : H_1 is TRUE / Effect size -0.4 / DIF negative
di "Scenario 7G - N=100"
@ -1528,7 +1528,7 @@ forvalues replication = 1/1000 {
export delimited using "`path'/scenario_8E_100.csv", replace
* Scenario 8F : H_1 is TRUE / Effect size 0.4 / DIF negative
* Scenario 8F : H_1 is TRUE / Effect size -0.2 / DIF negative
di "Scenario 8F - N=100"
@ -1585,7 +1585,7 @@ di "Scenario 8F - N=100"
export delimited using "`path'/scenario_8F_100.csv", replace
* Scenario 8G : H_1 is TRUE / Effect size 0.4 / DIF negative
* Scenario 8G : H_1 is TRUE / Effect size -0.4 / DIF negative
di "Scenario 8G - N=100"
@ -2040,7 +2040,7 @@ forvalues replication = 1/1000 {
* Scenario 9F : H_1 is TRUE / Effect size 0.4 / DIF negative x2
* Scenario 9F : H_1 is TRUE / Effect size -0.2 / DIF negative x2
di "Scenario 9F - N=100"
@ -2117,7 +2117,7 @@ forvalues replication = 1/1000 {
export delimited using "`path'/scenario_9F_100.csv", replace
* Scenario 9G : H_1 is TRUE / Effect size 0.4 / DIF negative x2
* Scenario 9G : H_1 is TRUE / Effect size -0.4 / DIF negative x2
di "Scenario 9G - N=100"
@ -2579,7 +2579,7 @@ forvalues replication = 1/1000 {
}
export delimited using "`path'/scenario_10E_100.csv", replace
* Scenario 10F : H_1 is TRUE / Effect size 0.4 / DIF negative x2
* Scenario 10F : H_1 is TRUE / Effect size -0.2 / DIF negative x2
di "Scenario 10F - N=100"
@ -2656,7 +2656,7 @@ forvalues replication = 1/1000 {
export delimited using "`path'/scenario_10F_100.csv", replace
* Scenario 10G : H_1 is TRUE / Effect size 0.4 / DIF negative x2
* Scenario 10G : H_1 is TRUE / Effect size -0.4 / DIF negative x2
di "Scenario 10G - N=100"
@ -3120,7 +3120,7 @@ forvalues replication = 1/1000 {
}
export delimited using "`path'/scenario_11E_100.csv", replace
* Scenario 11F : H_1 is TRUE / Effect size 0.4 / DIF negative x2
* Scenario 11F : H_1 is TRUE / Effect size -0.2 / DIF negative x2
di "Scenario 11F - N=100"
@ -3197,7 +3197,7 @@ forvalues replication = 1/1000 {
export delimited using "`path'/scenario_11F_100.csv", replace
* Scenario 11G : H_1 is TRUE / Effect size 0.4 / DIF negative x2
* Scenario 11G : H_1 is TRUE / Effect size -0.4 / DIF negative x2
di "Scenario 11G - N=100"
@ -3661,7 +3661,7 @@ forvalues replication = 1/1000 {
* Scenario 12F : H_1 is TRUE / Effect size 0.4 / DIF negative x2
* Scenario 12F : H_1 is TRUE / Effect size -0.2 / DIF negative x2
di "Scenario 12F - N=100"
@ -3738,7 +3738,7 @@ di "Scenario 12F - N=100"
export delimited using "`path'/scenario_12F_100.csv", replace
* Scenario 12G : H_1 is TRUE / Effect size 0.4 / DIF negative x2
* Scenario 12G : H_1 is TRUE / Effect size -0.4 / DIF negative x2
di "Scenario 12G - N=100"

24
Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_200.do
Unescape Escape View File

@ -1123,7 +1123,7 @@ forvalues replication = 1/1000 {
}
export delimited using "`path'/scenario_7E_200.csv", replace
* Scenario 7F : H_1 is TRUE / Effect size 0.4 / DIF negative
* Scenario 7F : H_1 is TRUE / Effect size -0.2 / DIF negative
di "Scenario 7F - N=200"
@ -1179,7 +1179,7 @@ di "Scenario 7F - N=200"
}
export delimited using "`path'/scenario_7F_200.csv", replace
* Scenario 7G : H_1 is TRUE / Effect size 0.4 / DIF negative
* Scenario 7G : H_1 is TRUE / Effect size -0.4 / DIF negative
di "Scenario 7G - N=200"
@ -1528,7 +1528,7 @@ forvalues replication = 1/1000 {
export delimited using "`path'/scenario_8E_200.csv", replace
* Scenario 8F : H_1 is TRUE / Effect size 0.4 / DIF negative
* Scenario 8F : H_1 is TRUE / Effect size -0.2 / DIF negative
di "Scenario 8F - N=200"
@ -1585,7 +1585,7 @@ di "Scenario 8F - N=200"
export delimited using "`path'/scenario_8F_200.csv", replace
* Scenario 8G : H_1 is TRUE / Effect size 0.4 / DIF negative
* Scenario 8G : H_1 is TRUE / Effect size -0.4 / DIF negative
di "Scenario 8G - N=200"
@ -2040,7 +2040,7 @@ forvalues replication = 1/1000 {
* Scenario 9F : H_1 is TRUE / Effect size 0.4 / DIF negative x2
* Scenario 9F : H_1 is TRUE / Effect size -0.2 / DIF negative x2
di "Scenario 9F - N=200"
@ -2117,7 +2117,7 @@ forvalues replication = 1/1000 {
export delimited using "`path'/scenario_9F_200.csv", replace
* Scenario 9G : H_1 is TRUE / Effect size 0.4 / DIF negative x2
* Scenario 9G : H_1 is TRUE / Effect size -0.4 / DIF negative x2
di "Scenario 9G - N=200"
@ -2579,7 +2579,7 @@ forvalues replication = 1/1000 {
}
export delimited using "`path'/scenario_10E_200.csv", replace
* Scenario 10F : H_1 is TRUE / Effect size 0.4 / DIF negative x2
* Scenario 10F : H_1 is TRUE / Effect size -0.2 / DIF negative x2
di "Scenario 10F - N=200"
@ -2656,7 +2656,7 @@ forvalues replication = 1/1000 {
export delimited using "`path'/scenario_10F_200.csv", replace
* Scenario 10G : H_1 is TRUE / Effect size 0.4 / DIF negative x2
* Scenario 10G : H_1 is TRUE / Effect size -0.4 / DIF negative x2
di "Scenario 10G - N=200"
@ -3120,7 +3120,7 @@ forvalues replication = 1/1000 {
}
export delimited using "`path'/scenario_11E_200.csv", replace
* Scenario 11F : H_1 is TRUE / Effect size 0.4 / DIF negative x2
* Scenario 11F : H_1 is TRUE / Effect size -0.2 / DIF negative x2
di "Scenario 11F - N=200"
@ -3197,7 +3197,7 @@ forvalues replication = 1/1000 {
export delimited using "`path'/scenario_11F_200.csv", replace
* Scenario 11G : H_1 is TRUE / Effect size 0.4 / DIF negative x2
* Scenario 11G : H_1 is TRUE / Effect size -0.4 / DIF negative x2
di "Scenario 11G - N=200"
@ -3661,7 +3661,7 @@ forvalues replication = 1/1000 {
* Scenario 12F : H_1 is TRUE / Effect size 0.4 / DIF negative x2
* Scenario 12F : H_1 is TRUE / Effect size -0.2 / DIF negative x2
di "Scenario 12F - N=200"
@ -3738,7 +3738,7 @@ di "Scenario 12F - N=200"
export delimited using "`path'/scenario_12F_200.csv", replace
* Scenario 12G : H_1 is TRUE / Effect size 0.4 / DIF negative x2
* Scenario 12G : H_1 is TRUE / Effect size -0.4 / DIF negative x2
di "Scenario 12G - N=200"

24
Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_300.do
Unescape Escape View File

@ -1123,7 +1123,7 @@ forvalues replication = 1/1000 {
}
export delimited using "`path'/scenario_7E_300.csv", replace
* Scenario 7F : H_1 is TRUE / Effect size 0.4 / DIF negative
* Scenario 7F : H_1 is TRUE / Effect size -0.2 / DIF negative
di "Scenario 7F - N=300"
@ -1179,7 +1179,7 @@ di "Scenario 7F - N=300"
}
export delimited using "`path'/scenario_7F_300.csv", replace
* Scenario 7G : H_1 is TRUE / Effect size 0.4 / DIF negative
* Scenario 7G : H_1 is TRUE / Effect size -0.4 / DIF negative
di "Scenario 7G - N=300"
@ -1528,7 +1528,7 @@ forvalues replication = 1/1000 {
export delimited using "`path'/scenario_8E_300.csv", replace
* Scenario 8F : H_1 is TRUE / Effect size 0.4 / DIF negative
* Scenario 8F : H_1 is TRUE / Effect size -0.2 / DIF negative
di "Scenario 8F - N=300"
@ -1585,7 +1585,7 @@ di "Scenario 8F - N=300"
export delimited using "`path'/scenario_8F_300.csv", replace
* Scenario 8G : H_1 is TRUE / Effect size 0.4 / DIF negative
* Scenario 8G : H_1 is TRUE / Effect size -0.4 / DIF negative
di "Scenario 8G - N=300"
@ -2040,7 +2040,7 @@ forvalues replication = 1/1000 {
* Scenario 9F : H_1 is TRUE / Effect size 0.4 / DIF negative x2
* Scenario 9F : H_1 is TRUE / Effect size -0.2 / DIF negative x2
di "Scenario 9F - N=300"
@ -2117,7 +2117,7 @@ forvalues replication = 1/1000 {
export delimited using "`path'/scenario_9F_300.csv", replace
* Scenario 9G : H_1 is TRUE / Effect size 0.4 / DIF negative x2
* Scenario 9G : H_1 is TRUE / Effect size -0.4 / DIF negative x2
di "Scenario 9G - N=300"
@ -2579,7 +2579,7 @@ forvalues replication = 1/1000 {
}
export delimited using "`path'/scenario_10E_300.csv", replace
* Scenario 10F : H_1 is TRUE / Effect size 0.4 / DIF negative x2
* Scenario 10F : H_1 is TRUE / Effect size -0.2 / DIF negative x2
di "Scenario 10F - N=300"
@ -2656,7 +2656,7 @@ forvalues replication = 1/1000 {
export delimited using "`path'/scenario_10F_300.csv", replace
* Scenario 10G : H_1 is TRUE / Effect size 0.4 / DIF negative x2
* Scenario 10G : H_1 is TRUE / Effect size -0.4 / DIF negative x2
di "Scenario 10G - N=300"
@ -3120,7 +3120,7 @@ forvalues replication = 1/1000 {
}
export delimited using "`path'/scenario_11E_300.csv", replace
* Scenario 11F : H_1 is TRUE / Effect size 0.4 / DIF negative x2
* Scenario 11F : H_1 is TRUE / Effect size -0.2 / DIF negative x2
di "Scenario 11F - N=300"
@ -3197,7 +3197,7 @@ forvalues replication = 1/1000 {
export delimited using "`path'/scenario_11F_300.csv", replace
* Scenario 11G : H_1 is TRUE / Effect size 0.4 / DIF negative x2
* Scenario 11G : H_1 is TRUE / Effect size -0.4 / DIF negative x2
di "Scenario 11G - N=300"
@ -3661,7 +3661,7 @@ forvalues replication = 1/1000 {
* Scenario 12F : H_1 is TRUE / Effect size 0.4 / DIF negative x2
* Scenario 12F : H_1 is TRUE / Effect size -0.2 / DIF negative x2
di "Scenario 12F - N=300"
@ -3738,7 +3738,7 @@ di "Scenario 12F - N=300"
export delimited using "`path'/scenario_12F_300.csv", replace
* Scenario 12G : H_1 is TRUE / Effect size 0.4 / DIF negative x2
* Scenario 12G : H_1 is TRUE / Effect size -0.4 / DIF negative x2
di "Scenario 12G - N=300"

3300
Scripts/Scenarios/DIF/scenarios_DIF_J7_baseline_100.do
View File

File diff suppressed because it is too large Load Diff

3300
Scripts/Scenarios/DIF/scenarios_DIF_J7_baseline_200.do
View File

File diff suppressed because it is too large Load Diff

3300
Scripts/Scenarios/DIF/scenarios_DIF_J7_baseline_300.do
View File

File diff suppressed because it is too large Load Diff

313
Scripts/Scenarios/NoDIF/scenarios_noDIF_baseline_100.do
Unescape Escape View File

@ -145,6 +145,84 @@ forvalues replication = 1/1000 {
}
export delimited using "`path'/scenario_1C_100.csv", replace
* Scenario 1D : H_0 is FALSE / Effect size = -0.2
di "SCENARIO 1D - 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_1D_100.csv", replace
* Scenario 1E : H_0 is FALSE / Effect size = -0.4
di "SCENARIO 1E - 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_1E_100.csv", replace
**------------------------------------------------------------------------------------**
@ -264,6 +342,88 @@ forvalues replication = 1/1000 {
export delimited using "`path'/scenario_2C_100.csv", replace
* Scenario 2D : H_0 is FALSE / Effect size = -0.2
di "SCENARIO 2D - 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_2D_100.csv", replace
* Scenario 2E : H_0 is FALSE / Effect size = -0.4
di "SCENARIO 2E - 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_2E_100.csv", replace
**------------------------------------------------------------------------------------**
@ -382,6 +542,84 @@ forvalues replication = 1/1000 {
export delimited using "`path'/scenario_3C_100.csv", replace
* Scenario 3D : H_0 is FALSE / Effect size = -0.2
di "SCENARIO 3D - 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_3D_100.csv", replace
* Scenario 3E : H_0 is FALSE / Effect size = -0.4
di "SCENARIO 3E - 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_3E_100.csv", replace
**------------------------------------------------------------------------------------**
@ -498,3 +736,78 @@ forvalues replication = 1/1000 {
}
}
export delimited using "`path'/scenario_4C_100.csv", replace
* Scenario 4D : H_0 is FALSE / Effect size = -0.2
di "SCENARIO 4D - 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_4D_100.csv", replace
* Scenario 4E : H_0 is FALSE / Effect size = -0.4
di "SCENARIO 4E - 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_4E_100.csv", replace

353
Scripts/Scenarios/NoDIF/scenarios_noDIF_baseline_200.do
Unescape Escape View File

@ -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=100/200/300
* outputs : scenario_1,scenario_2,scenario_3,scenario_4, 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'
@ -20,17 +20,21 @@ 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/N200"
local Nn = 200
** Scenario 1: J = 4 items / M = 2 modalities / N=100 per group / TT=treatment variable
** Scenario 1: J = 4 items / M = 2 modalities / N=200 per group / TT=treatment variable
* Scenario 1A : H_0 is TRUE
di "SCENARIO 1A - N=200"
* Scenario 1A : H_0 is TRUE
forvalues replication = 1/1000 {
di "SCENARIO 1A - N=200"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
@ -66,6 +70,7 @@ local Nn = 200
* Scenario 1B : H_0 is FALSE / Effect size = 0.2
di "SCENARIO 1B - N=200"
forvalues replication = 1/1000 {
@ -102,10 +107,11 @@ 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"
* Scenario 1C : H_0 is FALSE / Effect size = 0.4
forvalues replication = 1/1000 {
di "SCENARIO 1C - N=200"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
@ -139,11 +145,89 @@ local Nn = 200
}
export delimited using "`path'/scenario_1C_200.csv", replace
* Scenario 1D : H_0 is FALSE / Effect size = -0.2
di "SCENARIO 1D - 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_1D_200.csv", replace
* Scenario 1E : H_0 is FALSE / Effect size = -0.4
di "SCENARIO 1E - 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_1E_200.csv", replace
**------------------------------------------------------------------------------------**
** Scenario 2: J = 4 items / M = 4 modalities / N=100 per group / TT=treatment variable
** Scenario 2: J = 4 items / M = 4 modalities / N=200 per group / TT=treatment variable
* Scenario 2A : H_0 is TRUE
di "SCENARIO 2A - N=200"
@ -258,17 +342,99 @@ local Nn = 200
export delimited using "`path'/scenario_2C_200.csv", replace
**------------------------------------------------------------------------------------**
* Scenario 2D : H_0 is FALSE / Effect size = -0.2
di "SCENARIO 2D - 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
** Scenario 3: J = 7 items / M = 2 modalities / N=100 per group / TT=treatment variable
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_2D_200.csv", replace
* Scenario 3A : H_0 is TRUE
di "SCENARIO 3A - N=200"
* Scenario 2E : H_0 is FALSE / Effect size = -0.4
di "SCENARIO 2E - 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_2E_200.csv", replace
**------------------------------------------------------------------------------------**
** Scenario 3: J = 7 items / M = 2 modalities / N=200 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
@ -338,10 +504,10 @@ 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"
* Scenario 3C : H_0 is FALSE / Effect size = 0.4
di "SCENARIO 3C - N=200"
forvalues replication = 1/1000 {
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
@ -376,10 +542,88 @@ local Nn = 200
export delimited using "`path'/scenario_3C_200.csv", replace
* Scenario 3D : H_0 is FALSE / Effect size = -0.2
di "SCENARIO 3D - 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_3D_200.csv", replace
* Scenario 3E : H_0 is FALSE / Effect size = -0.4
di "SCENARIO 3E - 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_3E_200.csv", replace
**------------------------------------------------------------------------------------**
** Scenario 4: J = 7 items / M = 4 modalities / N=100 per group / TT=treatment variable
** Scenario 4: J = 7 items / M = 4 modalities / N=200 per group / TT=treatment variable
* Scenario 4A : H_0 is TRUE
di "SCENARIO 4A - N=200"
@ -492,3 +736,78 @@ local Nn = 200
}
}
export delimited using "`path'/scenario_4C_200.csv", replace
* Scenario 4D : H_0 is FALSE / Effect size = -0.2
di "SCENARIO 4D - 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_4D_200.csv", replace
* Scenario 4E : H_0 is FALSE / Effect size = -0.4
di "SCENARIO 4E - 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_4E_200.csv", replace

354
Scripts/Scenarios/NoDIF/scenarios_noDIF_baseline_300.do
Unescape Escape View File

@ -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=100/200/300
* outputs : scenario_1,scenario_2,scenario_3,scenario_4, for N=300/200/300
*
*
* Warning : To obtain reproduce the data obtained in the .csv files in this repository, use 'simirt_setseed.ado' instead of 'simirt.ado'
@ -20,18 +20,21 @@ 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
* Scenarios with : n = 100
*==========================
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 1: J = 4 items / M = 2 modalities / N=300 per group / TT=treatment variable
* Scenario 1A : H_0 is TRUE
di "SCENARIO 1A - N=300"
* Scenario 1A : H_0 is TRUE
forvalues replication = 1/1000 {
di "SCENARIO 1A - N=300"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
@ -67,6 +70,7 @@ local Nn = 300
* Scenario 1B : H_0 is FALSE / Effect size = 0.2
di "SCENARIO 1B - N=300"
forvalues replication = 1/1000 {
@ -103,10 +107,11 @@ 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"
* Scenario 1C : H_0 is FALSE / Effect size = 0.4
forvalues replication = 1/1000 {
di "SCENARIO 1C - N=300"
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
@ -140,11 +145,89 @@ local Nn = 300
}
export delimited using "`path'/scenario_1C_300.csv", replace
* Scenario 1D : H_0 is FALSE / Effect size = -0.2
di "SCENARIO 1D - 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_1D_300.csv", replace
* Scenario 1E : H_0 is FALSE / Effect size = -0.4
di "SCENARIO 1E - 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_1E_300.csv", replace
**------------------------------------------------------------------------------------**
** Scenario 2: J = 4 items / M = 4 modalities / N=100 per group / TT=treatment variable
** Scenario 2: J = 4 items / M = 4 modalities / N=300 per group / TT=treatment variable
* Scenario 2A : H_0 is TRUE
di "SCENARIO 2A - N=300"
@ -259,17 +342,99 @@ local Nn = 300
export delimited using "`path'/scenario_2C_300.csv", replace
**------------------------------------------------------------------------------------**
* Scenario 2D : H_0 is FALSE / Effect size = -0.2
di "SCENARIO 2D - 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
** Scenario 3: J = 7 items / M = 2 modalities / N=100 per group / TT=treatment variable
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_2D_300.csv", replace
* Scenario 3A : H_0 is TRUE
di "SCENARIO 3A - N=300"
* Scenario 2E : H_0 is FALSE / Effect size = -0.4
di "SCENARIO 2E - 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_2E_300.csv", replace
**------------------------------------------------------------------------------------**
** Scenario 3: J = 7 items / M = 2 modalities / N=300 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
@ -339,10 +504,10 @@ 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"
* Scenario 3C : H_0 is FALSE / Effect size = 0.4
di "SCENARIO 3C - N=300"
forvalues replication = 1/1000 {
forvalues replication = 1/1000 {
if mod(`replication',10)==0 {
di "replication = `replication'"
}
@ -377,10 +542,88 @@ local Nn = 300
export delimited using "`path'/scenario_3C_300.csv", replace
* Scenario 3D : H_0 is FALSE / Effect size = -0.2
di "SCENARIO 3D - 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_3D_300.csv", replace
* Scenario 3E : H_0 is FALSE / Effect size = -0.4
di "SCENARIO 3E - 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_3E_300.csv", replace
**------------------------------------------------------------------------------------**
** Scenario 4: J = 7 items / M = 4 modalities / N=100 per group / TT=treatment variable
** Scenario 4: J = 7 items / M = 4 modalities / N=300 per group / TT=treatment variable
* Scenario 4A : H_0 is TRUE
di "SCENARIO 4A - N=300"
@ -493,3 +736,78 @@ local Nn = 300
}
}
export delimited using "`path'/scenario_4C_300.csv", replace
* Scenario 4D : H_0 is FALSE / Effect size = -0.2
di "SCENARIO 4D - 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_4D_300.csv", replace
* Scenario 4E : H_0 is FALSE / Effect size = -0.4
di "SCENARIO 4E - 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_4E_300.csv", replace

Write Preview
Loading…
Cancel
Save