Beräkna en gemensam fördelning
Koden vi utvecklade i det senaste videoavsnittet visas här. Den definierar ett rutnät över de underliggande andelarna av klick (proportion_clicks) och möjliga utfall (n_visitors) i pars. Till detta läggs prior-sannolikheten för varje parameterkombination samt sannolikheten (likelihood) att varje proportion_clicks skulle generera motsvarande n_visitors.
Den här övningen är en del av kursen
Grunderna i bayesiansk dataanalys i R
Övningsinstruktioner
- Lägg till kolumnen
pars$probability: sannolikheten för varje kombination avproportion_clicksochn_visitors. Precis som i videon ska denna beräknas genom att multipliceralikelihoodmedprior. - Se till att kolumnen
pars$probabilitysummerar till1.0genom att normalisera den, det vill säga dividerapars$probabilitymed den totala summan avpars$probability.
Interaktiv övning med praktiskt arbete
Testa den här övningen genom att slutföra den här exempelkoden.
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 <- ___