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

Corrected pcmdif.do error

Browse Source
This commit is contained in:
Corentin Choisy
2024-02-05 12:07:20 +01:00
parent ef1fbfbd4f
commit d5baf77bc1
5 changed files with 762 additions and 754 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@ -296,6 +296,8 @@ putexcel A1=matrix(outmat), colnames
****** Scenarios with DIF on 2 items
* Load pcm.ado
adopath+"/home/corentin/Documents/These/Recherche/Simulations/Modules/"
** Scenario 9: J = 4 items / M = 2 modalities / DIF size 0.3 x2
local N = "100 200 300"
@ -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'" + "(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 `constrnt2'
qui `constrnt3'
@ -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'" + "(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 `constrnt2'
qui `constrnt3'
@ -666,12 +668,14 @@ putexcel A1=matrix(outmat), colnames
*==========================
* Scenarios with : J=4
* Scenarios with : J=7
*==========================
****** Scenarios with DIF on 2 items
* Load pcm.ado
adopath+"/home/corentin/Documents/These/Recherche/Simulations/Modules/"
** Scenario 13: J = 7 items / M = 2 modalities / DIF size 0.3 x2
local N = "100 200 300"
@ -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'" + "(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 `constrnt2'
qui `constrnt3'
@ -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'" + "(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 `constrnt2'
qui `constrnt3'
@ -1041,6 +1045,8 @@ putexcel A1=matrix(outmat), colnames
***** Scenarios with DIF on 3 items
* Load pcm.ado
adopath+"/home/corentin/Documents/These/Recherche/Simulations/Modules/"
** Scenario 17: J = 7 items / M = 2 modalities / DIF size 0.3 x3
local N = "100 200 300"
@ -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'" + "(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 `constrnt2'
qui `constrnt3'
@ -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'" + "(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 `constrnt2'
qui `constrnt3'