日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

Dirichlet分布与多项分布的共轭性

發布時間:2025/4/14 编程问答 42 豆豆
生活随笔 收集整理的這篇文章主要介紹了 Dirichlet分布与多项分布的共轭性 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

Dirichlet分布與多項分布的共軛性

  • 二項分布與Beta分布的共軛性
  • Dirichlet分布與多項分布的共軛性

關于多項分布與Dirichlet分布的基礎可以參考:
UA MATH564 概率論 多項分布
UA MATH564 概率論 Dirichlet分布

二項分布與Beta分布的共軛性

因為多項分布與Dirichlet分布分別是二項分布與Beta分布從一元到多元的推廣,所以在探討多項分布與Dirichlet分布的共軛性之前,我們有必要先看看二項分布與Beta分布的共軛性:

假設N~Binom(n,p)N\sim Binom(n,p)NBinom(n,p),其中p~Beta(α,β)p \sim Beta(\alpha,\beta)pBeta(α,β),則
π(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}π(pN=m)=m(N=m)π(N=mp)π(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)pN=mBeta(α+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)pDir(α),其中α=(α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=xp)=x1?!?xd?!n!?i=1d?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=1d?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={xRd:x1?+?+xd?=1,x0}

于是
π(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}π(pX=x)P(X=xp)π(p)i=1d?piαi?+xi??

因此
p∣X=x~Dir(α+x)p|X=x \sim Dir(\alpha+x)pX=xDir(α+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分布与多项分布的共轭性的全部內容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。