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