例による答えは次のとおりです。次の階層モデルがあるとします
for groups and observations aグループおよび既知の値および。
限界尤度は
積分の次元は
Y私グラム〜I N DN(θg、1 )θg〜I N DN(μ 、τ2)μ | τ2〜N(m 、τ2/ k)τ2〜IG (a 、b )
g= 1 、… 、Gi = 1 、… 、ngm 、k 、a 、by=(y1,1,…,yn1,1,y1,2,…,yn2,2,…,y1,G,…,ynG,G),
G + 2 Gp(y)=∫⋯∫∏g=1G[∏i=1ngN(yig;θg,1)]N(θg;μ,τ2)dθ1⋯dθGdμdτ2.
G+2およびが大きい場合、これは高次元の積分です。ほとんどの数値積分手法では、この積分の妥当な近似値を取得するために、極端な数のサンプルまたは反復が必要になります。
G
この積分はたまたま閉じた形で限界尤度を持っているので、数値積分手法が限界尤度をどれだけうまく推定できるかを評価できます。周辺尤度の計算が難しい理由を理解するには、たとえば、単一の観測値、単一のグループ、とが既知であるなど、単純なものから始めます。徐々に問題を徐々に難しくしていくことができます。数値積分手法が真実とどのように関連しているかをご覧ください。問題の次元、つまりが増加するにつれて、同じ精度を得るには、ますます多くのサンプルまたは反復が必要になることに気づくでしょう。最後に、σ 2 G Y I G I N D 〜 P O (E θ G)μσ2GYig∼indPo(eθg)そして今、あなたは閉じた形をしていないわずかな可能性を持っています。真実を知ったときの経験に基づいて、真実を知らない場合の数値の見積もりをどれだけ信じますか?数値の見積もりにあまり自信がないと思います。