这是我在这里的第一个问题,所以请原谅我(可能会或可能不会)犯的任何错误。
前提:
我得到了一个植被数据集,其中包含新旧观测的不同地块上的配对数据。我使用 'openxlsx'-package 来加载我的数据,并使用 'vegan'-package 来执行 NMDS,如下所示:
mydata <- read.xlsx(mydata)
mydataMDS <- metaMDS(mydata, k=2, trymax=500)
然后通过“envfit()”函数将结果用于模型,包括环境变量:
myenvdata <- read.xlsx(myenvdata)
mydataMDS_fit <- envfit(mydataMDS, myenvdata, perm=10000, na.rm=TRUE)
plot(mydataMDS, display="sites")
plot(mydataMDS_fit, p.max=0.01, axis=TRUE)
现在我有一个带有我的统计“mydataMDS”分析的图,包括由“mydataMDS_fit”R 计算产生的向量。
问题:
我想在这个情节中着色和连接某些点。由于“mydata”由同一图中不同时间的观察组成,我打算用一种颜色为旧观察的所有点着色,而用不同的颜色为所有新观察点着色。我已经阅读了有关添加列以对新旧观察结果进行分组的信息,但是由于我正在使用模型,因此没有列。如何编辑我的数据集(“mydata”、“mydataMDS”、“myenvdata”、“mydataMDS_fit”),以便以 2 种不同颜色(旧的一种颜色,新的一种颜色)显示新旧图,并连接成对的观察与线条?或者:是否有可能通过检查旧/新观察结果直接重新着色我的图形输出中的点?
(对不起,我觉得我的解释很复杂,但我还是希望有人能帮忙)