R中提供以下功能:
gamma 计算伽马函数digamma 计算log gamma函数的导数pgamma 计算不完整的伽马函数我想知道什么函数可以计算log不完整伽马函数的导数.我注意到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)
| 归档时间： |
|
| 查看次数： |
531 次 |
| 最近记录： |