Guest User

Untitled

a guest
May 23rd, 2018
71
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.61 KB | None | 0 0
  1. library(Nonpareil)
  2. library(tidyverse)
  3. f <- list.files(path="/scratch/antonio/nonpareil/out/", pattern = "npo")
  4. f <- list.files(path="~/tara_nonpareil/", pattern = "npo")
  5. sample.names <- sapply(strsplit(f, ".npo"), `[`, 1)
  6.  
  7. Nonpareil.curve(f[[1]])
  8.  
  9. results_np_tara <- lapply(file.path("~/tara_nonpareil/",f), Nonpareil.curve)
  10.  
  11. names(results_np_tara) <- sample.names
  12.  
  13. purrr::map_df(results_np_tara, "C") %>% gather() %>%
  14. ggplot(aes(1,value)) +
  15. geom_boxplot()
  16.  
  17. purrr::map_df(results_np_tara, "diversity") %>% gather() %>%
  18. ggplot(aes("1", value)) +
  19. geom_boxplot()
  20.  
  21. save(results_np_tara, file = "~/tara_nonpareil.Rdata")
Add Comment
Please, Sign In to add comment