R - 对数不完全 Gamma 函数的导数
R - Derivative of Log Incomplete Gamma Function
R 中可以使用以下函数：
gamma 计算伽马函数digamma 计算对数伽玛函数的导数pgamma 计算不完全伽马函数- ？计算对数不完全伽马函数的导数
我想知道什么函数可以计算对数不完全伽马函数的导数。我注意到 gsl 包有一个函数 gamma_inc 但不确定如何计算该函数对数的导数。
如果不存在函数，有没有一种简单的方法可以在R中逼近这个导数？
事实上，我的评论和 Ben Bolker 的评论一起给出了答案。
我建议使用链式求导法则：
log(u(x))' = u'(x) / u(x)
然后 Ben 指出 pgamma (CDF) 的导数就是 dgamma (PDF)。所以，我们有
dgamma(x) / pgamma(x)
正确定义的函数应该是
f <- function (x, shape, rate) dgamma(x, shape, rate) / pgamma(x, shape, rate)
R 中可以使用以下函数：
gamma计算伽马函数digamma计算对数伽玛函数的导数pgamma计算不完全伽马函数- ？计算对数不完全伽马函数的导数
我想知道什么函数可以计算对数不完全伽马函数的导数。我注意到 gsl 包有一个函数 gamma_inc 但不确定如何计算该函数对数的导数。
如果不存在函数，有没有一种简单的方法可以在R中逼近这个导数？
事实上，我的评论和 Ben Bolker 的评论一起给出了答案。
我建议使用链式求导法则：
log(u(x))' = u'(x) / u(x)
然后 Ben 指出 pgamma (CDF) 的导数就是 dgamma (PDF)。所以，我们有
dgamma(x) / pgamma(x)
正确定义的函数应该是
f <- function (x, shape, rate) dgamma(x, shape, rate) / pgamma(x, shape, rate)