diff --git a/src/bin/kmerclust_boot.R b/src/bin/kmerclust_boot.R
index eb2e344970e2376aaaf82b20b71452292962a92f..1857219b96640334ee69c8d459881ddf8e363d85 100644
--- a/src/bin/kmerclust_boot.R
+++ b/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"))
 
 
diff --git a/src/bin/kmerclust_plot.R b/src/bin/kmerclust_plot.R
index 5f50f46d9b2e84b0afce113c947848d652f851da..5a62bbe49c02f0883118f802023185be578cf483 100644
--- a/src/bin/kmerclust_plot.R
+++ b/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 %>%