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

Updated and created residuals analysis .do file

Browse Source
This commit is contained in:
Corentin Choisy
2024-04-17 12:00:36 +02:00
parent af45e435f9
commit 5cb862072c
2 changed files with 633 additions and 633 deletions
Download Patch File Download Diff File Expand all files Collapse all files

26
Scripts/Analysis/DIF-RESIDUS/pcm_dif_residus.do
View File

@ -16,12 +16,12 @@ local N = "50 100 200 300"
foreach s in `ss' { foreach s in `ss' {
foreach Nnn in `N' { foreach Nnn in `N' {
local Nn = `Nnn' local Nn = `Nnn'
local path_data = "/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N`Nn'" local path_data = "/home/corentin/Documents/These/Recherche/Simulations/Analysis/RESALI/Detection_data"
if (`s'<=4) { local path_res = "/home/corentin/Documents/These/Recherche/Simulations/Analysis/RESALI/Results/N`Nn'"
local path_data = "/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N`Nn'"
}
local path_res = "/home/corentin/Documents/These/Recherche/Simulations/Analysis/ROSALI-DIF/N`Nn'"
local scenarios = "A B C D E F G" local scenarios = "A B C D E F G"
if (`s' <= 4) {
local scenarios = "A B C D E"
}
foreach scen in `scenarios' { foreach scen in `scenarios' {
clear clear
import delim "`path_data'/scenario_`s'`scen'_`Nn'.csv", encoding(ISO-8859-2) case(preserve) clear import delim "`path_data'/scenario_`s'`scen'_`Nn'.csv", encoding(ISO-8859-2) case(preserve) clear
@ -40,12 +40,12 @@ local N = "50 100 200 300"
local nbitems=7 local nbitems=7
} }
if (mod(`s',2)==0) { if (mod(`s',2)==0) {
local nbmoda=3 local nbmoda=3
} }
else { else {
local nbmoda=1 local nbmoda=1
} }
if (`s'<=4) { if (`s'<=4) {
local nbdif=0 local nbdif=0
@ -75,7 +75,7 @@ else {
forvalues i=1/`nbdif' { forvalues i=1/`nbdif' {
local colna = "`colna'"+"real_dif_`i' " local colna = "`colna'"+"real_dif_`i' "
} }
local colna = "`colna'" + "beta " + "se_beta" local colna = "`colna'" + "beta " + "se_beta"
mat outmat = J(1000,`taillemat',.) mat outmat = J(1000,`taillemat',.)
@ -359,6 +359,6 @@ local colna = "`colna'" + "beta " + "se_beta"
} }
putexcel set "`path_res'/`s'`scen'_`Nn'_original.xls", sheet("outmat") replace putexcel set "`path_res'/`s'`scen'_`Nn'_original.xls", sheet("outmat") replace
putexcel A1=matrix(outmat), colnames putexcel A1=matrix(outmat), colnames
} }
} }
} }

16
Scripts/Analysis/DIF-ROSALI/pcm_dif_rosali.do
View File

@ -40,12 +40,12 @@ local N = "50 100 200 300"
local nbitems=7 local nbitems=7
} }
if (mod(`s',2)==0) { if (mod(`s',2)==0) {
local nbmoda=3 local nbmoda=3
} }
else { else {
local nbmoda=1 local nbmoda=1
} }
if (`s'<=4) { if (`s'<=4) {
local nbdif=0 local nbdif=0
@ -75,7 +75,7 @@ else {
forvalues i=1/`nbdif' { forvalues i=1/`nbdif' {
local colna = "`colna'"+"real_dif_`i' " local colna = "`colna'"+"real_dif_`i' "
} }
local colna = "`colna'" + "beta " + "se_beta" local colna = "`colna'" + "beta " + "se_beta"
mat outmat = J(1000,`taillemat',.) mat outmat = J(1000,`taillemat',.)
@ -359,6 +359,6 @@ local colna = "`colna'" + "beta " + "se_beta"
} }
putexcel set "`path_res'/`s'`scen'_`Nn'_original.xls", sheet("outmat") replace putexcel set "`path_res'/`s'`scen'_`Nn'_original.xls", sheet("outmat") replace
putexcel A1=matrix(outmat), colnames putexcel A1=matrix(outmat), colnames
} }
} }
} }