Bagaimana kita bisa mensimulasikan dari campuran geometris?

20

Jika f1,,fk kepadatan diketahui dari mana saya dapat mensimulasikan, yaitu, yang algoritma tersedia. dan jika produk

i=1kfi(x)αiα1,,αk>0
dapat diintegrasikan, apakah ada pendekatan umum untuk mensimulasikan dari kepadatan produk ini menggunakan simulator darifi ?
Xi'an
sumber
2
αi=1ϵ>0fiIifi10 < f i < ϵ I iI j = i j I if iPri(Ii)>1ϵ0<fi<ϵIiIj=ijIifisepertinya tidak ada hubungannya dengan .) Jadi, apa lagi yang bisa Anda ceritakan tentang ? f iIifi
whuber
1
(+10) Benar! Namun, menggunakan lebih kecil akan menyebabkan semua elemen rata dan karenanya mendukung tumpang tindih dukungan efektifnya ...αi
Xi'an
1
Seperti yang dikatakan whuber, sesak akan menjadi masalah, jadi saya akan mengambil transformasi (OR sampling istimewa) untuk membatalkan sesak sebelum membuat sampel acak. Ada satu pendekatan konstruktif yang saya pikir saya baca beberapa waktu lalu. Sec 10.7 dari link.springer.com/chapter/10.1007/978-1-4612-0209-7_10 Tidak yakin apakah diskritisasi juga dapat diterapkan di sini.
Henry.L

Jawaban:

3

Yah, tentu saja ada algoritma penerimaan-penolakan, yang akan saya terapkan untuk contoh Anda sebagai:

  1. (Inisialisasi) Untuk setiap , temukan . Sunting mencerminkan komentar Xi'an di bawah ini: Pilih distribusi yang sesuai dengan terkecil .A i = sup x { Π k j = 1 f j ( x ) α j / f i ( x ) } f i A iiAi=supx{Πj=1kfj(x)αj/fi(x)}fiAi
  2. Hasilkan dari .f ixfi
  3. Hitung .α=Πi=jkfj(x)αj/(Aifi(x))
  4. Hasilkan .uU(0,1)
  5. Jika , kembalikan , kalau tidak pergi ke 2.xuαx

Tergantung pada distribusinya, tentu saja, Anda mungkin memiliki tingkat penerimaan yang sangat rendah. Ketika itu terjadi, jumlah iterasi yang diharapkan sama dengan dipilih (dengan asumsi distribusi kontinu), jadi setidaknya Anda diperingatkan sebelumnya.Ai

Jbowman
sumber
3
(+1) Memang solusinya! Dengan asumsi batas memang ada untuk semua . Atau bahkan beberapa 's. Membandingkan [dengan asumsi mereka terbatas] juga dapat membantu dalam memilih paling efisien . i i A i f iAiiiAifi
Xi'an
1
Saya belum memikirkannya, tapi tentu saja Anda benar, sendiri sangat informatif, karena mereka juga sama dengan jumlah iterasi yang diharapkan untuk benar-benar menghasilkan angka acak jika Anda tetap menggunakan satu seluruh. Jadi, Anda ingin memilih distribusi dengan terkecil untuk digunakan sepanjang waktu. Saya akan mengedit jawaban sehingga poin Anda tidak hilang dalam komentar. i A iAiiAi
jbowman
Yaitu, dengan asumsi semua dinormalisasi dengan baik [untuk diintegrasikan ke satu], yang tidak harus merupakan kejadian standar. fi
Xi'an