kmerclust_boot.R

library(kmerclust)
args <- commandArgs(trailingOnly = TRUE)
print(args)
load(file = paste0(args[1], ".Rdata"))
res <- count %>%
  dplyr::select(count_m, count_f) %>%
  as.matrix() %>%
  compute_tpm() %>%
  log() %>%
  compare_models_constraint(count, nboot = as.numeric(args[2]) * 2, bootsize = nrow(count), core = as.numeric(args[2]))
save(res, file = paste0(args[1], "_boot_", args[3], ".Rdata"))