Berapa nilai yang diharapkan dari distribusi Dirichlet yang dimodifikasi? (masalah integrasi)

14

Sangat mudah untuk menghasilkan variabel acak dengan distribusi Dirichlet menggunakan variabel Gamma dengan parameter skala yang sama. Jika:

XiGamma(αi,β)

Kemudian:

(X1jXj,,XnjXj)Dirichlet(α1,,αn)

Masalah Apa yang terjadi jika parameter skala tidak sama?

XiGamma(αi,βi)

Lalu apa distribusi variabel ini?

(X1jXj,,XnjXj)?

Bagi saya itu akan cukup untuk mengetahui nilai yang diharapkan dari distribusi ini.
Saya membutuhkan rumus aljabar tertutup perkiraan yang dapat dievaluasi dengan sangat cepat oleh komputer.
Katakanlah perkiraan dengan ketelitian 0,01 sudah cukup.
Anda dapat mengasumsikan bahwa:

αi,βiN

Catatan Singkatnya, tugasnya adalah untuk menemukan perkiraan integral ini:

f(α,β)=R+nx1jxjjβjαjΓ(αj)xjαj1eβjxjdx1dxn

Łukasz Lew
sumber
1
@ Łukasz Bisakah Anda mengatakan sesuatu tentang parameter , , dan ? Dimungkinkan untuk mendapatkan ekspresi yang tepat untuk dan dengan demikian memperkirakan ekspektasi rasio, tetapi untuk kombinasi tertentu dari parameter seseorang dapat mengeksploitasi pendekatan Normal atau saddlepoint dengan sedikit kerja. Saya tidak berpikir akan ada metode pendekatan universal, yang mengapa pembatasan tambahan akan diterima. nαiβijXj
whuber
danj X j berkorelasi sehingga kita harus memperkirakan integral itu sendiri. α i sering berupa angka kecil seperti 1 atau 2 dan kadang-kadang sebesar 10.000. Demikian pula dengan β i tetapi biasanya 10 kali lebih besar dari α i . X1jXjαiβiαi
Łukasz Lew
Masalahnya dengan kecil i . Jika semua α i besar maka perkiraan baik dari keseluruhan integral adalah: α 1 / β 1αiαiα1/β1jαj/βj
Łukasz Lew
@ Łukasz Jika Anda perlu mengevaluasi ekspresi ekspektasi, mengapa Anda memerlukan formula aljabar? Saya berpikir dalam menerapkan beberapa trik numerik untuk mendapatkan harapan tetapi saya butuh umpan balik :)
deps_stats
Saya perlu mengevaluasi berkali-kali dalam program saya. Itu harus sangat cepat, yaitu tidak ada loop dan lebih disukai tidak terlalu banyak divisi.
Łukasz Lew

Jawaban:

2

Hanya komentar awal, jika Anda ingin kecepatan komputasi, Anda biasanya harus mengorbankan akurasi. "Lebih akurat" = "Lebih banyak waktu" secara umum. Anyways di sini adalah perkiraan urutan kedua, harus meningkatkan kira-kira "mentah" yang Anda sarankan dalam komentar Anda di atas:

E(XjiXi)E[Xj]E[iXi]cov[iXi,Xj]E[iXi]2+E[Xj]E[iXi]3Var[iXi]
=αjiβjβiαi×[11(iβjβiαi)+1(iαiβi)2(iαiβi2)]

EDIT Penjelasan untuk ekspansi di atas diminta. Jawaban singkatnya adalah wikipedia . Jawaban panjang diberikan di bawah ini.

f(x,y)=xyfXE(X)YE(Y)

2fx2=0
2fxy=1y2
2fy2=2xy3

Dan seri taylor hingga urutan kedua diberikan oleh:

xyμxμy+12(1μy22(xμx)(yμy)+2μxμy3(yμy)2)

Mengambil harapan hasil:

E[xy]μxμy1μy2E[(xμx)(yμy)]+μxμy3E[(yμy)2]

Itulah jawaban yang saya berikan. (meskipun saya awalnya lupa tanda minus pada periode kedua)

probabilityislogic
sumber
Ini persis seperti yang saya butuhkan. Bisakah Anda menjelaskan bagaimana Anda mendapatkan ekspansi ini? Saya mencoba dalam banyak hal dan tidak dapat melakukan itu ...
Lewukasz Lew