1

我使用 Rphyloseq函数纵坐标和plot_ordination一个phyloseq对象和一个先前计算的距离矩阵（未加权的 UniFrac 距离）作为输入，制作了微生物组数据的排序图。我想添加一些箭头来指示哪些物种的相对丰度主要驱动样本之间沿轴的距离。

使用 score 函数添加分数或specscores.dbrda无效。还有其他方法可以添加这些箭头/矢量吗？

这是我使用的代码：

MDS <- ordinate(physeqobject, "MDS", distance=unifracmatrix)
plot <- plot_ordination(physeqobject, MDS, color = "variable1")

speciesscores <- scores(MDS, display = "species")

第一部分工作正常，scores 函数返回此错误：

Error in scores.default(MDS, display = "species") : 
  cannot find scores

我也尝试过没有单独的距离矩阵，在纵坐标内计算（默认）Bray-Curtis 距离，但我仍然得到相同的错误：

MDS <-纵坐标（physeqobject，“MDS”）

4

0 回答 0