ROSALI and RESALI scripts updated

Scripts/Analysis/DIF-RESIDUS/pcm_dif_residus.do

@@ -342,7 +342,7 @@ local N = "50 100 200 300"
         }
         mat outmat[`k',colnumb(outmat,"dif_detect_`u'")] = `j'
         if (mod(`s',2)==0) {
-          mat outmat[`k',colnumb(outmat,"dif_detect_unif_`u'")] = `unif_`j''
+          mat outmat[`k',colnumb(outmat,"dif_detect_unif_`u'")] = `unif_`u''
         }
       }
     * Stocker les items de DIF originaux

Scripts/Analysis/DIF-ROSALI/pcm_dif_rosali.do

@@ -61,6 +61,9 @@ local N = "50 100 200 300"
   }
   * taillemat = Maximum J*M cases pour les items par et J*M cases pour les dif par + J cases pour les DIF detect + nbdif cases pour dif réel
   local taillemat=2*`nbitems'*`nbmoda'+`nbitems'+`nbdif'+2
+  if (mod(`s',2)==0) {
+    local taillemat=2*`nbitems'*`nbmoda'+`nbitems'+`nbitems'+`nbdif'+2
+  }
   local colna=""
   forvalues i=1/`nbitems' {
     forvalues z=1/`nbmoda' {
@@ -69,8 +72,13 @@ local N = "50 100 200 300"
     }
   }
   forvalues i=1/`nbitems' {
+    if (mod(`s',2)==1) {
       local colna = "`colna'"+"dif_detect_`i' "
     }
+    if (mod(`s',2)==0) {
+      local colna = "`colna'"+"dif_detect_`i' "+"dif_detect_unif_`i' "
+    }
+  }
 
   forvalues i=1/`nbdif' {
     local colna = "`colna'"+"real_dif_`i' "
@@ -213,7 +221,7 @@ local N = "50 100 200 300"
           local checker2: word 2 of `val'
           local checker3: word 3 of `val'
           local nummoda=r(r)
-          local nbmoda_`j'=`nummoda'
+          local nbmoda_`j'=`nummoda'-1
           if (`nummoda'==2) {
             qui recode item`j' (`checker'=0) (`checker2'=1)
           }
@@ -301,11 +309,14 @@ local N = "50 100 200 300"
           local mod = "`mod'"+"(2.item`v'<-THETA@2 tt)"
         }
         local w= 100+`u'
+        local unif_`u'=0
         if (resmat[1,`nbitems'+`v']!=testm[1,1] & `nbmoda_`v''==3) {
           local conformula = "`conformula'" + "`u' " + "`w' "
+          local unif_`u'=1
         }
         else if (resmat[1,`nbitems'+`v']!=testm[1,1] & `nbmoda_`v''==2) {
           local conformula = "`conformula'" + "`u' "
+          local unif_`u'=1
         }
     }
     if ("`conformula'" != "") {
@@ -330,7 +341,10 @@ local N = "50 100 200 300"
       local j=`difitems`u''
         forvalues z=1/`nbmoda_`j'' {
           mat outmat[`k',colnumb(outmat,"dif_`u'_`z'")] = W[1,colnumb(W,"`z'.item`j':tt")]
+        }
         mat outmat[`k',colnumb(outmat,"dif_detect_`u'")] = `j'
+        if (mod(`s',2)==0) {
+          mat outmat[`k',colnumb(outmat,"dif_detect_unif_`u'")] = `unif_`u''
         }
       }