R 带有数据子集向量的联合图
我正在使用R软件包vegan从社区数据生成NMDS排序图,并希望包括长度与所选物种重要性相对应的向量(即来自原点的箭头)。我如何才能将显示的箭头仅限于那些物种,比如说,在数据的前四分位?我可以计算每个向量的长度,但不知道如何将打印的箭头限制为符合所需标准的箭头。比如说,R 带有数据子集向量的联合图,r,plot,vegan,multivariate-testing,R,Plot,Vegan,Multivariate Testing,我正在使用R软件包vegan从社区数据生成NMDS排序图,并希望包括长度与所选物种重要性相对应的向量(即来自原点的箭头)。我如何才能将显示的箭头仅限于那些物种,比如说,在数据的前四分位?我可以计算每个向量的长度,但不知道如何将打印的箭头限制为符合所需标准的箭头。比如说, require(vegan) 数据(沙丘) mdsImportantSpecies行标识哪些点位于前四分位。然后,您可以使用它(如图所示)仅绘制这些点 ImportantSpecies = which(hyp > qua
require(vegan)
数据(沙丘)
mdsImportantSpecies
行标识哪些点位于前四分位。然后,您可以使用它(如图所示)仅绘制这些点
ImportantSpecies = which(hyp > quantile(hyp, 0.75)
plot(mds$points[ImportantSpecies,1], mds$point[ImportantSpecies,2])
arrows(0, 0, mds$species[ImportantSpecies,1],
mds$species[ImportantSpecies,2], col = "grey50")
(1) 在排序图中应始终使用
asp=1
。如果您使用纯素函数并发出plot(mds,dis=“site”)
,则会自动执行此操作,但如果您想使用自己更麻烦的代码,请添加asp=1
。(2) 使用hyp谢谢你,Jari。我真的很感谢你提供的信息,我要求并自愿提供这些信息。至于物种得分的箭头——你是对的:充其量只是误导。我试图让一位比我更不了解情况的同事满意。很高兴你让我想起来了。