这是我在将我的 OTU 表、元数据、分类和树作为一个对象导入 R 后当前使用的代码。我遇到的问题是在最后一行，我不断收到错误“dimnames(x) <- dn 中的错误：'dimnames' [1] 的长度不等于数组范围”
physeq = phyloseq(OTU, TAX, META, phy_tree)
library(phyloseq)
tps = transform_sample_counts(physeq, function(x) x / sum(x) )
ftps = filter_taxa(physeq, function(x) {sum(x >0)>1}, TRUE)
psfirm = subset_taxa(physeq, Phylum==“Firmicutes”)