问题标签 [gamma-function]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
1 回答
27 浏览

r - R中的多精度伽玛函数

我需要计算 R 中涉及每个元素的 gamma 函数的总和。当 gamma 的参数增加时,我得到 NaN,我怀疑问题是 gamma 评估的数值问题。我已经阅读了 Rmpfr 和 gmp 的文档,但我只发现整数的阶乘。我也在这里发布代码,也许您对错误的来源有更好的了解。

w1 是一个总和为 1 的正实数数组,a1 和 b1 是正值向量。当 a1(和 a1mat)变长并具有高值(~150)时,A 变为 NaN。