r - 按 r 中的列值标记 nmds 图

Question

我运行 metaMDS 并希望通过基于某些数据框字符的分组来绘制和颜色代码。在我的原始数据框中，df$yr 是年份，df$2 是网站。我想按年着色。

caltmds <- metaMDS(df[,3:12], k=3)
plot(caltmds, type = 'n')
cols <- c("red2", "mediumblue")
points(caltmds, col = cols[df$yr])

我也从这篇文章中尝试过：

scl <- 3
colvec <- c("red2", "mediumblue")
plot(caltmds, type = "n", scaling = scl) 
with(df, points(caltmds, display = "sites", col = colvec[yr], pch = 21, bg = colvec[yr]))
text(caltmds, display = "species", cex = 0.8, col = "darkcyan")
with(df, legend("topright", legend = levels(yr), bty = "n", col = colvec, pch = 21, pt.bg = colvec))

没有情节

Answer 1

#DATA
df1 = mtcars

mycolors = df1$cyl   #Identify the grouping vector

library(vegan)
m = metaMDS(df1)

x = scores(m)  #Extract co-ordinates

plot(x, col = as.numeric(as.factor(mycolors)))