我有几个月来在4个地点的物种数据。我已经使用vegan
R中的包成功创建了累积图,但是我想在一个图上绘制所有4个站点。
最初,我有一个包含所有站点和月份的数据表,但是当我绘制specaccum
结果时,无论位置如何,结果都是所有数据的累积曲线。
因此,我将每个站点拆分到一个单独的数据表中,然后将其加载到R中。在每个数据表中,第一行是物种名称,下面的每个其他行是一个月。
例如,我加载了我的网站“ FMR”之一的数据。然后,我执行以下操作:
FMR <-specaccum(FMRJ2F, "random")
plot(FMR)
我也同样对我的其他网站,PP
,DUX
,PM
。如何将所有4条线放在一个图上?
您可以只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
。
本文收集自互联网,转载请注明来源。
如有侵权,请联系[email protected] 删除。
我来说两句