Aktualizacja rozkładu posterior
Rozkład posterior twojego poparcia wyborczego \(p\) wynika zarówno z rozkładu prior, jak i z danych ankietowych \(X\). Uruchom skrypt po prawej stronie, aby przypomnieć sobie rozkład posterior, który powstał na podstawie pierwotnego rozkładu prior (Beta(45, 55)) i pierwotnych danych ankietowych (\(X = 6\) spośród \(n = 10\) ankietowanych wyborców popiera cię). Zdefiniowany model vote_model jest dostępny w twoim środowisku pracy.
W 3-etapowym ćwiczeniu sprawdzisz, jak użycie innego rozkładu prior lub zaobserwowanie nowych danych (albo kombinacji obu!) może wpłynąć na rozkład posterior.
To ćwiczenie jest częścią kursu
Modelowanie bayesowskie z RJAGS
Interaktywne ćwiczenie praktyczne
Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.
# COMPILE the model
vote_jags <- jags.model(textConnection(vote_model),
data = list(a = 45, b = 55, X = 6, n = 10),
inits = list(.RNG.name = "base::Wichmann-Hill", .RNG.seed = 100))
# SIMULATE the posterior
vote_sim <- coda.samples(model = vote_jags, variable.names = c("p"), n.iter = 10000)
# PLOT the posterior
plot(vote_sim, trace = FALSE, xlim = c(0,1), ylim = c(0,18))