Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- 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")
Add Comment
Please, Sign In to add comment