Kom igångKom igång gratis

Densitetsplotter för Markovkedjor

Medan en spårplot visar en Markovkedjas beteende över tid illustrerar en densitetsplot den slutliga fördelningen av kedjans värden. Densitetsplotten ger därmed en approximation av posteriorimodellen. Du ska konstruera och undersöka densitetsplotter för Markovkedjan \(m\) nedan. Objektet mcmc.list med namnet sleep_sim och dataramen sleep_chains finns i din arbetsmiljö:

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)

Den här övningen är en del av kursen

Bayesiansk modellering med RJAGS

Visa kurs

Övningsinstruktioner

  • Använd plot() på sleep_sim med trace = FALSE för att konstruera densitetsplotter för kedjorna \(m\) och \(s\).

  • Använd ggplot() på sleep_chains för att rekonstruera en densitetsplot för kedjan \(m\).

Interaktiv övning med praktiskt arbete

Testa den här övningen genom att slutföra den här exempelkoden.

# 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 = ___)) + 
    ___()
Redigera och kör kod