Estadística y estimación de parámetros: máxima verosimilitud, MAP e intervalos de confianza
Muestra, estimador, sesgo y por qué la varianza muestral divide entre n − 1, el método de máxima verosimilitud, la estimación bayesiana MAP y los intervalos de confianza, y por qué entrenar una red neuronal es estimar parámetros.
Serie Fundamentos de IA Parte 19 de 52
- Mapa de la inteligencia artificial: qué es cada cosa y cómo encaja todo
- Breve historia de la inteligencia artificial: de Turing a los modelos de lenguaje
- Datos, modelo, entrenamiento, inferencia: los conceptos básicos en un solo ejemplo
- Cómo leer la notación matemática: símbolos, subíndices, sumatorios y cuantificadores
- Conjuntos y relaciones: las piezas con las que se construyen las definiciones en IA
- Funciones y aplicaciones: por qué todo modelo de IA es una función
- Lógica y demostraciones: las reglas con las que se razona en IA
- Grafos y árboles: la estructura sobre la que se apoyan la búsqueda y las redes neuronales
- Álgebra lineal para aprendizaje automático: vectores, matrices y por qué la IA sobre todo multiplica
- Derivadas, gradiente y regla de la cadena: cómo sabe un modelo en qué dirección aprender
- Descenso de gradiente: cómo aprende un modelo
- Qué es la probabilidad: del lanzamiento de un dado a los axiomas de Kolmogórov
- Probabilidad condicionada y teorema de Bayes: cómo cambiar una estimación cuando te enteras de algo nuevo
- Variable aleatoria y función de distribución: una sola función que lo dice todo sobre la distribución
- Función de densidad y función de probabilidad: por qué una variable continua tiene probabilidad cero en cada punto
- Esperanza, varianza y covarianza: dónde está una distribución, qué anchura tiene y qué se relaciona con qué
- Resumen de las distribuciones de probabilidad: Bernoulli, binomial, Poisson, uniforme, normal y exponencial
- Ley de los grandes números y teorema central del límite: por qué las medias se comportan de forma tan predecible
- Estadística y estimación de parámetros: máxima verosimilitud, MAP e intervalos de confianza
- Teoría de la información: entropía, entropía cruzada y divergencia KL, el origen de las funciones de pérdida
- Búsqueda en el espacio de estados: BFS, DFS, UCS y A* con animaciones
- Juegos y adversarios: minimax, poda alfa-beta y Monte Carlo tree search
- Satisfacción de restricciones: backtracking, propagación de restricciones y AC-3 con un sudoku y un mapa de regiones checas
- Lógica y sistemas basados en conocimiento: encadenamiento hacia delante y hacia atrás, resolución y sistemas expertos
- Planificación: STRIPS, el mundo de bloques, heurísticas a partir de la relajación y la planificación en los agentes de IA
- Búsqueda local y recocido simulado: un viajante por las capitales regionales checas
- Algoritmos genéticos: población, cruce, mutación y selección en una mochila para la excursión
- Estrategias evolutivas y programación genética: de la regla de un quinto a la ley de Kepler
- Optimización por colonia de hormigas: cómo un rastro de feromona encuentra un camino corto
- Optimización por enjambre de partículas: inercia, experiencia propia y experiencia del enjambre
- Neuroevolución: evolución de los pesos y la forma de las redes neuronales, y el algoritmo NEAT
- Conjuntos difusos y funciones de pertenencia: cómo escribir «algo cálido»
- Lógica difusa y probabilidad: por qué «algo cálido» no significa «probablemente cálido»
- Lógica difusa: t-normas, implicaciones y reglas SI–ENTONCES
- Control difuso: Mamdani y Sugeno, del ventilador al péndulo invertido sobre un carro
- Agrupamiento difuso: el algoritmo fuzzy c-means y los puntos entre clusters
- Sistemas neuro-difusos (ANFIS): reglas difusas que aprenden de los datos
- Cadenas de Markov: matriz de transición, distribución estacionaria y PageRank
- Modelos ocultos de Márkov: el algoritmo de Viterbi y el casino deshonesto
- Redes bayesianas: por qué está mojada la hierba y qué significa eso para el aspersor
- El filtro de Kalman: cómo estimar a partir de mediciones imprecisas dónde está un coche y a qué velocidad va
- Tipos de aprendizaje automático: supervisado, no supervisado, semisupervisado y más
- Prompts, contexto y RAG: cómo darle a un modelo la información correcta
- Autocodificadores variacionales (VAE): imágenes nuevas a partir de un espacio latente continuo
- Redes generativas adversarias (GAN): el falsificador contra el detective
- Modelos de difusión: cómo el ruido puro se convierte en una imagen
- Procesos de decisión de Markov: estados, acciones, recompensas y ecuaciones de Bellman
- Q-learning y DQN: aprender el valor de las acciones, de la tabla a la red neuronal
- Gradiente de la política y PPO: aprender directamente la política
- Cómo funciona un agente de IA: un modelo de lenguaje en un bucle con herramientas
- Programación con agentes de IA: encargo, pruebas y revisión
- Evaluación y seguridad de los sistemas de IA: cómo medir y cómo vigilar
Temario completo de la serie 52 de 83 temas publicados
- Qué es un estimador, cuándo es insesgado y consistente, y por qué la varianza muestral divide entre n − 1.
- Cómo funcionan la máxima verosimilitud y la estimación MAP, y por qué son el entrenamiento y la regularización.
- Qué significa exactamente un intervalo de confianza del 95 % y cuál es el error más común al interpretarlo.
Conviene conocer antes: Ley de los grandes números y teorema central del límite: por qué las medias se comportan de forma tan predecible
La teoría de la probabilidad va del modelo a los datos: si sabes que la moneda es justa, calculas qué va a salir. La estadística hace el camino inverso, de los datos al modelo: han salido 6 cruces de 20, ¿cómo es esa moneda? Y eso es exactamente lo que hace el aprendizaje automático. El modelo tiene parámetros, los datos son lo que ha salido y el entrenamiento es una estimación de parámetros. Esta parte muestra qué es un buen estimador, de dónde sale el método de máxima verosimilitud, en el que se basa casi cualquier función de pérdida, y cómo decir con honestidad cuánta confianza podemos tener en una estimación.
Continúa la parte Ley de los grandes números y teorema central del límite.
1. Muestra, estadístico, estimador
Definición — muestra aleatoria
Una muestra aleatoria de tamaño n de una distribución con densidad (o función de probabilidad) f(x; θ) son n variables independientes e idénticamente distribuidas X1, …, Xn con esa distribución. El parámetro θ es desconocido.
Los datos de los que disponemos los entendemos como extraídos al azar de alguna distribución. El tipo de distribución lo elegimos nosotros (por ejemplo, una de Bernoulli: spam sí o no), pero su parámetro (qué parte del correo es spam) no lo conocemos y queremos averiguarlo a partir de los datos.
Definición — estadístico y estimador
Un estadístico es cualquier función de la muestra T = T(X1, …, Xn) que no depende del parámetro desconocido θ. Al estadístico que usamos para estimar θ lo llamamos estimador y lo denotamos θ̂.
Un estadístico es cualquier cálculo hecho con los datos: la media, el máximo, la mediana, la proporción de spam. Un estimador es un estadístico con una misión: dar con el parámetro desconocido. Como los datos son aleatorios, el estimador también lo es: con otros datos saldría distinto. Por eso nos preguntamos cómo se comportan los estimadores en promedio y a medida que hay más datos.
2. Un buen estimador: insesgado y consistente
Definición — estimador insesgado
Un estimador θ̂ es insesgado si E[θ̂] = θ para cualquier valor de θ. La diferencia E[θ̂] − θ se llama sesgo.
Un estimador insesgado se equivoca en cada intento hacia un lado o hacia el otro, pero en promedio da justo en el blanco. La media muestral es un estimador insesgado de la esperanza: lo vimos en la parte anterior.
Definición — estimador consistente
Un estimador θ̂n es consistente si para todo ε > 0 se cumple limn→∞ P(|θ̂n − θ| ≥ ε) = 0.
Con suficientes datos, el estimador acierta con tanta precisión como quieras. La media muestral es consistente: eso es justamente la ley de los grandes números.
Que un estimador sea insesgado no es algo que venga dado. Para estimar la varianza, lo natural parece tomar la desviación cuadrática media respecto a la media, (1/n) ∑ (Xi − X̄)². Pero ese estimador subestima la varianza de forma sistemática: medimos las desviaciones respecto a X̄, que se ha ajustado a los datos, en lugar de respecto a la verdadera μ.
Teorema — varianza muestral
Para una muestra aleatoria con varianza σ² se cumple E[(1/n) ∑ (Xi − X̄)²] = ((n − 1)/n) σ². Por eso el estimador insesgado de la varianza es la varianza muestral (o cuasivarianza)
s² = 1/(n − 1) · ∑i=1n (Xi − X̄)².
Hay que dividir entre n − 1, no entre n (corrección de Bessel). Un «grado de libertad» se gasta en estimar la media. En muestras pequeñas la diferencia se nota (con cinco valores, dividir entre cinco se queda corto en una quinta parte), pero con un millón de valores no importa.
Comprobémoslo: tiramos 2 000 veces cinco dados y cada vez calculamos las dos estimaciones de la varianza. La varianza real de un dado es 35/12 ≈ 2,92. La media de las estimaciones dividiendo entre n acabó en 2,38 (apunta a 4/5 · 2,92 ≈ 2,33); dividiendo entre n − 1, en 2,97:
3. Método de máxima verosimilitud
Definición — función de verosimilitud
Para los datos medidos x1, …, xn, la función de verosimilitud es la función del parámetro
L(θ) = ∏i=1n f(xi; θ).
Los datos están fijos; vamos cambiando el parámetro y nos preguntamos: «si el parámetro fuera θ, ¿qué probabilidad tendrían justo estos datos?». Para observaciones independientes, las probabilidades se multiplican. Ojo: L(θ) no es una probabilidad del parámetro; sumada sobre θ no tiene por qué dar 1.
Definición — estimador de máxima verosimilitud
El estimador de máxima verosimilitud (EMV, en inglés MLE) es θ̂ = argmaxθ L(θ). Como el logaritmo es creciente, es lo mismo que argmaxθ ∑i ln f(xi; θ).
Elige el parámetro con el que los datos medidos resultan más probables. El logaritmo convierte el producto en una suma: así es más cómodo calcular, y al ordenador no se le quedan en cero los números diminutos.
De 20 correos, 6 eran spam. Para una proporción de spam p se tiene L(p) = p6(1 − p)14. La derivada del logaritmo, 6/p − 14/(1 − p), se anula en p̂ = 6/20 = 0,3: exactamente la proporción de spam en los datos. En general, para k éxitos de n sale p̂ = k/n. Con 60 spams de 200 la estimación es la misma, pero el pico es mucho más estrecho: los datos dicen más sobre el parámetro.
Para la espera de un correo con distribución exponencial, el método de máxima verosimilitud da λ̂ = 1 / X̄. Con cinco esperas medidas de 1,2; 3,5; 0,4; 2,1 y 2,8 minutos, la media es de 2 minutos y λ̂ = 0,5 por minuto.
4. Estimación bayesiana MAP
Definición — estimador máximo a posteriori
Si π(θ) es la densidad a priori del parámetro, el estimador máximo a posteriori (MAP) es θ̂ = argmaxθ L(θ) · π(θ).
Es el teorema de Bayes de la parte sobre probabilidad condicionada aplicado a los parámetros: a posteriori ∝ verosimilitud × a priori. El MAP elige el parámetro más probable después de tener en cuenta lo que pensábamos de antemano. Si de entrada damos poco peso a los extremos (densidad a priori proporcional a p(1 − p)), la proporción de spam sale 7/22 ≈ 0,318: la estimación se desplaza un poco hacia el centro. A medida que llegan más datos, la influencia de la distribución a priori desaparece.
5. Intervalo de confianza
Definición — intervalo de confianza
El intervalo (T1, T2) con estadísticos T1 ≤ T2 es un intervalo de confianza para θ con nivel de confianza 1 − α si P(T1 < θ < T2) = 1 − α para todo θ.
En lugar de un solo número das un rango y, con él, una garantía: el procedimiento con el que calculas el rango acierta el parámetro real en el 95 % de los casos. Lo aleatorio son los extremos del intervalo, no el parámetro.
Para la esperanza y un n suficientemente grande, el teorema central del límite da un intervalo aproximado del 95 %: X̄ ± 1,96 · s/√n. Tiramos cuarenta veces 40 dados y con cada muestra calculamos un intervalo. La línea vertical es la esperanza real, 3,5:
El valor real está dentro de 38 de los 40 intervalos: el 95 %. En 20 000 repeticiones fue el 94,1 %, algo por debajo del 95 % prometido, porque σ lo estimamos a partir de los datos y n = 40 todavía no es infinito.
Un error frecuente: «con una probabilidad del 95 %, μ está en el intervalo (3,2; 3,9)». En un intervalo concreto, μ o está dentro o no lo está. Ese 95 % es una propiedad del procedimiento: cuántos intervalos aciertan a lo largo de muchas repeticiones. Una afirmación sobre la probabilidad del parámetro solo la da el enfoque bayesiano.
6. La estimación en la inteligencia artificial
- Entrenar es estimar por máxima verosimilitud. La entropía cruzada en clasificación es el logaritmo negativo de la verosimilitud dividido entre el número de ejemplos, y el error cuadrático medio es lo mismo para datos con ruido normal. Minimizar la pérdida = maximizar la verosimilitud.
- La regularización es MAP. La penalización de los pesos λ‖w‖² (L2, weight decay) corresponde a una distribución a priori normal sobre los pesos; la L1, a una de Laplace. La idea a priori de que «los pesos deben ser pequeños» frena el sobreajuste.
- Sesgo y varianza del modelo. Un modelo es un estimador; los modelos complejos tienen poco sesgo pero mucha varianza (sobreajuste), y los simples, al revés. La regularización acepta un poco más de sesgo a cambio de menos varianza.
- Conjunto de prueba e intervalo de confianza. La precisión en los datos de prueba es una estimación de la precisión real; sin un intervalo de confianza no sabrás si la diferencia entre dos modelos es real. La batch normalization estima la media y la varianza a partir de cada lote, y los modelos de lenguaje se entrenan por máxima verosimilitud: maximizan la verosimilitud de la siguiente palabra.
7. Resumen
- La estadística va de los datos al modelo: un estimador es un cálculo con los datos que debe dar con un parámetro desconocido.
- Un buen estimador es insesgado (en promedio apunta bien) y consistente (se vuelve más preciso con más datos); la varianza muestral divide entre n − 1.
- El método de máxima verosimilitud elige el parámetro con el que los datos son más probables; para una proporción, es k/n.
- El MAP añade una idea a priori; en el aprendizaje automático, eso es la regularización.
- Un intervalo de confianza del 95 % sale de un procedimiento que acierta el parámetro real en el 95 % de los casos.
¿Lo ves de otra manera? Escribe a support@coden-app.com.