wwwwww <- read.csv('/home/corentin/Documents/These/Recherche/Simulations/Analysis/RESALI/Detection_data/3A_200.csv')
table(wwwwww$dif.detect.1)
table(wwwwww$dif.detect.1)/400
sum(table(wwwwww$dif.detect.1)/400)
results <- c(sapply(1:4,function(x) paste0(x,c('A','B','C','D','E'))),sapply(5:9,function(x) paste0(x,c('A','B','C','D','E','F','G'))))
results2 <- c(sapply(10:20,function(x) paste0(x,c('A','B','C','D','E','F','G'))))
results <- c(sapply(c(50,100,200,300),function(x) paste0(results,'_',x)))
results2 <- c(sapply(c(50,100,200,300),function(x) paste0(results2,'_',x)))
results <- sort(results)
results2 <- sort(results2)
results <- c(results,results2)
## Importer l'analyse resali pour chaque scenario
for (r in results[473:length(results)]) {
cat('--------------------------------------------------------------------------\n')
cat(paste0(r,"\n"))
cat('--------------------------------------------------------------------------\n')
#### Importer les datas
scen <- as.numeric(gsub("[A,B,C,D,E,F,G,_]","",substr(r,0,3)))
if (substr(r,start=nchar(r)-1,stop=nchar(r))=="50") {
N <- 50
}
if (substr(r,start=nchar(r)-2,stop=nchar(r))=="100") {
N <- 100
}
if (substr(r,start=nchar(r)-2,stop=nchar(r))=="200") {
N <- 200
}
if (substr(r,start=nchar(r)-2,stop=nchar(r))=="300") {
N <- 300
}
if (scen<5) {
datt <- read.csv(paste0('/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N',N,'/scenario_',r,'.csv'))
}
if (scen>=5) {
datt <- read.csv(paste0('/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N',N,'/scenario_',r,'.csv'))
}
#### Importer l'analyse
analyse <- read.csv(paste0('/home/corentin/Documents/These/Recherche/Simulations/Analysis/RESALI/Detection/',r,".csv"))
#### Pour chaque replication
for (k in 1:1000) {
if (k%%100==0) {
cat(paste0("N = ",k," / 1000\n"))
}
datt[datt$replication==k,"dif.detect.1"] <- analyse[analyse$X==k,"dif.detect.1"]
datt[datt$replication==k,"dif.detect.2"] <- analyse[analyse$X==k,"dif.detect.2"]
datt[datt$replication==k,"dif.detect.3"] <- analyse[analyse$X==k,"dif.detect.3"]
datt[datt$replication==k,"dif.detect.4"] <- analyse[analyse$X==k,"dif.detect.4"]
datt[datt$replication==k,"dif.detect.unif.1"] <- analyse[analyse$X==k,"dif.detect.unif.1"]
datt[datt$replication==k,"dif.detect.unif.2"] <- analyse[analyse$X==k,"dif.detect.unif.2"]
datt[datt$replication==k,"dif.detect.unif.3"] <- analyse[analyse$X==k,"dif.detect.unif.3"]
datt[datt$replication==k,"dif.detect.unif.4"] <- analyse[analyse$X==k,"dif.detect.unif.4"]
if (scen==3 | scen==4 | scen>=13) {
datt[datt$replication==k,"dif.detect.5"] <- analyse[analyse$X==k,"dif.detect.5"]
datt[datt$replication==k,"dif.detect.6"] <- analyse[analyse$X==k,"dif.detect.6"]
datt[datt$replication==k,"dif.detect.7"] <- analyse[analyse$X==k,"dif.detect.7"]
datt[datt$replication==k,"dif.detect.unif.5"] <- analyse[analyse$X==k,"dif.detect.unif.5"]
datt[datt$replication==k,"dif.detect.unif.6"] <- analyse[analyse$X==k,"dif.detect.unif.6"]
datt[datt$replication==k,"dif.detect.unif.7"] <- analyse[analyse$X==k,"dif.detect.unif.7"]
}
}
datt[is.na(datt$dif.detect.1),"dif.detect.1"] <- ""
datt[is.na(datt$dif.detect.2),"dif.detect.2"] <- ""
datt[is.na(datt$dif.detect.3),"dif.detect.3"] <- ""
datt[is.na(datt$dif.detect.4),"dif.detect.4"] <- ""
datt[is.na(datt$dif.detect.unif.1),"dif.detect.unif.1"] <- ""
datt[is.na(datt$dif.detect.unif.2),"dif.detect.unif.2"] <- ""
datt[is.na(datt$dif.detect.unif.3),"dif.detect.unif.3"] <- ""
datt[is.na(datt$dif.detect.unif.4),"dif.detect.unif.4"] <- ""
if (scen==3 | scen==4 | scen>=13) {
datt[is.na(datt$dif.detect.5),"dif.detect.5"] <- ""
datt[is.na(datt$dif.detect.6),"dif.detect.6"] <- ""
datt[is.na(datt$dif.detect.7),"dif.detect.7"] <- ""
datt[is.na(datt$dif.detect.unif.5),"dif.detect.unif.5"] <- ""
datt[is.na(datt$dif.detect.unif.6),"dif.detect.unif.6"] <- ""
datt[is.na(datt$dif.detect.unif.7),"dif.detect.unif.7"] <- ""
}
write.csv(datt,paste0("/home/corentin/Documents/These/Recherche/Simulations/Analysis/RESALI/Detection_data/",r,".csv"))
}
results <- c(sapply(1:4,function(x) paste0(x,c('A','B','C','D','E'))),sapply(5:9,function(x) paste0(x,c('A','B','C','D','E','F','G'))))
results2 <- c(sapply(10:20,function(x) paste0(x,c('A','B','C','D','E','F','G'))))
results <- c(sapply(c(50,100,200,300),function(x) paste0(results,'_',x)))
results2 <- c(sapply(c(50,100,200,300),function(x) paste0(results2,'_',x)))
results <- sort(results)
results2 <- sort(results2)
results <- c(results,results2)
## Importer l'analyse resali pour chaque scenario
for (r in results) {
cat('--------------------------------------------------------------------------\n')
cat(paste0(r,"\n"))
cat('--------------------------------------------------------------------------\n')
#### Importer les datas
scen <- as.numeric(gsub("[A,B,C,D,E,F,G,_]","",substr(r,0,3)))
if (substr(r,start=nchar(r)-1,stop=nchar(r))=="50") {
N <- 50
}
if (substr(r,start=nchar(r)-2,stop=nchar(r))=="100") {
N <- 100
}
if (substr(r,start=nchar(r)-2,stop=nchar(r))=="200") {
N <- 200
}
if (substr(r,start=nchar(r)-2,stop=nchar(r))=="300") {
N <- 300
}
if (scen<5) {
datt <- read.csv(paste0('/home/corentin/Documents/These/Recherche/Simulations/Data/NoDIF/N',N,'/scenario_',r,'.csv'))
}
if (scen>=5) {
datt <- read.csv(paste0('/home/corentin/Documents/These/Recherche/Simulations/Data/DIF/N',N,'/scenario_',r,'.csv'))
}
#### Importer l'analyse
analyse <- read.csv(paste0('/home/corentin/Documents/These/Recherche/Simulations/Analysis/RESALI/Detection/',r,".csv"))
#### Pour chaque replication
for (k in 1:1000) {
if (k%%100==0) {
cat(paste0("N = ",k," / 1000\n"))
}
datt[datt$replication==k,"dif.detect.1"] <- analyse[analyse$X==k,"dif.detect.1"]
datt[datt$replication==k,"dif.detect.2"] <- analyse[analyse$X==k,"dif.detect.2"]
datt[datt$replication==k,"dif.detect.3"] <- analyse[analyse$X==k,"dif.detect.3"]
datt[datt$replication==k,"dif.detect.4"] <- analyse[analyse$X==k,"dif.detect.4"]
datt[datt$replication==k,"dif.detect.unif.1"] <- analyse[analyse$X==k,"dif.detect.unif.1"]
datt[datt$replication==k,"dif.detect.unif.2"] <- analyse[analyse$X==k,"dif.detect.unif.2"]
datt[datt$replication==k,"dif.detect.unif.3"] <- analyse[analyse$X==k,"dif.detect.unif.3"]
datt[datt$replication==k,"dif.detect.unif.4"] <- analyse[analyse$X==k,"dif.detect.unif.4"]
if (scen==3 | scen==4 | scen>=13) {
datt[datt$replication==k,"dif.detect.5"] <- analyse[analyse$X==k,"dif.detect.5"]
datt[datt$replication==k,"dif.detect.6"] <- analyse[analyse$X==k,"dif.detect.6"]
datt[datt$replication==k,"dif.detect.7"] <- analyse[analyse$X==k,"dif.detect.7"]
datt[datt$replication==k,"dif.detect.unif.5"] <- analyse[analyse$X==k,"dif.detect.unif.5"]
datt[datt$replication==k,"dif.detect.unif.6"] <- analyse[analyse$X==k,"dif.detect.unif.6"]
datt[datt$replication==k,"dif.detect.unif.7"] <- analyse[analyse$X==k,"dif.detect.unif.7"]
}
}
datt[is.na(datt$dif.detect.1),"dif.detect.1"] <- ""
datt[is.na(datt$dif.detect.2),"dif.detect.2"] <- ""
datt[is.na(datt$dif.detect.3),"dif.detect.3"] <- ""
datt[is.na(datt$dif.detect.4),"dif.detect.4"] <- ""
datt[is.na(datt$dif.detect.unif.1),"dif.detect.unif.1"] <- ""
datt[is.na(datt$dif.detect.unif.2),"dif.detect.unif.2"] <- ""
datt[is.na(datt$dif.detect.unif.3),"dif.detect.unif.3"] <- ""
datt[is.na(datt$dif.detect.unif.4),"dif.detect.unif.4"] <- ""
if (scen==3 | scen==4 | scen>=13) {
datt[is.na(datt$dif.detect.5),"dif.detect.5"] <- ""
datt[is.na(datt$dif.detect.6),"dif.detect.6"] <- ""
datt[is.na(datt$dif.detect.7),"dif.detect.7"] <- ""
datt[is.na(datt$dif.detect.unif.5),"dif.detect.unif.5"] <- ""
datt[is.na(datt$dif.detect.unif.6),"dif.detect.unif.6"] <- ""
datt[is.na(datt$dif.detect.unif.7),"dif.detect.unif.7"] <- ""
}
write.csv(datt,paste0("/home/corentin/Documents/These/Recherche/Simulations/Analysis/RESALI/Detection_data/",r,".csv"))
}
for (r in results) {
cat('--------------------------------------------------------------------------\n')
cat(paste0(r,"\n"))
cat('--------------------------------------------------------------------------\n')
#### Importer les datas
scen <- as.numeric(gsub("[A,B,C,D,E,F,G,_]","",substr(r,0,3)))
if (substr(r,start=nchar(r)-1,stop=nchar(r))=="50") {
N <- 50
}
if (substr(r,start=nchar(r)-2,stop=nchar(r))=="100") {
N <- 100
}
if (substr(r,start=nchar(r)-2,stop=nchar(r))=="200") {
N <- 200
}
if (substr(r,start=nchar(r)-2,stop=nchar(r))=="300") {
N <- 300
}
#### Importer l'analyse
analyse <- read.csv(paste0("/home/corentin/Documents/These/Recherche/Simulations/Analysis/RESALI/Detection_data/",r,".csv"))
names(analyse)[names(analyse)=="dif.detect.1"] <- "dif_detect_1"
names(analyse)[names(analyse)=="dif.detect.2"] <- "dif_detect_2"
names(analyse)[names(analyse)=="dif.detect.3"] <- "dif_detect_3"
names(analyse)[names(analyse)=="dif.detect.4"] <- "dif_detect_4"
names(analyse)[names(analyse)=="dif.detect.unif.1"] <- "dif_detect_unif_1"
names(analyse)[names(analyse)=="dif.detect.unif.2"] <- "dif_detect_unif_2"
names(analyse)[names(analyse)=="dif.detect.unif.3"] <- "dif_detect_unif_3"
names(analyse)[names(analyse)=="dif.detect.unif.4"] <- "dif_detect_unif_4"
analyse[analyse$dif_detect_1==999,"dif_detect_1"] <- ""
analyse[analyse$dif_detect_2==999,"dif_detect_2"] <- ""
analyse[analyse$dif_detect_3==999,"dif_detect_3"] <- ""
analyse[analyse$dif_detect_4==999,"dif_detect_4"] <- ""
analyse[analyse$dif_detect_unif_1==999,"dif_detect_unif_1"] <- ""
analyse[analyse$dif_detect_unif_2==999,"dif_detect_unif_2"] <- ""
analyse[analyse$dif_detect_unif_3==999,"dif_detect_unif_3"] <- ""
analyse[analyse$dif_detect_unif_4==999,"dif_detect_unif_4"] <- ""
if (scen==3 | scen==4 | scen>=13) {
names(analyse)[names(analyse)=="dif.detect.5"] <- "dif_detect_5"
names(analyse)[names(analyse)=="dif.detect.6"] <- "dif_detect_6"
names(analyse)[names(analyse)=="dif.detect.7"] <- "dif_detect_7"
names(analyse)[names(analyse)=="dif.detect.unif.5"] <- "dif_detect_unif_5"
names(analyse)[names(analyse)=="dif.detect.unif.6"] <- "dif_detect_unif_6"
names(analyse)[names(analyse)=="dif.detect.unif.7"] <- "dif_detect_unif_7"
analyse[analyse$dif_detect_5==999,"dif_detect_5"] <- ""
analyse[analyse$dif_detect_6==999,"dif_detect_6"] <- ""
analyse[analyse$dif_detect_7==999,"dif_detect_7"] <- ""
analyse[analyse$dif_detect_unif_5==999,"dif_detect_unif_5"] <- ""
analyse[analyse$dif_detect_unif_6==999,"dif_detect_unif_6"] <- ""
analyse[analyse$dif_detect_unif_7==999,"dif_detect_unif_7"] <- ""
}
analyse <- analyse[,!names(analyse) %in% c("X","X.1","X.2")]
write.csv(analyse,paste0("/home/corentin/Documents/These/Recherche/Simulations/Analysis/RESALI/Detection_data/",r,".csv"))
}
for (r in results[-c(1:4)]) {
cat('--------------------------------------------------------------------------\n')
cat(paste0(r,"\n"))
cat('--------------------------------------------------------------------------\n')
#### Importer les datas
scen <- as.numeric(gsub("[A,B,C,D,E,F,G,_]","",substr(r,0,3)))
if (substr(r,start=nchar(r)-1,stop=nchar(r))=="50") {
N <- 50
}
if (substr(r,start=nchar(r)-2,stop=nchar(r))=="100") {
N <- 100
}
if (substr(r,start=nchar(r)-2,stop=nchar(r))=="200") {
N <- 200
}
if (substr(r,start=nchar(r)-2,stop=nchar(r))=="300") {
N <- 300
}
#### Importer l'analyse
analyse <- read.csv(paste0("/home/corentin/Documents/These/Recherche/Simulations/Analysis/RESALI/Detection_data/",r,".csv"))
names(analyse)[names(analyse)=="dif.detect.1"] <- "dif_detect_1"
names(analyse)[names(analyse)=="dif.detect.2"] <- "dif_detect_2"
names(analyse)[names(analyse)=="dif.detect.3"] <- "dif_detect_3"
names(analyse)[names(analyse)=="dif.detect.4"] <- "dif_detect_4"
names(analyse)[names(analyse)=="dif.detect.unif.1"] <- "dif_detect_unif_1"
names(analyse)[names(analyse)=="dif.detect.unif.2"] <- "dif_detect_unif_2"
names(analyse)[names(analyse)=="dif.detect.unif.3"] <- "dif_detect_unif_3"
names(analyse)[names(analyse)=="dif.detect.unif.4"] <- "dif_detect_unif_4"
analyse[analyse$dif_detect_1==999,"dif_detect_1"] <- ""
analyse[analyse$dif_detect_2==999,"dif_detect_2"] <- ""
analyse[analyse$dif_detect_3==999,"dif_detect_3"] <- ""
analyse[analyse$dif_detect_4==999,"dif_detect_4"] <- ""
analyse[analyse$dif_detect_unif_1==999,"dif_detect_unif_1"] <- ""
analyse[analyse$dif_detect_unif_2==999,"dif_detect_unif_2"] <- ""
analyse[analyse$dif_detect_unif_3==999,"dif_detect_unif_3"] <- ""
analyse[analyse$dif_detect_unif_4==999,"dif_detect_unif_4"] <- ""
if (scen==3 | scen==4 | scen>=13) {
names(analyse)[names(analyse)=="dif.detect.5"] <- "dif_detect_5"
names(analyse)[names(analyse)=="dif.detect.6"] <- "dif_detect_6"
names(analyse)[names(analyse)=="dif.detect.7"] <- "dif_detect_7"
names(analyse)[names(analyse)=="dif.detect.unif.5"] <- "dif_detect_unif_5"
names(analyse)[names(analyse)=="dif.detect.unif.6"] <- "dif_detect_unif_6"
names(analyse)[names(analyse)=="dif.detect.unif.7"] <- "dif_detect_unif_7"
analyse[analyse$dif_detect_5==999,"dif_detect_5"] <- ""
analyse[analyse$dif_detect_6==999,"dif_detect_6"] <- ""
analyse[analyse$dif_detect_7==999,"dif_detect_7"] <- ""
analyse[analyse$dif_detect_unif_5==999,"dif_detect_unif_5"] <- ""
analyse[analyse$dif_detect_unif_6==999,"dif_detect_unif_6"] <- ""
analyse[analyse$dif_detect_unif_7==999,"dif_detect_unif_7"] <- ""
}
analyse <- analyse[,!names(analyse) %in% c("X","X.1","X.2")]
write.csv(analyse,paste0("/home/corentin/Documents/These/Recherche/Simulations/Analysis/RESALI/Detection_data/",r,".csv"))
}
for (r in results) {
cat('--------------------------------------------------------------------------\n')
cat(paste0(r,"\n"))
cat('--------------------------------------------------------------------------\n')
#### Importer les datas
scen <- as.numeric(gsub("[A,B,C,D,E,F,G,_]","",substr(r,0,3)))
if (substr(r,start=nchar(r)-1,stop=nchar(r))=="50") {
N <- 50
}
if (substr(r,start=nchar(r)-2,stop=nchar(r))=="100") {
N <- 100
}
if (substr(r,start=nchar(r)-2,stop=nchar(r))=="200") {
N <- 200
}
if (substr(r,start=nchar(r)-2,stop=nchar(r))=="300") {
N <- 300
}
#### Importer l'analyse
analyse <- read.csv(paste0("/home/corentin/Documents/These/Recherche/Simulations/Analysis/RESALI/Detection_data/",r,".csv"))
names(analyse)[names(analyse)=="dif.detect.1"] <- "dif_detect_1"
names(analyse)[names(analyse)=="dif.detect.2"] <- "dif_detect_2"
names(analyse)[names(analyse)=="dif.detect.3"] <- "dif_detect_3"
names(analyse)[names(analyse)=="dif.detect.4"] <- "dif_detect_4"
names(analyse)[names(analyse)=="dif.detect.unif.1"] <- "dif_detect_unif_1"
names(analyse)[names(analyse)=="dif.detect.unif.2"] <- "dif_detect_unif_2"
names(analyse)[names(analyse)=="dif.detect.unif.3"] <- "dif_detect_unif_3"
names(analyse)[names(analyse)=="dif.detect.unif.4"] <- "dif_detect_unif_4"
if (scen==3 | scen==4 | scen>=13) {
names(analyse)[names(analyse)=="dif.detect.5"] <- "dif_detect_5"
names(analyse)[names(analyse)=="dif.detect.6"] <- "dif_detect_6"
names(analyse)[names(analyse)=="dif.detect.7"] <- "dif_detect_7"
names(analyse)[names(analyse)=="dif.detect.unif.5"] <- "dif_detect_unif_5"
names(analyse)[names(analyse)=="dif.detect.unif.6"] <- "dif_detect_unif_6"
names(analyse)[names(analyse)=="dif.detect.unif.7"] <- "dif_detect_unif_7"
#analyse[analyse$dif_detect_5==999,"dif_detect_5"] <- ""
#analyse[analyse$dif_detect_6==999,"dif_detect_6"] <- ""
#analyse[analyse$dif_detect_7==999,"dif_detect_7"] <- ""
#analyse[analyse$dif_detect_unif_5==999,"dif_detect_unif_5"] <- ""
#analyse[analyse$dif_detect_unif_6==999,"dif_detect_unif_6"] <- ""
#analyse[analyse$dif_detect_unif_7==999,"dif_detect_unif_7"] <- ""
}
analyse <- analyse[,!names(analyse) %in% c("X","X.1","X.2")]
write.csv(analyse,paste0("/home/corentin/Documents/These/Recherche/Simulations/Analysis/RESALI/Detection_data/",r,".csv"))
}
for (r in results[-c(1:4)]) {
cat('--------------------------------------------------------------------------\n')
cat(paste0(r,"\n"))
cat('--------------------------------------------------------------------------\n')
#### Importer les datas
scen <- as.numeric(gsub("[A,B,C,D,E,F,G,_]","",substr(r,0,3)))
if (substr(r,start=nchar(r)-1,stop=nchar(r))=="50") {
N <- 50
}
if (substr(r,start=nchar(r)-2,stop=nchar(r))=="100") {
N <- 100
}
if (substr(r,start=nchar(r)-2,stop=nchar(r))=="200") {
N <- 200
}
if (substr(r,start=nchar(r)-2,stop=nchar(r))=="300") {
N <- 300
}
#### Importer l'analyse
analyse <- read.csv(paste0("/home/corentin/Documents/These/Recherche/Simulations/Analysis/RESALI/Detection_data/",r,".csv"))
names(analyse)[names(analyse)=="dif.detect.1"] <- "dif_detect_1"
names(analyse)[names(analyse)=="dif.detect.2"] <- "dif_detect_2"
names(analyse)[names(analyse)=="dif.detect.3"] <- "dif_detect_3"
names(analyse)[names(analyse)=="dif.detect.4"] <- "dif_detect_4"
names(analyse)[names(analyse)=="dif.detect.unif.1"] <- "dif_detect_unif_1"
names(analyse)[names(analyse)=="dif.detect.unif.2"] <- "dif_detect_unif_2"
names(analyse)[names(analyse)=="dif.detect.unif.3"] <- "dif_detect_unif_3"
names(analyse)[names(analyse)=="dif.detect.unif.4"] <- "dif_detect_unif_4"
analyse[analyse$dif_detect_1==999,"dif_detect_1"] <- ""
analyse[analyse$dif_detect_2==999,"dif_detect_2"] <- ""
analyse[analyse$dif_detect_3==999,"dif_detect_3"] <- ""
analyse[analyse$dif_detect_4==999,"dif_detect_4"] <- ""
analyse[analyse$dif_detect_unif_1==999,"dif_detect_unif_1"] <- ""
analyse[analyse$dif_detect_unif_2==999,"dif_detect_unif_2"] <- ""
analyse[analyse$dif_detect_unif_3==999,"dif_detect_unif_3"] <- ""
analyse[analyse$dif_detect_unif_4==999,"dif_detect_unif_4"] <- ""
if (scen==3 | scen==4 | scen>=13) {
names(analyse)[names(analyse)=="dif.detect.5"] <- "dif_detect_5"
names(analyse)[names(analyse)=="dif.detect.6"] <- "dif_detect_6"
names(analyse)[names(analyse)=="dif.detect.7"] <- "dif_detect_7"
names(analyse)[names(analyse)=="dif.detect.unif.5"] <- "dif_detect_unif_5"
names(analyse)[names(analyse)=="dif.detect.unif.6"] <- "dif_detect_unif_6"
names(analyse)[names(analyse)=="dif.detect.unif.7"] <- "dif_detect_unif_7"
analyse[analyse$dif_detect_5==999,"dif_detect_5"] <- ""
analyse[analyse$dif_detect_6==999,"dif_detect_6"] <- ""
analyse[analyse$dif_detect_7==999,"dif_detect_7"] <- ""
analyse[analyse$dif_detect_unif_5==999,"dif_detect_unif_5"] <- ""
analyse[analyse$dif_detect_unif_6==999,"dif_detect_unif_6"] <- ""
analyse[analyse$dif_detect_unif_7==999,"dif_detect_unif_7"] <- ""
}
analyse <- analyse[,!names(analyse) %in% c("X","X.1","X.2")]
write.csv(analyse,paste0("/home/corentin/Documents/These/Recherche/Simulations/Analysis/RESALI/Detection_data/",r,".csv"))
}
for (r in results) {
cat('--------------------------------------------------------------------------\n')
cat(paste0(r,"\n"))
cat('--------------------------------------------------------------------------\n')
#### Importer les datas
scen <- as.numeric(gsub("[A,B,C,D,E,F,G,_]","",substr(r,0,3)))
if (substr(r,start=nchar(r)-1,stop=nchar(r))=="50") {
N <- 50
}
if (substr(r,start=nchar(r)-2,stop=nchar(r))=="100") {
N <- 100
}
if (substr(r,start=nchar(r)-2,stop=nchar(r))=="200") {
N <- 200
}
if (substr(r,start=nchar(r)-2,stop=nchar(r))=="300") {
N <- 300
}
#### Importer l'analyse
analyse <- read.csv(paste0("/home/corentin/Documents/These/Recherche/Simulations/Analysis/RESALI/Detection_data/",r,".csv"))
names(analyse)[names(analyse)=="dif.detect.1"] <- "dif_detect_1"
names(analyse)[names(analyse)=="dif.detect.2"] <- "dif_detect_2"
names(analyse)[names(analyse)=="dif.detect.3"] <- "dif_detect_3"
names(analyse)[names(analyse)=="dif.detect.4"] <- "dif_detect_4"
names(analyse)[names(analyse)=="dif.detect.unif.1"] <- "dif_detect_unif_1"
names(analyse)[names(analyse)=="dif.detect.unif.2"] <- "dif_detect_unif_2"
names(analyse)[names(analyse)=="dif.detect.unif.3"] <- "dif_detect_unif_3"
names(analyse)[names(analyse)=="dif.detect.unif.4"] <- "dif_detect_unif_4"
analyse[analyse$dif_detect_1==999,"dif_detect_1"] <- ""
analyse[analyse$dif_detect_2==999,"dif_detect_2"] <- ""
analyse[analyse$dif_detect_3==999,"dif_detect_3"] <- ""
analyse[analyse$dif_detect_4==999,"dif_detect_4"] <- ""
analyse[analyse$dif_detect_unif_1==999,"dif_detect_unif_1"] <- ""
analyse[analyse$dif_detect_unif_2==999,"dif_detect_unif_2"] <- ""
analyse[analyse$dif_detect_unif_3==999,"dif_detect_unif_3"] <- ""
analyse[analyse$dif_detect_unif_4==999,"dif_detect_unif_4"] <- ""
if (scen==3 | scen==4 | scen>=13) {
names(analyse)[names(analyse)=="dif.detect.5"] <- "dif_detect_5"
names(analyse)[names(analyse)=="dif.detect.6"] <- "dif_detect_6"
names(analyse)[names(analyse)=="dif.detect.7"] <- "dif_detect_7"
names(analyse)[names(analyse)=="dif.detect.unif.5"] <- "dif_detect_unif_5"
names(analyse)[names(analyse)=="dif.detect.unif.6"] <- "dif_detect_unif_6"
names(analyse)[names(analyse)=="dif.detect.unif.7"] <- "dif_detect_unif_7"
analyse[analyse$dif_detect_5==999,"dif_detect_5"] <- ""
analyse[analyse$dif_detect_6==999,"dif_detect_6"] <- ""
analyse[analyse$dif_detect_7==999,"dif_detect_7"] <- ""
analyse[analyse$dif_detect_unif_5==999,"dif_detect_unif_5"] <- ""
analyse[analyse$dif_detect_unif_6==999,"dif_detect_unif_6"] <- ""
analyse[analyse$dif_detect_unif_7==999,"dif_detect_unif_7"] <- ""
}
analyse[is.na(analyse)] <- ""
analyse <- analyse[,!names(analyse) %in% c("X","X.1","X.2")]
write.csv(analyse,paste0("/home/corentin/Documents/These/Recherche/Simulations/Analysis/RESALI/Detection_data/",r,".csv"))
}
for (r in results) {
cat('--------------------------------------------------------------------------\n')
cat(paste0(r,"\n"))
cat('--------------------------------------------------------------------------\n')
#### Importer les datas
scen <- as.numeric(gsub("[A,B,C,D,E,F,G,_]","",substr(r,0,3)))
if (substr(r,start=nchar(r)-1,stop=nchar(r))=="50") {
N <- 50
}
if (substr(r,start=nchar(r)-2,stop=nchar(r))=="100") {
N <- 100
}
if (substr(r,start=nchar(r)-2,stop=nchar(r))=="200") {
N <- 200
}
if (substr(r,start=nchar(r)-2,stop=nchar(r))=="300") {
N <- 300
}
#### Importer l'analyse
analyse <- read.csv(paste0("/home/corentin/Documents/These/Recherche/Simulations/Analysis/RESALI/Detection_data/",r,".csv"))
analyse[is.na(analyse)] <- ""
names(analyse)[names(analyse)=="dif.detect.1"] <- "dif_detect_1"
names(analyse)[names(analyse)=="dif.detect.2"] <- "dif_detect_2"
names(analyse)[names(analyse)=="dif.detect.3"] <- "dif_detect_3"
names(analyse)[names(analyse)=="dif.detect.4"] <- "dif_detect_4"
names(analyse)[names(analyse)=="dif.detect.unif.1"] <- "dif_detect_unif_1"
names(analyse)[names(analyse)=="dif.detect.unif.2"] <- "dif_detect_unif_2"
names(analyse)[names(analyse)=="dif.detect.unif.3"] <- "dif_detect_unif_3"
names(analyse)[names(analyse)=="dif.detect.unif.4"] <- "dif_detect_unif_4"
analyse[analyse$dif_detect_1==999,"dif_detect_1"] <- ""
analyse[analyse$dif_detect_2==999,"dif_detect_2"] <- ""
analyse[analyse$dif_detect_3==999,"dif_detect_3"] <- ""
analyse[analyse$dif_detect_4==999,"dif_detect_4"] <- ""
analyse[analyse$dif_detect_unif_1==999,"dif_detect_unif_1"] <- ""
analyse[analyse$dif_detect_unif_2==999,"dif_detect_unif_2"] <- ""
analyse[analyse$dif_detect_unif_3==999,"dif_detect_unif_3"] <- ""
analyse[analyse$dif_detect_unif_4==999,"dif_detect_unif_4"] <- ""
if (scen==3 | scen==4 | scen>=13) {
names(analyse)[names(analyse)=="dif.detect.5"] <- "dif_detect_5"
names(analyse)[names(analyse)=="dif.detect.6"] <- "dif_detect_6"
names(analyse)[names(analyse)=="dif.detect.7"] <- "dif_detect_7"
names(analyse)[names(analyse)=="dif.detect.unif.5"] <- "dif_detect_unif_5"
names(analyse)[names(analyse)=="dif.detect.unif.6"] <- "dif_detect_unif_6"
names(analyse)[names(analyse)=="dif.detect.unif.7"] <- "dif_detect_unif_7"
}
analyse <- analyse[,!names(analyse) %in% c("X","X.1","X.2")]
write.csv(analyse,paste0("/home/corentin/Documents/These/Recherche/Simulations/Analysis/RESALI/Detection_data/",r,".csv"))
}
for (r in results) {
cat('--------------------------------------------------------------------------\n')
cat(paste0(r,"\n"))
cat('--------------------------------------------------------------------------\n')
#### Importer les datas
scen <- as.numeric(gsub("[A,B,C,D,E,F,G,_]","",substr(r,0,3)))
if (substr(r,start=nchar(r)-1,stop=nchar(r))=="50") {
N <- 50
}
if (substr(r,start=nchar(r)-2,stop=nchar(r))=="100") {
N <- 100
}
if (substr(r,start=nchar(r)-2,stop=nchar(r))=="200") {
N <- 200
}
if (substr(r,start=nchar(r)-2,stop=nchar(r))=="300") {
N <- 300
}
#### Importer l'analyse
analyse <- read.csv(paste0("/home/corentin/Documents/These/Recherche/Simulations/Analysis/RESALI/Detection_data/",r,".csv"))
analyse[is.na(analyse)] <- ""
names(analyse)[names(analyse)=="dif.detect.1"] <- "dif_detect_1"
names(analyse)[names(analyse)=="dif.detect.2"] <- "dif_detect_2"
names(analyse)[names(analyse)=="dif.detect.3"] <- "dif_detect_3"
names(analyse)[names(analyse)=="dif.detect.4"] <- "dif_detect_4"
names(analyse)[names(analyse)=="dif.detect.unif.1"] <- "dif_detect_unif_1"
names(analyse)[names(analyse)=="dif.detect.unif.2"] <- "dif_detect_unif_2"
names(analyse)[names(analyse)=="dif.detect.unif.3"] <- "dif_detect_unif_3"
names(analyse)[names(analyse)=="dif.detect.unif.4"] <- "dif_detect_unif_4"
analyse[analyse$dif_detect_1==999,"dif_detect_1"] <- ""
analyse[analyse$dif_detect_2==999,"dif_detect_2"] <- ""
analyse[analyse$dif_detect_3==999,"dif_detect_3"] <- ""
analyse[analyse$dif_detect_4==999,"dif_detect_4"] <- ""
analyse[analyse$dif_detect_unif_1==999,"dif_detect_unif_1"] <- ""
analyse[analyse$dif_detect_unif_2==999,"dif_detect_unif_2"] <- ""
analyse[analyse$dif_detect_unif_3==999,"dif_detect_unif_3"] <- ""
analyse[analyse$dif_detect_unif_4==999,"dif_detect_unif_4"] <- ""
if (scen==3 | scen==4 | scen>=13) {
names(analyse)[names(analyse)=="dif.detect.5"] <- "dif_detect_5"
names(analyse)[names(analyse)=="dif.detect.6"] <- "dif_detect_6"
names(analyse)[names(analyse)=="dif.detect.7"] <- "dif_detect_7"
names(analyse)[names(analyse)=="dif.detect.unif.5"] <- "dif_detect_unif_5"
names(analyse)[names(analyse)=="dif.detect.unif.6"] <- "dif_detect_unif_6"
names(analyse)[names(analyse)=="dif.detect.unif.7"] <- "dif_detect_unif_7"
analyse[analyse$dif_detect_5==999,"dif_detect_5"] <- ""
analyse[analyse$dif_detect_6==999,"dif_detect_6"] <- ""
analyse[analyse$dif_detect_7==999,"dif_detect_7"] <- ""
analyse[analyse$dif_detect_unif_5==999,"dif_detect_unif_5"] <- ""
analyse[analyse$dif_detect_unif_6==999,"dif_detect_unif_6"] <- ""
analyse[analyse$dif_detect_unif_7==999,"dif_detect_unif_7"] <- ""
}
analyse <- analyse[,!names(analyse) %in% c("X","X.1","X.2")]
write.csv(analyse,paste0("/home/corentin/Documents/These/Recherche/Simulations/Analysis/RESALI/Detection_data/",r,".csv"))
}
