r - R中的累积曲线

Question

我有几个月内 4 个地点的物种数据。我已经使用包 vegan 成功创建了累积图在 R 中，但我想在一张图上绘制所有 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)

此代码片段仅加载 vegan 中的内置 BSI 数据集, 并创建一个包含 5 specaccum 的列表运行对象 specaccum(...)在 BCI 中的列的子集上。您不需要这样做，因为您已经拥有 specaccum 对象。

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