Atualizando com três moedas
Suponha que, em vez de a moeda ser apenas justa ou viciada, haja três possibilidades: a moeda é justa (50% de caras), baixa (25% de caras) ou alta (75% de caras). Há 80% de chance de ela ser justa, 10% de chance de ser viciada para baixo e 10% de chance de ser viciada para cima.
Você observa que 14 de 20 lançamentos deram cara. Qual é a probabilidade de a moeda ser justa?
Este exercício faz parte do curso
Fundamentos de Probabilidade em R
Instruções do exercício
- Use a função
rbinom()para simular 80.000 amostras da moeda justa, 10.000 amostras da moeda alta e 10.000 amostras da moeda baixa, com cada amostra contendo 20 lançamentos. Salve-as comoflips_fair,flips_higheflips_low, respectivamente. - Para cada um desses tipos, calcule quantas moedas resultaram em 14. Salve-as como
fair_14,high_14elow_14, respectivamente. - Encontre a probabilidade a posteriori de a moeda ser justa, dividindo o número de moedas justas que resultaram em 14 pelo total de moedas que resultaram em 14.
Exercício interativo prático
Experimente este exercício completando este código de exemplo.
# Simulate 80,000 draws from fair coin, 10,000 from each of high and low coins
flips_fair <-
flips_high <-
flips_low <-
# Compute the number of coins that resulted in 14 heads from each of these piles
fair_14 <-
high_14 <-
low_14 <-
# Compute the posterior probability that the coin was fair