82.1. Una distribución cambia cuando conocemos algo

C81 introdujo distribuciones conjuntas, marginalización e independencia. Una distribución conjunta describe cómo se reparte probabilidad entre configuraciones posibles antes de saber cuál de ellas ocurrió. El siguiente paso es preguntar qué ocurre cuando obtenemos información parcial.

Supongamos que conocemos la distribución conjunta de dos variables XX e YY, pero observamos que Y=1.Y=1.

Las configuraciones incompatibles con esa observación dejan de ser posibles para la pregunta actual. Las probabilidades restantes deben renormalizarse para volver a sumar uno. Esa es la idea central del condicionamiento.

La probabilidad condicional cambia la distribución relevante una vez conocido un evento, sin que ello introduzca necesariamente una nueva dinámica ni una nueva causalidad. Puede representar una actualización temporal, pero también una selección, una restricción lógica o simplemente una pregunta diferente sobre la misma distribución conjunta.

C82 desarrollará cuatro piezas conectadas, a saber, la probabilidad condicional, la ley de probabilidad total, el teorema de Bayes y la verosimilitud. El objetivo es mantenerlas separadas. La notación es compacta, mientras que los objetos matemáticos no son intercambiables.

82.2. Definición de probabilidad condicional

Sean AA y BB eventos con P(B)>0.P(B)>0.

La probabilidad de AA condicionada a BB se define por

P(A∣B)=P(A∩B)P(B).P(A\mid B) = \frac{P(A\cap B)}{P(B)}.

El numerador contiene los resultados que satisfacen simultáneamente ambos eventos. El denominador restringe el universo de referencia a aquellos resultados donde BB ocurre.

La definición puede reorganizarse como

P(A∩B)=P(A∣B)P(B).P(A\cap B) = P(A\mid B)P(B).

Intercambiando AA y BB,

P(A∩B)=P(B∣A)P(A).P(A\cap B) = P(B\mid A)P(A).

Estas dos factorizaciones de la misma probabilidad conjunta producirán el teorema de Bayes.

La división cumple una función precisa, pues normaliza la masa que permanece después de restringirnos a BB.

82.3. Retomar la tabla de C81

C81 estudió dos componentes binarios con

P(X=1)=0.3,P(Y=1)=0.4,P(X=1)=0.3, \qquad P(Y=1)=0.4,

y una distribución dependiente donde P(X=1,Y=1)=0.25.P(X=1,Y=1)=0.25.

La tabla conjunta era

Y=0Y=1X=00.550.15X=10.050.25.\begin{array}{c|cc} & Y=0 & Y=1\\ \hline X=0 & 0.55 & 0.15\\ X=1 & 0.05 & 0.25 \end{array}.

Si observamos Y=1Y=1, solo importa la segunda columna. Su masa total es P(Y=1)=0.40.P(Y=1)=0.40.

Por tanto,

P(X=1∣Y=1)=0.250.40=0.625.P(X=1\mid Y=1) = \frac{0.25}{0.40} = 0.625.

Análogamente,

P(X=0∣Y=1)=0.150.40=0.375.P(X=0\mid Y=1) = \frac{0.15}{0.40} = 0.375.

Ambas cantidades suman uno.

Condicionamiento de una tabla conjunta.
Condicionamiento de una tabla conjunta Ampliar figura ↗

El condicionamiento construye una nueva distribución sobre XX dentro del subconjunto de realizaciones compatibles con Y=1Y=1, sin alterar la tabla original.

82.4. El orden de la barra importa

En general,

P(A∣B)≠P(B∣A).P(A\mid B) \neq P(B\mid A).

En el ejemplo,

P(X=1∣Y=1)=0.625,P(X=1\mid Y=1)=0.625,

mientras

P(Y=1∣X=1)=0.250.30=56≈0.833333.P(Y=1\mid X=1) = \frac{0.25}{0.30} = \frac56 \approx0.833333.

La intersección P(X=1,Y=1)=0.25P(X=1,Y=1)=0.25 es la misma. Lo que cambia es el conjunto utilizado para normalizar.

Esta diferencia es responsable de numerosos errores de interpretación. Saber que un estado produce con frecuencia cierta señal no equivale a saber con qué frecuencia ese estado está presente entre quienes muestran la señal. Bayes es la identidad que conecta ambas direcciones.

82.5. Independencia como invariancia frente al condicionamiento

Si AA y BB son independientes,

P(A∩B)=P(A)P(B).P(A\cap B)=P(A)P(B).

Cuando P(B)>0P(B)>0,

P(A∣B)=P(A)P(B)P(B)=P(A).P(A\mid B) = \frac{P(A)P(B)}{P(B)} = P(A).

Por tanto, una forma equivalente de expresar independencia es

