Modelo de Poisson para la Expectativa de Goles: Análisis Predictivo en Fútbol

1. Fundamentos Matemáticos del Proceso de Poisson en el Fútbol

El fútbol profesional presenta una particularidad estadística fundamental que lo diferencia de la inmensa mayoría de disciplinas deportivas: los goles son sucesos cuantitativamente escasos, discretos e integrados a lo largo de un intervalo continuo fijo de 90 minutos reglamentarios. En un partido estándar de las principales competiciones europeas o latinoamericanas se anotan entre 2,5 y 2,8 goles combinando ambos conjuntos, lo que equivale a una frecuencia media de apenas un gol cada 32 a 36 minutos de juego efectivo.

Debido a que las anotaciones se producen de manera discontinua y relativamente independiente en una línea temporal continua, la dinámica estocástica del marcador en un encuentro de fútbol se ajusta con rigor a los axiomas de un Proceso de Poisson. En la teoría formal de probabilidades, un proceso de Poisson se cimenta sobre tres propiedades esenciales:

1. Discretización y No Negatividad: La variable aleatoria $X$, que representa la cantidad de goles anotados por un equipo, adopta exclusivamente valores enteros no negativos pertenecientes al conjunto: $k\in \{0,1,2,3,\dots \}$. Es físicamente imposible que un club marque una fracción de gol o registre un saldo negativo.
2. Independencia Temporal (Ausencia de Memoria): La probabilidad de que se marque un gol en un lapso infinitesimal $[t,t+\mathrm{\Delta }t]$ es estadísticamente independiente de los eventos acaecidos con anterioridad al instante $t$. El hecho de haber anotado en el minuto 15 no altera intrínsecamente la probabilidad objetiva de anotar en el minuto 16, bajo una calidad táctica constante.
3. Proporcionalidad y No Simultaneidad: La probabilidad de que acontezca exactamente un gol durante un breve incremento $\mathrm{\Delta }t$ es directamente proporcional a la amplitud del intervalo ($\lambda \mathrm{\Delta }t+o(\mathrm{\Delta }t)$), mientras que la probabilidad de que se produzcan dos o más goles en el mismo instante infinitesimal es asintóticamente despreciable ($o(\mathrm{\Delta }t)$).

Matemáticamente, la distribución de Poisson surge como el límite analítico de una distribución Binomial $B(n,p)$ cuando el número de ensayos $n$ tiende a infinito y la probabilidad de éxito $p$ se aproxima a cero, de modo que el producto $np=\lambda$ permanece constante. La Función de Masa de Probabilidad (FMP) que determina la probabilidad de que un equipo anote con exactitud $k$ goles en función del parámetro de intensidad media $\lambda >0$ se define como:

Donde $e\approx 2.71828$ denota la constante de Euler y $k!$ representa el factorial de $k$ (con $0!=1$). Un atributo sobresaliente de la distribución de Poisson es la equidispersión: la esperanza matemática y la varianza son idénticas al parámetro de intensidad $\lambda$:

Poisson tiene igual media y varianza por definición. La dispersión muestral cerca de uno puede ser diagnóstica, pero no prueba distribución completa, independencia ni calibración. Defina población, ventana e incertidumbre antes de interpretar la razón; no se establece rango universal 0.98–1.04.

2. Descomposición de la Expectativa de Goles: Ataque, Defensa y Factor Local

Para emplear el modelo de Poisson con fines analíticos y predictivos, el analista no puede asignar valores arbitrarios a $\lambda$. En su lugar, las tasas de goles esperados para el conjunto local (${\lambda }_{local}$) y el conjunto visitante (${\mu }_{visitante}$) se deducen a partir de la descomposición de datos históricos consolidados:

Donde cada término representa una dimensión de rendimiento deportivo:

• ${\mu }_{liga}$ (Media Histórica de la Liga): Promedio de goles anotados por equipo por partido en el campeonato global (generalmente entre 1.30 y 1.45 goles por equipo).
• ${\alpha }_i$ (Potencia Ofensiva del Equipo $i$): Mide la capacidad relativa de ataque. Un registro de $\alpha =1.25$ denota que el equipo genera un 25% más de goles que la media del torneo frente a un rival promedio.
• ${\beta }_j$ (Vulnerabilidad Defensiva del Equipo $j$): Evalúa la fragilidad defensiva relativa. Un coeficiente de $\beta =0.80$ indica que la plantilla encaja un 20% menos de goles que el estándar de la liga (un sistema defensivo de élite).
• $\gamma$ (Ventaja de Localía): Multiplicador que cuantifica el impacto positivo de disputar el encuentro en el estadio propio, propiciado por la presión ambiental del público, la familiaridad con el terreno de juego y la ausencia de desgaste por desplazamientos. En el fútbol de alta competición, $\gamma$ oscila entre $1.15$ y $1.25$.

