r - 自定义 metaMDS() 图

Question

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

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

现在我有了这个 NMDS，我的观点得到了很好的传播。但是，我想添加情节的图形。我想给图中的点一个不同的颜色，这取决于我的数据集中的一个分类变量（变量称为“区域”），它有两个值（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