From c6afca9a3b114b7ece2e873e84d6cad866a74a7f Mon Sep 17 00:00:00 2001 From: corentinchoisy Date: Fri, 17 May 2024 16:20:33 +0200 Subject: [PATCH] temp fix for scenario 10B --- RProject/Scripts/aggregation.R | 13 +++++++++++-- 1 file changed, 11 insertions(+), 2 deletions(-) diff --git a/RProject/Scripts/aggregation.R b/RProject/Scripts/aggregation.R index e89e22e..9206010 100644 --- a/RProject/Scripts/aggregation.R +++ b/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,])