Gráficos de densidade da cadeia de Markov
Enquanto um gráfico de trajetória (trace plot) captura o comportamento longitudinal de uma cadeia de Markov, um gráfico de densidade mostra a distribuição final dos valores da cadeia. Por sua vez, o gráfico de densidade fornece uma aproximação do modelo a posteriori. Você vai construir e analisar gráficos de densidade da cadeia \(m\) abaixo. Os objetos mcmc.list sleep_sim e o data frame sleep_chains estão no seu ambiente de trabalho:
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)
Este exercicio faz parte do curso
Modelagem Bayesiana com RJAGS
Instruções do exercicio
Aplique
plot()asleep_simcomtrace = FALSEpara construir gráficos de densidade para as cadeias \(m\) e \(s\).Aplique
ggplot()asleep_chainspara reconstruir um gráfico de densidade da cadeia \(m\).
exercicio interativo prático
Tente este exercicio completando este código de exemplo.
# 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 = ___)) +
___()