我正在尝试使用R-distance从横断面数据中获得丰富的鸟类估计值。我非常有信心我的代码是正确的,但我得到了“Inf”的丰度估计,我假设它意味着无限。谁能告诉我为什么? 这是丰富的代码:
dfunc <- F.dfunc.estim(M.ob, likelihood="halfnorm", w.hi=40)
plot(dfunc) # construct distance function
##
abun<- F.abund.estim(dfunc, detection.data=M.ob,
transect.data=M.t,
area=100, R=100, ci=0.95, plot.bs=TRUE)
M.ob是perp的观察结果。距横断面的距离,M.t是横断面的细节(例如长度)。
感谢您的帮助。如果我需要有关我的数据的更多详细信息,请告诉我。
编辑: 这些是“R-distance”包中的功能。
以下是我的数据:
> head(M.ob)
siteID groupsize dist sample
1 L1 1 5 1
2 L1 1 5 1
3 L1 1 5 1
4 L1 1 5 1
5 L1 1 5 1
6 L1 1 5 1
> str(M.ob)
'data.frame': 85 obs. of 4 variables:
$ siteID : Factor w/ 3 levels "L1","L2","L3": 1 1 1 1 1 1 1 1 1 1 ...
$ groupsize: num 1 1 1 1 1 1 1 1 1 1 ...
$ dist : num 5 5 5 5 5 5 5 5 10 10 ...
$ sample : num 1 1 1 1 1 1 1 1 1 1 ...
和M.t:
> head(M.t)
siteID Length lead.
1 L1 500 low
2 L2 500 low
3 L3 500 low
> str(M.t)
'data.frame': 3 obs. of 3 variables:
$ siteID: Factor w/ 3 levels "L1","L2","L3": 1 2 3
$ Length: num 500 500 500
$ lead. : Factor w/ 1 level "low": 1 1 1