2
Lấy mẫu từ phân phối bivariate với mật độ đã biết bằng MCMC
Tôi đã cố gắng mô phỏng từ mật độ bivariate bằng thuật toán Metropolis trong R và không gặp may. Mật độ có thể được biểu thị là , trong đó là phân phối Singh-Maddalap(x,y)p(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) p(x)=aqxa−1ba(1+(xb)a)1+qp(x)=aqxa−1ba(1+(xb)a)1+qp(x)=\dfrac{aq x^{a-1}}{b^a (1 + (\frac{x}{b})^a)^{1+q}} với các tham số , , và là log-normal với …