P(A∣B)=P(A).P(A\mid B)=P(A).

Conocer BB no modifica la probabilidad de AA.

En términos de variables, si XX e YY son independientes, la distribución de XX condicionada a conocer YY coincide con su distribución marginal. Si cambia, existe dependencia probabilística.

Esta formulación también aclara que independencia y ausencia de información están relacionadas respecto de la distribución considerada, no respecto de cualquier descripción posible del sistema.

82.6. Condicionar puede crear dependencia

C81 terminó con un ejemplo que permite mostrar una sutileza adicional.

Sean UU y VV bits independientes y equiprobables. Definamos

W=U⊕V.W=U\oplus V.

Marginalmente, UU y VV son independientes.

Ahora condicionemos en W=0.W=0.

Los únicos estados compatibles son (U,V)=(0,0)(U,V)=(0,0) y (U,V)=(1,1),(U,V)=(1,1), cada uno con probabilidad condicional 1/21/2.

Así,

P(U=1∣W=0)=12,P(U=1\mid W=0)=\frac12,
P(V=1∣W=0)=12,P(V=1\mid W=0)=\frac12,

pero

P(U=1,V=1∣W=0)=12.P(U=1,V=1\mid W=0)=\frac12.

Si fueran independientes después de condicionar, esperaríamos

12⋅12=14.\frac12\cdot\frac12=\frac14.

La igualdad no se cumple.

La independencia marginal puede desaparecer al restringir el espacio de resultados. El condicionamiento modifica la estructura de dependencia, además de reducir el conjunto de realizaciones consideradas.

82.7. Probabilidad condicional en variables continuas

La definición

P(A∣B)=P(A∩B)P(B)P(A\mid B) = \frac{P(A\cap B)}{P(B)}

requiere P(B)>0P(B)>0.

Si YY es una variable absolutamente continua,

P(Y=y)=0P(Y=y)=0

para cualquier valor individual. Por tanto, no podemos definir ingenuamente

P(X∈A∣Y=y)P(X\in A\mid Y=y)

dividiendo por P(Y=y)P(Y=y).

Cuando existe una densidad conjunta, utilizamos la densidad condicional

fX∣Y(x∣y)=fX,Y(x,y)fY(y),f_{X\mid Y}(x\mid y) = \frac{f_{X,Y}(x,y)}{f_Y(y)},

para valores donde fY(y)>0.f_Y(y)>0.

La teoría general utiliza distribuciones condicionales regulares. Para este nivel basta con retener que condicionar en un valor continuo exacto difiere de condicionar en un evento discreto de probabilidad positiva.

82.8. Particiones y ley de probabilidad total

Supongamos que

H1,…,HmH_1,\ldots,H_m

forman una partición del espacio muestral:

Hi∩Hj=∅(i≠j),H_i\cap H_j=\varnothing \quad (i\neq j),

y

⋃i=1mHi=Ω.\bigcup_{i=1}^m H_i=\Omega.

Para cualquier evento EE,

E=⋃i=1m(E∩Hi).E = \bigcup_{i=1}^m(E\cap H_i).

Las intersecciones son disjuntas, de modo que

P(E)=∑i=1mP(E∩Hi).P(E) = \sum_{i=1}^mP(E\cap H_i).

Usando

P(E∩Hi)=P(E∣Hi)P(Hi),P(E\cap H_i) = P(E\mid H_i)P(H_i),

obtenemos

P(E)=∑i=1mP(E∣Hi)P(Hi).P(E) = \sum_{i=1}^m P(E\mid H_i)P(H_i).

Esta es la ley de probabilidad total, que permite calcular la probabilidad marginal de una observación combinando sus probabilidades bajo distintos estados posibles.

82.9. Un régimen oculto y un sensor imperfecto

Consideremos un sistema que puede encontrarse en un régimen oculto HH con probabilidad previa P(H)=0.02.P(H)=0.02.

El régimen complementario es HcH^c, con probabilidad 0.980.98.

Disponemos de un sensor binario que emite una señal positiva ++. Si el sistema está en HH,

P(+∣H)=0.95.P(+\mid H)=0.95.

Si no está en HH,

P(−∣Hc)=0.90.P(-\mid H^c)=0.90.

Por tanto,

P(+∣Hc)=0.10.P(+\mid H^c)=0.10.

Antes de invertir la inferencia necesitamos la probabilidad marginal de una señal positiva. Por probabilidad total,

P(+)=P(+∣H)P(H)+P(+∣Hc)P(Hc).P(+) = P(+\mid H)P(H) + P(+\mid H^c)P(H^c).

Sustituyendo,

P(+)=0.95(0.02)+0.10(0.98)=0.117.P(+) = 0.95(0.02)+0.10(0.98) = 0.117.

