R 平均物种累积曲线?

R 平均物种累积曲线?,r,vegan,R,Vegan,我感兴趣的是在多个调查中生成每分钟调查的平均鱼种累积曲线 我理解随着更多的地点被取样,specaccum素食主义者通常是如何看待物种的积累的 这个特殊的样本设计是研究当使用调查工具摄像机调查鱼类时,物种是如何随着时间累积的。这个想法是,相机在水中浸泡的时间越长,就应该观察到更多的鱼类在水中游动。然后我多次重复这个过程,最后得到了几个类似于BCI数据集的表,每一行都是一分钟的时间 一次15分钟的试验可能如下所示: library(vegan) counts = c(0,0,0,0,0,

我感兴趣的是在多个调查中生成每分钟调查的平均鱼种累积曲线

我理解随着更多的地点被取样,specaccum素食主义者通常是如何看待物种的积累的

这个特殊的样本设计是研究当使用调查工具摄像机调查鱼类时,物种是如何随着时间累积的。这个想法是,相机在水中浸泡的时间越长,就应该观察到更多的鱼类在水中游动。然后我多次重复这个过程,最后得到了几个类似于BCI数据集的表,每一行都是一分钟的时间

一次15分钟的试验可能如下所示:

library(vegan)

counts = c(0,0,0,0,0,
       0,0,0,0,0,
       1,2,3,4,5)



#Survey over 15 minutes
fish = data.frame(Time = seq(1,15,1),
              S.mys = sample(counts,15,replace = T),
              S.ros = sample(counts,15,replace = T),
              S.pau = sample(counts,15,replace = T),
              S.ens = sample(counts,15,replace = T),
              S.con = sample(counts,15,replace = T),
              S.sax = sample(counts,15,replace = T),
              S.sim = sample(counts,15,replace = T))

#taken from the help file example:
sp1 = specaccum(fish)
sp2 = specaccum(fish, 'random')
plot(sp1, ci.type = 'poly',col = 'blue',lwd = 2, ci.lty = 0, ci.col = 'lightblue')
boxplot(sp2,add = T, col = 'peachpuff')

现在,如果我已经进行了多次试验,几个15分钟的时间段,有没有办法平均物种累积结果?

如果你想得到每项研究的平均值,你可以将所有研究数据放在一个data.frame中,并使用data.table这样聚合itre


谢谢@CCurtis,我认为这是一个很好的开始。我仍在弄清楚这是否是我应该采取的方法,根据几次调查计算一条平均物种累积曲线——或者是否应该计算单个曲线,然后在事后平均。无论如何,我需要学习应用函数。如果你平均,你就积累。这意味着:平均值比单个观测值有更多的物种。如果物种积累要有意义,就必须通过简单的观察,然后对曲线进行平均。rowMeanscbindsp1$richness,sp2$richness是你想要的吗?嗨@JariOksanen,我相信这一点和你下面的其他有用的评论都回答了我问题的编码和理论部分
require(data.table)
fish<-cbind(fish,Study=rep(1,15))   #Add study column
setDT(fish)                         #Add class data.table
fish[, lapply(.SD, mean), by=Study] #Aggregate data by study and take mean

   Study Time    S.mys S.ros    S.pau     S.ens S.con     S.sax
1:     1    8 1.266667   0.8 1.666667 0.2666667   1.6 0.6666667
      S.sim 
1: 1.133333