Arcadia-Science / sourmashconsumr

Working with the outputs of sourmash in R
https://arcadia-science.github.io/sourmashconsumr/
Other
25 stars 3 forks source link

example of making rarefaction curves from signatures representing fastq files #25

Closed taylorreiter closed 2 years ago

taylorreiter commented 2 years ago
remotes::install_github("Arcadia-Science/sourmashconsumr")
library(sourmashconsumr)
library(dplyr)
library(ggplot2)
library(purrr)

sigs <- Sys.glob("*100k.sig") %>%
  map_dfr(read_signature) %>%
  filter(ksize == 21)

rarefaction_df <- from_signatures_to_rarefaction_df(sigs)
plot_signatures_rarefaction(rarefaction_df) # +
  # theme_minimal() +
  # geom_point(aes(color = name))

image

uncomment lines to get colored curves and no grey background.

taylorreiter commented 2 years ago

closed by #46