R -

Question

다음 기능 R에서 사용할 수있는 로그인 불완전 감마 함수의 유도체

digamma 로그 감마 함수의 도함수를 계산하는

• gamma 계산할 감마 함수
• pgamma 불완전 감마 함수를 계산
• ? 로그 불완전 감마 함수의 미분을 계산하는 방법

로그 불완전 감마 함수의 미분을 계산할 수있는 함수는 무엇입니까? gsl 패키지에는 gamma_inc 기능이 있지만이 기능 로그의 파생물을 계산하는 방법을 알지 못합니다.

기능이없는 경우 R에서이 미분을 근사 할 수있는 간단한 방법이 있습니까?

Answer 1

사실 내 의견과 벤 볼커가 함께 답을 알려줍니다.

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

다음 벤 pgamma (CDF)의 유도체는 단지 dgamma (PDF)이라고 지적 :

나는 유도의 체인 규칙의 사용을 제안했다. 그래서, 우리는

dgamma(x)/pgamma(x) 

적절히 정의 함수

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