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

Transférer les fichiers vers 'Scripts/Analysis/DIF'

Browse Source
This commit is contained in:
Corentin Choisy
2024-02-05 09:37:30 +00:00
parent ef1fbfbd4f
commit e8e7f555d6
1 changed files with 46 additions and 40 deletions
Download Patch File Download Diff File Expand all files Collapse all files

86
Scripts/Analysis/DIF/pcm_dif.do
View File

@ -12,7 +12,7 @@
*================================================================================================================================================ *================================================================================================================================================
* Load pcm.ado * Load pcm.ado
adopath+"/home/corentin/Documents/These/Recherche/Simulations/Modules/" adopath+"/home/sphere/Corentin/These/Recherche/Simulations/Modules/"
*========================== *==========================
* Scenarios with : J=4 * Scenarios with : J=4
@ -24,8 +24,8 @@
local N = "100 200 300" local N = "100 200 300"
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/sphere/Corentin/These/Recherche/Simulations/Data/DIF/N`Nn'"
local path_res = "/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N`Nn'" local path_res = "/home/sphere/Corentin/These/Recherche/Simulations/Analysis/DIF/N`Nn'"
local scenarios = "A B C D E F G" local scenarios = "A B C D E F G"
foreach scen in `scenarios' { foreach scen in `scenarios' {
clear clear
@ -83,8 +83,8 @@
local N = "100 200 300" local N = "100 200 300"
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/sphere/Corentin/These/Recherche/Simulations/Data/DIF/N`Nn'"
local path_res = "/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N`Nn'" local path_res = "/home/sphere/Corentin/These/Recherche/Simulations/Analysis/DIF/N`Nn'"
local scenarios = "A B C D E F G" local scenarios = "A B C D E F G"
foreach scen in `scenarios' { foreach scen in `scenarios' {
clear clear
@ -161,8 +161,8 @@ putexcel A1=matrix(outmat), colnames
local N = "100 200 300" local N = "100 200 300"
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/sphere/Corentin/These/Recherche/Simulations/Data/DIF/N`Nn'"
local path_res = "/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N`Nn'" local path_res = "/home/sphere/Corentin/These/Recherche/Simulations/Analysis/DIF/N`Nn'"
local scenarios = "A B C D E F G" local scenarios = "A B C D E F G"
foreach scen in `scenarios' { foreach scen in `scenarios' {
clear clear
@ -220,8 +220,8 @@ putexcel A1=matrix(outmat), colnames
local N = "100 200 300" local N = "100 200 300"
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/sphere/Corentin/These/Recherche/Simulations/Data/DIF/N`Nn'"
local path_res = "/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N`Nn'" local path_res = "/home/sphere/Corentin/These/Recherche/Simulations/Analysis/DIF/N`Nn'"
local scenarios = "A B C D E F G" local scenarios = "A B C D E F G"
foreach scen in `scenarios' { foreach scen in `scenarios' {
clear clear
@ -296,13 +296,15 @@ putexcel A1=matrix(outmat), colnames
****** Scenarios with DIF on 2 items ****** Scenarios with DIF on 2 items
* Load pcm.ado
adopath+"/home/sphere/Corentin/These/Recherche/Simulations/Modules/"
** Scenario 9: J = 4 items / M = 2 modalities / DIF size 0.3 x2 ** Scenario 9: J = 4 items / M = 2 modalities / DIF size 0.3 x2
local N = "100 200 300" local N = "100 200 300"
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/sphere/Corentin/These/Recherche/Simulations/Data/DIF/N`Nn'"
local path_res = "/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N`Nn'" local path_res = "/home/sphere/Corentin/These/Recherche/Simulations/Analysis/DIF/N`Nn'"
local scenarios = "A B C D E F G" local scenarios = "A B C D E F G"
foreach scen in `scenarios' { foreach scen in `scenarios' {
clear clear
@ -375,8 +377,8 @@ putexcel A1=matrix(outmat), colnames
local N = "100 200 300" local N = "100 200 300"
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/sphere/Corentin/These/Recherche/Simulations/Data/DIF/N`Nn'"
local path_res = "/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N`Nn'" local path_res = "/home/sphere/Corentin/These/Recherche/Simulations/Analysis/DIF/N`Nn'"
local scenarios = "A B C D E F G" local scenarios = "A B C D E F G"
foreach scen in `scenarios' { foreach scen in `scenarios' {
clear clear
@ -427,7 +429,7 @@ forvalues k=1/1000 {
local mod = "`mod'"+"(1.item`i'<-THETA@1)(2.item`i'<-THETA@2)(3.item`i'<-THETA@3)" local mod = "`mod'"+"(1.item`i'<-THETA@1)(2.item`i'<-THETA@2)(3.item`i'<-THETA@3)"
} }
} }
local mod = "`mod'" + "(THETA<-tt), mlogit tol(0.01) iterate(500) latent(THETA) nocapslatent constraint(1 2)" local mod = "`mod'" + "(THETA<-tt), mlogit tol(0.01) iterate(500) latent(THETA) nocapslatent constraint(1 2 3 4)"
qui `constrnt' qui `constrnt'
qui `constrnt2' qui `constrnt2'
qui `constrnt3' qui `constrnt3'
@ -485,8 +487,8 @@ putexcel A1=matrix(outmat), colnames
local N = "100 200 300" local N = "100 200 300"
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/sphere/Corentin/These/Recherche/Simulations/Data/DIF/N`Nn'"
local path_res = "/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N`Nn'" local path_res = "/home/sphere/Corentin/These/Recherche/Simulations/Analysis/DIF/N`Nn'"
local scenarios = "A B C D E F G" local scenarios = "A B C D E F G"
foreach scen in `scenarios' { foreach scen in `scenarios' {
clear clear
@ -559,8 +561,8 @@ putexcel A1=matrix(outmat), colnames
local N = "100 200 300" local N = "100 200 300"
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/sphere/Corentin/These/Recherche/Simulations/Data/DIF/N`Nn'"
local path_res = "/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N`Nn'" local path_res = "/home/sphere/Corentin/These/Recherche/Simulations/Analysis/DIF/N`Nn'"
local scenarios = "A B C D E F G" local scenarios = "A B C D E F G"
foreach scen in `scenarios' { foreach scen in `scenarios' {
clear clear
@ -611,7 +613,7 @@ forvalues k=1/1000 {
local mod = "`mod'"+"(1.item`i'<-THETA@1)(2.item`i'<-THETA@2)(3.item`i'<-THETA@3)" local mod = "`mod'"+"(1.item`i'<-THETA@1)(2.item`i'<-THETA@2)(3.item`i'<-THETA@3)"
} }
} }
local mod = "`mod'" + "(THETA<-tt), mlogit tol(0.01) iterate(500) latent(THETA) nocapslatent constraint(1 2)" local mod = "`mod'" + "(THETA<-tt), mlogit tol(0.01) iterate(500) latent(THETA) nocapslatent constraint(1 2 3 4)"
qui `constrnt' qui `constrnt'
qui `constrnt2' qui `constrnt2'
qui `constrnt3' qui `constrnt3'
@ -666,19 +668,21 @@ putexcel A1=matrix(outmat), colnames
*========================== *==========================
* Scenarios with : J=4 * Scenarios with : J=7
*========================== *==========================
****** Scenarios with DIF on 2 items ****** Scenarios with DIF on 2 items
* Load pcm.ado
adopath+"/home/sphere/Corentin/These/Recherche/Simulations/Modules/"
** Scenario 13: J = 7 items / M = 2 modalities / DIF size 0.3 x2 ** Scenario 13: J = 7 items / M = 2 modalities / DIF size 0.3 x2
local N = "100 200 300" local N = "100 200 300"
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/sphere/Corentin/These/Recherche/Simulations/Data/DIF/N`Nn'"
local path_res = "/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N`Nn'" local path_res = "/home/sphere/Corentin/These/Recherche/Simulations/Analysis/DIF/N`Nn'"
local scenarios = "A B C D E F G" local scenarios = "A B C D E F G"
foreach scen in `scenarios' { foreach scen in `scenarios' {
clear clear
@ -751,8 +755,8 @@ putexcel A1=matrix(outmat), colnames
local N = "100 200 300" local N = "100 200 300"
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/sphere/Corentin/These/Recherche/Simulations/Data/DIF/N`Nn'"
local path_res = "/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N`Nn'" local path_res = "/home/sphere/Corentin/These/Recherche/Simulations/Analysis/DIF/N`Nn'"
local scenarios = "A B C D E F G" local scenarios = "A B C D E F G"
foreach scen in `scenarios' { foreach scen in `scenarios' {
clear clear
@ -803,7 +807,7 @@ forvalues k=1/1000 {
local mod = "`mod'"+"(1.item`i'<-THETA@1)(2.item`i'<-THETA@2)(3.item`i'<-THETA@3)" local mod = "`mod'"+"(1.item`i'<-THETA@1)(2.item`i'<-THETA@2)(3.item`i'<-THETA@3)"
} }
} }
local mod = "`mod'" + "(THETA<-tt), mlogit tol(0.01) iterate(500) latent(THETA) nocapslatent constraint(1 2)" local mod = "`mod'" + "(THETA<-tt), mlogit tol(0.01) iterate(500) latent(THETA) nocapslatent constraint(1 2 3 4)"
qui `constrnt' qui `constrnt'
qui `constrnt2' qui `constrnt2'
qui `constrnt3' qui `constrnt3'
@ -861,8 +865,8 @@ putexcel A1=matrix(outmat), colnames
local N = "100 200 300" local N = "100 200 300"
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/sphere/Corentin/These/Recherche/Simulations/Data/DIF/N`Nn'"
local path_res = "/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N`Nn'" local path_res = "/home/sphere/Corentin/These/Recherche/Simulations/Analysis/DIF/N`Nn'"
local scenarios = "A B C D E F G" local scenarios = "A B C D E F G"
foreach scen in `scenarios' { foreach scen in `scenarios' {
clear clear
@ -935,8 +939,8 @@ putexcel A1=matrix(outmat), colnames
local N = "100 200 300" local N = "100 200 300"
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/sphere/Corentin/These/Recherche/Simulations/Data/DIF/N`Nn'"
local path_res = "/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N`Nn'" local path_res = "/home/sphere/Corentin/These/Recherche/Simulations/Analysis/DIF/N`Nn'"
local scenarios = "A B C D E F G" local scenarios = "A B C D E F G"
foreach scen in `scenarios' { foreach scen in `scenarios' {
clear clear
@ -987,7 +991,7 @@ forvalues k=1/1000 {
local mod = "`mod'"+"(1.item`i'<-THETA@1)(2.item`i'<-THETA@2)(3.item`i'<-THETA@3)" local mod = "`mod'"+"(1.item`i'<-THETA@1)(2.item`i'<-THETA@2)(3.item`i'<-THETA@3)"
} }
} }
local mod = "`mod'" + "(THETA<-tt), mlogit tol(0.01) iterate(500) latent(THETA) nocapslatent constraint(1 2)" local mod = "`mod'" + "(THETA<-tt), mlogit tol(0.01) iterate(500) latent(THETA) nocapslatent constraint(1 2 3 4)"
qui `constrnt' qui `constrnt'
qui `constrnt2' qui `constrnt2'
qui `constrnt3' qui `constrnt3'
@ -1041,13 +1045,15 @@ putexcel A1=matrix(outmat), colnames
***** Scenarios with DIF on 3 items ***** Scenarios with DIF on 3 items
* Load pcm.ado
adopath+"/home/sphere/Corentin/These/Recherche/Simulations/Modules/"
** Scenario 17: J = 7 items / M = 2 modalities / DIF size 0.3 x3 ** Scenario 17: J = 7 items / M = 2 modalities / DIF size 0.3 x3
local N = "100 200 300" local N = "100 200 300"
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/sphere/Corentin/These/Recherche/Simulations/Data/DIF/N`Nn'"
local path_res = "/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N`Nn'" local path_res = "/home/sphere/Corentin/These/Recherche/Simulations/Analysis/DIF/N`Nn'"
local scenarios = "A B C D E F G" local scenarios = "A B C D E F G"
foreach scen in `scenarios' { foreach scen in `scenarios' {
clear clear
@ -1137,8 +1143,8 @@ putexcel A1=matrix(outmat), colnames
local N = "100 200 300" local N = "100 200 300"
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/sphere/Corentin/These/Recherche/Simulations/Data/DIF/N`Nn'"
local path_res = "/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N`Nn'" local path_res = "/home/sphere/Corentin/These/Recherche/Simulations/Analysis/DIF/N`Nn'"
local scenarios = "A B C D E F G" local scenarios = "A B C D E F G"
foreach scen in `scenarios' { foreach scen in `scenarios' {
clear clear
@ -1207,7 +1213,7 @@ forvalues k=1/1000 {
local mod = "`mod'"+"(1.item`i'<-THETA@1)(2.item`i'<-THETA@2)(3.item`i'<-THETA@3)" local mod = "`mod'"+"(1.item`i'<-THETA@1)(2.item`i'<-THETA@2)(3.item`i'<-THETA@3)"
} }
} }
local mod = "`mod'" + "(THETA<-tt), mlogit tol(0.01) iterate(500) latent(THETA) nocapslatent constraint(1 2)" local mod = "`mod'" + "(THETA<-tt), mlogit tol(0.01) iterate(500) latent(THETA) nocapslatent constraint(1 2 3 4 5 6)"
qui `constrnt' qui `constrnt'
qui `constrnt2' qui `constrnt2'
qui `constrnt3' qui `constrnt3'
@ -1279,8 +1285,8 @@ putexcel A1=matrix(outmat), colnames
local N = "100 200 300" local N = "100 200 300"
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/sphere/Corentin/These/Recherche/Simulations/Data/DIF/N`Nn'"
local path_res = "/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N`Nn'" local path_res = "/home/sphere/Corentin/These/Recherche/Simulations/Analysis/DIF/N`Nn'"
local scenarios = "A B C D E F G" local scenarios = "A B C D E F G"
foreach scen in `scenarios' { foreach scen in `scenarios' {
clear clear
@ -1370,8 +1376,8 @@ putexcel A1=matrix(outmat), colnames
local N = "100 200 300" local N = "100 200 300"
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/sphere/Corentin/These/Recherche/Simulations/Data/DIF/N`Nn'"
local path_res = "/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N`Nn'" local path_res = "/home/sphere/Corentin/These/Recherche/Simulations/Analysis/DIF/N`Nn'"
local scenarios = "A B C D E F G" local scenarios = "A B C D E F G"
foreach scen in `scenarios' { foreach scen in `scenarios' {
clear clear
@ -1440,7 +1446,7 @@ forvalues k=1/1000 {
local mod = "`mod'"+"(1.item`i'<-THETA@1)(2.item`i'<-THETA@2)(3.item`i'<-THETA@3)" local mod = "`mod'"+"(1.item`i'<-THETA@1)(2.item`i'<-THETA@2)(3.item`i'<-THETA@3)"
} }
} }
local mod = "`mod'" + "(THETA<-tt), mlogit tol(0.01) iterate(500) latent(THETA) nocapslatent constraint(1 2)" local mod = "`mod'" + "(THETA<-tt), mlogit tol(0.01) iterate(500) latent(THETA) nocapslatent constraint(1 2 3 4 5 6)"
qui `constrnt' qui `constrnt'
qui `constrnt2' qui `constrnt2'
qui `constrnt3' qui `constrnt3'