Mulai sekarangMulai gratis

Plot kerapatan rantai Markov

Jika plot jejak menangkap perilaku longitudinal suatu rantai Markov, plot kerapatan menggambarkan sebaran akhir nilai rantai. Pada gilirannya, plot kerapatan memberikan aproksimasi dari model posterior. Anda akan membangun dan menelaah plot kerapatan untuk rantai Markov \(m\) di bawah ini. Objek mcmc.list sleep_sim dan data frame sleep_chains sudah tersedia di workspace Anda:

sleep_sim <- coda.samples(model = sleep_jags, variable.names = c("m", "s"), n.iter = 10000)
sleep_chains <- data.frame(sleep_sim[[1]], iter = 1:10000)

Latihan ini merupakan bagian dari kursus

Pemodelan Bayesian dengan RJAGS

Lihat Kursus

Instruksi latihan

  • Terapkan plot() pada sleep_sim dengan trace = FALSE untuk membangun plot kerapatan bagi rantai \(m\) dan \(s\).

  • Terapkan ggplot() pada sleep_chains untuk membangun kembali plot kerapatan dari rantai \(m\).

Latihan interaktif langsung praktik

Cobalah latihan ini dengan melengkapi kode contoh ini.

# Use plot() to construct density plots of the m and s chains


# Use ggplot() to construct a density plot of the m chain
ggplot(___, aes(x = ___)) + 
    ___()
Edit dan Jalankan Kode