From 100ce17ead4ce2b01deb8dc7563a0aaeaba5f48a Mon Sep 17 00:00:00 2001 From: corentinchoisy Date: Thu, 29 Feb 2024 18:18:07 +0100 Subject: [PATCH] Corrected bug in S18 and S20, n=50 --- RProject/Scripts/pcm.R | 502 +++++++++++++------------- Scripts/Analysis/DIF/pcm_dif_50_J7.do | 20 +- 2 files changed, 262 insertions(+), 260 deletions(-) diff --git a/RProject/Scripts/pcm.R b/RProject/Scripts/pcm.R index 2c5deb5..74cce92 100644 --- a/RProject/Scripts/pcm.R +++ b/RProject/Scripts/pcm.R @@ -2132,340 +2132,340 @@ write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysi ###################################### Scenario 5 -dat1 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_5A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_5B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_5C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_5D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_5E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_5F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_5G_50.csv') - -res_nodif <- 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/NoDIF/N50/scenario_5A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_5B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_5C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_5D_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_5F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_5G_50.csv') +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/NoDIF/N50/scenario_6A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_6B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_6C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_6D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_6E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_6F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_6G_50.csv') +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_nodif <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 1,difsize = 0.3)) +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/NoDIF/N50/scenario_6A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_6B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_6C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_6D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_6E_50.csv') -write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_6F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_6G_50.csv') +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/NoDIF/N50/scenario_7A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_7B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_7C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_7D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_7E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_7F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_7G_50.csv') +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_nodif <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 1,difsize = 0.5)) +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/NoDIF/N50/scenario_7A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_7B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_7C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_7D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_7E_50.csv') -write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_7F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_7G_50.csv') +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/NoDIF/N50/scenario_8A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_8B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_8C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_8D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_8E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_8F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_8G_50.csv') +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_nodif <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 1,difsize = 0.5)) +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/NoDIF/N50/scenario_8A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_8B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_8C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_8D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_8E_50.csv') -write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_8F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_8G_50.csv') +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/NoDIF/N50/scenario_9A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_9B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_9C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_9D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_9E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_9F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_9G_50.csv') +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_nodif <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 2,difsize = 0.3)) +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/NoDIF/N50/scenario_9A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_9B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_9C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_9D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_9E_50.csv') -write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_9F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_9G_50.csv') +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/NoDIF/N50/scenario_10A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_10B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_10C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_10D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_10E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_10F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_10G_50.csv') +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_nodif <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 2,difsize = 0.3)) +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/NoDIF/N50/scenario_10A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_10B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_10C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_10D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_10E_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_10G_50.csv') +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/NoDIF/N50/scenario_11A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_11B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_11C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_11D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_11E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_11F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_11G_50.csv') +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_nodif <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 2,difsize = 0.5)) +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/NoDIF/N50/scenario_11A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_11B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_11C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_11D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_11E_50.csv') -write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_11F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_11G_50.csv') +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/NoDIF/N50/scenario_12A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_12B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_12C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_12D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_12E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_12F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_12G_50.csv') +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_nodif <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 2,difsize = 0.5)) +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/NoDIF/N50/scenario_12A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_12B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_12C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_12D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_12E_50.csv') -write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_12F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_12G_50.csv') +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/NoDIF/N50/scenario_13A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_13B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_13C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_13D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_13E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_13F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_13G_50.csv') +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') -res_nodif <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 2,difsize = 0.3)) +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/NoDIF/N50/scenario_13A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_13B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_13C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_13D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_13E_50.csv') -write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_13F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_13G_50.csv') +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/NoDIF/N50/scenario_14A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_14B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_14C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_14D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_14E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_14F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_14G_50.csv') +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_nodif <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 2,difsize = 0.3)) +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/NoDIF/N50/scenario_14A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_14B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_14C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_14D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_14E_50.csv') -write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_14F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_14G_50.csv') +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/NoDIF/N50/scenario_15A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_15B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_15C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_15D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_15E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_15F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_15G_50.csv') +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_nodif <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 2,difsize = 0.5)) +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/NoDIF/N50/scenario_15A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_15B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_15C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_15D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_15E_50.csv') -write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_15F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_15G_50.csv') +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/NoDIF/N50/scenario_16A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_16B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_16C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_16D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_16E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_16F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_16G_50.csv') +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_nodif <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 2,difsize = 0.5)) +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/NoDIF/N50/scenario_16A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_16B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_16C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_16D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_16E_50.csv') -write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_16F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_16G_50.csv') +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/NoDIF/N50/scenario_17A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_17B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_17C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_17D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_17E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_17F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_17G_50.csv') +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_nodif <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 2,difsize = 0.3)) +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/NoDIF/N50/scenario_17A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_17B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_17C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_17D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_17E_50.csv') -write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_17F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_17G_50.csv') +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/NoDIF/N50/scenario_18A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_18B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_18C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_18D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_18E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_18F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_18G_50.csv') +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_nodif <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 2,difsize = 0.3)) +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/NoDIF/N50/scenario_18A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_18B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_18C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_18D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_18E_50.csv') -write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_18F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_18G_50.csv') +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/NoDIF/N50/scenario_19A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_19B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_19C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_19D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_19E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_19F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_19G_50.csv') +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_nodif <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 2,difsize = 0.5)) +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/NoDIF/N50/scenario_19A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_19B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_19C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_19D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_19E_50.csv') -write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_19F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_19G_50.csv') +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/NoDIF/N50/scenario_20A_50.csv') -dat2 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_20B_50.csv') -dat3 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_20C_50.csv') -dat4 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_20D_50.csv') -dat5 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_20E_50.csv') -dat6 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_20F_50.csv') -dat7 <- read.csv(file = '/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N50/scenario_20G_50.csv') +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_nodif <- pbmclapply(c('dat1','dat2','dat3','dat4','dat5','dat6','dat7'),function(x) replicate_pcm_analysis(get(x),nbdif = 2,difsize = 0.5)) +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/NoDIF/N50/scenario_20A_50.csv') -write.csv(res[[2]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_20B_50.csv') -write.csv(res[[3]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_20C_50.csv') -write.csv(res[[4]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_20D_50.csv') -write.csv(res[[5]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_20E_50.csv') -write.csv(res[[6]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_20F_50.csv') -write.csv(res[[7]],'/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_20G_50.csv') +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') ############################################################################## #----------------------------------------------------------------------------# @@ -2493,11 +2493,11 @@ results <- c(results,results2) compile_simulation <- function(scenario) { name <- as.numeric(gsub("[^0-9.-]", "", substr(scenario,start=0,stop=2))) - if (substr(scenario,start=nchar(scenario)-2,stop=nchar(scenario))=="50" & name<=4) { + if (substr(scenario,start=nchar(scenario)-1,stop=nchar(scenario))=="50" & name<=4) { s <- read.csv(paste0('/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_',scenario,'.csv')) } - if (substr(scenario,start=nchar(scenario)-2,stop=nchar(scenario))=="50" & name>4) { - s <- read.csv(paste0('/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_',scenario,'_nodif.csv')) + if (substr(scenario,start=nchar(scenario)-1,stop=nchar(scenario))=="50" & name>4) { + s <- read.csv(paste0('/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N50/scenario_',scenario,'.csv')) } if (substr(scenario,start=nchar(scenario)-2,stop=nchar(scenario))=="100" & name<=4) { s <- read.csv(paste0('/home/corentin/Documents/These/Recherche/Simulations/Analysis/NoDIF/N100/scenario_',scenario,'.csv')) @@ -2542,10 +2542,11 @@ compile_simulation <- function(scenario) { ) } } - zz <- substr(scenario,start=0,stop=nchar(scenario)-4) + N <- ifelse(substr(scenario,start=nchar(scenario)-1,stop=nchar(scenario))=="50","50",substr(scenario,start=nchar(scenario)-2,stop=nchar(scenario))) + zz <- ifelse(N=="50",substr(scenario,start=0,stop=nchar(scenario)-3),substr(scenario,start=0,stop=nchar(scenario)-4)) b <- data.frame(scenario=zz, scenario.type=substr(zz,start=nchar(zz),stop=nchar(zz)), - N=substr(scenario,start=nchar(scenario)-2,stop=nchar(scenario)), + N=N, J=unique(s$J), M=unique(s$M), eff.size=unique(s$eff.size), @@ -2611,7 +2612,7 @@ results <- c(results,results2)[81:528] compile_simulation2 <- function(scenario) { name <- as.numeric(gsub("[^0-9.-]", "", substr(scenario,start=0,stop=2))) - if (substr(scenario,start=nchar(scenario)-2,stop=nchar(scenario))=="50" & name>4) { + if (substr(scenario,start=nchar(scenario)-1,stop=nchar(scenario))=="50" & name>4) { s <- read_excel(paste0('/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N50/',scenario,'.xls')) } if (substr(scenario,start=nchar(scenario)-2,stop=nchar(scenario))=="100" & name>4) { @@ -2623,6 +2624,7 @@ compile_simulation2 <- function(scenario) { if (substr(scenario,start=nchar(scenario)-2,stop=nchar(scenario))=="300" & name>4) { s <- read_excel(paste0('/home/corentin/Documents/These/Recherche/Simulations/Analysis/DIF/N300/',scenario,'.xls')) } + s <- s[!is.na(s$beta),] J <- max(which(sapply(1:7,function(x) paste0('item',x) %in% colnames(s) | paste0('item',x,'_1') %in% colnames(s)))) M <- 1+sum(sapply(1:3,function(x) paste0('item1_',x) %in% colnames(s) )) if (M==1) {M <- 2} @@ -2652,8 +2654,8 @@ compile_simulation2 <- function(scenario) { ) } } - zz <- substr(scenario,start=0,stop=nchar(scenario)-4) - N <- substr(scenario,start=nchar(scenario)-2,stop=nchar(scenario)) + N <- ifelse(substr(scenario,start=nchar(scenario)-1,stop=nchar(scenario))=="50","50",substr(scenario,start=nchar(scenario)-2,stop=nchar(scenario))) + zz <- ifelse(N=="50",substr(scenario,start=0,stop=nchar(scenario)-3),substr(scenario,start=0,stop=nchar(scenario)-4)) eff.size <- unique(res.dat[res.dat$scenario==zz & res.dat$N==N,'eff.size']) dif.size <- unique(res.dat[res.dat$scenario==zz & res.dat$N==N,'dif.size']) b <- data.frame(scenario=zz, diff --git a/Scripts/Analysis/DIF/pcm_dif_50_J7.do b/Scripts/Analysis/DIF/pcm_dif_50_J7.do index 6ed4a2f..bd73ae7 100644 --- a/Scripts/Analysis/DIF/pcm_dif_50_J7.do +++ b/Scripts/Analysis/DIF/pcm_dif_50_J7.do @@ -86,7 +86,7 @@ adopath+"/home/corentin/Documents/These/Recherche/Simulations/Modules/" mat outmat[`k',`nbitems'+6] = `difitemsmax' // numéro item de dif restore } - putexcel set "`path_res'/out/13`scen'_`Nn'.xls", sheet("outmat") replace + putexcel set "`path_res'/13`scen'_`Nn'.xls", sheet("outmat") replace putexcel A1=matrix(outmat), colnames } } @@ -440,7 +440,7 @@ forvalues j=1/`nbitems' { restore } } -putexcel set "`path_res'/out/14`scen'_`Nn'.xls", sheet("outmat") replace +putexcel set "`path_res'/14`scen'_`Nn'.xls", sheet("outmat") replace putexcel A1=matrix(outmat), colnames } } @@ -515,7 +515,7 @@ putexcel A1=matrix(outmat), colnames mat outmat[`k',`nbitems'+6] = `difitemsmax' // numéro item de dif restore } - putexcel set "`path_res'/out/15`scen'_`Nn'.xls", sheet("outmat") replace + putexcel set "`path_res'/15`scen'_`Nn'.xls", sheet("outmat") replace putexcel A1=matrix(outmat), colnames } } @@ -870,7 +870,7 @@ forvalues j=1/`nbitems' { restore } } -putexcel set "`path_res'/out/16`scen'_`Nn'.xls", sheet("outmat") replace +putexcel set "`path_res'/16`scen'_`Nn'.xls", sheet("outmat") replace putexcel A1=matrix(outmat), colnames } } @@ -965,7 +965,7 @@ adopath+"/home/corentin/Documents/These/Recherche/Simulations/Modules/" mat outmat[`k',`nbitems'+8] = `difitemsmax' // numéro item de dif restore } - putexcel set "`path_res'/out/17`scen'_`Nn'.xls", sheet("outmat") replace + putexcel set "`path_res'/17`scen'_`Nn'.xls", sheet("outmat") replace putexcel A1=matrix(outmat), colnames } } @@ -1254,7 +1254,7 @@ if (`difcheck'==1) { } } } -else if (`difcheck'==1) { +else if (`difcheck'==0) { if (`difcheck2'==1) { if (`difcheck3'==1) { local mod = "`mod'" + "(THETA<-tt), mlogit tol(0.01) iterate(500) latent(THETA) nocapslatent constraint(1 2 4 5 6)" @@ -1483,7 +1483,7 @@ else { restore } } -putexcel set "`path_res'/out/18`scen'_`Nn'.xls", sheet("outmat") replace +putexcel set "`path_res'/18`scen'_`Nn'.xls", sheet("outmat") replace putexcel A1=matrix(outmat), colnames } } @@ -1573,7 +1573,7 @@ putexcel A1=matrix(outmat), colnames mat outmat[`k',`nbitems'+8] = `difitemsmax' // numéro item de dif restore } - putexcel set "`path_res'/out/19`scen'_`Nn'.xls", sheet("outmat") replace + putexcel set "`path_res'/19`scen'_`Nn'.xls", sheet("outmat") replace putexcel A1=matrix(outmat), colnames } } @@ -1862,7 +1862,7 @@ if (`difcheck'==1) { } } } -else if (`difcheck'==1) { +else if (`difcheck'==0) { if (`difcheck2'==1) { if (`difcheck3'==1) { local mod = "`mod'" + "(THETA<-tt), mlogit tol(0.01) iterate(500) latent(THETA) nocapslatent constraint(1 2 4 5 6)" @@ -2091,7 +2091,7 @@ else { restore } } -putexcel set "`path_res'/out/20`scen'_`Nn'.xls", sheet("outmat") replace +putexcel set "`path_res'/20`scen'_`Nn'.xls", sheet("outmat") replace putexcel A1=matrix(outmat), colnames } }