R:我如何在给定其他值的情况下找到beta函数的参数?

时间:2018-12-13 22:49:14

标签: r function math parameters probability

我有一个具有beta分布的变量Y:Beta(alpha,1/3)。

我必须找到alpha的值,例如: P(Y <= 0.416)= 0.2

有关更多详细信息,请在此处阅读我在数学stackexchange中的问题: https://math.stackexchange.com/questions/3038125/beta-distribution-find-the-parameter-alpha-of-mathcalbe-alpha-frac1

我写了这个函数(我假设x是alpha,函数的根):

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 solution in R中读到,该方法需要更强有力的假设以确保根的存在:f(下)* f(上)<0;但是我 有一个积极的功能,然后我不能使用此功能! R中是否存在替代功能? 谁能给我建议代码,如何用R查找alpha? 感谢您的提前帮助!

0 个答案:

没有答案