From 413b10ae29a9a6c11fc6f2fafd72cdaaa5b4fe5d Mon Sep 17 00:00:00 2001
From: corentinchoisy <corentin.choisy@proton.me>
Date: Fri, 17 May 2024 17:09:21 +0200
Subject: [PATCH] Corrected data generation bug for scenario 10B

---
 RProject/Scripts/aggregation.R                       | 12 ++++++------
 .../Scenarios/DIF/scenarios_DIF_J4_baseline_100.do   |  2 +-
 .../Scenarios/DIF/scenarios_DIF_J4_baseline_200.do   |  2 +-
 .../Scenarios/DIF/scenarios_DIF_J4_baseline_300.do   |  2 +-
 .../Scenarios/DIF/scenarios_DIF_J4_baseline_50.do    |  2 +-
 5 files changed, 10 insertions(+), 10 deletions(-)

diff --git a/RProject/Scripts/aggregation.R b/RProject/Scripts/aggregation.R
index 9206010..f1fb458 100644
--- a/RProject/Scripts/aggregation.R
+++ b/RProject/Scripts/aggregation.R
@@ -1202,12 +1202,12 @@ res.dat[res.dat$dif.size!=0 & res.dat$nb.dif==0,]$nb.dif <- 1
 res.dat.dif <- res.dat.dif %>%
   relocate(method, .after = theoretical.power)
 
-aa <- res.dat.dif[res.dat.dif$scenario=="10B",]
-bb <- res.dat[res.dat$scenario=="10B",]
-res.dat.dif[res.dat.dif$scenario=="10B",] <- bb
-res.dat[res.dat$scenario=="10B",] <- aa
-res.dat[res.dat$method=="PERFECT",]$method <- "NONE"
-res.dat.dif[res.dat.dif$method=="NONE",]$method <- "PERFECT"
+#aa <- res.dat.dif[res.dat.dif$scenario=="10B",]
+#bb <- res.dat[res.dat$scenario=="10B",]
+#res.dat.dif[res.dat.dif$scenario=="10B",] <- bb
+#res.dat[res.dat$scenario=="10B",] <- aa
+#res.dat[res.dat$method=="PERFECT",]$method <- "NONE"
+#res.dat.dif[res.dat.dif$method=="NONE",]$method <- "PERFECT"
 
 # Items dichotomiques
 
diff --git a/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_100.do b/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_100.do
index 8385e51..dc1bd12 100644
--- a/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_100.do
+++ b/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_100.do
@@ -2291,7 +2291,7 @@ forvalues replication = 1/1000 {
     qui local difi :  word 1 of `ItemsDIF'
     qui local difj :  word 2 of `ItemsDIF'
 
-      mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
+      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
diff --git a/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_200.do b/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_200.do
index 789cf71..b0e48a0 100644
--- a/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_200.do
+++ b/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_200.do
@@ -2291,7 +2291,7 @@ forvalues replication = 1/1000 {
     qui local difi :  word 1 of `ItemsDIF'
     qui local difj :  word 2 of `ItemsDIF'
 
-      mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
+      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
diff --git a/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_300.do b/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_300.do
index fefb4fe..fe2d7a6 100644
--- a/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_300.do
+++ b/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_300.do
@@ -2291,7 +2291,7 @@ forvalues replication = 1/1000 {
     qui local difi :  word 1 of `ItemsDIF'
     qui local difj :  word 2 of `ItemsDIF'
 
-      mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
+      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
diff --git a/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_50.do b/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_50.do
index 09d2152..a83c506 100644
--- a/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_50.do
+++ b/Scripts/Scenarios/DIF/scenarios_DIF_J4_baseline_50.do
@@ -2291,7 +2291,7 @@ forvalues replication = 1/1000 {
     qui local difi :  word 1 of `ItemsDIF'
     qui local difj :  word 2 of `ItemsDIF'
 
-      mat B= (-1.84,-0.84,0.16 \ -1.25,-0.25,0.75 \ -0.75,0.25,1.25 \ 0.16,0.84,1.84)
+      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