Tidak benar bahwa MCMC yang memenuhi keseimbangan terinci selalu menghasilkan distribusi stasioner. Anda juga membutuhkan proses untuk menjadi ergodik . Mari kita lihat mengapa:
Anggap sebagai keadaan himpunan semua keadaan yang mungkin, dan identifikasikan dengan indeks . Dalam proses markov, distribusi berkembang sesuai dengani p t ( i )xsayahalt( i )
halt( i ) = ∑jΩj → ihalt - 1( j )
di mana adalah matriks yang menunjukkan probabilitas transisi ( Anda ).Ωj → iq( x | y)
Jadi, kita memilikinya
halt( i ) = ∑j( Ωj → i)thal0( j )
Fakta bahwa adalah probabilitas transisi menyiratkan bahwa nilai eigennya harus termasuk dalam interval [0,1].Ωj → i
Untuk memastikan bahwa setiap distribusi awal menyatu dengan asimptotik, Anda harus memastikan bahwahal0( j )
- 1 Hanya ada satu nilai eigen dengan nilai 1 dan memiliki vektor eigen bukan nol yang unik.Ω
Untuk memastikan bahwa adalah distribusi asimptotik, Anda perlu memastikannyaπ
- 2 Vektor eigen yang terkait dengan nilai eigen 1 adalah .π
Ergodisitas menyiratkan 1., keseimbangan terperinci menyiratkan 2., dan itulah sebabnya keduanya membentuk kondisi konvergensi asimptotik yang perlu dan cukup.
Mengapa saldo terperinci menyiratkan 2:
Mulai dari
p ( i ) Ωsaya j= Ωj ip ( j )
dan menjumlahkan di kedua sisi, kita dapatkanj
p ( i ) = ∑jΩj ip ( j )
karena , karena Anda selalu transit ke suatu tempat.∑jΩsaya j= 1
Persamaan di atas adalah definisi nilai eigen 1, (lebih mudah untuk melihat apakah Anda menulisnya dalam bentuk vektor :)
1. v = Ω ⋅ v
Saya pikir memang demikian, karena untuk MC yang tidak dapat direduksi jika keseimbangan terinci terpenuhi maka ia memiliki distribusi stasioner yang unik, tetapi agar independen dari distribusi awal juga harus aperiodik.
Dalam kasus MCMC kita mulai dari titik data dan kemudian mengusulkan titik baru. Kami mungkin atau mungkin tidak pindah ke titik yang diusulkan yaitu kami memiliki loop diri yang membuat MC aperiodik tak tereduksi.
Sekarang karena memuaskan DB itu juga memiliki keadaan berulang positif, yaitu waktu pengembalian rata-rata ke negara terbatas. Jadi rantai yang kami buat di MCMC tidak dapat direduksi, aperiodik, dan berulang positif, yang artinya merupakan rantai ergodik.
Kita tahu bahwa untuk rantai ergodik tak tereduksi distribusi stasioner ada yang unik dan independen dari distribusi awal.
sumber