La mayoría de las oportunidades de producir falsos positivos procede del estado frecuente HcH^c. Ese dato pesará al interpretar la señal.

82.10. Derivación elemental del teorema de Bayes

Partimos de

P(H∩+)=P(+∣H)P(H).P(H\cap +) = P(+\mid H)P(H).

También,

P(H∩+)=P(H∣+)P(+).P(H\cap +) = P(H\mid +)P(+).

Igualando,

P(H∣+)P(+)=P(+∣H)P(H).P(H\mid +)P(+) = P(+\mid H)P(H).

Si P(+)>0P(+)>0,

P(H∣+)=P(+∣H)P(H)P(+).P(H\mid +) = \frac{P(+\mid H)P(H)}{P(+)}.

Sustituyendo la probabilidad total en el denominador,

P(H∣+)=P(+∣H)P(H)P(+∣H)P(H)+P(+∣Hc)P(Hc).P(H\mid +) = \frac{ P(+\mid H)P(H) }{ P(+\mid H)P(H) + P(+\mid H^c)P(H^c) }.

En nuestro ejemplo,

P(H∣+)=0.95(0.02)0.117≈0.162393.P(H\mid +) = \frac{0.95(0.02)}{0.117} \approx 0.162393.

Una señal positiva aumenta la probabilidad de HH desde 0.02 hasta aproximadamente 0.162, lejos del 0.95 que caracteriza al sensor.

82.11. La misma cuenta mediante frecuencias naturales

La aritmética puede visualizarse sin porcentajes abstractos.

Imaginemos 10.000 realizaciones. Aproximadamente 200200 están en HH, y 98009800 en HcH^c.

Entre las 200 realizaciones en HH, el sensor produce aproximadamente 0.95(200)=1900.95(200)=190 positivos.

Entre las 9800 realizaciones en HcH^c, produce 0.10(9800)=9800.10(9800)=980 falsos positivos.

Por tanto, existen aproximadamente 190+980=1170190+980=1170 señales positivas. De ellas, solo 190 corresponden a HH.

Así,

1901170=0.162393.\frac{190}{1170} = 0.162393.

Este razonamiento es Bayes expresado mediante conteos esperados, con la ventaja de dejar visible el efecto de la tasa previa.

82.12. La tasa previa forma parte del resultado

El número

P(+∣H)=0.95P(+\mid H)=0.95

describe el comportamiento del sensor bajo un estado conocido, cantidad distinta de la probabilidad posterior del estado después de observar la señal.

La figura siguiente mantiene fija la sensibilidad en 0.95 y la tasa de falsos positivos en 0.10, pero varía P(H).P(H).

Posterior según probabilidad previa.
Posterior según probabilidad previa Ampliar figura ↗

Con el mismo sensor, una señal positiva tiene implicaciones distintas si HH era extremadamente raro o relativamente frecuente.

Ignorar la tasa previa equivale a descartar parte de la distribución necesaria para invertir la condición.

82.13. Una señal negativa

La misma estructura funciona para una señal negativa.

Tenemos

P(−∣H)=0.05,P(-\mid H)=0.05,
P(−∣Hc)=0.90.P(-\mid H^c)=0.90.

La probabilidad total es

P(−)=0.05(0.02)+0.90(0.98)=0.883.P(-) = 0.05(0.02)+0.90(0.98) = 0.883.

Por Bayes,

P(H∣−)=0.05(0.02)0.883≈0.001133.P(H\mid -) = \frac{0.05(0.02)}{0.883} \approx 0.001133.

Una señal negativa reduce la probabilidad del régimen oculto.

Bayes es una transformación coherente de una distribución previa en una posterior según qué observación ocurrió, con independencia de que el resultado suba o baje la confianza en la hipótesis.

82.14. Forma de odds y razón de verosimilitudes

Para dos hipótesis complementarias podemos trabajar con odds.

Las odds previas de HH son

O(H)=P(H)P(Hc).O(H) = \frac{P(H)}{P(H^c)}.

Después de observar EE,

O(H∣E)=P(H∣E)P(Hc∣E).O(H\mid E) = \frac{P(H\mid E)}{P(H^c\mid E)}.

Aplicando Bayes,

O(H∣E)=O(H)P(E∣H)P(E∣Hc).O(H\mid E) = O(H) \frac{P(E\mid H)}{P(E\mid H^c)}.

El factor

Λ(E)=P(E∣H)P(E∣Hc)\Lambda(E) = \frac{P(E\mid H)}{P(E\mid H^c)}

es una razón de verosimilitudes entre las dos hipótesis simples.

Para una señal positiva,

Λ(+)=0.950.10=9.5.\Lambda(+) = \frac{0.95}{0.10} = 9.5.

Las odds previas son

0.020.98≈0.0204.\frac{0.02}{0.98}\approx0.0204.

