Biarkan menjadi fungsi non-negatif. Saya tertarik menemukan sedemikian rupa sehingga Peringatan itu : yang bisa saya lakukan adalah sampel pada titik di . Aku bisa, bagaimanapun, memilih lokasi di mana saya sampel secara acak, jika saya menginginkannya.
Pertanyaan:
- Apakah mungkin untuk mendapatkan estimasi tidak bias setelah banyak sampel? Jika demikian, berapakah varian terkecil yang dapat dimiliki dari estimasi tersebut setelah sampel ?
- Jika tidak, prosedur apa yang tersedia untuk memperkirakan , dan apa waktu konvergensi yang terkait.
Seperti yang ditunjukkan oleh Douglas Zare dalam komentar, ini bisa sangat sulit dilakukan jika fungsinya mendekati nol atau sangat besar. Untungnya, fungsi yang perlu saya gunakan untuk ini dibatasi dari atas dan bawah, jadi anggaplah . Selain itu, kita juga dapat mengasumsikan bahwa adalah Lipschitz atau bahkan dapat dibedakan jika itu membantu.
sampling
monte-carlo
quantiles
quasi-monte-carlo
Robinson
sumber
sumber
Jawaban:
Seperti yang ditunjukkan kardinal dalam komentarnya, pertanyaan Anda dapat disajikan kembali sebagai berikut.
Dengan aljabar sederhana, persamaan integral dapat ditulis ulang sebagai
MembiarkanX menjadi variabel acak dengan kerapatan g . Menurut definisi,P{X≤z}=∫z0g(x)dx , jadi persamaan integral Anda setara dengan
"MembiarkanX menjadi variabel acak dengan kerapatan g . Temukan medianX . "
Untuk memperkirakan medianX , gunakan metode simulasi apa pun untuk menggambar sampel nilai X dan ambil sebagai estimasi Anda median sampel.
Salah satu kemungkinan adalah menggunakan algoritma Metropolis-Hastings untuk mendapatkan sampel poin dengan distribusi yang diinginkan. Karena ekspresi probabilitas penerimaan dalam algoritma Metropolis-Hastings, kita tidak perlu mengetahui nilai konstanta normalisasi∫10f(t)dt kepadatan g . Jadi, kita tidak perlu melakukan integrasi ini.
Kode di bawah ini menggunakan bentuk sederhana dari algoritma Metropolis-Hastings yang dikenal sebagai Indepence Sampler, yang menggunakan proposal yang distribusinya tidak bergantung pada nilai rantai saat ini. Saya telah menggunakan proposal seragam independen. Sebagai perbandingan, skrip menampilkan minimum Monte Carlo dan hasilnya ditemukan dengan optimasi standar. Titik sampel disimpan dalam vektor10000 titik yang membentuk apa yang disebut periode "terbakar" dalam simulasi.
chain
, tetapi kami membuang yang pertamaBerikut ini beberapa hasilnya:
Kode ini dimaksudkan hanya sebagai titik awal untuk apa yang benar-benar Anda butuhkan. Karenanya, gunakan dengan hati-hati.
sumber
Kualitas perkiraan integral, setidaknya dalam kasus sesederhana 1D, diberikan oleh (Teorema 2.10 dalam Niederreiter (1992) ):
Jadi jelas untuk meminimalkan kesalahan dalam perkiraan integral, setidaknya dalam RHS persamaan Anda, Anda perlu mengambilxn=(2n−1)/2N . Sekrup evaluasi acak, mereka berisiko memiliki kesenjangan acak pada fitur penting dari fungsi.
Kerugian besar dari pendekatan ini adalah Anda harus berkomitmen pada suatu nilaiN untuk menghasilkan urutan yang didistribusikan secara seragam ini. Jika Anda tidak puas dengan kualitas perkiraan yang diberikannya, yang dapat Anda lakukan hanyalah menggandakan nilaiN dan tekan semua titik tengah dari interval yang dibuat sebelumnya.
Jika Anda ingin memiliki solusi di mana Anda dapat meningkatkan jumlah poin secara lebih bertahap, Anda dapat terus membaca buku itu, dan belajar tentang urutan van der Corput dan invers radikal. Lihat urutan perbedaan rendah di Wikipedia, ini memberikan semua detail.
Perbarui: untuk dipecahkanz , tentukan jumlah parsial
sumber