Distribución estable

En teoría de la probabilidad, una distribución se denomina estable (o una variable aleatoria se denomina estable) si es una combinación lineal de dos o más copias independientes de una muestra aleatoria que tiene la misma distribución de probabilidad, salvo por quizá algún parámetro de localización o factor de escala.

Distribuciones estables
Distribuciones estables simétricas
Distribución α-estable simétrica con factor de escala unitario
Distribuciones estables asimétricas centradas.
Distribuciones estables asimétricas centradas con factor de escala unitario
Función de densidad de probabilidad
Funciones de distribución para las distribuciones α-estables simétricas
Función de distribución para las distribuciones α-estables simétricas
Función de distribución para distribuciones de Lévy asimétricas
Función de distribución para distribuciones estables asimétricas
Función de distribución de probabilidad
Parámetros

α ∈ (0, 2] — parámetro de estabilidad
β ∈ [−1, 1] — parámetro de asimetría (nótese que la asimetría no está definida)
c ∈ (0, ∞) — parámetro de escala

μ ∈ (−∞, ∞) — parámetro de localización
Dominio xR, o x ∈ [μ, +∞) si α < 1 y β = 1, o x ∈ (-∞, μ] si α < 1 y β = −1
Función de densidad (pdf) no hay forma analítica explícita, excepto para algunos valores de los parámetros
Función de distribución (cdf) no hay forma analítica explícita, excepto para algunos valores de los parámetros
Media μ cuando α > 1 y no definida en el resto de casos
Mediana μ cuando β = 0 y no definida en el resto de casos
Moda μ cuando β = 0 y no definida en el resto de casos
Varianza 2c2 cuando α = 2, para otros casos no es finita
Coeficiente de simetría 0 cuando α = 2, para otros casos no es finita
Curtosis 0 when α = 2, para otros casos no es finita
Entropía no hay forma analítica explícita, excepto para algunos valores de los parámetros
Función generadora de momentos (mgf) no definida
Función característica
donde

La familia de distribuciones estables a veces se denomina también distribución α-estable de Lévy, en honor a Paul Lévy, el primer en estudiar este tipo de distribuciones.[1][2]

De los cuatro parámetros que definen una distribución estable, el más significativo es el parámetro de estabilidad, α (ver ficha lateral). Las distribuciones estables satisfacen que 0 < α ≤ 2, correspondiendo el valor máximo con una distribución normal (que es el caso más sencillo de distribución estable). El valor α = 1 corresponde a la distribución de Cauchy. Las distribuciones estables no tienen una varianza finita si α < 2, más aún si α ≤ 1 ni siquiera tienen media finita. La importancia práctica de las distribuciones estables es que son "atractores" para la distribución de sumas de variables aleatorias independientes e idénticamente distribuidas (que además pertenecen a espacios Lp). La distribución normal define una subfamilia de distribuciones estables. Por el clásico teorema del límite central la suma de un conjunto de variables con idéntica distribución e independientes y que además tenga varianza finita, tenderá a una distribución normal a medida que el número de variables que interviene en la suma crece. Sin la restricción de varianza finita, el teorema del límite central no será aplicable, pero la suma de esas variables tenderá hacia una distribución estable (α < 2).

Mandelbrot denominó a las distribuciones estables con α < 2 como "distribuciones estables paretianas",[3][4][5]​ por Vilfredo Pareto. Mandelbrot usó el término para distribuciones estables "positivas" (es decir, máximamente asimétricas hacia la dirección positiva) con 1<α<2 como "distribuciones de Pareto-Lévy".[1]​ Además consideró que estas últimas eran relevantes para describir los precios de acciones y productos de consumo.

DefiniciónEditar

Una distribución no degenerada es estable si satisface la siguiente propiedad:

Sean X1 y X2 dos copias de una variable aleatoria X (es decir, dos instancias de variables aleatorias independientes, de variables que tienen la misma distribución que X). Entonces X se denomina estable si existen dos constantes a > 0 y b > 0 tales que la nueva variable aleatoria aX1 + bX2 tenga la misma distribución que cX + d para otras dos constantes c > 0 y d. La distribución se llama estictamente estable si esto sigue siendo cierto aún con d = 0.[6]