Multiplicarlas por 9.5 produce las odds posteriores. Convertidas de nuevo a probabilidad recuperan exactamente 0.1623930.162393.

82.15. Qué es una verosimilitud

La palabra verosimilitud cambia cuál variable consideramos fija.

Supongamos una familia de modelos parametrizados por θ\theta, con distribución

p(x∣θ).p(x\mid\theta).

Antes de observar los datos, para un θ\theta fijo,

x↦p(x∣θ)x\mapsto p(x\mid\theta)

es una distribución de probabilidad sobre resultados posibles.

Después de observar

x=xobs,x=x_{\mathrm{obs}},

definimos la función de verosimilitud

L(θ;xobs)=p(xobs∣θ).L(\theta;x_{\mathrm{obs}}) = p(x_{\mathrm{obs}}\mid\theta).

Ahora los datos permanecen fijos y variamos θ\theta.

La expresión algebraica puede ser la misma, mientras el objeto matemático cambia. Por sí sola, la verosimilitud carece de la normalización que exigiría una distribución de probabilidad sobre θ\theta.

82.16. Siete éxitos en diez ensayos de Bernoulli

Supongamos diez ensayos de Bernoulli independientes con probabilidad de éxito θ\theta y observamos k=7k=7 éxitos.

La probabilidad de observar exactamente siete éxitos, para un θ\theta fijado, es

P(K=7∣θ)=(107)θ7(1−θ)3.P(K=7\mid\theta) = \binom{10}{7} \theta^7(1-\theta)^3.

Una vez observado K=7K=7, la verosimilitud es

L(θ)=(107)θ7(1−θ)3.L(\theta) = \binom{10}{7} \theta^7(1-\theta)^3.
Verosimilitud binomial.
Verosimilitud binomial Ampliar figura ↗

La función alcanza su máximo en

θ^MV=710=0.7.\widehat\theta_{\mathrm{MV}} = \frac{7}{10} = 0.7.

Esto no significa

P(θ=0.7∣K=7)=1.P(\theta=0.7\mid K=7)=1.

Ni siquiera hemos especificado una distribución de probabilidad para θ\theta.

82.17. La verosimilitud no necesita integrar uno

Una densidad de probabilidad sobre θ\theta tendría que integrar uno.

La verosimilitud anterior satisface

∫01L(θ) dθ=(107)B(8,4)=111≈0.0909.\int_0^1 L(\theta)\,d\theta = \binom{10}{7} B(8,4) = \frac{1}{11} \approx0.0909.

Ese valor no plantea ninguna dificultad, porque la verosimilitud no está obligada a normalizar respecto de θ\theta.

Si multiplicamos L(θ)L(\theta) por una constante positiva cc que no depende de θ\theta,

L~(θ)=cL(θ),\widetilde L(\theta)=cL(\theta),

las comparaciones relativas entre valores de θ\theta no cambian.

Por eso a menudo se grafica

L(θ)max⁡θL(θ).\frac{L(\theta)}{\max_\theta L(\theta)}.

La escala vertical absoluta no debe interpretarse como una probabilidad del parámetro.

82.18. Verosimilitud, prior y posterior

Para convertir información sobre los datos en una distribución posterior sobre θ\theta necesitamos una distribución previa π(θ).\pi(\theta).

La forma continua de Bayes puede escribirse

π(θ∣x)=p(x∣θ)π(θ)∫p(x∣u)π(u) du.\pi(\theta\mid x) = \frac{ p(x\mid\theta)\pi(\theta) }{ \int p(x\mid u)\pi(u)\,du }.

El numerador es

verosimilitud×prior.\text{verosimilitud}\times\text{prior}.

El denominador normaliza el producto para que la posterior integre uno. También es la probabilidad marginal o densidad predictiva de los datos bajo el modelo previo completo.

Si elegimos una prior uniforme en [0,1][0,1] para el ejemplo binomial, la posterior es proporcional a θ7(1−θ)3.\theta^7(1-\theta)^3.

Ahora sí, después de normalizar, obtenemos una distribución de probabilidad sobre θ\theta.

La distinción debe conservarse:

L(θ;x)≠π(θ∣x).L(\theta;x) \neq \pi(\theta\mid x).

Bayes conecta ambos objetos sin identificarlos. Mantener esta separación evita leer una curva de verosimilitud como si ya expresara grados posteriores de creencia o frecuencias poblacionales de los parámetros.

82.19. La evidencia en el denominador

En la forma discreta,

P(Hi∣E)=P(E∣Hi)P(Hi)P(E).P(H_i\mid E) = \frac{ P(E\mid H_i)P(H_i) }{ P(E) }.

El denominador

P(E)=∑jP(E∣Hj)P(Hj)P(E) = \sum_j P(E\mid H_j)P(H_j)

