switch to poisson clustering

src/bin/kmerclust_boot.R

@@ -5,9 +5,7 @@ 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]))
+  poiss_compare_models(count, nboot = as.numeric(args[2]) * 2, bootsize = 1, core = as.numeric(args[2]))
 save(res, file = paste0(args[1], "_boot_", args[3], ".Rdata"))
 
 

src/bin/kmerclust_load.R

@@ -2,7 +2,6 @@ library(kmerclust)
 library(vroom)
 
 args <- commandArgs(trailingOnly = TRUE)
-
 data <- vroom::vroom(args[2], show_col_types = FALSE)
 format(object.size(data), units = "Mb")
 annotation <- parse_annotation(args[1])