首页 > 解决方案 > R:在给定其他值的情况下,我如何找到 beta 函数的参数?

问题描述

我有一个带有 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

我在这里读到 R 中的 Uniroot 解决方案,这种方法需要更强的假设来确保根的存在:f(lower) * f(upper) < 0; 但我有一个积极的功能,然后我不能使用这个功能!R中是否存在替代功能?任何人都可以建议我的代码,如何用 R 找到 alpha?我在这里先向您的帮助表示感谢!!

标签: rfunctionmathparametersprobability

解决方案


推荐阅读