Altere o modelo para usar um prior informativo
O código exibido é o modelo antigo que você desenvolveu do zero no capítulo 2.
Este exercicio faz parte do curso
Fundamentos da Análise de Dados Bayesiana em R
Instruções do exercicio
- Altere este modelo para usar o novo prior informativo para
proportion_clicksque você acabou de selecionar:
rbeta(n_draws, shape1 = 5, shape2 = 95)
exercicio interativo prático
Tente este exercicio completando este código de exemplo.
n_draws <- 100000
n_ads_shown <- 100
# Change the prior on proportion_clicks
proportion_clicks <-
runif(n_draws, min = 0.0, max = 0.2)
n_visitors <-
rbinom(n_draws, size = n_ads_shown,
prob = proportion_clicks)
prior <-
data.frame(proportion_clicks, n_visitors)
posterior <-
prior[prior$n_visitors == 13, ]
# This plots the prior and the posterior in the same plot
par(mfcol = c(2, 1))
hist(prior$proportion_clicks,
xlim = c(0, 0.25))
hist(posterior$proportion_clicks,
xlim = c(0, 0.25))