Para asegurar que los parámetros del modelo se encuentren matemáticamente definidos en un torneo compuesto por $N$ escuadras, se fijan condiciones de normalización:

El vector paramétrico $\overrightarrow{\theta }=(\overrightarrow{\alpha },\overrightarrow{\beta },\gamma )$ se estima rigurosamente aplicando el método de Máxima Verosimilitud (MLE) sobre un histórico de $M$ partidos:

3. Distribución Bivariada de Poisson y la Matriz de Marcadores Correctos

Una vez determinados los parámetros de intensidad media $\lambda$ (local) y $\mu$ (visitante), la probabilidad de cualquier resultado exacto —sea un triunfo por 2-1, una igualdad 0-0 o una sorpresa visitante— se deriva a través de la Distribución Bivariada de Poisson.

Asumiendo de forma preliminar la independencia entre ambas fuerzas, la probabilidad conjunta de que el equipo local consiga $x$ goles y el visitante obtenga $y$ goles corresponde al producto de sus probabilidades marginales:

A modo de ejemplo ilustrativo, consideremos un duelo donde el modelo proyecta $\lambda =1.68$ goles esperados para el equipo anfitrión y $\mu =1.15$ goles para la escuadra foránea. Estructuramos la matriz completa de probabilidades para los marcadores comprendidos entre 0 y 5 goles por bando:

La cuadrícula muestra marcadores de 0 a 5 como ilustración. Las probabilidades de los mercados usan el modelo completo, calculado hasta 45 goles con masa omitida despreciable a estas intensidades. Cada celda visible se redondea por separado; sumar la cuadrícula redondeada no reproduce todas las probabilidades del mercado completo.

Cada celda refleja la probabilidad modelada redondeada de un marcador concreto. El resultado individual más probable es el empate 1-1 (11.40%), seguido de cerca por la victoria local 1-0 (9.91%) y el 2-1 (9.58%).

4. Fijación de Cuotas en Mercados Derivados

La ventaja operativa suprema de este modelo estriba en que la matriz bivariada permite calcular de forma unificada las cuotas justas de los mercados principales.

Mercado de Resultado Final (1X2)

• Victoria Local ($P_H$): Suma de casillas situadas bajo la diagonal ($x>y$):
  $$P_H=\sum _{x>y}P(X=x,Y=y)=49.754448\%⟹O_H=\frac{1}{0.49754448}=2.0099$$
• Empate ($P_D$): Suma de las casillas de la diagonal ($x=y$):
  $$P_D=\sum _{x=y}P(X=x,Y=y)=24.145909\%⟹O_D=\frac{1}{0.24145909}=4.1415$$
• Victoria Visitante ($P_A$): Suma de casillas sobre la diagonal ($x<y$):
  $$P_A=\sum _{x<y}P(X=x,Y=y)=26.099644\%⟹O_A=\frac{1}{0.26099644}=3.8315$$

Mercado de Más / Menos de 2.5 Goles (Over / Under)

Sumando los marcadores correspondientes a dos o menos goles combinados:

Ambos Equipos Anotan (BTTS)

5. La Corrección Bivariada de Dixon y Coles (1997)

A pesar de su claridad formal, el modelo simple adolece de un fallo econométrico sistemático: la generación de goles entre ambos contrincantes no es enteramente independiente.

Para corregir este desajuste, Mark J. Dixon y Stuart G. Coles propusieron un factor de acoplamiento $\rho$ y un modificador $\tau (x,y,\lambda ,\mu ,\rho )$ que calibra de forma selectiva los cuatro marcadores de baja anotación:

Donde la función de ajuste $\tau$ se define como:

