kmerclust_plot.R

library(kmerclust)
library(tidyverse)
library(ggplot2)

args <- commandArgs(trailingOnly = TRUE)
print(args)
load(file = paste0(args[1], "_boot.Rdata"))

print(res)
res %>%
  ggplot(aes(x = name, y = BIC)) +
  geom_violin() +
  theme_bw()
ggsave(paste0(args[1], "_BIC.pdf"))

load(file = paste0(args[1], ".Rdata"))
count %>%
  sample_frac(0.1) %>%
  mutate(
    count_m = log1p(count_m),
    count_f = log1p(count_f),
  )
  ggplot(aes(x = count_m, y = count_f)) +
  geom_point() +
  theme_bw()
ggsave(paste0(args[1], "_scatter.pdf"))