add kmerclust_plot.R

src/bin/kmerclust_plot.R

+library(kmerclust)
+library(ggplot2)
+
+args <- commandArgs(trailingOnly = TRUE)
+print(args)
+load(file = paste0(args[3], ".Rdata"))
+load(file = paste0(args[3], "_boot.Rdata"))
+
+res %>%
+  ggplot(aes(x = name, y = BIC)) +
+  geom_violin() %>%
+  theme_bw()
+ggsave(paste0(args[3], "_BIC.pdf"))
+res %>%
+  ggplot(aes(x = name, y = WSS_f / BSS)) +
+  geom_violin() %>%
+  theme_bw()
+ggsave(paste0(args[3], "_WSS_BSS.pdf"))
+count %>%
+  ggplot(aes(x = count_m, y = count_f)) %>%
+  theme_bw()
+ggsave(paste0(args[3], "_scatter.pdf"))