Корни функции 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 есть числовой пакет (возможно, rootolve?), Который может помочь. Не уверен, что мне здесь не хватает - любые советы будут оценены. Благодаря!

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