2
Pengambilan sampel dari distribusi bivariat dengan kerapatan yang diketahui menggunakan MCMC
Saya mencoba mensimulasikan dari kepadatan bivariat menggunakan algoritma Metropolis di R dan tidak beruntung. Kepadatan dapat dinyatakan sebagai , di mana adalah distribusi Singh-Maddalap ( x , y)hal(x,y)p(x,y)p(y|x)p(x)p(y|x)p(x)p(y|x)p(x)p(x)p(x)p(x) hal(x)=aqxa−1ba(1+(xb)a)1 + qp(x)=aqxSebuah-1bSebuah(1+(xb)Sebuah)1+qp(x)=\dfrac{aq x^{a-1}}{b^a (1 + (\frac{x}{b})^a)^{1+q}} dengan parameter , , , dan adalah log-normal dengan log-mean sebagai fraksi , dan log-sd …