Een gezamenlijke verdeling berekenen
De delen van de code die we in de laatste video hebben ontwikkeld, staan hieronder. Ze definiëren in pars een rooster over de onderliggende verhoudingen van clicks (proportion_clicks) en mogelijke uitkomsten (n_visitors). Vervolgens voegen we de prior-kans van elke parametercombinatie toe en de likelihood dat elke proportion_clicks het bijbehorende n_visitors zou opleveren.
Deze oefening maakt deel uit van de cursus
Basisprincipes van Bayesiaanse data-analyse in R
Oefeninstructies
- Voeg de kolom
pars$probabilitytoe: de kans van elke combinatie vanproportion_clicksenn_visitors. Net als in de video bereken je dit door delikelihoodte vermenigvuldigen met deprior. - Zorg dat de kolom
pars$probabilityoptelt tot1.0door te normaliseren: deelpars$probabilitydoor de totale som vanpars$probability.
Interactieve oefening met praktijkervaring
Probeer deze oefening door deze voorbeeldcode aan te vullen.
n_ads_shown <- 100
proportion_clicks <- seq(0, 1, by = 0.01)
n_visitors <- seq(0, 100, by = 1)
pars <- expand.grid(proportion_clicks = proportion_clicks,
n_visitors = n_visitors)
pars$prior <- dunif(pars$proportion_clicks, min = 0, max = 0.2)
pars$likelihood <- dbinom(pars$n_visitors,
size = n_ads_shown, prob = pars$proportion_clicks)
# Add the column pars$probability and normalize it
pars$probability <- ___