Está a punto de acabar un partido de fútbol entre la Real Sociedad y el Barcelona y me temo que esta jornada vuelven a ganar los equipos que más dinero ganan por los derechos televisivos de la liga española. El Real Madrid está a 7 puntos del Barcelona. Está claro que el Barcelona ganará al Real en su campo. Si asumimos que los dos equipos tienen un 80% de posibilidades de ganar el partido, un 10% de empatarlo y un 10% de perderlo, ¿cuál es la probabilidad de que gane alguno de los dos equipos?
Pido disculpas por la “poca profesionalidad” del código, pero últimamente pico poco código en R; de todas formas es muy sencillo de entender y no voy a entrar en profundidad:
library(memisc)
final_madrid <- data.frame()
for (i in 1:1000) {
madrid = runif(16)
madrid <- recode(madrid,
0 <- range(0, 0.1),
1 <- range(0.1, 0.2),
3 <- range(0.2, 1))
final_madrid = rbind(final_madrid, sum(madrid, 52))
}
final_barcelona <- data.frame()
for (i in 1:1000) {
barcelona = runif(16)
barcelona <- recode(barcelona,
0 <- range(0, 0.1),
1 <- range(0.1, 0.2),
3 <- range(0.2, 1))
final_barcelona = rbind(final_barcelona, sum(barcelona, 45, 3))
}
resultado_liga = cbind(final_madrid, final_barcelona)
names(resultado_liga) = c("madrid", "barcelona")
resultado_liga$campeon = ifelse(resultado_liga$madrid > resultado_liga$barcelona, "Madrid", "Barcelona")
table(resultado_liga$campeon)
He lanzado 1000 simulaciones y he obtenido unas probabilidades de 26% a favor del Barcelona y un 74% a favor del Madrid. ¿En las casas de apuestas está el resultado final 4 a 1? A ver qué pasa; en unos meses volveré sobre esta entrada y espero decir que al final ganó la liga otro equipo, pero eso sí que es improbable. Saludos.