自定义metaMDS（）图

Question

我使用'vegan'包创建了一个NMDS图，如下所示：

y=metaMDS(data,type="p").
plot(y)

现在我的NMDS有很好的分数。但是，我想添加图表的图形。我想在图中给出不同颜色的点，这取决于我的数据集中的分类变量（变量称为'regio'），它有两个值（1或2）。

这可能吗？如果是这样，怎么样？

最佳， 柯恩

Answer 1

最简单的方法是使用分组变量regio索引到要绘制的颜色矢量。例如，（未经测试，因为我没有data ...）

colvec <- c("red","blue")
plot(y, type = "n")
points(y, display = "sites", col = colvec[data$regio])
## or
text(y, display = "sites", col = colvec[data$regio])
## depending on how you want to represent the sample scores