Bisakah para ahli Monte Carlo kami menjelaskan ekspektasi "tak terduga" di akhir jawaban ini ? Ringkasan ex post facto dari pertanyaan / jawaban lain: jika adalah variabel acak IID dan harapan ada, maka argumen simetri sederhana menunjukkan bahwa , tetapi percobaan Monte Carlo dengan tampaknya...
Simulasi ilmu semi-komputer terkait masalah di sini. Saya punya distribusi di mana P (x) =(eb−1)eb(n−x)ebn+b−1(eb−1)eb(n−x)ebn+b−1\frac{(e^b-1) e^{b (n-x)}}{e^{b n+b}-1} untuk beberapa konstanta b dan n, dan x adalah bilangan bulat sehingga .0≤x≤n0≤x≤n0\leq x \leq n Sekarang, saya perlu sampel...
Semangat pertanyaan ini berasal dari "Monte Carlo Biasa", yang juga dikenal sebagai "Monte Carlo kuno yang baik" Misalkan saya punya variabel acak XXX, dengan μ : = E[ X]σ2: = Va r [ X]μ: =E[X]σ2: =VSebuahr[X]\mu := E[X]\\ \sigma^2:=Var[X] Keduanya adalah nilai yang tidak diketahui, karena...
Saya telah mencoba untuk merasakan masalah yang berbeda dalam pengaturan sering di mana MCMC digunakan. Saya tahu bahwa MCMC (atau Monte Carlo) digunakan dalam pemasangan GLMM dan mungkin juga algoritma Monte Carlo EM. Apakah ada masalah yang lebih sering terjadi di mana MCMC
Saya membaca pertanyaan ini , dan berpikir untuk mensimulasikan jumlah yang diperlukan. Masalahnya adalah sebagai berikut: JikaAAA dan BBB apakah standar normal, apa E(A2|A+B)E(A2|A+B)E(A^2|A+B)? Jadi saya ingin mensimulasikanE(A2|A+B)E(A2|A+B)E(A^2|A+B). (untuk nilai yang dipilih...
Saat ini saya menggunakan Latin Hypercube Sampling (LHS) untuk menghasilkan angka acak seragam spasi baik untuk prosedur Monte Carlo. Walaupun reduksi varians yang saya dapatkan dari LHS sangat bagus untuk 1 dimensi, sepertinya tidak efektif dalam 2 dimensi atau lebih. Melihat bagaimana LHS adalah...
Saya memiliki dua kelompok yang terdiri dari 10 peserta yang dinilai tiga kali selama percobaan. Untuk menguji perbedaan antara kelompok dan di tiga penilaian, saya menjalankan ANOVA desain campuran 2x3 dengan group(kontrol, eksperimental), time(pertama, kedua, tiga), dan group x time. Keduanya...
Saya mencoba menerapkan HMC dengan matriks massa non-diagonal, tetapi saya tersandung oleh beberapa terminologi. Menurut BDA3 dan ulasan Neal, istilah energi kinetik (yang saya kira selalu digunakan karena kenyamanan) adalah K( p ) =halTM.- 1hal2.K(p)=pTM−1p2. K(p) = \frac{p^T M^{-1} p}{2}...
Saya telah melihat pertanyaan augmentasi data MCMC; bentuk umum dari pertanyaan adalah sebagai berikut: Misalkan data yang dikumpulkan pada suatu proses menyarankan dan prior untuk parameter rate disarankan sebagai . Data dicatat dan disajikan dalam bentuk yang khas (yaitu jumlah kemunculan dari...
Dalam rantai Metropolis-Hastings Markov Monte Carlo, distribusi proposal dapat berupa apa saja termasuk Gaussian (menurut Wikipedia). T: Apa motivasi untuk menggunakan selain Gaussian? Gaussian bekerja, mudah untuk mengevaluasi, cepat dan semua orang memahaminya. Mengapa saya mempertimbangkan hal...
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. fffz∈[0,1]z∈[0,1]z \in...