diff --git a/RProject/Scripts/aggregation.R b/RProject/Scripts/aggregation.R index fcd9141..857eba7 100644 --- a/RProject/Scripts/aggregation.R +++ b/RProject/Scripts/aggregation.R @@ -264,6 +264,14 @@ res.dat.simple <- res.dat[,c(1:8,13,16:18)] res.dat.simple$m.beta <- round(res.dat.simple$m.beta,3) res.dat.simple +is.nan.data.frame <- function(x) + do.call(cbind, lapply(x, is.nan)) + + +res.dat[is.nan(res.dat)] <- NA +res.dat$bias <- res.dat$eff.size-res.dat$m.beta + + ############################################################################## #----------------------------------------------------------------------------# ########################### AGGREGATION DIF MATRICES ######################### diff --git a/RProject/Scripts/pcm_nodif.R b/RProject/Scripts/pcm_nodif.R index 8fb36e9..9cd3e30 100644 --- a/RProject/Scripts/pcm_nodif.R +++ b/RProject/Scripts/pcm_nodif.R @@ -43,7 +43,6 @@ replicate_pcm_analysis_m4 <- function(df=NULL,treatment='TT',irtmodel='PCM2',met truebeta <- eff.size if (method=='MML') { n <- max(df[,sequence]) - print(n) tam1 <- lapply(seq(1,n), function(x) pcm_analysis(df=df[df[,sequence]==x,],treatment=treatment,irtmodel=irtmodel) ) @@ -91,7 +90,6 @@ replicate_pcm_analysis_m2 <- function(df=NULL,treatment='TT',irtmodel='PCM2',met resp <- df[,sapply(seq(1,nbitems),function(x) paste0('item',x))] if (method=='MML') { n <- max(df[,sequence]) - print(n) tam1 <- lapply(seq(1,n), function(x) pcm_analysis(df=df[df[,sequence]==x,],treatment=treatment,irtmodel=irtmodel) ) @@ -2057,412 +2055,322 @@ write.csv(res_nodif[[3]],'/home/corentin/Documents/These/Recherche/Simulations/A registerDoMC(5) -###################################### Scenario 1 - - -dat1 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_1A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_1B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_1C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_1D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_1E_50.csv') - -res <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5'),function(x) replicate_pcm_analysis(get(x))) +###################################### Scenario A + + +dat1 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_1A_50.csv') +dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_2A_50.csv') +dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_3A_50.csv') +dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_4A_50.csv') +dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_5A_50.csv') +dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_6A_50.csv') +dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_7A_50.csv') +dat8 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_8A_50.csv') +dat9 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_9A_50.csv') +dat10 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_10A_50.csv') +dat11 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_11A_50.csv') +dat12 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_12A_50.csv') +dat13 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_13A_50.csv') +dat14 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_14A_50.csv') +dat15 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_15A_50.csv') +dat16 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_16A_50.csv') +dat17 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_17A_50.csv') +dat18 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_18A_50.csv') +dat19 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_19A_50.csv') +dat20 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_20A_50.csv') + +res <- pbmclapply(paste0("dat",1:20),function(x) replicate_pcm_analysis(get(x),eff.size = 0)) write.csv(res[[1]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_1A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_1B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_1C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_1D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_1E_50.csv') - - -###################################### Scenario 2 - - -dat1 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_2A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_2B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_2C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_2D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_2E_50.csv') - -res <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5'),function(x) replicate_pcm_analysis(get(x))) - -write.csv(res[[1]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_2A_50.csv') +write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_2A_50.csv') +write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_3A_50.csv') +write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_4A_50.csv') +write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_5A_50.csv') +write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_6A_50.csv') +write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_7A_50.csv') +write.csv(res[[8]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_8A_50.csv') +write.csv(res[[9]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_9A_50.csv') +write.csv(res[[10]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_10A_50.csv') +write.csv(res[[11]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_11A_50.csv') +write.csv(res[[12]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_12A_50.csv') +write.csv(res[[13]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_13A_50.csv') +write.csv(res[[14]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_14A_50.csv') +write.csv(res[[15]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_15A_50.csv') +write.csv(res[[16]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_16A_50.csv') +write.csv(res[[17]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_17A_50.csv') +write.csv(res[[18]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_18A_50.csv') +write.csv(res[[19]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_19A_50.csv') +write.csv(res[[20]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_20A_50.csv') + + +###################################### Scenario B + + +dat1 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_1B_50.csv') +dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_2B_50.csv') +dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_3B_50.csv') +dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_4B_50.csv') +dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_5B_50.csv') +dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_6B_50.csv') +dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_7B_50.csv') +dat8 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_8B_50.csv') +dat9 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_9B_50.csv') +dat10 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_10B_50.csv') +dat11 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_11B_50.csv') +dat12 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_12B_50.csv') +dat13 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_13B_50.csv') +dat14 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_14B_50.csv') +dat15 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_15B_50.csv') +dat16 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_16B_50.csv') +dat17 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_17B_50.csv') +dat18 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_18B_50.csv') +dat19 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_19B_50.csv') +dat20 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_20B_50.csv') + +res <- pbmclapply(paste0("dat",1:20),function(x) replicate_pcm_analysis(get(x),eff.size = 0.2)) + +write.csv(res[[1]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_1B_50.csv') write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_2B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_2C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_2D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_2E_50.csv') - - -###################################### Scenario 3 - - -dat1 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_3A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_3B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_3C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_3D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_3E_50.csv') - -res <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5'),function(x) replicate_pcm_analysis(get(x))) - -write.csv(res[[1]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_3A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_3B_50.csv') +write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_3B_50.csv') +write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_4B_50.csv') +write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_5B_50.csv') +write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_6B_50.csv') +write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_7B_50.csv') +write.csv(res[[8]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_8B_50.csv') +write.csv(res[[9]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_9B_50.csv') +write.csv(res[[10]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_10B_50.csv') +write.csv(res[[11]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_11B_50.csv') +write.csv(res[[12]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_12B_50.csv') +write.csv(res[[13]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_13B_50.csv') +write.csv(res[[14]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_14B_50.csv') +write.csv(res[[15]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_15B_50.csv') +write.csv(res[[16]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_16B_50.csv') +write.csv(res[[17]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_17B_50.csv') +write.csv(res[[18]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_18B_50.csv') +write.csv(res[[19]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_19B_50.csv') +write.csv(res[[20]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_20B_50.csv') + + +###################################### Scenario C + + +dat1 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_1C_50.csv') +dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_2C_50.csv') +dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_3C_50.csv') +dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_4C_50.csv') +dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_5C_50.csv') +dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_6C_50.csv') +dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_7C_50.csv') +dat8 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_8C_50.csv') +dat9 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_9C_50.csv') +dat10 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_10C_50.csv') +dat11 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_11C_50.csv') +dat12 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_12C_50.csv') +dat13 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_13C_50.csv') +dat14 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_14C_50.csv') +dat15 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_15C_50.csv') +dat16 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_16C_50.csv') +dat17 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_17C_50.csv') +dat18 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_18C_50.csv') +dat19 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_19C_50.csv') +dat20 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_20C_50.csv') + +res <- pbmclapply(1:20,function(x) replicate_pcm_analysis(get(paste0("dat",x)),eff.size = ifelse(x>4,0.2,0.4))) + +write.csv(res[[1]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_1C_50.csv') +write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_2C_50.csv') write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_3C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_3D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_3E_50.csv') - -###################################### Scenario 4 - - -dat1 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_4A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_4B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_4C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_4D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_4E_50.csv') - -res <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5'),function(x) replicate_pcm_analysis(get(x))) - -write.csv(res[[1]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_4A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_4B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_4C_50.csv') +write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_4C_50.csv') +write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_5C_50.csv') +write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_6C_50.csv') +write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_7C_50.csv') +write.csv(res[[8]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_8C_50.csv') +write.csv(res[[9]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_9C_50.csv') +write.csv(res[[10]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_10C_50.csv') +write.csv(res[[11]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_11C_50.csv') +write.csv(res[[12]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_12C_50.csv') +write.csv(res[[13]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_13C_50.csv') +write.csv(res[[14]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_14C_50.csv') +write.csv(res[[15]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_15C_50.csv') +write.csv(res[[16]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_16C_50.csv') +write.csv(res[[17]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_17C_50.csv') +write.csv(res[[18]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_18C_50.csv') +write.csv(res[[19]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_19C_50.csv') +write.csv(res[[20]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_20C_50.csv') + + + + +###################################### Scenario D + + +dat1 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_1D_50.csv') +dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_2D_50.csv') +dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_3D_50.csv') +dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_4D_50.csv') +dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_5D_50.csv') +dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_6D_50.csv') +dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_7D_50.csv') +dat8 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_8D_50.csv') +dat9 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_9D_50.csv') +dat10 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_10D_50.csv') +dat11 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_11D_50.csv') +dat12 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_12D_50.csv') +dat13 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_13D_50.csv') +dat14 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_14D_50.csv') +dat15 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_15D_50.csv') +dat16 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_16D_50.csv') +dat17 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_17D_50.csv') +dat18 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_18D_50.csv') +dat19 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_19D_50.csv') +dat20 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_20D_50.csv') + +res <- pbmclapply(1:20,function(x) replicate_pcm_analysis(get(paste0("dat",x)),eff.size = ifelse(x>4,0.4,-0.2))) + +write.csv(res[[1]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_1D_50.csv') +write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_2D_50.csv') +write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_3D_50.csv') write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_4D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_4E_50.csv') - - - -###################################### Scenario 5 - - -dat1 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_5A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_5B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_5C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_5D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_5E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_5F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_5G_50.csv') - -res <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 1,difsize = 0.3)) - -write.csv(res[[1]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_5A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_5B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_5C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_5D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_5E_50.csv') -write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_5F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_5G_50.csv') - -###################################### Scenario 6 - - -dat1 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_6A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_6B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_6C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_6D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_6E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_6F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_6G_50.csv') - -res <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 1,difsize = 0.3)) - -write.csv(res[[1]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_6A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_6B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_6C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_6D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_6E_50.csv') -write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_6F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_6G_50.csv') - -###################################### Scenario 7 - - -dat1 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_7A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_7B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_7C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_7D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_7E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_7F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_7G_50.csv') - -res <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 1,difsize = 0.5)) - -write.csv(res[[1]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_7A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_7B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_7C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_7D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_7E_50.csv') -write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_7F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_7G_50.csv') - -###################################### Scenario 8 - - -dat1 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_8A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_8B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_8C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_8D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_8E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_8F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_8G_50.csv') - -res <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 1,difsize = 0.5)) - -write.csv(res[[1]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_8A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_8B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_8C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_8D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_8E_50.csv') -write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_8F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_8G_50.csv') - - -###################################### Scenario 9 - - -dat1 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_9A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_9B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_9C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_9D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_9E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_9F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_9G_50.csv') - -res <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 2,difsize = 0.3)) - -write.csv(res[[1]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_9A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_9B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_9C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_9D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_9E_50.csv') -write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_9F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_9G_50.csv') - -###################################### Scenario 10 - - -dat1 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_10A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_10B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_10C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_10D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_10E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_10F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_10G_50.csv') - -res <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 2,difsize = 0.3)) - -write.csv(res[[1]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_10A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_10B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_10C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_10D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_10E_50.csv') -write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_10F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_10G_50.csv') - -###################################### Scenario 11 - - -dat1 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_11A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_11B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_11C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_11D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_11E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_11F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_11G_50.csv') - -res <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 2,difsize = 0.5)) - -write.csv(res[[1]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_11A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_11B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_11C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_11D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_11E_50.csv') -write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_11F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_11G_50.csv') - -###################################### Scenario 12 - - -dat1 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_12A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_12B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_12C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_12D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_12E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_12F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_12G_50.csv') - -res <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 2,difsize = 0.5)) - -write.csv(res[[1]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_12A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_12B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_12C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_12D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_12E_50.csv') -write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_12F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_12G_50.csv') - - -###################################### Scenario 13 - - -dat1 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_13A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_13B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_13C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_13D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_13E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_13F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_13G_50.csv') +write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_5D_50.csv') +write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_6D_50.csv') +write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_7D_50.csv') +write.csv(res[[8]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_8D_50.csv') +write.csv(res[[9]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_9D_50.csv') +write.csv(res[[10]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_10D_50.csv') +write.csv(res[[11]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_11D_50.csv') +write.csv(res[[12]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_12D_50.csv') +write.csv(res[[13]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_13D_50.csv') +write.csv(res[[14]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_14D_50.csv') +write.csv(res[[15]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_15D_50.csv') +write.csv(res[[16]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_16D_50.csv') +write.csv(res[[17]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_17D_50.csv') +write.csv(res[[18]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_18D_50.csv') +write.csv(res[[19]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_19D_50.csv') +write.csv(res[[20]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_20D_50.csv') + + + +###################################### Scenario E + + +dat1 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_1E_50.csv') +dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_2E_50.csv') +dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_3E_50.csv') +dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_4E_50.csv') +dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_5E_50.csv') +dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_6E_50.csv') +dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_7E_50.csv') +dat8 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_8E_50.csv') +dat9 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_9E_50.csv') +dat10 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_10E_50.csv') +dat11 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_11E_50.csv') +dat12 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_12E_50.csv') +dat13 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_13E_50.csv') +dat14 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_14E_50.csv') +dat15 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_15E_50.csv') +dat16 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_16E_50.csv') +dat17 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_17E_50.csv') +dat18 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_18E_50.csv') +dat19 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_19E_50.csv') +dat20 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_20E_50.csv') + +res <- pbmclapply(1:20,function(x) replicate_pcm_analysis(get(paste0("dat",x)),eff.size = ifelse(x>4,0.4,-0.4))) + +write.csv(res[[1]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_1E_50.csv') +write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_2E_50.csv') +write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_3E_50.csv') +write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_4E_50.csv') +write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_5E_50.csv') +write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_6E_50.csv') +write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_7E_50.csv') +write.csv(res[[8]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_8E_50.csv') +write.csv(res[[9]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_9E_50.csv') +write.csv(res[[10]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_10E_50.csv') +write.csv(res[[11]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_11E_50.csv') +write.csv(res[[12]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_12E_50.csv') +write.csv(res[[13]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_13E_50.csv') +write.csv(res[[14]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_14E_50.csv') +write.csv(res[[15]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_15E_50.csv') +write.csv(res[[16]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_16E_50.csv') +write.csv(res[[17]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_17E_50.csv') +write.csv(res[[18]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_18E_50.csv') +write.csv(res[[19]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_19E_50.csv') +write.csv(res[[20]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_20E_50.csv') + +###################################### Scenario F + +dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_5F_50.csv') +dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_6F_50.csv') +dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_7F_50.csv') +dat8 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_8F_50.csv') +dat9 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_9F_50.csv') +dat10 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_10F_50.csv') +dat11 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_11F_50.csv') +dat12 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_12F_50.csv') +dat13 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_13F_50.csv') +dat14 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_14F_50.csv') +dat15 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_15F_50.csv') +dat16 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_16F_50.csv') +dat17 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_17F_50.csv') +dat18 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_18F_50.csv') +dat19 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_19F_50.csv') +dat20 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_20F_50.csv') + +res <- pbmclapply(5:20,function(x) replicate_pcm_analysis(get(paste0("dat",x)),eff.size = -0.2)) + +write.csv(res[[1]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_5F_50.csv') +write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_6F_50.csv') +write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_7F_50.csv') +write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_8F_50.csv') +write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_9F_50.csv') +write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_10F_50.csv') +write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_11F_50.csv') +write.csv(res[[8]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_12F_50.csv') +write.csv(res[[9]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_13F_50.csv') +write.csv(res[[10]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_14F_50.csv') +write.csv(res[[11]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_15F_50.csv') +write.csv(res[[12]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_16F_50.csv') +write.csv(res[[13]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_17F_50.csv') +write.csv(res[[14]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_18F_50.csv') +write.csv(res[[15]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_19F_50.csv') +write.csv(res[[16]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_20F_50.csv') + +###################################### Scenario G + +dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_5G_50.csv') +dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_6G_50.csv') +dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_7G_50.csv') +dat8 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_8G_50.csv') +dat9 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_9G_50.csv') +dat10 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_10G_50.csv') +dat11 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_11G_50.csv') +dat12 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_12G_50.csv') +dat13 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_13G_50.csv') +dat14 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_14G_50.csv') +dat15 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_15G_50.csv') +dat16 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_16G_50.csv') +dat17 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_17G_50.csv') +dat18 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_18G_50.csv') +dat19 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_19G_50.csv') +dat20 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/N50/scenario_20G_50.csv') + +res <- pbmclapply(5:20,function(x) replicate_pcm_analysis(get(paste0("dat",x)),eff.size = -0.4)) + +write.csv(res[[1]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_5G_50.csv') +write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_6G_50.csv') +write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_7G_50.csv') +write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_8G_50.csv') +write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_9G_50.csv') +write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_10G_50.csv') +write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_11G_50.csv') +write.csv(res[[8]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_12G_50.csv') +write.csv(res[[9]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_13G_50.csv') +write.csv(res[[10]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_14G_50.csv') +write.csv(res[[11]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_15G_50.csv') +write.csv(res[[12]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_16G_50.csv') +write.csv(res[[13]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_17G_50.csv') +write.csv(res[[14]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_18G_50.csv') +write.csv(res[[15]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_19G_50.csv') +write.csv(res[[16]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_20G_50.csv') -res <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 2,difsize = 0.3)) -write.csv(res[[1]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_13A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_13B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_13C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_13D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_13E_50.csv') -write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_13F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_13G_50.csv') - -###################################### Scenario 14 - - -dat1 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_14A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_14B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_14C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_14D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_14E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_14F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_14G_50.csv') - -res <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 2,difsize = 0.3)) - -write.csv(res[[1]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_14A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_14B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_14C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_14D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_14E_50.csv') -write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_14F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_14G_50.csv') - -###################################### Scenario 15 - - -dat1 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_15A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_15B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_15C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_15D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_15E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_15F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_15G_50.csv') - -res <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 2,difsize = 0.5)) - -write.csv(res[[1]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_15A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_15B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_15C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_15D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_15E_50.csv') -write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_15F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_15G_50.csv') - -###################################### Scenario 16 - - -dat1 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_16A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_16B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_16C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_16D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_16E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_16F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_16G_50.csv') - -res <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 2,difsize = 0.5)) - -write.csv(res[[1]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_16A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_16B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_16C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_16D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_16E_50.csv') -write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_16F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_16G_50.csv') - -###################################### Scenario 17 - - -dat1 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_17A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_17B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_17C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_17D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_17E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_17F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_17G_50.csv') - -res <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 3,difsize = 0.3)) - -write.csv(res[[1]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_17A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_17B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_17C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_17D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_17E_50.csv') -write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_17F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_17G_50.csv') - -###################################### Scenario 18 - - -dat1 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_18A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_18B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_18C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_18D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_18E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_18F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_18G_50.csv') - -res <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 3,difsize = 0.3)) - -write.csv(res[[1]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_18A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_18B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_18C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_18D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_18E_50.csv') -write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_18F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_18G_50.csv') - -###################################### Scenario 19 - - -dat1 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_19A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_19B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_19C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_19D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_19E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_19F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_19G_50.csv') - -res <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 3,difsize = 0.5)) - -write.csv(res[[1]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_19A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_19B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_19C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_19D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_19E_50.csv') -write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_19F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_19G_50.csv') - -###################################### Scenario 20 - - -dat1 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_20A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_20B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_20C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_20D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_20E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_20F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N50/scenario_20G_50.csv') - -res <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 3,difsize = 0.5)) - -write.csv(res[[1]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_20A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_20B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_20C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_20D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_20E_50.csv') -write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_20F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/scenario_20G_50.csv') diff --git a/RProject/shinyapps/results-dashboard/app.R b/RProject/shinyapps/results-dashboard/app.R new file mode 100644 index 0000000..8c20d10 --- /dev/null +++ b/RProject/shinyapps/results-dashboard/app.R @@ -0,0 +1,265 @@ +library(shiny) +library(shinydashboard) +library(shinyBS) +library(shinyWidgets) + +conditionalPanel2 <- function(condition, ..., ns = NS(NULL), inline = FALSE, container = if (inline) span else div) { + container(`data-display-if` = condition, `data-ns-prefix` = ns(""), ...) +} +conditionalMenuItem <- function(label,tabName,iconName,condition) { + t <- menuItem( + label, + tabName = tabName, + icon = icon(iconName) + ) + t$attribs <- list('data-display-if'=condition,'data-ns-prefix'='') + t +} + +# Define UI for application that draws a histogram +ui <- + dashboardPage( + dashboardHeader(title = "ResCenter", + dropdownMenu(type = "notifications", + notificationItem( + text = "Return to services list", + icon = icon("right-from-bracket"), + href = 'https://service.corentinchoisy.xyz', + status = 'danger'), + badgeStatus = NULL, + icon = icon('right-from-bracket'), + headerText = '' + ), + dropdownMenu(type='messages', + messageItem(from = 'View source on Gitea', + message = 'Browse and download', + icon = icon('git'), + href = 'https://gitlab.com/corentinchoisy/kids'), + badgeStatus = NULL, + icon = icon('circle-info'), + headerText = "App info") + ), + + dashboardSidebar( + collapsed=TRUE, + sidebarMenu( + menuItem("Tables", tabName = "tables", icon = icon("table")) + #,conditionalMenuItem(condition='input.predictions == 1',label = "Autres métriques", tabName = "metriques", iconName = "square-poll-vertical") + ) + ), + + dashboardBody( + tabItems(tabItem('tables', + fluidRow( + box(title = "Input",status="success",width=12, + pickerInput( + inputId = "typeres", + label = "DIF recovery", + options=list(style='btn'), + choices = c("None", "Perfect", "ROSALI", "Residuals")), + conditionalPanel2(condition='input.typeres == "ROSALI" | input.typeres == "Residuals"',checkboxInput(inputId = "weak",label = "Show only weak detection scenarios (p.moreflexible<0.15)") + ,conditionalPanel2(condition="!input.dif",checkboxInput(inputId = "causal",label = "Hide DIF detection data")), + conditionalPanel2(condition='!input.causal',checkboxInput(inputId = "dif",label = "Hide causal inference data")) + ) + ) + ), + fluidRow( + box(title = "Compiled results",status="primary",width = 12,solidHeader = T,DT::dataTableOutput("res")) + ) + ) + ) + + + + + ,tags$head(tags$style(HTML(' + .content-wrapper { overflow: auto; } + + #options .box-header{ display: none} + + button.btn.dropdown-toggle.btn-default { + background-color: #d32927; + border-color: #d32927; + border-radius: 0px; + } + + .box.box-danger{ + background:#f5eceb + } + + .box.box-warning{ + background:#f7f4ed + } + + .box.box-danger .form-control{ + background-color:#f5eceb !important + } + + .box.box-success{ + background:#e1fff1 + } + + .box.box-success .form-control{ + background-color:#e1fff1 !important + } + + .box.box-primary{ + background:#ebf4f9 + } + + + .box.box-solid{ + background:#ffffff + } + + .navbar-nav > .messages-menu > .dropdown-menu > li .menu, .navbar-nav > .notifications-menu > .dropdown-menu > li .menu, .navbar-nav > .tasks-menu > .dropdown-menu > li .menu { + max-height: 500px; overflow-y: hidden; + } + + .content { + height: auto; overflow-y: auto; + } + + /* navbar (rest of the header) */ + .skin-blue .main-header .navbar { + background-color: #d32927; + } + + /* navbar (rest of the header) */ + .skin-blue .main-header .logo { + background-color: #d32927; + } + + /* logo when hovered */ + .skin-blue .main-header .logo:hover { + background-color: #931c1b; + } + /* toggle button when hovered */ + .skin-blue .main-header .navbar .sidebar-toggle:hover{ + background-color: #931c1b; + } + + /* main sidebar */ + .skin-blue .main-sidebar { + background-color: #4e4c4c; + } + + /* active selected tab in the sidebarmenu */ + .skin-blue .main-sidebar .sidebar .sidebar-menu .active a{ + background-color: #282828; + border-left: 3px solid #d32927; + } + + /* other links in the sidebarmenu when hovered */ + .skin-blue .main-sidebar .sidebar .sidebar-menu a:hover{ + background-color: #282828; + border-left: 3px solid #d32927; + } + ')))) + ) + + + + + + + + + + + + + + +# Define server logic required to draw a histogram +server <- function(input, output) { + res.dat <- read.csv("res_dat.csv") + res.dat.dif <- read.csv("res_dat_dif.csv") + res.dat.dif.resali <- read.csv("res_dat_dif_resali.csv") + res.dat.dif.rosali <- read.csv("res_dat_dif_rosali.csv") + + df_res <- reactive({ + if(input$typeres=="None") { + dd <- res.dat + } + if(input$typeres=="Perfect") { + dd <- res.dat.dif + } + if(input$typeres=="ROSALI") { + dd <- res.dat.dif.rosali + } + if(input$typeres=="Residuals") { + dd <- res.dat.dif.resali + } + dd + }) + + + + df_res_r <- reactive({ + rs <- df_res() + if(input$typeres=="None") { + rs <- rs[,c('scenario','N','J','M','eff.size',"nb.dif","dif.size","m.beta","se.empirical.beta","se.analytical.beta", + "m.low.ci.beta","m.high.ci.beta","true.value.in.ci.p","h0.rejected.p","beta.same.sign.truebeta.p", + "beta.same.sign.truebeta.signif.p","bias")] + rs <- dplyr::mutate_if(rs,is.numeric,round,3) + rs <- dplyr::mutate_all(rs,as.character) + } + if(input$typeres=="Perfect") { + rs <- rs[,c('scenario','N','J','M','eff.size',"nb.dif","dif.size","m.beta","se.empirical.beta","se.analytical.beta", + "m.low.ci.beta","m.high.ci.beta","true.value.in.ci.p","h0.rejected.p","beta.same.sign.truebeta.p", + "beta.same.sign.truebeta.signif.p","bias")] + rs <- dplyr::mutate_if(rs,is.numeric,round,3) + rs <- dplyr::mutate_all(rs,as.character) + } + if(input$typeres=="ROSALI") { + rs <- rs[,c('scenario','N','J','M','eff.size',"nb.dif","dif.size","m.beta","se.empirical.beta","se.analytical.beta", + "m.low.ci.beta","m.high.ci.beta","true.value.in.ci.p","h0.rejected.p","beta.same.sign.truebeta.p", + "beta.same.sign.truebeta.signif.p","bias","dif.detected","prop.perfect","flexible.detect","moreflexible.detect","any.detect")] + if (input$causal) { + rs <- rs[,c('scenario','N','J','M','eff.size',"nb.dif","dif.size","m.beta","se.empirical.beta","se.analytical.beta", + "m.low.ci.beta","m.high.ci.beta","true.value.in.ci.p","h0.rejected.p","beta.same.sign.truebeta.p", + "beta.same.sign.truebeta.signif.p","bias")] + } + if (input$dif) { + rs <- rs[,c('scenario','N','J','M','eff.size',"nb.dif","dif.size","dif.detected","prop.perfect","flexible.detect","moreflexible.detect","any.detect")] + } + if (input$weak) { + rs <- rs[rs$moreflexible.detect<0.15 & !is.na(rs$moreflexible.detect),] + } + + rs <- dplyr::mutate_if(rs,is.numeric,round,3) + rs <- dplyr::mutate_all(rs,as.character) + } + if(input$typeres=="Residuals") { + rs <- rs[,c('scenario','N','J','M','eff.size',"nb.dif","dif.size","m.beta","se.empirical.beta","se.analytical.beta", + "m.low.ci.beta","m.high.ci.beta","true.value.in.ci.p","h0.rejected.p","beta.same.sign.truebeta.p", + "beta.same.sign.truebeta.signif.p","bias","dif.detected","prop.perfect","flexible.detect","moreflexible.detect","any.detect")] + if (input$causal) { + rs <- rs[,c('scenario','N','J','M','eff.size',"nb.dif","dif.size","m.beta","se.empirical.beta","se.analytical.beta", + "m.low.ci.beta","m.high.ci.beta","true.value.in.ci.p","h0.rejected.p","beta.same.sign.truebeta.p", + "beta.same.sign.truebeta.signif.p","bias")] + } + if (input$dif) { + rs <- rs[,c('scenario','N','J','M','eff.size',"nb.dif","dif.size","dif.detected","prop.perfect","flexible.detect","moreflexible.detect","any.detect")] + } + if (input$weak) { + rs <- rs[rs$moreflexible.detect<0.15 & !is.na(rs$moreflexible.detect),] + } + rs <- dplyr::mutate_if(rs,is.numeric,round,3) + rs <- dplyr::mutate_all(rs,as.character) + } + rs + }) + + + output$res <- DT::renderDataTable(df_res_r(), + options = list(columnDefs = list(list(className = 'dt-center',width=ifelse(input$dif & !is.na(input$dif),"101px","50px"), targets = "_all")), + pageLength=100,autoWidth=T,scrollY=400,scrollX = T,server=F,buttons=c("csv","excel")), + rownames = FALSE, + selection="single", + filter=list(plain=F,clear=F,position="top")) +} + +# Run the application +shinyApp(ui = ui, server = server)