recibe a menudo el nombre de evidencia o probabilidad marginal de los datos.

Su función inmediata es normalizar las probabilidades posteriores:

∑iP(Hi∣E)=1.\sum_i P(H_i\mid E)=1.

El denominador tiene además significado probabilístico propio. Responde a cuán probable era observar EE bajo la mezcla completa de hipótesis, antes de saber cuál de ellas produjo los datos.

En el ejemplo del sensor, P(+)=0.117.P(+)=0.117.

Ese valor depende tanto de las características del sensor como de la distribución previa de los estados, de modo que cambia con el contexto.

82.20. Actualización secuencial

Bayes puede aplicarse más de una vez.

Supongamos que recibimos dos señales positivas

E1=+,E2=+.E_1=+, \qquad E_2=+.

Si asumimos que son condicionalmente independientes dado HH y dado HcH^c,

P(E1,E2∣H)=0.952,P(E_1,E_2\mid H) = 0.95^2,

mientras

P(E1,E2∣Hc)=0.102.P(E_1,E_2\mid H^c) = 0.10^2.

Entonces

P(H∣+,+)=0.02(0.95)20.02(0.95)2+0.98(0.10)2≈0.648115.P(H\mid +,+) = \frac{ 0.02(0.95)^2 }{ 0.02(0.95)^2+0.98(0.10)^2 } \approx 0.648115.
Actualización secuencial.
Actualización secuencial Ampliar figura ↗

El segundo positivo aporta información adicional porque hemos supuesto que, una vez fijado el estado oculto, las dos mediciones aportan evidencia independiente.

82.21. Repetir evidencia no siempre multiplica información

La independencia condicional del ejemplo anterior es un supuesto.

Si dos sensores comparten la misma fuente de error, o si la segunda lectura reutiliza datos de la primera, entonces

P(E1,E2∣H)≠P(E1∣H)P(E2∣H)P(E_1,E_2\mid H) \neq P(E_1\mid H)P(E_2\mid H)

en general.

Multiplicar dos razones de verosimilitud como si provinieran de observaciones independientes puede contar dos veces la misma información.

Esta es una conexión directa con C81. La dependencia importa para describir el sistema y también para saber cuánto aprendizaje aporta una colección de observaciones.

Una secuencia de cien mediciones casi idénticas y fuertemente correlacionadas no equivale automáticamente a cien piezas independientes de evidencia.

82.22. Independencia condicional no es independencia marginal

También puede ocurrir lo contrario del ejemplo XOR.

Supongamos un estado oculto HH que modifica simultáneamente dos sensores AA y BB. Dentro de cada estado, los errores de ambos sensores pueden ser independientes:

P(A,B∣H)=P(A∣H)P(B∣H),P(A,B\mid H) = P(A\mid H)P(B\mid H),

y

P(A,B∣Hc)=P(A∣Hc)P(B∣Hc).P(A,B\mid H^c) = P(A\mid H^c)P(B\mid H^c).

Sin embargo, al ignorar HH, AA y BB pueden aparecer dependientes porque comparten una causa estadística común dentro del modelo.

Condicionar puede entonces eliminar dependencia o crearla, según la arquitectura de la distribución conjunta.

No existe una regla general según la cual «condicionar siempre reduce incertidumbre entre variables» o «si dos variables son independientes lo seguirán siendo bajo cualquier condición».

82.23. Condicional no significa causal

La expresión

P(Y∣X)P(Y\mid X)

no afirma que XX cause YY.

Puede describir una relación temporal, una asociación simultánea, una selección o una restricción puramente matemática.

En XOR, condicionar en WW crea dependencia entre UU y VV, aunque WW fue definido a partir de ambos.

La probabilidad condicional pertenece al lenguaje de distribuciones. La causalidad requiere estructura adicional sobre qué variables pueden intervenirse, qué relaciones se mantienen bajo intervención y qué supuestos conectan observación con mecanismo.

Esta separación reaparecerá cuando el currículo llegue a inferencia causal. En C82 basta con no introducir causalidad mediante la barra vertical.

82.24. La inversión accidental de la condición

Supongamos que un patrón EE aparece frecuentemente cuando un sistema está en régimen HH:

P(E∣H)=0.95.P(E\mid H)=0.95.

No podemos concluir

P(H∣E)=0.95.P(H\mid E)=0.95.

Para calcular la segunda cantidad necesitamos al menos P(H)P(H) y la frecuencia de EE cuando HH no está presente.

El error surge porque ambas expresiones comparten los mismos símbolos, pero normalizan sobre conjuntos diferentes.

Una forma práctica de detectarlo es traducir la notación:

P(E∣H)P(E\mid H)

significa «probabilidad de la evidencia si HH fuera el estado».

