scripts_r/rscript_rjmcmc.R

In alexanderbuchholz/distbayesianmc: Distributed Bayesian Model Choice

#!/usr/bin/env Rscript
library(distbayesianmc)
library(R2BGLiMS)
#task_id_string <- Sys.getenv("SLURM_ARRAY_TASK_ID")
#task_id <- as.numeric(task_id_string)

mrep <- 20
n.mil <- 10
n.cores <- 1
dataset <-  "sim2"
nobsdata <- 4000
typesplit <-  "random"
savestring <-  paste(dataset, typesplit, sep = "_")


list_params_model_onesplit <- list(scale =  1,
                                   ssplits =  1,
                                   typesplit =  typesplit,
                                   dataset = dataset, 
                                   mrep = mrep,
                                   n.mil = n.mil, 
                                   n.cores = n.cores, 
                                   savestring = savestring, 
                                   nobsdata = nobsdata)

list_params_model_multisplits1 <- list(scale =  1,
                                       ssplits =  2,
                                       typesplit =  typesplit,
                                       dataset = dataset,#"",
                                       mrep = mrep, 
                                       n.mil = n.mil,
                                       n.cores = n.cores,
                                       savestring = savestring,
                                       nobsdata = nobsdata)

list_params_model_multisplits2 <- list(scale =  1,
                                       ssplits =  3,
                                       typesplit =  typesplit,
                                       dataset = dataset,#"",
                                       mrep = mrep,
                                       n.mil = n.mil,
                                       n.cores = n.cores,
                                       savestring = savestring,
                                       nobsdata = nobsdata)

list_params_model_multisplits3 <- list(scale =  1,
                                       ssplits =  5,
                                       typesplit =  typesplit,
                                       dataset = dataset,#"",
                                       mrep = mrep,
                                       n.mil = n.mil,
                                       n.cores = n.cores,
                                       savestring = savestring,
                                       nobsdata = nobsdata)

list_params_model_multisplits4 <- list(scale =  1,
                                       ssplits =  10,
                                       typesplit =  typesplit,
                                       dataset = dataset,#"",
                                       mrep = mrep,
                                       n.mil = n.mil,
                                       n.cores = n.cores,
                                       savestring = savestring,
                                       nobsdata = nobsdata)


list_params_model <- list(list_params_model_onesplit, 
                          list_params_model_multisplits1, 
                          list_params_model_multisplits2, 
                          list_params_model_multisplits3, 
                          list_params_model_multisplits4)



# run only if run by itself
if (sys.nframe() == 0){
  args = commandArgs(trailingOnly=TRUE)
  print(args)
  sim_id <- as.numeric(args[1])
  
  setwd("/scratch/alexander/distbayesianmc_rjmcmc/")
  #library(doParallel)
  #registerDoParallel(cores=6)
  for(params_model in list_params_model){
    #set.seed(sim_id)
    #x_wait <- rexp(1,1)
    #Sys.sleep(x_wait)
    #foreach(sim_id = 1:mrep) %dopar% {
      res_sim <- f_single_run_rep_rjmcmc(params_model, sim_id)
    #}
    #save(res_sim, file = paste("res_split_id", sim_id, "_splits_", params_model$ssplits, ".RData", sep = ""))
  }
}