Porovnání Beta prior rozdělení

Rozdělení Beta(\(a\),\(b\)) je definováno na intervalu od 0 do 1, a proto přirozeně a flexibilně modeluje míru volební podpory \(p\). Tvarové parametry \(a\) a \(b\) lze ladit tak, aby vznikaly různé prior modely. Níže porovnáš původní prior Beta(45,55) se dvěma alternativami: Beta(1, 1) a Beta(100, 100). V prostředí máš k dispozici původních 10 000 vzorků prior_A vytvořených z Beta(45,55).

Toto cvičení je součástí kurzu

Bayesovské modelování s RJAGS

Zobrazit kurz

Pokyny k cvičení

Vyber 10 000 vzorků z prior rozdělení Beta(1,1). Výsledek ulož do prior_B.
Vyber 10 000 vzorků z prior rozdělení Beta(100,100). Výsledek ulož do prior_C.
Datový rámec prior_sim kombinuje vzorky prior_A, prior_B a prior_C s příslušným identifikátorem priors. Pro vykreslení grafu hustot pomocí ggplot() pro tato 3 oddělená prior rozdělení ve stejném grafu zadej fill = priors do aes().

Interaktivní cvičení na vyzkoušení si v praxi

Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.

# Sample 10000 draws from the Beta(1,1) prior
prior_B <- rbeta(n = ___, shape1 = ___, shape2 = ___)    

# Sample 10000 draws from the Beta(100,100) prior
prior_C <- rbeta(n = ___, shape1 = ___, shape2 = ___)

# Combine the results in a single data frame
prior_sim <- data.frame(samples = c(prior_A, prior_B, prior_C),
        priors = rep(c("A","B","C"), each = 10000))

# Plot the 3 priors
ggplot(___, aes(x = ___, fill = ___)) + 
    geom_density(alpha = 0.5)

Upravit a spustit kód