En cambio,

P(H∣E)P(H\mid E)

significa «probabilidad del estado dentro de los casos compatibles con la evidencia».

La traducción verbal no sustituye la matemática, pero hace visible una inversión accidental.

82.25. Bayes no elimina la dependencia del modelo

El teorema de Bayes es una identidad matemática. La calidad de una aplicación depende de las probabilidades que introducimos.

Si la sensibilidad del sensor está mal estimada, si la tasa previa corresponde a otra población o si asumimos independencia condicional donde no existe, la posterior puede estar mal calibrada aunque la aritmética sea impecable.

Bayes organiza coherentemente supuestos y datos, sin garantizar que esos supuestos representen bien el sistema.

La transparencia del cálculo es una ventaja porque permite localizar dónde entra cada supuesto:

posterior∝verosimilitud×prior.\text{posterior} \propto \text{verosimilitud} \times \text{prior}.

Cada término puede discutirse por separado.

82.26. Qué significa una prior

Una distribución previa no tiene una única interpretación obligatoria.

Puede representar información acumulada antes de la observación actual, describir una frecuencia poblacional o servir para regularizar parámetros. En otros contextos formaliza incertidumbre epistemológica sobre valores desconocidos.

Lo relevante es que su papel matemático sea explícito.

En el ejemplo del sensor, P(H)=0.02P(H)=0.02 es una probabilidad previa de estados. Si esa tasa cambia entre contextos, la posterior cambia aunque el mecanismo de medición permanezca igual.

En problemas con parámetros continuos, elegir una prior exige más decisiones sobre soporte, escala, concentración y sensibilidad de la posterior frente a elecciones alternativas. Este capítulo no desarrollará aún inferencia bayesiana completa. El objetivo es entender la arquitectura elemental de la actualización.

82.27. Bayes con más de dos estados posibles

La forma binaria HH frente a HcH^c es útil para derivar el teorema, pero muchos sistemas admiten varios estados mutuamente excluyentes.

Supongamos una partición

H1,H2,…,HmH_1,H_2,\ldots,H_m

con probabilidades previas

pi=P(Hi),∑i=1mpi=1.p_i=P(H_i), \qquad \sum_{i=1}^m p_i=1.

Después observamos un evento EE. Para cada estado conocemos o postulamos la probabilidad

ℓi=P(E∣Hi).\ell_i=P(E\mid H_i).

La probabilidad total de la evidencia es

P(E)=∑j=1mℓjpj.P(E) = \sum_{j=1}^m \ell_j p_j.

Bayes produce entonces

P(Hi∣E)=ℓipi∑j=1mℓjpj.P(H_i\mid E) = \frac{\ell_i p_i} {\sum_{j=1}^m \ell_jp_j}.

Esta expresión puede interpretarse como dos operaciones. Primero asignamos a cada estado un peso no normalizado

wi=ℓipi.w_i=\ell_i p_i.

Después normalizamos:

P(Hi∣E)=wi∑jwj.P(H_i\mid E) = \frac{w_i}{\sum_jw_j}.

La evidencia redistribuye masa entre las posibilidades según dos factores, el peso que tenían antes y la compatibilidad de la observación con cada una, sin que ninguna quede seleccionada automáticamente.

Esto importa cuando varias explicaciones producen la misma señal con probabilidades parecidas. Una observación puede aumentar el peso de una familia de hipótesis sin concentrar la posterior en una sola.

También muestra por qué una hipótesis con verosimilitud máxima no tiene necesariamente la probabilidad posterior máxima. Si dos hipótesis H1H_1 y H2H_2 tienen

P(E∣H1)>P(E∣H2),P(E\mid H_1)>P(E\mid H_2),

pero H1H_1 tenía una probabilidad previa mucho menor, puede ocurrir

P(H1∣E)<P(H2∣E).P(H_1\mid E)<P(H_2\mid E).

La inferencia por máxima verosimilitud y la actualización bayesiana responden, por tanto, a preguntas relacionadas pero distintas.

82.28. Ceros, soporte y actualización

Bayes también obliga a prestar atención a los ceros.

Si

P(Hi)=0,P(H_i)=0,

entonces, siempre que la formulación ordinaria sea aplicable,

P(Hi∣E)=0.P(H_i\mid E)=0.

Ninguna verosimilitud finita puede recuperar masa posterior desde una hipótesis a la que el modelo asignó probabilidad previa exactamente nula.

Un cero previo es, por tanto, una afirmación estructural fuerte que excluye esa posibilidad del soporte del modelo, aunque de ello no se sigue que toda prior deba asignar probabilidad positiva a cualquier posibilidad imaginable.

Existe un segundo tipo de cero. Si

P(E∣Hi)=0,P(E\mid H_i)=0,

