Question

我正在计算某些指标的贝叶斯后验分布（Layman，CA，Arrington，DA，Montana，CG，＆amp; Post，DM（2007）Can稳定同位素比率能否提供社区范围的营养结构测量？生态学， 88,42-48。）有一个分组变量G所以我希望获得每个G级别的后验分布

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

require(siar)
data(geese2demo)

gee <-data.frame(geese2demo)

summary(gee)
unique(gee$Group)

# There are 8 groups

me <- siber.hull.metrics(gee$d13CPl,gee$d15NPl,gee$Group,R=10^3)

# I expect 8000 rows, 1000 per group

nrow(me)

siber.hull.metrics

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

Answer 1

我认为你误解了构成SIBER的两种方法。

计算凸包区域，dN_range，到质心距离等的“基于船体的”方法只能适用于整个社区。也就是说，在您的情况下，它将向所有八个社区成员估算TA，CD等。您无法使用此方法在社区的各个社区成员上计算这些指标。

如果您希望将8个社区成员相互比较，那么您将需要使用siber.ellipses（）函数。

这是混淆的常见原因。我有两个播客解释这个，以及一些示例数据文件：http://www.tcd.ie/Zoology/research/research/theoretical/Rpodcasts.php#siber

感谢在Twitter上向我发誓，并再次感谢bug提示。我已将此记录为正在进行的工作的一部分，以推出大幅更新的软件包。

欢呼声

安德鲁杰克逊 @yodacomplex

从siar包中了解siber.hull.metrics函数

1 个答案: