genomewalker
5/18/2018 - 11:42 AM

nonpareil parser

library(Nonpareil)
library(tidyverse)
f <- list.files(path="/scratch/antonio/nonpareil/out/", pattern = "npo")
f <- list.files(path="~/tara_nonpareil/", pattern = "npo")
sample.names <- sapply(strsplit(f, ".npo"), `[`, 1)

Nonpareil.curve(f[[1]])

results_np_tara <- lapply(file.path("~/tara_nonpareil/",f), Nonpareil.curve)

names(results_np_tara) <- sample.names

purrr::map_df(results_np_tara, "C") %>% gather() %>% 
  ggplot(aes(1,value)) +
  geom_boxplot()

purrr::map_df(results_np_tara, "diversity") %>% gather() %>% 
  ggplot(aes("1", value)) +
  geom_boxplot()

  save(results_np_tara, file = "~/tara_nonpareil.Rdata")