r语言 siar 代码_从siar包中了解siber.hull.metrics函数(Understanding siber.hull.metrics function from siar packag...

从siar包中了解siber.hull.metrics函数(Understanding siber.hull.metrics function from siar package)

我正在计算一些指标的贝叶斯后验分布(Layman，CA，Arrington，DA，Montana，CG，＆Post，DM(2007)Can稳定同位素比率可以提供社区范围的营养结构测量吗？生态学，88,42- 48.)使用分组变量G，因此我希望获得每个G级别的后验分布

我正在使用siar包中的函数siber.hull.metrics函数，

require(siar)

data(geese2demo)

gee

summary(gee)

unique(gee$Group)

# There are 8 groups

me

# I expect 8000 rows, 1000 per group

nrow(me)

siber.hull.metrics

我查看了函数的代码，但仍然不明白对组有什么影响，除了有一个bug：重复修复为10 ^ 4。

I am calculating the bayesian posterior distribution of some metrics (Layman, C.A., Arrington, D.A., Montana, C.G., & Post, D.M. (2007) Can stable isotope ratios provide for community-wide measures of trophic structure? Ecology, 88, 42-48.) with a grouping variable