Reproducerbarhet
Nu när du har slutfört (och klarat!) några Markov-kedjediagnostiker är du redo att färdigställa din RJAGS-simulering. För detta är reproducerbarhet avgörande. För att få reproducerbara simuleringsresultat måste du ange ett startvärde för RJAGS slumptalsgenerator. Det fungerar annorlunda än i bas-R. I stället för att använda set.seed() anger du ett startvärde med inits = list(.RNG.name = "base::Wichmann-Hill", .RNG.seed = ___) när du kompilerar din modell.
Den här övningen är en del av kursen
Bayesiansk modellering med RJAGS
Övningsinstruktioner
Kör den angivna koden några gånger. Lägg märke till att
summary()-statistiken ändras varje gång.För reproducerbara resultat anger du slumptalsgeneratorns
initstilljags.model(). Ange startvärdet 1989.Kör den nya koden några gånger. Lägg märke till att
summary()-statistiken INTE ändras!
Interaktiv övning med praktiskt arbete
Testa den här övningen genom att slutföra den här exempelkoden.
# COMPILE the model
sleep_jags <- jags.model(textConnection(sleep_model), data = list(Y = sleep_study$diff_3))
# SIMULATE the posterior
sleep_sim <- coda.samples(model = sleep_jags, variable.names = c("m", "s"), n.iter = 10000)
# Summarize the m and s chains of sleep_sim
summary(sleep_sim)