1. Učit se
    2. /
  2. Kurzy
    3. /
  3. Bayesovské modelování s RJAGS
Connected

cvičení

Ukládání Markovových řetězců

Nechť \(m\) je průměrná změna reakční doby po 3 dnech spánkové deprivace. V předchozím cvičení jsi získal/a přibližný vzorek 10 000 hodnot z posteriorního modelu \(m\). Výsledný objekt mcmc.list jsi uložil/a jako sleep_sim, který je načten v tvém pracovním prostředí:

sleep_sim <- coda.samples(model = sleep_jags, variable.names = c("m", "s"), n.iter = 10000)

Vzorek hodnot \(m\) uložený v sleep_sim je závislý Markovův řetězec, jehož rozdělení konverguje k posterioru. Prozkoumáš obsah sleep_sim a pro lepší kontrolu nad analýzou ho uložíš do datového rámce.

Pokyny

100 XP

  • Podívej se na výstup funkce head() aplikované na seznam sleep_sim.

  • První položka seznamu sleep_sim obsahuje řetězce \(m\) a \(s\). Ulož je do datového rámce s názvem sleep_chains. Přidej proměnnou iter, která zaznamenává příslušné číslo iterace, 1:10000, pro každý prvek řetězce.

  • Prohlédni si prvních 6 řádků sleep_chains.