temp fix for scenario 10B

RProject/Scripts/aggregation.R

@@ -1198,7 +1198,16 @@ res.dat[res.dat$N==50,]$dif.size <- sapply(which(res.dat$N==50),function(k) res.
 res.dat.dif[res.dat.dif$N==50,]$dif.size <- sapply(which(res.dat.dif$N==50),function(k) res.dat.dif[k-1,]$dif.size)
 res.dat.dif.rosali[res.dat.dif.rosali$N==50,]$dif.size <- sapply(which(res.dat.dif.rosali$N==50),function(k) res.dat.dif.rosali[k-1,]$dif.size)
 res.dat.dif.resali[res.dat.dif.resali$N==50,]$dif.size <- sapply(which(res.dat.dif.resali$N==50),function(k) res.dat.dif.resali[k-1,]$dif.size)
+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"
 
 # Items dichotomiques
 
@@ -1210,11 +1219,11 @@ res.dat.dif.resali$method <- "RESIDUS"
 res.dat.dicho <- res.dat[res.dat$M==2,]
 res.dat.dicho <- rbind(res.dat.dicho,res.dat.dif[res.dat.dif$M==2,])
 res.dat.dicho <- rbind.fill(res.dat.dicho,res.dat.dif.rosali[res.dat.dif.rosali$M==2,])
-res.dat.dicho <- rbind(res.dat.dicho,res.dat.dif.resali[res.dat.dif.resali$M==2,])
+res.dat.dicho <- rbind.fill(res.dat.dicho,res.dat.dif.resali[res.dat.dif.resali$M==2,])
 
 # Items polytomiques
 
 res.dat.poly <- res.dat[res.dat$M==4,]
 res.dat.poly <- rbind(res.dat.poly,res.dat.dif[res.dat.dif$M==4,])
 res.dat.poly <- rbind.fill(res.dat.poly,res.dat.dif.rosali[res.dat.dif.rosali$M==4,])
-res.dat.poly <- rbind(res.dat.poly,res.dat.dif.resali[res.dat.dif.resali$M==4,])
+res.dat.poly <- rbind.fill(res.dat.poly,res.dat.dif.resali[res.dat.dif.resali$M==4,])