Análisis completo de Lotería Nacional

12 pruebas estadísticas con 6.131 sorteos y 18.393 números. 1 de 12 dieron señal (p < 0,01). Con 12 pruebas, que alguna salga por casualidad pasa en ~11 % de los casos.

Pruebas estadísticas de Lotería Nacional
PruebaEstadísticoValor pResultado
Frecuencia de los 100 númerosχ² = 69,1 · gl 990,99Normal
Frecuencia en 1raχ² = 85,8 · gl 990,83Normal
Frecuencia en 2daχ² = 88,6 · gl 990,76Normal
Frecuencia en 3raχ² = 82,9 · gl 990,88Normal
Bombo de las decenasχ² = 5,7 · gl 90,77Normal
Bombo de las unidadesχ² = 7,1 · gl 90,63Normal
Decenas y unidades, independientesχ² = 56,5 · gl 810,98Normal
Memoria (¿importa cuánto lleva sin salir?)χ² = 6,8 · gl 70,45Normal
Repetición del sorteo anteriorz = -0,110,91Normal
Tendencia de los últimos 90 díasΣz² = 94,7 · gl 990,60Normal
Pares que salen juntos (palé)χ² = 4.889 · gl 4.9490,73Normal
Rachas de pares e impares en 1raz = 3,230,001Señal
¿Qué tan buenas son estas probabilidades?

Pusimos a competir 6 modelos en los últimos 1.500 sorteos. Antes de cada sorteo, cada modelo dio la probabilidad de los 100 números usando solo lo anterior, y se calificó con lo que pasó. Ningún modelo mejoró al azar: las probabilidades de todos los números son, en la práctica, iguales.

Competencia de modelos
ModeloPérdida logarítmicaMejora (log)Mejora (Brier)
Azar puro (2,97 % cada número)0,13347——
Bayes con todo el historial0,13347+0,000 %+0,000 %
Bayes por posición (el de los números del día)0,13347+0,000 %+0,000 %
Bayes con lo reciente (vida media 30 días)0,13347+0,000 %+0,000 %
Bayes con lo reciente (vida media 90 días)0,13347+0,000 %+0,000 %
Bayes con lo reciente (vida media 1 año)0,13347+0,000 %+0,000 %
Calibración del modelo de los números del díaCuando el modelo dijo una probabilidad (eje de abajo), ¿con qué frecuencia salió de verdad el número (eje de la izquierda)? Si los puntos caen sobre la diagonal, las probabilidades son honestas.
2,8 %2,8 %3,0 %3,0 %3,1 %3,1 %

L = −(1⁄n) Σ [y ln p + (1 − y) ln(1 − p)]  ·  B = (1⁄n) Σ (p − y)²  ·  mejora = 1 − Lmodelo ⁄ Lazar

¿Qué tan lejos del azar está?Curva de lo que da el azar puro (χ² con 99 grados de libertad) y lo observado en los 100 números. Si la línea roja cae en la zona normal, no hay sesgo. Valor p = 0,99. Entropía de Shannon: 6,641 de 6,644 bits (99,96 % del máximo).
observado: 69,199 (lo típico)43169

χ² = Σ (Oᵢ − Eᵢ)² ⁄ Eᵢ  ·  H = −Σ pᵢ log₂ pᵢ

¿Importa cuánto lleva sin salir?Probabilidad de que un número salga según los sorteos que lleva sin salir. Si los "atrasados" fueran reales, los puntos subirían hacia la derecha. Línea: 2,97 % (azar). Valor p = 0,45.
3,0 %012–34–78–1516–3132–6364–127128+

h(g) = salidas con g sorteos sin salir ⁄ veces que estuvo g sorteos sin salir  ·  azar: P(G = g) = (1 − p)g−1 p

Los dos bombosPrimer dígito (decenas, p = 0,77) y último dígito (unidades, p = 0,63). Independencia entre ellos: p = 0,98.

Decenas (primer dígito)

Por azar, cada decena sale unas 1839 veces (línea punteada). Dentro de la franja gris es normal: todas están dentro.

esperado: 1839Del 00 al 091842Del 10 al 191847Del 20 al 291845Del 30 al 391830Del 40 al 491868Del 50 al 591759Del 60 al 691871Del 70 al 791811Del 80 al 891872Del 90 al 991848

Unidades (último dígito)

Por azar, cada terminación sale unas 1839 veces (línea punteada). Dentro de la franja gris es normal: todas están dentro.

esperado: 1839Termina en 01861Termina en 11896Termina en 21848Termina en 31793Termina en 41853Termina en 51865Termina en 61772Termina en 71848Termina en 81802Termina en 91855

χ²ind = Σ (Odu − RdCu⁄N)² ⁄ (RdCu⁄N), gl = 81

¿Se repite el sorteo anterior?En cuántos sorteos salió algún número del sorteo anterior del mismo tipo, frente a lo que da el azar. Valor p = 0,91.
8,6 %0ObservadoAzar

E = Σ [1 − (1 − |Sₜ₋₁|⁄100)³]  ·  z = (O − E) ⁄ √Σ q(1 − q)

¿Hay números "en racha"?Para cada número, cuánto cambió su frecuencia en los últimos 90 días frente al resto (en desviaciones estándar). Si hubiera rachas reales, las barras se saldrían de la curva normal. Valor p = 0,60.
-3-2-10+1+2+3

z = (r₉₀ − rresto) ⁄ √(p̄(1 − p̄)(1⁄n₉₀ + 1⁄nresto))  ·  Σz² ~ χ²₉₉

¿Hay palés que salen más?Cuántos de los 4.950 pares posibles salieron juntos 0, 1, 2… veces (barras), frente a la distribución de Poisson del azar (línea), con λ = 3,63. Valor p = 0,73.
0246810

P(k) = λᵏ e−λ ⁄ k!  ·  λ = pares ⁄ 4.950

Rachas de pares e imparesVeces que la 1ra cambió de par a impar (o al revés) frente a lo esperado por azar (prueba de Wald-Wolfowitz). Valor p = 0,001.
3.0660ObservadoAzar

μ = 2n₁n₂⁄n + 1  ·  σ² = 2n₁n₂(2n₁n₂ − n) ⁄ (n²(n − 1))  ·  z = (R − μ) ⁄ σ

Cada prueba compara lo observado con lo que daría un sorteo justo. Un valor p alto (cerca de 1) significa "normal"; uno muy bajo (menos de 0,01) es una señal que conviene vigilar. Ninguna prueba dice qué número saldrá.