0

这是我如何分离数据集中前 10 个订单以绘制相对丰度的代码。我怎样才能将前 10 个订单 +第 11 组称为其他人的剩余数据包含在内,以使我的堆叠条填充到 100%?

order.sum = tapply(taxa_sums(relabun), tax_table(relabun)[, "Order"], sum, na.rm=TRUE) top10orders = names(sort(order.sum, TRUE))[1:10] ps2 = prune_taxa((tax_table(relabun)[, "Order"] %in% top10orders), relabun) ps_top10orders <- tax_glom(ps2,taxrank = "Order") ps_top10orders<-transform_sample_counts(ps_top10orders,function(x) x / sum(x)) ps_top10orders<-psmelt(ps_top10orders) ps_top10orders<-arrange(ps_top10orders,Order)

4

0 回答 0