元の投稿には、間違った記号やマトリックスの削除など、いくつかの数学のタイプミスがあるため、最初から始めます。V
あなたは前に指定したと可能性:P (Y | β )= N(Bのβ 、σ 2 I )。p(β)=N(0,σ2V)p(y|β)=N(Bβ,σ2I)
これらのそれぞれを純粋にβに依存する内の項の式として記述し、βに関連しないすべての項を単一の定数にグループ化できます。expββ
ログp (β)+ const = − 12つのσ2βTV− 1β
ログp (y| β)+ const = − 12つのσ2(βTBTBの β- 2 YTBの β)(なお、常に)yTBの β= βTBTy
これらをログスペースに追加し、同様の用語を収集すると、正規化されていないログが得られます
ログp (β| y)+ const = − 12つのσ2(βT(V− 1+ BTB )β- 2 YTBの β) (1)
バツTA x + xTCx = xT(A + C)xバツA 、C
β
logp(β|y)+const=(β−μp)TΛp(β−μp)=βTΛpβ−2μTpΛpβ+ μTpΛpμp
μp、Λp
まあ、検査式によって。(1)設定すると、このフォームによく似ています
Λp= V− 1+ BTB そして
μp= Λ− 1pBTy
詳細には、この置換により(1)から必要な各項が作成されることを示すことができます。
βTΛpβ= βT(V− 1+ BTB )β
μTpΛpβ= (Λ− 1pBTy)TΛpβ= yTB Λ− 1pΛpβ= yTBの β
(A B )T= BTあT(Λ− 1p)T= Λ− 1pΛp
μTpΛpμpμp、Λp
logp(β|y)+const=−12σ2[(β−μp)TΛp(β−μp)−μpΛpμp]
β