diff --git a/src/bin/kmerclust_boot.R b/src/bin/kmerclust_boot.R index 89a70070fce509fa89daec5cb4a06aa4e073432e..eb2e344970e2376aaaf82b20b71452292962a92f 100644 --- a/src/bin/kmerclust_boot.R +++ b/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")) diff --git a/src/bin/kmerclust_load.R b/src/bin/kmerclust_load.R index 235cab928de032a2d8098b34e7c428adc761e965..d8afb0bbd294d4cb582913e3623ed2bf9430ade7 100644 --- a/src/bin/kmerclust_load.R +++ b/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])