MathLabs
定理証明済み

事後分布に対するベイズの定理

内容

θ\theta が事前密度 π(θ)\pi(\theta) を持ち、データ xx が尤度 L(x∣θ)L(x\mid\theta) を持つとき、xx を与えたときの θ\theta の事後密度は π(θ∣x)=L(x∣θ)π(θ)m(x)\pi(\theta\mid x)=\dfrac{L(x\mid\theta)\pi(\theta)}{m(x)} である。ここで m(x)=∫L(x∣θ′)π(θ′) dθ′m(x)=\int L(x\mid\theta')\pi(\theta')\,d\theta'。同値に π(θ∣x)∝L(x∣θ)π(θ)\pi(\theta\mid x)\propto L(x\mid\theta)\pi(\theta) とも書ける。

なぜ正しいのか?

これは θ\theta と xx に同時に適用した条件付き密度の定義に他ならない:事後分布は (θ,x)(\theta,x) の同時密度を xx だけの周辺密度で割ったものであり、同時密度は事前分布と尤度の積に分解される。

証明の概略

(θ,x)(\theta,x) の同時密度は二通りに分解できる:π(θ)L(x∣θ)\pi(\theta)L(x\mid\theta)(尤度を θ\theta を与えたときの xx の条件付き密度と定義することによる、事前分布と尤度の積)と、π(θ∣x)m(x)\pi(\theta\mid x)m(x)(事後分布を xx を与えたときの θ\theta の条件付き密度と定義することによる、事後分布と xx の周辺密度の積)である。両方の式は同じ同時密度に等しいので、π(θ)L(x∣θ)=π(θ∣x)m(x)\pi(\theta)L(x\mid\theta)=\pi(\theta\mid x)m(x)。

π(θ∣x)\pi(\theta\mid x) について解くと、m(x)>0m(x)>0 である限り π(θ∣x)=L(x∣θ)π(θ)m(x)\pi(\theta\mid x)=\dfrac{L(x\mid\theta)\pi(\theta)}{m(x)} が得られる。

残るは m(x)=∫L(x∣θ′)π(θ′) dθ′m(x)=\int L(x\mid\theta')\pi(\theta')\,d\theta' がまさに正しい正規化定数であることの確認である:分解 π(θ)L(x∣θ)=π(θ∣x)m(x)\pi(\theta)L(x\mid\theta)=\pi(\theta\mid x)m(x) の両辺を θ\theta について積分すると、左辺は定義により ∫π(θ)L(x∣θ) dθ=m(x)\int \pi(\theta)L(x\mid\theta)\,d\theta=m(x) となり、右辺は π(⋅∣x)\pi(\cdot\mid x) が確率密度であることから m(x)∫π(θ∣x) dθ=m(x)×1m(x)\int\pi(\theta\mid x)\,d\theta=m(x)\times1 となる。両辺が一致することから m(x)m(x) の整合性が確かめられ、π(θ∣x)\pi(\theta\mid x) が密度として要求される通り積分して 11 になることが確認される。

この定理を使うトピック

ステップごとの証明

この定理のステップごとの証明はまだありません。

参考文献

  1. Andrew Gelman, John B. Carlin, Hal S. Stern, David B. Dunson, Aki Vehtari, Donald B. Rubin (2013). Bayesian Data Analysis (3rd ed.)
  2. Matthew D. Hoffman, Andrew Gelman (2014). The No-U-Turn Sampler: Adaptively Setting Path Lengths in Hamiltonian Monte Carlo · arXiv:1111.4246