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"))