update boot script to use 1e-3 of the reads for the fit

src/bin/kmerclust_boot.R

@@ -5,7 +5,7 @@ load(file = paste0(args[1], ".Rdata"))
 res <- count %>%
   dplyr::select(count_m, count_f) %>%
   as.matrix() %>%
-  poiss_compare_models(count, nboot = as.numeric(args[2]) * 2, bootsize = 1, core = as.numeric(args[2]))
+  poiss_compare_models(count, nboot = as.numeric(args[2]) * 2, frac = 1e-3 bootsize = 1, core = as.numeric(args[2]))
 save(res, file = paste0(args[1], "_boot_", args[3], ".Rdata"))
 
 

src/bin/kmerclust_plot.R

@@ -7,11 +7,13 @@ print(args)
 load(file = paste0(args[1], "_boot.Rdata"))
 
 print(res)
-res %>%
-  ggplot(aes(x = name, y = BIC)) +
+res %>% 
+  pivot_longer(cols = -c("name"), names_to = "metric") %>% 
+  ggplot(aes(x =  name, y = value)) +
   geom_violin() +
+  facet_wrap(~metric) +
   theme_bw()
-ggsave(paste0(args[1], "_BIC.pdf"))
+ggsave(paste0(args[1], "_BIC_loglik.pdf"))
 
 load(file = paste0(args[1], ".Rdata"))
 count %>%