У меня есть переменная Y с бета-распределением: бета (альфа, 1/3).
Мне нужно найти значение альфа, например: P (Y <= 0,416) = 0,2 </p>
для получения более подробной информации прочитайте мой вопрос в математике stackexchange здесь: https://math.stackexchange.com/questions/3038125/beta-distribution-find-the-parameter-alpha-of-mathcalbe-alpha-frac1
Я написал эту функцию (я предположил, что x это альфа, корень функции):
f=function(x){
pbeta(0.416,x,1/3)
}
и я попытался использовать uniroot:
uniroot(f,interval=c(0,5),tol=1e-5)
Я не понимаю это сообщение: Error in uniroot(f, interval = c(0, 5), tol = 1e-05) : f() values at end points not of opposite sign.
Я прочитал здесь Решение Uniroot в R , что этоМетод method нуждается в более строгом допущении для обеспечения существования корня: f (нижний) * f (верхний) <0;но у меня есть положительная функция, и тогда я не мог использовать эту функцию!Существует ли альтернативная функция в R?Может кто-нибудь предложить мне код, как найти альфа с R?Заранее спасибо за помощь !! </p>