我有一个具有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? 感谢您的提前帮助!