Dichteplots der Markow-Kette

Während ein Trace-Plot das zeitliche Verhalten einer Markow-Kette zeigt, veranschaulicht ein Dichteplot die resultierende Verteilung der Kettenwerte. Der Dichteplot liefert damit eine Approximation des posterioren Modells. Du erstellst und untersuchst Dichteplots der Markow-Kette für \(m\) unten. Das mcmc.list-Objekt sleep_sim und der Data Frame sleep_chains befinden sich in deinem Workspace:

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)

Diese Übung ist Teil des Kurses

<Kurs>Bayesianisches Modellieren mit RJAGS</Kurs>

Kurs ansehen

Übungsanweisungen

Wende plot() auf sleep_sim mit trace = FALSE an, um Dichteplots für die \(m\)- und \(s\)-Ketten zu erstellen.
Wende ggplot() auf sleep_chains an, um einen Dichteplot der \(m\)-Kette erneut zu erstellen.

Interaktive praktische Übung

Versuche dich an dieser Übung, indem du diesen Beispielcode vervollständigst.

# 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 = ___)) + 
    ___()

Code bearbeiten und ausführen