Skip to content
Snippets Groups Projects
Commit f6625b45 authored by Simeon's avatar Simeon
Browse files

minor bug in clustered_seqs var call

parent 8f0af997
No related branches found
No related tags found
No related merge requests found
...@@ -45,7 +45,9 @@ export(translate_and_count_stops) ...@@ -45,7 +45,9 @@ export(translate_and_count_stops)
export(variant_classifier) export(variant_classifier)
export(veganify_asvcounts) export(veganify_asvcounts)
export(veganify_generic_wide_tbl) export(veganify_generic_wide_tbl)
import(Biostrings)
import(DECIPHER) import(DECIPHER)
import(DNAStringSetList)
import(dplyr) import(dplyr)
import(forcats) import(forcats)
import(ggplot2) import(ggplot2)
......
...@@ -8,6 +8,7 @@ ...@@ -8,6 +8,7 @@
#' #'
#' @return A modified master table with variant classifications #' @return A modified master table with variant classifications
#' @import dplyr purrr stringr tibble tidyr #' @import dplyr purrr stringr tibble tidyr
#' @import Biostrings DNAStringSetList
#' @export #' @export
variant_classifier <- function( variant_classifier <- function(
...@@ -26,8 +27,8 @@ variant_classifier <- function( ...@@ -26,8 +27,8 @@ variant_classifier <- function(
# Create a table for clustered sequences and filter out clusters with only one sequence # Create a table for clustered sequences and filter out clusters with only one sequence
clustab_tbl <- cluster_tbl_named(clustered_sequences) %>% clustab_tbl <- cluster_tbl_named(clustered_sequences) %>%
left_join(tibble(seqnames = unlist(map(clustered_seqs, names)), left_join(tibble(seqnames = unlist(map(clustered_sequences, names)),
seqs = as.character(unlist(clustered_sequences))), seqs = as.character(Biostrings:: unlist(clustered_sequences))),
by = 'seqnames') %>% by = 'seqnames') %>%
dplyr::filter(clus_size > 1) dplyr::filter(clus_size > 1)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment