inst/analysis/MGA.analysis_SI20/scripts/CIMseqSwarm.R

In EngeLab/CIMseq.publication: The CIMseq.publication package reproduces the figures associated with the CIMseq publication

#ARGS
args <- commandArgs(TRUE)

#PACKAGES
packages <- c("tidyverse", "future")
purrr::walk(packages, library, character.only = TRUE)
rm(packages)

currPath <- getwd()

#load data
if(file.exists(file.path(currPath, 'data/CIMseqData.rda'))) {
  load(file.path(currPath, 'data/CIMseqData.rda'))
}

if(!is.na(args[1])) {
  library(CIMseq, lib.loc = "/home/jason/R/x86_64-redhat-linux-gnu-library/3.5")
  library(CIMseq.data, lib.loc = "/home/jason/R/x86_64-redhat-linux-gnu-library/3.5")
  library(CIMseq.testing, lib.loc = "/home/jason/R/x86_64-redhat-linux-gnu-library/3.5")
  out <- if(is.na(args[2])) {out <- 'tmp'} else {out <- as.character(args[2])}
  
  #gives 620 swarm members
  baseSeed <- 43892
  init <- map(1:2, function(i) {
    cbind(
      swarmInit(cObjSng, 2, null.weight = 1, seed = baseSeed + i), 
      swarmInit(cObjSng, 3, null.weight = 1, seed = baseSeed + i)
    )
  }) %>% do.call(cbind, .)

  options(future.globals.maxSize = Inf)
  runSwarmUppmax(
    cObjSng, cObjMul, swarmInit = init, maxiter = 100,
    nSyntheticMultiplets = 2000, eps.stagnate = 1, maxit.stagnate = 5,
    currPath = out, args = args
  )
  
} else {
  library(CIMseq)
  library(CIMseq.data)
  library(CIMseq.testing)
  
  #gives 620 swarm members
  baseSeed <- 43892
  init <- map(1:2, function(i) {
    cbind(
      swarmInit(cObjSng, 2, null.weight = 1, seed = baseSeed + i), 
      swarmInit(cObjSng, 3, null.weight = 1, seed = baseSeed + i)
    )
  }) %>% do.call(cbind, .)

  options(future.globals.maxSize = Inf)  
  runSwarmMultiprocess(
    cObjSng, cObjMul, swarmInit = init, maxiter = 100,
    nSyntheticMultiplets = 2000, eps.stagnate = 1, maxit.stagnate = 5,
    currPath = getwd()
  )
}

writeLines(capture.output(sessionInfo()), file.path(currPath, "logs/sessionInfo_CIMseqSwarm.txt"))
print("finished")