From 12621d49749ac0e373048b081561dadcc203fb1e Mon Sep 17 00:00:00 2001 From: Simeon <51403284+simeross@users.noreply.github.com> Date: Wed, 8 Nov 2023 08:56:34 +0100 Subject: [PATCH] bug fix previous commit --- R/per_cluster_abundance.R | 5 ++++- R/read_and_write_cluster_abundance.R | 6 ++---- 2 files changed, 6 insertions(+), 5 deletions(-) diff --git a/R/per_cluster_abundance.R b/R/per_cluster_abundance.R index 2348497..d189373 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 255b5fe..dc6361c 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]] -- GitLab