beta 分布是否有一个共轭先验?

机器算法验证 贝塔分布 共轭先验
2022-01-30 13:00:44

我知道 beta 分布与二项式是共轭的。但是 beta 的共轭先验是什么?谢谢你。

4个回答

是的,它在指数族中有一个共轭先验。考虑三个参数族 对于某些值,这是可积的,尽管我还没有完全弄清楚哪个(我相信应该起作用 -对应于独立的指数分布所以这绝对有效,并且共轭更新涉及增加所以这表明也有效)。

π(α,βa,b,p){Γ(α+β)Γ(α)Γ(β)}pexp(aα+bβ).
(a,b,p)p0a<0,b<0p=0pp>0

问题,并且至少部分原因没有人使用它,是 即归一化常数没有封闭形式。

00{Γ(α+β)Γ(α)Γ(β)}pexp(aα+bβ)=?

看来你已经放弃了共轭。只是为了记录,我见过人们做的一件事(但不记得确切的位置,抱歉)是这样的重新参数化。如果是有条件的 iid,给定,使得,请记住 因此,您可以重新参数化可能性并用作先验 X1,,Xnα,βXiα,βBeta(α,β)

E[Xiα,β]=αα+β=:μ
Var[Xiα,β]=αβ(α+β)2(α+β+1)=:σ2.
μσ2
σ2μU[0,μ(1μ)]μU[0,1].
现在您已准备好计算后验并通过您最喜欢的计算方法对其进行探索。

理论上β分布应该有一个共轭先验。这是因为

然而推导看起来很困难,引用A Bouchard-Cote 的指数族和共轭先验

一个重要的观察是,这个配方并不总是产生一个计算上易于处理的共轭先验。

与此一致,在D Fink 的 A Compendium of Conjugate Priors中没有 Beta 分布的先验

Robert 和 Casella (RC) 碰巧在他们的书的示例 3.6(第 71 - 75 页)中描述了 beta 分布的共轭先验族,在 R 中介绍蒙特卡罗方法,Springer,2010。但是,他们引用了结果而没有引用来源。

为响应 gung 的详细要求而添加。RC 声明对于分布,共轭先验是“...的形式B(α,β)

π(α,β){Γ(α+β)Γ(α)Γ(β)}λx0αy0β

其中是超参数,因为后验等于{λ,x0,y0}

π(α,β|x){Γ(α+β)Γ(α)Γ(β)}λ(xx0)α((1x)y0)β."

该示例的其余部分涉及从进行重要性采样,以计算的边际似然。π(α,β|x)x