la observación EE descarta HiH_i dentro del modelo, porque esa hipótesis declaraba imposible la evidencia observada.

Un tercer caso es más problemático. Si

P(E)=0,P(E)=0,

la expresión elemental

P(Hi∣E)=P(E∣Hi)P(Hi)P(E)P(H_i\mid E) = \frac{P(E\mid H_i)P(H_i)}{P(E)}

no está definida. En modelos continuos, como vimos al condicionar en Y=yY=y, este problema requiere trabajar con densidades o distribuciones condicionales regulares en lugar de dividir probabilidades puntuales.

Los ceros exactos deben distinguirse de números muy pequeños. Una verosimilitud de 10−1210^{-12} no declara imposible un resultado. Puede volverse relevante si otras hipótesis asignan probabilidades aún menores.

En implementación numérica, esta diferencia también importa. Redondear a cero probabilidades extremadamente pequeñas puede convertir una afirmación de «muy improbable» en «imposible» y alterar actualizaciones posteriores de forma irreversible.

82.29. Verosimilitud relativa y log-verosimilitud

La verosimilitud se interpreta principalmente mediante comparaciones.

Para dos valores de parámetro θ1\theta_1 y θ2\theta_2,

L(θ1;x)L(θ2;x)\frac{L(\theta_1;x)} {L(\theta_2;x)}

indica cuánto más compatibles son los datos observados con uno que con otro dentro de la familia especificada.

Si multiplicamos toda la función por una constante c>0c>0 independiente de θ\theta,

cL(θ1;x)cL(θ2;x)=L(θ1;x)L(θ2;x).\frac{cL(\theta_1;x)} {cL(\theta_2;x)} = \frac{L(\theta_1;x)} {L(\theta_2;x)}.

Por eso los factores que no dependen del parámetro pueden omitirse al buscar el máximo de verosimilitud.

En muestras con muchas observaciones independientes, la verosimilitud suele ser un producto:

L(θ)=∏i=1np(xi∣θ).L(\theta) = \prod_{i=1}^{n} p(x_i\mid\theta).

Los productos de muchos números menores que uno pueden volverse extremadamente pequeños en aritmética de precisión finita. Se utiliza entonces la log-verosimilitud

ℓ(θ)=log⁡L(θ)=∑i=1nlog⁡p(xi∣θ).\ell(\theta) = \log L(\theta) = \sum_{i=1}^{n} \log p(x_i\mid\theta).

Como el logaritmo es estrictamente creciente,

arg max⁡θL(θ)=arg max⁡θℓ(θ).\operatorname*{arg\,max}_{\theta}L(\theta) = \operatorname*{arg\,max}_{\theta}\ell(\theta).

El máximo no cambia, pero el cálculo se vuelve más estable y los productos se transforman en sumas.

Esta operación tampoco convierte la verosimilitud en probabilidad. La log-verosimilitud puede ser negativa, positiva o desplazarse por una constante según la representación utilizada. Su función es comparar parámetros respecto de los mismos datos.

82.30. Protocolo elemental para una actualización bayesiana

Un cálculo de Bayes puede ser aritméticamente corto y conceptualmente incompleto. Antes de utilizar una posterior conviene explicitar la arquitectura del problema.

Primero se define el conjunto de hipótesis o estados. En el ejemplo fueron HH y HcH^c. Esa elección determina qué alternativas pueden recibir probabilidad. Si falta una explicación relevante, ninguna actualización dentro del conjunto cerrado puede descubrirla por sí sola.

Después se especifica la distribución previa. Debe quedar claro qué población, estado inicial o información previa representa. Una tasa previa transportada desde otro contexto puede ser matemáticamente válida y empíricamente inadecuada.

El tercer componente es el modelo de observación:

P(E∣Hi).P(E\mid H_i).

Aquí se concentran supuestos sobre sensibilidad, errores, dependencia entre mediciones y estabilidad del instrumento. Si estas probabilidades cambian con el contexto, el modelo debe representarlo.

El cuarto paso es calcular la evidencia

P(E)=∑iP(E∣Hi)P(Hi),P(E) = \sum_iP(E\mid H_i)P(H_i),

y después la posterior.

Finalmente debe examinarse sensibilidad. Si pequeñas modificaciones plausibles de la prior o del modelo de observación producen cambios grandes en la posterior, esa fragilidad forma parte de la conclusión.

Este protocolo evita presentar Bayes como una calculadora que recibe dos porcentajes y entrega una verdad actualizada. El teorema garantiza coherencia entre probabilidades dentro del modelo. La adecuación del modelo sigue siendo una cuestión científica.

82.31. Una posterior alta no demuestra que el conjunto de modelos sea adecuado

Existe una limitación que conviene introducir desde el comienzo. Bayes redistribuye probabilidad dentro del conjunto de hipótesis que el modelo contiene.

