polarolouis/mia-stage-2024
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

🌱 Adding the seed to the clustering

Browse source
This commit is contained in:
Louis Lacoste 2024-06-04 11:43:47 +02:00
parent 2a40e89a88
commit 7706600bac
2 changed files with 21 additions and 7 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

19
code/applications/dore/01_dore_clusterize.R
View file

@ -7,6 +7,7 @@ library(here)
arg <- commandArgs(trailingOnly = TRUE)
model <- "iid"
seed <- 1234L
if (length(arg) == 0L) {
message("No arguments provided, using default.")
@ -16,17 +17,26 @@ if (length(arg) == 0L) {
} else {
message("No model provided, defaulting to iid.")
}
if ("--seed" %in% arg) {
seed <- try(as.integer(arg[(which(arg == "--seed") + 1L)]))
} else {
message("No seed provided, defaulting to 1234.")
}
}
#  Arguments checks
allowed_model <- c("iid", "pi", "rho", "pirho")
stopifnot(
"Unknown model, should be : iid, pi, rho or pirho" = (model %in% allowed_model)
"Unknown model, should be : iid, pi, rho or pirho" = (model %in% allowed_model),
"Seed isn't castable to integer" = (is.integer(seed))
)
message("Début du clustering des données Doré et al. avec le modèle ", model)
message(
"Début du clustering des données Doré et al. avec le modèle ", model,
" et la seed ", seed
)
set.seed(1234, "L'Ecuyer-CMRG")
set.seed(seed, "L'Ecuyer-CMRG")
base_data_folder <- file.path(here(), "code", "data", "dore")
save_folder <- file.path(here(), "code", "results", "applications", "dore")
@ -97,7 +107,8 @@ list_collection <- clusterize_bipartite_networks(
save_file <- file.path(
save_folder, paste0(
"dore_collection_",
model, "_", start_time, ".Rds"
model, "_seed_", seed,
"_", start_time, ".Rds"
)
)

9
code/scripts/migale_application_dore.sh
View file

@ -3,7 +3,7 @@
#$ -cwd
#$ -N Dore_app
#$ -m besa
#$ -t 1:4
#$ -t 1:16
#$ -q long.q
#$ -pe thread 64
#$ -M louis.lacoste+migale@agroparistech.fr
@ -20,12 +20,15 @@ fi
# Constant data
MODELARRAY=("iid" "pi" "rho" "pirho")
ID=$((SGE_TASK_ID - 1))
MODEL=${MODELARRAY[$((ID % 4))]}
SEED=$(($((ID + $((ID / 4)))) % 4))
MODEL=${MODELARRAY[$(($((SGE_TASK_ID - 1)) % 4))]}
# Finding directory
APPLICATIONS_DIR=$(echo "$BASE_DIR/code/applications")
echo $APPLICATIONS_DIR
Rscript "${APPLICATIONS_DIR}/dore/01_dore_clusterize.R" --model $MODEL
Rscript "${APPLICATIONS_DIR}/dore/01_dore_clusterize.R" --model $MODEL --seed $SEED