我使用'vegan'包创建了一个NMDS图，如下所示：
y=metaMDS(data,type="p").
plot(y)
现在我的NMDS有很好的分数。但是，我想添加图表的图形。我想在图中给出不同颜色的点，这取决于我的数据集中的分类变量（变量称为'regio'），它有两个值（1或2）。
这可能吗？如果是这样，怎么样？
最佳， 柯恩
答案 0 :(得分：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