Distribución estable
Distribuciones estables | ||
---|---|---|
Distribución α-estable simétrica con factor de escala unitario Distribuciones estables asimétricas centradas con factor de escala unitario Función de densidad de probabilidad | ||
Función de distribución para las distribuciones α-estables simé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 | |
Dominio | x ∈ R, 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 | |
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 escala.
La familia de distribuciones estables a veces se denomina también distribución α-estable de Lévy, en honor a Paul Lévy, el primero 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ón
[editar]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.
Referencias
[editar]- ↑ a b B. Mandelbrot, "The Pareto-Lévy Law and the Distribution of Income", International Economic Review 1960
- ↑ Paul Lévy, Calcul des probabilités, 1925.
- ↑ B.Mandelbrot, Stable Paretian Random Functions and the Multiplicative Variation of Income, Econometrica 1961 http://www.jstor.org/stable/pdfplus/1911802.pdf
- ↑ B. Mandelbrot, The variation of certain Speculative Prices, The Journal of Business 1963 [1]
- ↑ Eugene F. Fama, Mandelbrot and the Stable Paretian Hypothesis, The Journal of Business 1963
- ↑ a b Nolan, John P. «Stable Distributions - Models for Heavy Tailed Data». Archivado desde el original el 17 de julio de 2011. Consultado el 21 de febrero de 2009.
- ↑ Voit, Johannes (2005). The Statistical Mechanics of Financial Markets - Springer. Springer. doi:10.1007/b137351.
Enlaces externos
[editar]- The STABLE program for Windows is available from John Nolan's stable webpage: http://academic2.american.edu/~jpnolan/stable/stable.html Archivado el 30 de octubre de 2006 en Wayback Machine.. 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