0

我正在尝试对相对丰度 >= 35% 并且属于分组内 >= 70% 的样本的分类群进行子集化(或过滤?)(在我的情况下,它是我的数据中的“集群”的数量)。

是否有关于如何执行此操作的简单代码行?我已经开始用这行代码做到这一点

subset <- filter_taxa(phyloseq_object, function (x) sum (x) > 0.35, TRUE)

我无法弄清楚如何应用此过滤步骤来查看这些分类群是否属于 >= 70% 的样本。也许在分组中使用这样的东西?

prevalenceThreshold = 0.70 * nsamples(phyloseq_object)

我的 sample_data 文件中有一个名为“cluster”的列,其中包含 8 个已定义的集群。也许我想多了,但我希望这是有道理的。谢谢!

编辑** 我在 R 版本 4.1.2 中使用 MacOS Monterey 12.2.1 执行此操作。这是我的数据的一个子集:

这是我的数据的一个子集。我无法在图像中放置“集群”列。每个样本属于 1 个集群。

4

0 回答 0