Supongamos tres hipótesis H1,H2,H3H_1,H_2,H_3. Después de observar EE, obtenemos

P(H1∣E)=0.90.P(H_1\mid E)=0.90.

Esto significa que H1H_1 concentra el 90 % de la masa posterior entre esas tres alternativas, dadas las probabilidades especificadas. No implica que H1H_1 reproduzca bien todos los rasgos de los datos ni que no exista una cuarta explicación ausente del modelo.

Puede ocurrir que la evidencia observada sea muy improbable bajo todas las hipótesis:

P(E∣H1)=10−6,P(E∣H2)=10−8,P(E∣H3)=10−9.P(E\mid H_1)=10^{-6}, \qquad P(E\mid H_2)=10^{-8}, \qquad P(E\mid H_3)=10^{-9}.

H1H_1 puede dominar la posterior relativa porque explica EE mejor que sus competidoras, aunque EE siga siendo extraordinariamente inesperado bajo H1H_1. La inferencia relativa y la adecuación absoluta son preguntas diferentes.

Este punto importa en sistemas complejos, donde el conjunto de modelos suele ser una simplificación deliberada. Una posterior concentrada puede justificar preferir provisionalmente una explicación dentro de la familia estudiada, pero no certifica que esa familia abarque el mecanismo real.

Por ello, además de actualizar probabilidades, conviene preguntar si el modelo genera observaciones comparables con las que realmente aparecen. Más adelante el currículo estudiará comprobaciones predictivas y comparación formal de modelos. Aquí basta conservar la advertencia conceptual:

mejor entre los modelos considerados\text{mejor entre los modelos considerados}

no equivale necesariamente a

modelo adecuado del fenoˊmeno.\text{modelo adecuado del fenómeno}.

Bayes resuelve cómo actualizar dentro de una arquitectura probabilística especificada. La suficiencia de esa arquitectura queda fuera de su alcance.

82.32. Verificación computacional

El paquete reproducible ejecuta controles analíticos y numéricos.

Primero reconstruye la tabla conjunta dependiente de C81 y verifica

P(X=1∣Y=1)=0.625P(X=1\mid Y=1)=0.625

y

P(Y=1∣X=1)=56.P(Y=1\mid X=1)=\frac56.

Después comprueba la ley de probabilidad total para el sensor, obteniendo P(+)=0.117,P(+)=0.117, y calcula las posteriores para señales positiva y negativa. La forma en odds reproduce exactamente los mismos resultados.

El código verifica además que dos señales positivas elevan la posterior a aproximadamente 0.648115.0.648115.

Esa comprobación utiliza explícitamente la hipótesis de independencia condicional.

Para la verosimilitud binomial con siete éxitos en diez ensayos, localiza el máximo cerca de

θ=0.7,\theta=0.7,

compara la verosimilitud con θ=0.5\theta=0.5 y comprueba que

∫01L(θ) dθ=111,\int_0^1L(\theta)\,d\theta=\frac1{11},

mostrando que la verosimilitud no está normalizada como densidad sobre θ\theta. Al multiplicarla por una prior uniforme y normalizar, la posterior sí integra aproximadamente uno.

Finalmente se reutiliza XOR para verificar que UU y VV, independientes marginalmente, dejan de ser independientes al condicionar en W=0W=0.

Se ejecutaron 29 comprobaciones automáticas y las 29 resultaron satisfactorias.

82.33. Qué queda establecido

La probabilidad condicional restringe una distribución a un subconjunto de resultados compatible con información conocida y renormaliza la masa restante.

La ley de probabilidad total reconstruye una probabilidad marginal a partir de una partición. Bayes invierte el orden de una condición utilizando esa marginal como normalización.

La verosimilitud responde a otra operación. Fija los datos observados y compara cómo distintos parámetros o hipótesis los harían probables, sin ser por sí sola una probabilidad del parámetro.

Estas distinciones permiten modelar sistemas parcialmente observados sin convertir toda incertidumbre en una sola cantidad. El estado oculto, la observación, la tasa previa y el mecanismo de medición ocupan lugares diferentes en la distribución.

El próximo capítulo, C83, introducirá la esperanza, la varianza y la covarianza. Allí pasaremos de describir distribuciones completas a construir resúmenes cuantitativos de localización, dispersión y dependencia lineal, con una cautela ya establecida. Resumir una distribución no equivale a conservar toda su estructura.


Referencias


Recursos reproducibles

Descargar recursos del capítulo (ZIP, 270 kB).

4 archivos CSV · 4 figuras · 1 programa Python.

El programa se ha ejecutado en una carpeta vacía y genera los archivos de datos y figuras incluidos en el paquete.

Ver instrucciones y contenido del paquete.