1

R中提供以下功能：

  • gamma计算伽玛函数
  • digamma计算对数伽玛函数的导数
  • pgamma计算不完全伽马函数
  • ? 计算对数不完全伽马函数的导数

我想知道什么函数可以计算对数不完全伽马函数的导数。我注意到这个gsl包有一个函数gamma_inc，但不知道如何计算这个函数对数的导数。

如果不存在函数，是否有一种简单的方法来近似这个导数R？

4

1 回答 1

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)
于 2016-11-19T18:28:34.150 回答