将站点添加到R中的DCA

时间:2013-08-01 13:01:25

标签: r vegan

我使用下面的脚本创建了一个显示物种的Detrended对应分析(DCA),但我想要将网站添加到其中。我知道我可以使用简单的脚本plot(vare.dca)来做到这一点,但是一旦我改进了图形,我就不确定如何做到这一点。我是R的新手,如果这是一个非常简单的问题,请道歉。非常感谢您的帮助!

bugs<-read.csv(file= "bugs.csv", row.names='Sites', sep=",", header=TRUE)
env<-read.csv(file= "envir.csv", row.names='Sites', sep=",", header=TRUE) 
library(vegan)

shnam <- make.cepnames(names(bugs)) # abbreviate Latin names
identify(pl, "sp", labels=shnam)
plot(vare.dca, dis="sp", type="n")
sel <- orditorp (vare.dca, dis="species", lab=shnam, pcol = "grey", pch="+")

&#39;错误的数据&#39;和&#39; env&#39;在这里:https://gist.github.com/anonymous/dcf0fad859eb879014b2

1 个答案:

答案 0 :(得分:1)

一般方法是使用points()方法,如

points(vare.dca, display = "sites", col = "red")

如果您需要网站/样本标签,请使用text()方法。

也可以使用

orditorp(),只需将参数更改为display = "sites",您可以通过阅读?orditorp确定参数。

还有许多其他选项,尤其是ordipointlabel(),它们可以更加努力地阻止各种标签的重叠。

我在素食主义者中撰写了一些关于这些功能的博客文章,您可以参考这些帖子获取更多解释和示例

  1. Decluttering ordination plots in vegan part 1: ordilabel()
  2. Decluttering ordination plots in vegan part 2: orditorp()
  3. Decluttering ordination plots part 3: ordipointlabel()