Question

[在此输入图片说明] [1] 我从纯素包中遇到metaMDS和ordiplot的问题。我是新手，所以我无法想出一个超华丽的例子。

使用ordiplot我想根据类型（A，B，C，D，E组）标记图。我有一个数据框＆＃34; df＆＃34;带有一列绘图名称和一列类型。我不再关心这种类型是用颜色还是用不同类型的点来说明的。我的基本例子：

data(dune)
data(dune.env)
dist.dune<-distShell(dune, distSimRatio)
ord<-metaMDS(dist.dune)
metaMDS(comm=dune)
plot(ord, type="n")
colfactor<-factor(dune.env$Management)
points(ord, display="sites", cex=2, pch=21, col=c(1,2,3,4,5)[colfactor])
text(ord, display="sites", cex=2, col=c(1,2,3,4,5)[colfactor])

更新：通过朋友脚本检查并找到了略有不同的解决方案。这适用于沙丘数据并且易于验证。但是，相同的解决方案不适用于我的数据。我可以非常轻松地检查我的NMDS结果（我尝试着色的属性是我的点名称的一部分）并且它们的着色毫无意义，但是在正确数量的绘图以一种颜色着色的意义上它是正确的。它也总是被分成一种颜色的相同点，似乎是随机的。

更新：警告消息已消失。现在我只得到物种分数不可用的警告，这是明确的，因为我使用距离矩阵进行计算

更新：[看下图，下划线后面的字母是情节类型，与我在colfactor中使用的类型相同。显然结果是关闭的：/] [2]

更新：[标签的最后一个字母指定了我的表格中的类型，因此strR应该是所有三角形所在的位置，因为两个标签命令都相同] [1]

标记R中的NMDS图

0 个答案: