3

我正在尝试使用嵌套因子制作数据的 nMDS 图。我希望 nMDS 通过使用符号和颜色在一个图上显示这两个因素。

在这个可重现的示例中，如果use嵌套在 中moisture，我希望该图显示Moisture为不同的符号，然后显示Use为不同的颜色。

到目前为止，我已经弄清楚了这一点：

library("vegan")
library("BiodiversityR")

data(dune, dune.env)
MDS <- metaMDS(dune, distance="bray", strata=dune.env$Moisture)
MDS

plot(MDS$points[,2], MDS$points[,1], type="n", main="Communities by Use", 
     xlab="NMDS Axis 1", ylab="NMDS Axis 2", xlim=c(-1.5,1.5), ylim=c(-1.5,1.5))
ordisymbol(MDS, dune.env, factor="Use", cex=1.25, rainbow=T, legend=T)

这给了我不同的用途作为不同的符号和颜色，但没有告诉我水分。是否可以让它显示不同的因素？我假设它可能在MDS$points[,]争论中的某个地方，但我不确定它们到底在做什么。

4

1 回答 1

2

通过修改这个问题的答案来解决这个问题：Plot points of metaMDS

data(dune, dune.env)
dune.MDS <- metaMDS(dune, distance = "bray", strata=dune.env$Moisture)
dune.MDS

pchs<- c(0:5)
gr.moi <- factor(dune.env$Moisture)
gr.use <- factor(dune.env$Use)
col.gr <- c("red", "blue", "purple")


plot(dune.MDS, type = "n", display = "sites")
orditorp(dune.MDS,display="species",col="dark grey",air=0.01)
points(dune.MDS, display = "sites", pch = pchs[gr.moi], col = col.gr[gr.use])
legend("topright", legend=levels(gr.moi), bty = "n", col= c("black"), pch = pchs)
legend("bottomright", legend = levels(gr.use), bty = "n", col = col.gr, pch=c(20),)

它会产生一个带有符号和颜色的可爱情节，正是我想要的:) 显示这两个因素的嵌套 nMDS

于 2015-04-13T00:22:31.327 回答