diff --git a/R/per_cluster_abundance.R b/R/per_cluster_abundance.R index 234849754ae49989b6d720dfc48a9e7d3618525b..d1893738d48aee71c1af0a0a838d05d4f5743b97 100644 --- a/R/per_cluster_abundance.R +++ b/R/per_cluster_abundance.R @@ -12,10 +12,13 @@ per_cluster_abundance <- function(seqs_of_one_cluster = DNAStringSet, seqtab_nochim = stab, reference_seqs = NULL){ + stab <- seqtab_nochim %>% + t() %>% as.data.frame() %>% + rownames_to_column(var = 'seqs') seq_tbl <- tibble(seqs = as.data.frame(seqs_of_one_cluster)[[1]], ID = names(seqs_of_one_cluster)) - named_stab <- left_join(seq_tbl, seqtab_nochim, by = 'seqs') + named_stab <- left_join(seq_tbl, stab, by = 'seqs') if(!is.null(reference_seqs)){ named_stab <- named_stab %>% diff --git a/R/read_and_write_cluster_abundance.R b/R/read_and_write_cluster_abundance.R index 255b5feba15c5972fc72ce0bf907b7a03a9fb1f3..dc6361c87ead3dcdad7de60de1f4a3bbf0fdbdd4 100644 --- a/R/read_and_write_cluster_abundance.R +++ b/R/read_and_write_cluster_abundance.R @@ -22,12 +22,10 @@ read_and_write_cluster_abundance <- function( seqtab_nochim = 'seqtab_nochim.rds', outpath = path) { - stab <- readRDS(seqtab_nochim) %>% - t() %>% as.data.frame() %>% - rownames_to_column(var = 'seqs') + stab <- readRDS(seqtab_nochim) out <- lapply(cluster_sequence_list, per_cluster_abundance, - seqtab_nochim = seqtab_nochim, reference_seqs = reference_seqs) + seqtab_nochim = stab, reference_seqs = reference_seqs) write_cluster_abundance_tbl <- function(clus_name){ stab_tbl <- out[[clus_name]]