diff --git a/R/pcrslayer.R b/R/pcrslayer.R index c8b8387..52d6fbc 100644 --- a/R/pcrslayer.R +++ b/R/pcrslayer.R @@ -423,8 +423,8 @@ check_pcr_repl <- function(metabarlist, funcpcr <- ifelse(funcpcr, "0ok", "dyspcr") } - reads_stdt <- reads/rowSums(reads) + replicates <- paste0("bar_", as.character(replicates)) bar <- rowsum(reads_stdt, replicates)/as.vector(table(replicates)) all <- rbind(reads_stdt, bar) all.d <- vegdist(all, "bray")