NetLogo计算测试版分布/函数

Question

我正在使用 NetLogo 中的统计扩展来计算 beta 函数:

设置 z (统计数据:beta (H + 1) (T + 1))

其中 H 和 T 是抛硬币中正面和反面的数量。stats 扩展的使用不是必需的，我在使用 beta 函数的阶乘表达式时也遇到同样的问题。

问题是，当 H + T > 168 时，NetLogo 报告 z = 0 并且有一些程序我无法执行[特别是 beta 发行版]

是否有任何方法可以近似 Netlogo 中的 beta 函数(或分布)，这样就不会遇到此问题？

Answer 1

基于 Chris 的回答，统计扩展确实有一个 Gamma 函数的对数，stats:loggamma。它处理的参数远远超过 171，因此 a 和 b 相加也可以远远超过 172。

observer> show stats:beta 85 86
observer: 1.2864854397253604E-52
observer> show exp (stats:loggamma 85 + stats:loggamma 86 - stats:loggamma (85 + 86))
observer: 1.2864854397251408E-52
observer> show stats:beta 86 86
observer: 0
observer> show exp (stats:loggamma 86 + stats:loggamma 86 - stats:loggamma (86 + 86))
observer: 6.394810665301235E-53
observer> show exp (stats:loggamma 200 + stats:loggamma 200 - stats:loggamma (200 + 200))
observer: 9.713217247613997E-122

新版本的统计扩展 (v1.4.0) 已发布，其中包含使用日志的“bigBeta”函数。