问题描述
我正在尝试为整数参数实现 Gamma PDF 函数。首先，我用指数项实现了它，但存在溢出，所以我试图计算它的对数，但随着 alpha 和 beta 的增加，对数变为 -inf。
有没有想法如何从头开始实现 Gamma PDF？
这是我的实现：
log_gamma = a * np.log(b) + (a-1) * np.log(lam) - b * lam - np.log(factorial(a-1))
解决方法
这仅仅是因为计算 'a' 的阶乘，然后应用对数。相反，我们可以将阶乘的对数计算为数字 1 到 a 的对数之和。
