Mulai sekarangMulai gratis

Probabilitas posterior

Anda telah menggunakan keluaran RJAGS untuk menelusuri dan mengkuantifikasi tren dan ketidakpastian posterior \(b\). Anda juga dapat menggunakan keluaran RJAGS untuk menilai hipotesis tertentu. Misalnya: Berapa probabilitas posterior bahwa, rata-rata, berat meningkat lebih dari 1,1 kg untuk setiap kenaikan 1 cm pada tinggi badan? Artinya, berapa probabilitas posterior bahwa \(b > 1.1\)?

Anda akan mengaproksimasi probabilitas ini dengan proporsi nilai rantai Markov \(b\) yang melebihi 1.1. Data frame weight_chains dengan keluaran rantai Markov 100.000 iterasi ada di ruang kerja Anda.

Latihan ini merupakan bagian dari kursus

Pemodelan Bayesian dengan RJAGS

Lihat Kursus

Instruksi latihan

  • Bangun plot kerapatan dari nilai rantai Markov \(b\) dan gunakan geom_vline() untuk menimpa garis vertikal pada 1.1.
  • Gunakan table() untuk merangkum jumlah nilai rantai Markov \(b\) yang melebihi 1.1.
  • Gunakan mean() untuk menghitung proporsi nilai rantai Markov \(b\) yang melebihi 1.1.

Latihan interaktif langsung praktik

Cobalah latihan ini dengan melengkapi kode contoh ini.

# Mark 1.1 on a posterior density plot for b
ggplot(___, aes(x = ___)) + 
    geom_density() + 
    geom_vline(xintercept = ___, color = "red")

# Summarize the number of b chain values that exceed 1.1


# Calculate the proportion of b chain values that exceed 1.1
Edit dan Jalankan Kode