Puesto que la distribución normal, la distribución de Cauchy y la distribución de Lévy satisfacen esta propiedad, son casos particulares de distribuciones estables. Más aún, Lévy demostró que el conjunto de todas las distribuciones estables puede representarse como una familia dada por cuatro parámetros de distribuciones continuas. Dos de los parámetros representan el parámetro de localización μ y el factor de escala c, mientras que los otros dos β y α, se corresponden con el grado de asimetría y concentración (ver gráficas de la ficha lateral).

Aunque la densidad de probabilidad de estas distribuciones estables no admite una fórmula analítica cerrada (expresable en términos de funciones elementales), sin embargo, su función característica si admite una fórmula analítica cerrada. Cualquier distribución de probaiblidad dada por la transformada de Fourier de una función característica φ(t) del tipo:

 

Una variable aleatoria X es estable si su función característica puede escribirse como[6][7]

 

donde sgn(t) es simplemente la función signo de t y Φ viene dada por

 

para todo α, excepto α = 1, en cuyo caso:

 

μ ∈ R es un parámetro de localización y β ∈ [−1, 1] se denomina parámetro de asimetría. Nótese que en este contexto la asimetría "usual" no está bien definida, cuando α < 2 ya que la distribución no admite moments de segundo orden ni superiores (pero la asimetría usual usa en su definición el tercer momento central).

La razón de que una distribución que satisface las condiciones anteriores sea estable es que la función característica para la suma de dos variables aleatorias es igual al producto de las correspondientes funciones características. Al sumar dos variables aleatorias que tienen una distribución estable se obtiene otra variable con los mismos valores de α y β, pero posiblemente valores diferentes de μ y c (cambia la localización y la escala, pero no la el parámetro de estabilidad ni el parámetro de asimetría).

No cualquier función es la función característica de una distribución de probabilidad legítima, ya que una función de distribución real varía entre 0 y 1 son decrecer, pero la función caractecterística dada anteriormente podría no ser adecuada si los pámetros no están dentro del rango adecuado. El valor de la función característica para cierto valor t es el complejo conjugado de su valor en −t como debería ser para que una función de distribución sea real. En el caso más sencillo β = 0, la función característica es simplemente una función exponencial estirada, la distribución es simétrica alrededor de μ y en ese caso se denomina distribución α-estable simétrica de Lévy, frecuentemente abreviada como SαS. Cuando α < 1 y β = 1, la distribución tiene un soporte en [μ, ∞). El parámetro c > 0 es un factor de escala que mide el ancho típico de la distribución mientras que α es el exponente o índice de la distribución que especifica el comportamiento asintótico de dicha distribución.


ReferenciasEditar

  1. a b B. Mandelbrot, "The Pareto-Lévy Law and the Distribution of Income", International Economic Review 1960
  2. Paul Lévy, Calcul des probabilités, 1925.
  3. B.Mandelbrot, Stable Paretian Random Functions and the Multiplicative Variation of Income, Econometrica 1961 http://www.jstor.org/stable/pdfplus/1911802.pdf
  4. B. Mandelbrot, The variation of certain Speculative Prices, The Journal of Business 1963 [1]
  5. Eugene F. Fama, Mandelbrot and the Stable Paretian Hypothesis, The Journal of Business 1963
  6. a b Nolan, John P. «Stable Distributions - Models for Heavy Tailed Data». Consultado el 21 de febrero de 2009. 
  7. Voit, Johannes (2005). The Statistical Mechanics of Financial Markets - Springer. Springer. doi:10.1007/b137351. 

Enlaces externosEditar

  • The STABLE program for Windows is available from John Nolan's stable webpage: http://academic2.american.edu/~jpnolan/stable/stable.html. It calculates the density (pdf), cumulative distribution function (cdf) and quantiles for a general stable distribution, and performs maximum likelihood estimation of stable parameters and some exploratory data analysis techniques for assessing the fit of a data set.
  • Matlab codes for simulation of stable variables and estimation of stable parameters are available from RePEc: https://ideas.repec.org/e/pwe42.html#software