Dirichlet分布与多项分布的共轭性
Dirichlet分布與多項分布的共軛性
- 二項分布與Beta分布的共軛性
- Dirichlet分布與多項分布的共軛性
關于多項分布與Dirichlet分布的基礎可以參考:
UA MATH564 概率論 多項分布
UA MATH564 概率論 Dirichlet分布
二項分布與Beta分布的共軛性
因為多項分布與Dirichlet分布分別是二項分布與Beta分布從一元到多元的推廣,所以在探討多項分布與Dirichlet分布的共軛性之前,我們有必要先看看二項分布與Beta分布的共軛性:
假設N~Binom(n,p)N\sim Binom(n,p)N~Binom(n,p),其中p~Beta(α,β)p \sim Beta(\alpha,\beta)p~Beta(α,β),則
π(p∣N=m)=π(N=m∣p)π(p)m(N=m)=Cnm(1?p)n?mpmΓ(α+β)pα?1(1?p)β?1Γ(α)Γ(β)∫01Γ(α+β)Γ(α)Γ(β)Cnm(1?p)n?mpmpα?1(1?p)β?1dp∝pm+α?1(1?p)n?m+β?1\pi(p|N=m)=\frac{\pi(N=m|p)\pi(p)}{m(N=m)} \\ =\frac{C_n^m(1-p)^{n-m}p^m\Gamma(\alpha+\beta)p^{\alpha-1}(1-p)^{\beta-1}}{\Gamma(\alpha)\Gamma(\beta)\int_{0}^1 \frac{\Gamma(\alpha+\beta)}{\Gamma(\alpha)\Gamma(\beta)}C_n^m(1-p)^{n-m}p^mp^{\alpha-1}(1-p)^{\beta-1}dp} \\ \propto p^{m+\alpha-1}(1-p)^{n-m+\beta-1}π(p∣N=m)=m(N=m)π(N=m∣p)π(p)?=Γ(α)Γ(β)∫01?Γ(α)Γ(β)Γ(α+β)?Cnm?(1?p)n?mpmpα?1(1?p)β?1dpCnm?(1?p)n?mpmΓ(α+β)pα?1(1?p)β?1?∝pm+α?1(1?p)n?m+β?1
這是Beta分布的kernel,因此
p∣N=m~Beta(α+m,β+n?m)p|N=m \sim Beta(\alpha+m,\beta+n-m)p∣N=m~Beta(α+m,β+n?m)
Dirichlet分布與多項分布的共軛性
假設隨機向量X=(X1,?,Xd)X=(X_1,\cdots,X_d)X=(X1?,?,Xd?)服從多項分布MN(n,p)MN(n,p)MN(n,p),其中p=(p1,?,pd)Tp=(p_1,\cdots,p_d)^Tp=(p1?,?,pd?)T,并且p~Dir(α)p \sim Dir(\alpha)p~Dir(α),其中α=(α1,?,αd)T\alpha=(\alpha_1,\cdots,\alpha_d)^Tα=(α1?,?,αd?)T;
P(X=x∣p)=n!x1!?xd!∏i=1dpixi,x∈Δd?1P(X=x|p)=\frac{n!}{x_1!\cdots x_d!}\prod_{i=1}^d p_i^{x_i},x \in \Delta^{d-1}P(X=x∣p)=x1?!?xd?!n!?i=1∏d?pixi??,x∈Δd?1
其中x1+?+xd=nx_1+\cdots+x_d=nx1?+?+xd?=n,且所有的xix_ixi?都是非負整數;
π(p)=Γ(∑i=1dαi)∏i=1dΓ(αi)∏i=1dpiαi\pi(p)=\frac{\Gamma(\sum_{i=1}^d \alpha_i)}{\prod_{i=1}^d \Gamma(\alpha_i)}\prod_{i=1}^d p_i^{\alpha_i}π(p)=∏i=1d?Γ(αi?)Γ(∑i=1d?αi?)?i=1∏d?piαi??
且p∈Δd?1p \in \Delta^{d-1}p∈Δd?1,
Δd?1={x∈Rd:x1+?+xd=1,x≥0}\Delta^{d-1}=\{x \in \mathbb{R}^d:x_1+\cdots+x_d=1,x \ge 0\}Δd?1={x∈Rd:x1?+?+xd?=1,x≥0}
于是
π(p∣X=x)∝P(X=x∣p)π(p)∝∏i=1dpiαi+xi\pi(p|X=x) \propto P(X=x|p)\pi(p) \propto \prod_{i=1}^d p_i^{\alpha_i+x_i}π(p∣X=x)∝P(X=x∣p)π(p)∝i=1∏d?piαi?+xi??
因此
p∣X=x~Dir(α+x)p|X=x \sim Dir(\alpha+x)p∣X=x~Dir(α+x)
后驗均值為
E[pi∣X]=αi+xin+∑i=1dαi,i=1,?,dE[p_i|X]=\frac{\alpha_i+x_i}{n+\sum_{i=1}^d \alpha_i},i=1,\cdots,dE[pi?∣X]=n+∑i=1d?αi?αi?+xi??,i=1,?,d
總結
以上是生活随笔為你收集整理的Dirichlet分布与多项分布的共轭性的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 物理光学8 多波束干涉
- 下一篇: 马尔可夫“折棍子”过程 Markovia