R中的累积曲线

Question

我有几个月的4个地点的物种数据.我已经使用veganR中的包成功创建了累积图,但我想在一个图上绘制所有4个站点.

起初我有一个包含所有站点和月份的数据表,但是当我绘制specaccum结果时,所有数据的累积曲线都与站点无关.

因此,我将每个站点拆分为一个单独的数据表,我将其加载到R中.在每个数据表中,第一行是物种名称,下面的每一行都是一个月.

例如,我加载了我的一个网站"FMR"的数据.然后我做了以下事情:

FMR <-specaccum(FMRJ2F, "random")
plot(FMR)

我也同样对我的其他网站,PP,DUX,PM.如何将所有4条线放在一个图上？

Answer 1

你可以使用add=T参数plot.specaccum(...)

library(vegan)
data(BCI) 
df <- lapply(c(1,21,41,61,81),function(i)specaccum(BCI[,seq(i,i+19)], method="random"))
plot(df[[1]])
for (i in 2:5) plot(df[[i]],add=T, col=i)

此代码段仅加载内置的BSI数据集vegan,并specaccum通过specaccum(...)在BCI中的列的子集上运行来创建5个对象的列表.您不需要这个,因为您已经拥有了specaccum对象.

然后,我们创建第一个图,并添加每个新曲线add=T.