实现Gamma PDF函数

Question

我正在尝试为整数参数实现Gamma PDF函数。首先，我用指数项实现了它，但有溢出，所以我试图计算它的对数，但随着α和β的增加，对数变成了-inf。有没有办法从头开始实现Gamma PDF？下面是我的实现：

log_gamma = a * np.log(b) + (a-1) * np.log(lam) - b * lam - np.log(factorial(a-1))

Answer 1

这很简单，因为要计算'a‘的阶乘，然后应用对数。相反，我们可以计算阶乘的对数作为数字1到a的对数之和。