如何将物种箭头/向量添加到排序图

Question

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

使用功能分数或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”）