R中提供以下功能：
gamma计算gamma功能digamma计算log gamma函数的导数pgamma计算不完整的伽玛函数我想知道什么函数可以计算log不完整伽玛函数的导数。我注意到gsl包有一个函数gamma_inc但不知道如何计算这个函数的日志的导数。
如果不存在任何函数，是否有一种简单的方法可以在R中近似这个导数？
答案 0 :(得分：8)
事实上，我的评论和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)