0

使用 计算相对丰度phyloseq

我想为每个样本创建一个图而不使用merge_samples.

这是我的代码:

 carbom1 <- subset_taxa(carbom, Phylum %in% c(" p__Firmicutes", " p__Bacteroidota", " p__Proteobacteria",
               " p__Fusobacteriota", " p__Actinobacteriota"))
                   
ps <- tax_glom(carbom1, "Phylum")

ps0 <- transform_sample_counts(ps, function(x) x / sum(x))

ps0

phyloseq-class 实验级对象

otu_table() OTU 表:[5 个分类单元和 18 个样本]

sample_data() 样本数据:[ 2 个样本变量的 18 个样本]

tax_table() 分类表:[ 7 个分类等级的 5 个分类群]

如何绘制每个样本的图ps0

4

0 回答 0