Perebutan dan korelasi dalam urutan perbedaan rendah (Halton / Sobol)

14

Saat ini saya sedang mengerjakan proyek di mana saya menghasilkan nilai acak menggunakan set titik diskrepansi / kuasi-acak yang rendah , seperti set titik Halton dan Sobol. Ini pada dasarnya adalah vektor dimensi yang meniru variabel seragam d- dimensi (0,1), tetapi memiliki penyebaran yang lebih baik. Secara teori, mereka seharusnya membantu mengurangi varian perkiraan saya di bagian lain dari proyek.dd

Sayangnya, saya telah mengalami masalah dalam bekerja dengan mereka dan banyak literatur tentang mereka padat. Karena itu saya berharap untuk mendapatkan wawasan dari seseorang yang memiliki pengalaman dengan mereka, atau setidaknya mencari cara untuk menilai secara empiris apa yang sedang terjadi:

Jika Anda pernah bekerja dengan mereka:

  • Apa sebenarnya pengacakan? Dan apa efeknya pada aliran poin yang dihasilkan? Secara khusus, apakah ada efek ketika dimensi poin yang dihasilkan meningkat?

  • Mengapa jika saya menghasilkan dua aliran poin Sobol dengan pengacakan MatousekAffineOwen, saya mendapatkan dua aliran poin yang berbeda. Mengapa ini tidak terjadi ketika saya menggunakan reverse-radix scrambling dengan poin Halton? Apakah ada metode pengacakan lainnya yang ada untuk set poin ini - dan jika demikian, apakah ada implementasi MATLAB dari mereka?

Jika Anda belum bekerja dengan mereka:

  • Katakanlah saya memiliki urutan S 1 , S 2 , ... , S n dari angka yang seharusnya acak, jenis statistik apa yang harus saya gunakan untuk menunjukkan bahwa mereka tidak berkorelasi satu sama lain? Dan apa jumlah n akan saya butuhkan untuk membuktikan bahwa hasil saya adalah signifikan secara statistik? Juga, bagaimana saya bisa melakukan hal yang sama jika saya memiliki n sekuens S 1 , S 2 , ... , S n dari d- dimensional dimensional [ 0 , 1 ] vektor?nS1,S2,,SnnnS1,S2,,Snd[0,1]

Pertanyaan Tindak Lanjut tentang Jawaban Kardinal

  1. Secara teoritis, dapatkah kita memasangkan metode pengacakan dengan urutan perbedaan rendah? MATLAB hanya memungkinkan saya untuk menerapkan scrambling reverse-radix pada urutan Halton, dan saya bertanya-tanya apakah itu hanya masalah implementasi atau masalah kompabilitas.

  2. Saya mencari cara yang akan memungkinkan saya untuk menghasilkan dua (t, m, s) jaring yang tidak berkorelasi satu sama lain. Apakah MatouseAffineOwen mengizinkan saya melakukan ini? Bagaimana jika saya menggunakan algoritma pengacakan deterministik dan hanya memutuskan untuk memilih setiap nilai 'kth' di mana k adalah prima?

Berk U.
sumber
apa yang Anda maksud dengan dua jaring yang tidak berkorelasi? Secara khusus, apa artinya ini ketika Anda mengatakan "kami [algoritma ing deterministic]"? Banyak algoritma pengacakan dapat diterapkan ke jaring ( t , m , s ) yang berubah-ubah ; Sejujurnya saya tidak tahu apakah semua skema bisa. Saya akan membayangkan jawabannya mungkin "tidak". (Yaitu, orang dapat membuat scrambling yang cukup khusus sehingga mempertahankan properti penutupan untuk urutan tertentu , tetapi tidak secara umum. Saya tidak tahu tentang itu begitu saja.)(t,m,s)(t,m,s)
kardinal
(t,m,s)PQ{pi}1100{qi}1100{pi}1100{qi}1100{pi}1200{qi}1200
(t,m,s)

Jawaban:

10

(t,m,s)bb=2b

d=2

masukkan deskripsi gambar di sini

bn

(t,m,s)(t,m,s)(t,m,s) bersih di tempat pertama, ini sangat diinginkan.

Mengenai jenis pengacakan, perebutan reverse-radix adalah pengacakan deterministik . Algoritma pengacakan Matousek adalah pengacakan acak , dilakukan, sekali lagi, dengan cara untuk mempertahankan properti penutupan. Jika Anda mengatur seed acak sebelum Anda melakukan panggilan ke fungsi scrambling, Anda harus selalu mendapatkan net yang sama.

Anda mungkin juga tertarik dengan proyek MinT .

kardinal
sumber
Terima kasih banyak atas ini. Saya punya beberapa pertanyaan lanjutan jika Anda tidak keberatan. Karena kotak komentar tidak memungkinkan saya untuk mencantumkannya dengan jelas, saya memasukkannya dalam pos.
Berk U.