3
からサンプリングする方法は
密度f (a )∝ c a d a − 1に従ってサンプリングしたい f(a )∝ cada − 1Γ (a )1(1 、∞ )(a)f(a)∝cada−1Γ(a)1(1,∞)(a) f(a) \propto \frac{c^a d^{a-1}}{\Gamma(a)} 1_{(1,\infty)}(a) ここで、cccとdddは厳密に正です。(動機:これは、ガンマ密度の形状パラメーターが均一な事前分布を持つ場合のギブスサンプリングに役立ちます。) 誰でもこの密度から簡単にサンプリングする方法を知っていますか?たぶんそれは標準的なもので、私が知らないことなのでしょうか? 私は、多かれ少なかれ仕事(モードを見つけるでしょう愚かな拒絶sampliingアルゴリズムと考えることができます*のF、サンプル(、U )大きな箱に均一からを[ 0 、10 * ] × [ 0 、F (A ∗)]およびu > f (a ))の場合は拒否しますが、(i)それはまったく効率的ではなく、(ii)f (a ∗)a∗a∗a^*fff(a,u)(a,u)(a,u)[0,10a∗]×[0,f(a∗)][0,10a∗]×[0,f(a∗)][0,10a^*]\times [0,f(a^*)]u>f(a)u>f(a)u>f(a)f(a∗)f(a∗)f(a^*)コンピュータが大きすぎて、適度に大きいおよびdでも簡単に処理できません。(大きなcとdのモードはおよそa = c dであることに注意してください。)cccdddcccddda=cda=cda=cd 助けてくれてありがとう!