OpenBetaMatemática

← Temas

Funciones de masa y densidad

La CDF te da la probabilidad acumulada hasta un punto. Pero para describir el comportamiento local —dónde están las modas, la asimetría, el grosor de las colas— necesitás su ritmo de cambio. Esa es la densidad: la tasa local de probabilidad, cuánta probabilidad hay concentrada cerca de cada punto. Es la derivada de la CDF en el caso continuo y el salto en el discreto.

Dos definiciones, una idea

PMF (caso discreto). Si XX es discreta, su función de masa es pX(x)=P(X=x)=FX(x)FX(x).p_X(x) = P(X = x) = F_X(x) - F_X(x^-). Es, literalmente, el salto de FXF_X en xx.

PDF (caso continuo). Si XX es continua, su densidad fX0f_X \ge 0 es la única función tal que FX(c)=cfX(t)dtpara todo c.F_X(c) = \int_{-\infty}^{c} f_X(t)\, dt \quad \text{para todo } c. Por el teorema fundamental del cálculo, fX=FXf_X = F_X' donde FXF_X es derivable.

Son el salto y la pendiente de la misma CDF. Movelo y miralo:

Densidad = derivada de la CDF · área = probabilidad — interactivo
1F(c) = P(X ≤ c)pendiente = f(c)-3-2-10123f(x) = F′(x)

El área sombreada bajo la densidad hasta c es exactamente F(c), y la pendiente de la CDF en c es la densidad f(c). Dos caras de lo mismo: f = F′, F = ∫ f. La probabilidad es área, no altura —un punto solo tiene ancho cero, por eso P(X = c) = 0.

Las propiedades que la caracterizan

No hacen falta axiomas nuevos: son la traducción de los de Kolmogorov (no negatividad y normalización) al lenguaje de densidades. Una función es PMF/PDF válida si y solo si:

  1. No negatividad: f(x)0f(x) \ge 0 para todo xx.
  2. Normalización: xpX(x)=1\displaystyle\sum_x p_X(x) = 1 (discreto) o fX(x)dx=1\displaystyle\int_{-\infty}^{\infty} f_X(x)\, dx = 1 (continuo).

Una densidad no es una probabilidad

Es el malentendido más común, y vale la pena clavarlo. En el caso continuo, fX(x)f_X(x) no es P(X=x)P(X = x) —de hecho P(X=x)=0P(X = x) = 0 para todo xx—. La densidad es probabilidad por unidad de longitud; lo que es una probabilidad es el área:

P(XA)=AfX(x)dx(continuo),P(XA)=xApX(x)(discreto).P(X \in A) = \int_A f_X(x)\, dx \quad (\text{continuo}), \qquad P(X \in A) = \sum_{x \in A} p_X(x) \quad (\text{discreto}).
PMF pX(x)p_X(x)PDF fX(x)f_X(x)
Interpretaciónprobabilidad de X=xX = xdensidad por unidad de longitud
Rango[0,1][0, 1][0,)[0, \infty)puede pasar de 1
P(X=x)P(X = x)=pX(x)= p_X(x)=0= 0 siempre
Probabilidad de AAxApX(x)\sum_{x \in A} p_X(x)AfX(x)dx\int_A f_X(x)\, dx

Que fX(x)=5,7f_X(x) = 5{,}7 es perfectamente válido. El caso extremo: XUniforme(0,0,1)X \sim \text{Uniforme}(0,\,0{,}1) tiene fX(x)=10f_X(x) = 10 sobre [0,0,1][0,\,0{,}1] —y aun así integra a 11, porque el área es 10×0,110 \times 0{,}1. El catálogo de formas habituales:

Catálogo de densidades — interactivo
-4-3-2-101234

La gaussiana es el caballo de batalla (CLT, error de medición, log-precios en Black–Scholes). Su pico vale 1/(σ√2π): con σ chico la densidad pasa de 1 —y está bien, porque lo que integra a 1 es el área, no la altura.

Operativa con densidades

Casi todo lo que se calcula con una variable pasa por su densidad:

  • Probabilidades: el área (o la suma) sobre el conjunto, como arriba.
  • Esperanza (anticipo, LOTUS): E[g(X)]=Rg(x)fX(x)dx\displaystyle E[g(X)] = \int_{\mathbb R} g(x)\, f_X(x)\, dx (continuo) o xg(x)pX(x)\sum_x g(x)\, p_X(x) (discreto).
  • Transformación (anticipo): para Y=g(X)Y = g(X) con gg monótona y derivable, fY(y)=fX(g1(y))ddyg1(y)\displaystyle f_Y(y) = f_X\bigl(g^{-1}(y)\bigr)\,\Bigl|\tfrac{d}{dy} g^{-1}(y)\Bigr| —el jacobiano corrige el estiramiento del eje (es la derivada de la identidad de CDFs FY=FXg1F_Y = F_X \circ g^{-1} que vimos antes).
  • Soporte: supp(X)={x:fX(x)>0}\operatorname{supp}(X) = \{x : f_X(x) > 0\}, la región donde la variable "vive".
  • Cuantiles: qα=FX1(α)=inf{x:FX(x)α}q_\alpha = F_X^{-1}(\alpha) = \inf\{x : F_X(x) \ge \alpha\} —la mediana es q0,5q_{0{,}5}, y el VaR al 5%5\% es un cuantil.

La unificación: derivada de Radon–Nikodym

PMF y PDF parecen objetos distintos, pero bajo la lente de la teoría de la medida son el mismo: ambos son la derivada de Radon–Nikodym de la ley PXP_X respecto de una medida de referencia μ\mu,

P(XA)=AdPXdμ(x)dμ(x).P(X \in A) = \int_A \frac{dP_X}{d\mu}(x)\, d\mu(x).

Lo único que cambia es qué μ\mu usás:

CasoMedida de referencia μ\mudPX/dμdP_X/d\mu
Discretomedida de conteoPMF pXp_X
Continuomedida de LebesguePDF fXf_X

Por eso se puede usar una sola notación para los dos casos: el cálculo es idéntico, solo cambia la medida. Y es también la razón por la que un cambio de medida dQ=ZdPdQ = Z\, dP es simplemente cambiar la densidad respecto de una misma referencia: valuar bajo Q\mathbb Q es, operativamente, tomar una esperanza bajo P\mathbb P con otra densidad.

Lo que sigue

Con la densidad ya tenemos el objeto que se manipula a diario. El próximo paso es exprimirle números: esperanza, varianza y momentos —el resumen que convierte una distribución entera en unas pocas cantidades— y de ahí a los estimadores.