有没有一种方法可以只在NMDS ordiplot上显示选定的物种？

Question

我有一个NMDS图，使用plot()函数在其上覆盖了物种负载。物种负荷几乎都是重要的，但我只想显示影响最大的少数物种。我不确定如何制作一个可重现的示例，因为大多数数据集所包含的物种数量并不多。有没有一种方法可以从envfit()加载函数中选择某些物种，然后调用plot()？

下面的代码：

ordiplot(inverts_test_log10_10.nmds, choices = c(1, 2), type="n", display="sites", xlab="Axis 1", ylab="Axis 2")

points(inverts_test_log10_10.nmds, col=as.vector(SW.sitegroup$color), bg=as.vector(SW.sitegroup$color),pch=21, cex=1)

plot(sp.load, p.max=.001, col="blue")

带有载荷的NMDS图：