Reprodutibilidade
Agora que você concluiu (e passou!) alguns diagnósticos de cadeias de Markov, está pronto para finalizar sua simulação com RJAGS. Para isso, a reprodutibilidade é essencial. Para obter resultados de simulação reproduzíveis, você precisa definir a semente do gerador de números aleatórios do RJAGS. Isso funciona de forma diferente do R base. Em vez de usar set.seed(), você vai indicar uma semente inicial usando inits = list(.RNG.name = "base::Wichmann-Hill", .RNG.seed = ___) quando você compilar seu modelo.
Este exercicio faz parte do curso
Modelagem Bayesiana com RJAGS
Instruções do exercicio
Execute o código fornecido algumas vezes. Observe que as estatísticas de
summary()mudam a cada execução.Para resultados reproduzíveis, forneça o gerador de números aleatórios
initsparajags.model(). Especifique a semente inicial como 1989.Execute o novo código algumas vezes. Observe que as estatísticas de
summary()NÃO mudam!
exercicio interativo prático
Tente este exercicio completando este código de exemplo.
# 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)