是这样的：

我是一个贝叶斯的学生(我的问题不涉及贝叶斯）。 老师用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