紧急求助，想知道，原分布的系数

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

是这样的：

我是一个贝叶斯的学生(我的问题不涉及贝叶斯）。 老师用R 交给我们如何来算一个DISTRIBUTION， 比如BETA（2,10）的  HPD的interval（Highest Probability Density)。 CODE 如下：

alpha<-2
beta<-10

# HPD interval
f <- function(a)
{
b<-qbeta(0.95+pbeta(a,alpha,beta),alpha,beta)
return((dbeta(a,alpha,beta)-dbeta(b,alpha,beta))^2)
}
hpdmin <- optimize(f,lower=0,upper=qbeta(0.05,alpha,beta))$minimum
hpdmax <- qbeta(pbeta(hpdmin,alpha,beta)+0.95,alpha,beta)
c(hpdmin,hpdmax)



但是我现在需要做的，是把上述过程逆过来。  也就是我现在知道HPD的INTERVAL， 但是要推算出来， 是从哪个BETA分布来的。 也就是说BETA(多少，多少）。   有哪位高人能给指点一下， 现在这里谢谢啦。

我目前手头这个问题， HPD 的INTERVAL 是 0.82 , 0.9,   mode（也就是density分布图最高那点）是0.85.                  根据BETA的性质：mode<-(alpha-1)/(alpha+beta-2)根据已知MODE（0.85）我自己能退出来ALPHA 和BETA之间的关系： beta <- (3*alpha +14)/17

求高手帮忙啦。谢谢 如果可能的话，加我QQ 也OK， 84739375

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：紧急求助 急求助 distribution Probability Interval 求助 系数

试试递归，或许可以。

interval应该是alpha和beta的函数，然后可以解出另外一个关系，解方程组就可以了。