文章/答案/技术大牛

发布

社区首页 >问答首页 >对数不完全Gamma函数的r-导数

问对数不完全Gamma函数的r-导数
EN

Stack Overflow用户

提问于 2016-11-19 18:10:26

回答 1查看 697关注 0票数 1

在R中有下列功能：

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

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

回答已采纳

发布于 2016-11-19 18:28:34

事实上，我和Bolker的评论一起给出了答案。

我建议使用链式衍生规则：

log(u(x))' = u'(x) / u(x)

然后本指出pgamma的导数就是dgamma (PDF)。所以，我们有

dgamma(x) / pgamma(x)

正确定义的函数将是

f <- function (x, shape, rate) dgamma(x, shape, rate) / pgamma(x, shape, rate)

票数 8

页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持

原文链接：

https://stackoverflow.com/questions/40696292

复制

相似问题

问对数不完全Gamma函数的r-导数EN