digamma 함수의 근본 원인은 R

Question

입니다. 최대 우도 추정량에 대해 연구 중이며 매개 변수 중 하나는 digamma 함수를 사용하여 추정됩니다. 나는 방정식을 풀기 위해 uniroot를 사용하려하지만 그렇게 할 수 없다.

이

dig = function(alpha){ 
    digamma(2 + alpha) - digamma(alpha) - (1/(2+alpha)) + (2/(2+alpha)) 
} 

curve(dig, from = 0, to = 10) 
uniroot(dig, lower = 0, upper = 10) 

이 다음과 같은 오류가 발생합니다 : 여기에 내 코드의

이

Error in uniroot(dig, lower = 0, upper = 10) : f.lower = f(lower) is NA 
In addition: Warning messages: 
1: In digamma(alpha) : NaNs produced 
2: In digamma(alpha) : NaNs produced 

첫 번째 오류는 종류의 곡선에 따라 의미가 있지만, 두 번째는 나에게 붙어있다. digamma 함수의 근원을 찾는 방법을 오해하고 있거나 R에 숫자로 된 패키지 (아마도 rootsolve?)가있을 수 있다는 것이 완전히 가능할 수 있습니다. 내가 무엇을 여기에서 놓치고 있는지 명확히하지 않고있다 - 어떤 팁이라도 인정 될 것이다. 감사!

Answer 1

는 다음

curve(dig, from = 0.01, to = 10) 
uniroot(dig, lower = 0.01, upper = 10) 

Error in uniroot(dig, lower = 0.01, upper = 10) : f() values at end points not of opposite sign