Rho es parámetro de corrección de marcadores bajos, no correlación de Pearson ni covarianza. Estímelo con factores positivos. Para tasas positivas el dominio estricto es $-\min (1/\lambda ,1/\mu )<\rho <\min (\mathrm{1,1}/(\lambda \mu ))$. Rho negativo factible aumenta 0:0 y1:1; positivo invierte la dirección. En soporte completo conserva masa y marginales; una cuadrícula finita tiene cola.

Ponderación Temporal Exponencial (Time-Decay)

Dixon y Coles integraron asimismo una función de ponderación temporal decreciente para otorgar mayor peso a los enfrentamientos recientes en comparación con partidos disputados meses atrás:

Calibrada usualmente con una vida media de 180 a 240 días, asegurando una rápida adaptación a los cambios en el estado de forma de los planteles.

6. Modernización del Modelo con Expected Goals (xG)

En el panorama cuantitativo actual, los sindicatos de apuestas ya no alimentan sus modelos exclusivamente con el registro de goles históricos. El motivo reside en la reducida muestra de partidos.

A lo largo de 38 jornadas ligueras se anotan muy pocos goles, interviniendo factores puramente fortuitos (rebotes involuntarios, decisiones arbitrales discutibles o errores atípicos de guardametas). Utilizar goles brutos introduce ruido estadístico en las proyecciones.

Las entradas xG estiman probabilidades de convertir remates. Pueden aportar información adicional a los goles, pero dependen de variables, proveedor, entrenamiento y calibración. Compare modelos de goles y xG en los mismos partidos fuera de muestra, con métrica e incertidumbre definidas. No se establece reducción fija del 35% en varianza.

7. Caso Práctico Detallado: Real Madrid vs Atlético de Madrid

Las fuerzas de los equipos, la media de la liga y la cuota siguiente son datos ilustrativos del cálculo. No constituyen una medición fechada de los equipos citados. Las probabilidades conservan toda la precisión de esos datos.

Analicemos un caso práctico completo: Real Madrid (Local) vs Atlético de Madrid (Visitante).

• Media liguera: ${\mu }_{liga}=1.32$ | Ventaja de Localía: $\gamma =1.18$
• Ataque Real Madrid: ${\alpha }_H=1.40$ | Defensa Real Madrid: ${\beta }_H=0.75$
• Ataque Atlético: ${\alpha }_A=1.15$ | Defensa Atlético: ${\beta }_A=0.80$

Determinación de las tasas de goles:

Calculando probabilidades bajo Dixon-Coles ($\rho =-0.10$):

• Victoria local: $P_H\approx 50.450218\%,O_H\approx 1.9822$
• Empate: $P_D\approx 25.836350\%,O_D\approx 3.8705$
• Victoria visitante: $P_A\approx 23.713432\%,O_A\approx 4.2170$
• Más de 2.5: $P_{Over}\approx 55.009591\%,O_{Over}\approx 1.8179$
• Menos de 2.5: $P_{Under}\approx 44.990409\%,O_{Under}\approx 2.2227$
• Ambos marcan: $P_{BTTS}\approx 57.205214\%,O_{BTTS}\approx 1.7481$

Con la cuota ilustrativa de 2.1 para la victoria local, este modelo da un EV estimado de aproximadamente 5.9455%. El resultado depende de las probabilidades y de reglas de liquidación compatibles; evaluar el modelo con datos fuera de la estimación es un paso distinto.

8. Consideraciones Operativas y Factores Críticos

El analista cuantitativo debe integrar consideraciones dinámicas que trascienden el marco estático de Poisson:

1. Dinámica del Marcador en Vivo (Game State): Un equipo con ventaja tiende a proteger el resultado replegando sus líneas, mientras que el rezagado asume mayores riesgos ofensivos, alterando los parámetros originales.
2. Las expulsiones cambian el estado. El efecto depende de tiempo restante, marcador, fuerza y modelo ajustado. Estime y valide el ajuste por estado, no multiplicadores universales del 50% y30%.
3. Rotaciones y Bajas Clave: La ausencia de futbolistas determinantes exige ajustar manualmente los coeficientes $\alpha$ y $\beta$.

9. Conclusión y Herramientas Prácticas

El Modelo de Poisson enriquecido con la corrección de Dixon-Coles y la integración de métricas xG constituye la base analítica de referencia en el análisis cuantitativo de fútbol.

Comprueba estas formulaciones en tiempo real con nuestra Calculadora de Modelo de Poisson gratuita.