Distribuciones de Pearson

De Wikipedia, la enciclopedia libre
Saltar a: navegación, búsqueda
Diagrama del Sistema de Pearson, mostrando distribuciones de los tipos I, III, VI, V, y IV en términos de β1 (asimetría cuadrada) y β2 (curtosis tradicional)

La distribución de Pearson en una familia de distribuciones probabilísticas continuas. Fue publicada por primera vez por Karl Pearson en 1895 y subsecuentemente extendida por él en 1901 y 1916 en una serie de artículos de bioestadística.

Historia[editar]

El sistema Pearson fue originalmente ideado en un esfuerzo para modelar observaciones visiblemente asimétricas. Era bien conocido en aquel tiempo cómo ajustar un modelo teórico para acomodar los primeros dos cumulantes o los momentos de observados datos: Cualquier distribución de probabilidad puede estar extendida directamente para formar una familia de escala de posición. Excepto en los casos patológicos, una familia de escala de posición puede estar hecha para acomodar la media (primer cumulante) y la varianza (segundo cumulante) arbitrariamente bien. Sin embargo, no era conocido cómo construir distribuciones de probabilidad en las cuales la asimetría (tercer cumulante estándar) y la curtosis (cuarto cumulante estándar) pudieron estar ajustados igualmente. Esta necesidad surgió al intentar acomodar modelos teóricos conocidos a datos observados que exhibieron asimetría. Los ejemplos de Pearson incluyen datos de supervivencia, cuáles son usualmente asimétricos. En su escrito original, Pearson (1895, p. 360) identificó cuatro tipos de distribuciones (numeradas del I al IV) además de la distribución normal (la cual era originalmente conocida como tipo V). La clasificación dependió en si las distribuciones estaban definidas en un intervalo definido, en una semirrecta, o en los reales y si estaban potencialmente asimétricas o necesariamente simétricas. Un segundo escrito (Pearson 1901) arregló dos omisiones: Redefinió la distribución de tipo V (originalmente incluía la distribución normal, ahora incorporaba la distribución gamma inversa) e introdujo la distribución de tipo VI. Conjuntamente los primeros dos documentos de identificación cubren los cinco tipos principales del sistema Pearson (I, III, VI, V y IV). En un tercer escrito, Pearson (1916) introdujo aún más casos especiales y subtipos (del VII al XII).

Rhind (1909, pp. 430–432) ideó una forma sencilla de visualizar el espacio de parámetros del sistema Pearson, el cual fue adoptado por Pearson (1916, plate 1 and pp. 430ff., 448ff.). Los tipos de Pearson son caracterizados por dos cantidades, comúnmente referidas como β1 y β2. El primero es el cuadrado de la asimetría: \beta_1 = \gamma_1^2 donde γ1 es la asimetría o el tercer momento estandarizado. El segundo es el curtosis tradicional o cuarto momento estandarizado: β2 = γ2 + 3. Tratamientos modernos definen kurtosis γ2 en términos de cumunlant en vez de momentos, por lo tanto una distribución normal tenemos γ2 = 0 y β2 = 3. Aquí seguimos el precedente histórico y usamos β2. EL diagrama a la derecha muestra dada una distribución concreta a qué tipo de Pearson pertenece (identificado por el punto (β1, β2)). Muchas de las distribuciones asimétricas y no mesocúrtica que hoy nos son familiares, no eran conocidas a principios de 1890. Lo que hoy se conoce como distribución beta había sido usada por Thomas Bayes como la Probabilidad a posteriori del parámetro de la distribución de Bernoulli en su trabajo de 1763 sobre la probabilidad inversa. La distribución beta ganó prominencia debido a su pertenencia al sistema Pearson y era conocida hasta los años 1940 como la distribución Pearson tipo I. [1] (La distribución de Pearson tipo II es un caso especial derivada del tipo I, pero ya no es tratada por separado.) La distribución gamma originada como resultado del trabajo de Pearson (Pearson 1893, p. 331; Pearson 1895, pp. 357, 360, 373–376) y era conocida como la distribución de Pearson tipo III, antes de adquirir su nombre moderno en 1930s y 1940s. .[2] El artículo de Pearson escrito en 1895 introdujo la distribución de tipo IV, la cual contiene la distribución t-Student como caso especial, precediendo por varios años a William Sealy Gosset. En su artículo de 1901 introdujo la distribución gamma inversa (tipo V) y la distribución beta prima (tipo VI).

Definición[editar]

Una función de densidad de Pearson, p, está definida para ser una solución válida a una ecuación diferencial (cf. Pearson 1895, p. 381)

\frac{p'(x)}{p(x)} + \frac{a+x-\lambda}{b_2 (x-\lambda)^2 + b_1 (x-\lambda) + b_0} = 0. \qquad (1) \!
Con:
 b_0=\frac{4 \beta_2-3 \beta_1}{10 \beta_2 -12\beta_1 -18} \mu_2 ,
 a=b_1=\sqrt{\mu_2} \sqrt{\beta_1}\frac{\beta_2+3}{10 \beta_2-12\beta_1 -18},
b_2=\frac{2 \beta_2-3 \beta_1 -6}{10 \beta_2-12\beta_1 -18} .

Según Ord,[3] Pearson ideó la forma subyacente de la ecuación (1), con base, primeramente, en la fórmula para la derivada del logaritmo de la función de densidad de la distribución normal (la cual da una función lineal) y, en segundo lugar , de una relación de recurrencia para los valores en la función de probabilidad de la masa de la distribución hipergeométrica (que produce la función lineal dividida por una estructura cuadrática).

En la ecuación (1), el parámetro a determina un punto estacionario, y por lo tanto bajo ciertas condiciones un moda de la distribución, ya que

p'(\lambda-a) = 0 \!

sale directamente de la ecuación diferencial.

Dado que nos enfrentamos a una ecuación diferencial lineal de primer orden con coeficientes variables, su solución es directa:

p(x) \propto \exp\left( -\!\int\!\!\frac{x-a}{b_2 x^2 + b_1 x + b_0} \,\mathrm{d}x \right).

La integral en esta solución simplifica considerablemente cuando ciertos casos especiales de integrando son considerados. Pearson (1895, p. 367) distingue dos casos principales, determinados por el signo del discriminante (y por tanto el número de raíces reales) de la función cuadrática.

f(x) = b_2\,x^2 + b_1\,x + b_0. \qquad (2)\!

Tipos particulares de distribución[editar]

Caso 1, discriminante negativo. La distribución de Pearson tipo IV[editar]

Si el discriminante de la función cuadrática (2) es negativo (b_1^2 - 4 b_2 b_0 < 0) no tiene raíces reales. Luego se define

 y = x + \frac{b_1}{2\,b_2} \!   y
\alpha = \frac{\sqrt{4\,b_2\,b_0 - b_1^2\,}}{2\,b_2}. \!

Observe que α es un número real bien definido y α ≠ 0, porque por suposición 4 b_2 b_0 - b_1^2 > 0 y por tanto b2 ≠ 0. Aplicando estas tres sustituciones, la función cuadrática (2) es transformada en

f(x) = b_2\,(y^2 + \alpha^2). \!

La ausencia de raíces reales es obvio en esta formulación ya que α2 es necesariamente positiva.

Ahora expresamos la solución de la ecuación diferencial (1) en función de y:

p(y) \propto \exp\left(- \frac{1}{b_2}\, \int\frac{y - \frac{b_1}{2\,b_2} - a}{y^2 + \alpha^2} \,\mathrm{d}y  \right). \!

Pearson (1895, p. 362) lo llamo el "caso trigonométrico" , debido a la integral

\int\frac{y-\frac{2\,b_2\,a + b_1}{2\,b_2}}{y^2 + \alpha^2} \,\mathrm{d}y = \frac{1}{2} \ln(y^2 + \alpha^2) - \frac{2\,b_2\,a + b_1}{2\,b_2\,\alpha}\arctan\left(\frac{y}{\alpha}\right) + C_0

Involucra la función trigonométrica inversa arcotangente. Entonces

p(y) \propto \exp\left[ -\frac{1}{2\,b_2} \ln\!\left(1+\frac{y^2}{\alpha^2}\right) -\frac{\ln\alpha}{b_2} +\frac{2\,b_2\,a + b_1}{2\,b_2^2\,\alpha} \arctan\left(\frac{y}{\alpha}\right) + C_1 \right]

Finalmente sea

m = \frac{1}{2\,b_2} \!   y
\nu = -\frac{2\,b_2\,a + b_1}{2\,b_2^2\,\alpha} \!

Aplicando estas sustituciones, obtenemos la función paramétrica :

p(y) \propto  \left[1 + \frac{y^2}{\alpha^2}\right]^{-m} \exp\left[-\nu \arctan\left(\frac{y}{\alpha}\right)\right]

Esta función de densidad sin normalizar tiene soporte en toda la línea real. Depende del parámetro de escala α > 0 y el parámetro de forma m>1/2 y v. Un parámetro se perdió cuando preferimos encontrar la solución a la ecuación diferencial(1) como una función de y o de x. Por lo tanto volvemos a introducir un cuarto parámetro, llamado parámetro de posición λ. Así hemos derivado la función densidad de la distribución de tipo Pearson IV:

p(x) = \frac{\left|\frac{\Gamma\!\left(m+\frac{\nu}{2}i\right)}{\Gamma(m)}\right|^2}{\alpha\,\mathrm{\Beta}\!\left(m-\frac12, \frac12\right)}
\left[1 + \left(\frac{x-\lambda}{\alpha}\right)^{\!2\,} \right]^{-m} \exp\left[-\nu \arctan\left(\frac{x-\lambda}{\alpha}\right)\right].

La normalización de las constantes involucra función gamma compleja (Γ) y la función beta (B).

Distribución de Pearson tipo VII[editar]

Diagrama de densidades de distribuciones Pearson tipo VII con λ = 0, σ = 1, y: γ2 = ∞ (rojo); γ2 = 4 (azul); y γ2 = 0 (negro)

El parámetro de la forma ν de la distribución de Pearson tipo IV controla su asimetría. Si fijamos su valor a cero, obtenemos una familia simétrica de tres parámetros. Este caso especial es conocido como Distribución de Pearson tipo VII (cf. Pearson 1916, p. 450). Su función de densidad es

p(x) = \frac{1}{\alpha\,\mathrm{\Beta}\!\left(m-\frac12, \frac12\right)} \left[1 + \left(\frac{x-\lambda}{\alpha}\right)^{\!2\,} \right]^{-m},

donde B denota la función Beta.

Una parametrización alternativa (y una ligera especialización) de la distribución tipo VII es obtenida permitiendo

\alpha = \sigma\,\sqrt{2\,m-3}, \!

Lo cual requiere m>3/2. Esto conlleva una pérdida menor de generalidad pero asegura que la varianza de la distribución existe y es igual a σ2. Ahora el parámetro m solo controla la curtosis de la distribución. Si m tiende a infinito como λ y σ se mantiene constante, la distribución normal emerge como un caso especial:

\lim_{m\to\infty}\frac{1}{\sigma\,\sqrt{2\,m-3}\,\mathrm{\Beta}\!\left(m-\frac12, \frac12\right)} \left[1 + \left(\frac{x-\lambda}{\sigma\,\sqrt{2\,m-3}}\right)^{\!2\,} \right]^{-m}
= \frac{1}{\sigma\,\sqrt{2}\,\Gamma\!\left(\frac12\right)} \times \lim_{m\to\infty} \frac{\Gamma(m)}{\Gamma\!\left(m-\frac12\right) \sqrt{m-\frac32}} \times \lim_{m\to\infty} \left[1 + \frac{\left(\frac{x-\lambda}{\sigma}\right)^2}{2\,m-3} \right]^{-m}
= \frac{1}{\sigma\sqrt{2\,\pi}} \times 1 \times \exp\!\left[-\frac12 \left(\frac{x-\lambda}{\sigma}\right)^{\!2\,} \right]

Esta es la función de densidad de la distribución normal con media λ y desviación estándar σ.

Es conveniente exigir que m > 5/2 y dejar que:

m = \frac52 + \frac{3}{\gamma_2}. \!

Esta es otra especialización, y garantiza que los primeros cuatro momentos de la distribución existan. Más aun, la distribución de Pearson tipo VII parametrizada en términos de (λ, σ, γ2) tiene como media λ, como desviación estándar σ, asimetría cero y curtosis exceso es γ2).


Distribución t-Student[editar]

La distribución de Pearson tipo VII es equivalente a la distribución t-Student no estandarizada con parámetros ν > 0, μ, σ2 aplicando las siguientes sustituciones a su parametrización original.

 \lambda = \mu, \!
 \alpha = \sqrt{\nu\sigma^2}, \!  y
 m = \frac{\nu+1}{2}, \!

Observe que la restricción m > ½ se satisface.

La función de densidad resultante es:

p(x|\mu,\sigma^2,\nu) = \frac{1}{\sqrt{\nu\sigma^2}\,\mathrm{\Beta}\!\left(\frac{\nu}{2}, \frac12\right)} \left(1+\frac{1}{\nu}\frac{(x-\mu)^2}{\sigma^2}\right)^{-\frac{\nu+1}{2}},

La cual es más conocida como la densidad de distribución t-student.

Note además que esto implica que la Distribución de Pearson tipo VII subsume la distribución t-Student estándar y también la distribución de Cauchy estándar. En particular, la distribución t-Student estándar emerge como un subcaso cuando μ = 0 y σ2 = 1, equivalente a las siguientes sustituciones.

\lambda = 0, \!
\alpha = \sqrt{\nu}, \!   y
m = \frac{\nu+1}{2}, \!

La densidad de esta restringida familia de un solo parámetro es una t-student estándar:

p(x) = \frac{1}{\sqrt{\nu}\,\mathrm{\Beta}\!\left(\frac{\nu}{2}, \frac12\right)} \left(1 + \frac{x^2}{\nu} \right)^{-\frac{\nu+1}{2}},

Caso 2, discriminante no negativo[editar]

Si la función cuadrática (2) tiene discriminante no negativo ((b_1^2 - 4 b_2 b_0 \geq 0), tiene como raíces reales a1 y a2 (no necesariamente distintas):

a_1 = \frac{-b_1 - \sqrt{b_1^2 - 4 b_2 b_0}}{2 b_2}, \!
a_2 = \frac{-b_1 + \sqrt{b_1^2 - 4 b_2 b_0}}{2 b_2}, \!

En presencia de raíces reales, la función cuadrática (2) puede ser escrita como

f(x) = b_2\,(x-a_1)(x-a_2), \!

y por lo tanto la solución de la ecuación diferencial es:

p(x) \propto \exp\left( -\frac{1}{b_2} \int\!\!\frac{x-a}{(x - a_1) (x - a_2)} \,\mathrm{d}x \right). \!

Pearson (1895, p. 362) la llamó el "caso logarítmico", debido a la integral

\int\!\!\frac{x-a}{(x - a_1) (x - a_2)} \,\mathrm{d}x = \frac{(a_1-a)\ln(x-a_1) - (a_2-a)\ln(x-a_2)}{a_1-a_2} + C

involucra solo la función logarítmica, y no la función arcotangente como en el caso anterior.

Usando la sustitución

\nu = \frac{1}{b_2\,(a_1-a_2)} \!

obtenemos la siguiente solución a la ecuación diferencial (1):

p(x) \propto (x-a_1)^{-\nu (a_1-a)} (x-a_2)^{\nu (a_2-a)}.

Dado que esta densidad es sólo sabida hasta una constante escondida de proporcionalidad, esa constante puede variarse y la densidad puede escrita como sigue:

p(x) \propto \left(1-\frac{x}{a_1}\right)^{-\nu (a_1-a)} \left(1-\frac{x}{a_2}\right)^{ \nu (a_2-a)}

Distribución de Pearson tipo I[editar]

La Distribución de Pearson tipo I (una generalización de la distribución beta surge cuando las raíces de la ecuación cuadrática (2) son de signos opuestos, eso es , a_1 < 0 < a_2. Luego la solución p es soportada en intervalo (a_1, a_2). Aplicando la sustitución

x = a_1 + y (a_2 - a_1) \qquad \mbox{where}\ 0<y<1, \!

la cuál produce una solución en términos de y que esta soportada en el intervalo (0, 1):

p(y) \propto \left(\frac{a_1-a_2}{a_1}\;y\right)^{(-a_1+a)\nu} \left(\frac{a_2-a_1}{a_2}\;(1-y)\right)^{(a_2-a)\nu}.

Uno puede definir

 m_1=\frac{a-a_1}{b_2 (a_1-a_2)} \!
 m_2=\frac{a-a_2}{b_2 (a_2-a_1)}\!

Reagrupando las constantes y parámetros, esto se simplifica a:

p(y) \propto y^{m_1} (1-y)^{m_2}, \!

Así \frac{x-\lambda-a_1}{a_2-a_1} sigue a \Beta(m_1+1,m_2+1) con \lambda=\mu_1-(a_2-a_1) \frac{m_1+1}{m_1+m_2+2}-a_1.

Resulta que m1, m2 > −1 es necesario y suficiente para que p sea una función de densidad de probabilidades.

Distribución de Pearson tipo II[editar]

La distribución de Pearson de tipo II es un caso especial de la familia de Pearson de tipo I restringida a distribuciones simétricas.

Para la curva de Pearson de tipo II,[4]

y = y_{0}\left(1-\frac{x^2}{a^2}\right)^m

Donde

x = \sum d^2/2 -(n^3-n)/12

La ordenada, y, es la frecuencia de \sum d^2. La curva de Pearson de tipo II es usada en computar la tabla de coeficientes de correlación significativos para el coeficiente de correlación de Spearman cuando el número de elementos en una serie es menor a 100(o 30 dependiendo en algunas fuentes). Luego, la distribución imita una distribución t- student estándar. Para la tabla de valores, ciertos valores son usados como constantes en la ecuación previa:

m     = \frac{5\beta_{2}-9}{2(3-\beta_{2})}
a^2   = \frac{2\mu_{2}\beta_{2}}{3-\beta_{2}}
y_{0} = \frac{N[\Gamma(2m+2)]}{a[2^{2m+1}][\Gamma(m+1)]}

Los momentos de x usada son

\mu_{2} = (n-1)[(n^2+n)/12]^2
\beta_{2}=\frac{3(25n^4-13n^3-73n^2+37n+72)}{25n(n+1)^2(n-1)}

Distribución de Pearson tipo III[editar]

\lambda= \mu_1 + \frac{b_0}{b_1} - (m+1) b_1\!
b_0+b_1 (x-\lambda)\! es \mathrm{Gamma}(m+1,b_1^2)\!

La distribución de Pearson tipo III es una distribución gamma o una distribución chi-cuadrado.

Distribución de Pearson tipo V[editar]

Definiendo nuevos parámetros:

 C_1=\frac{b_1}{2 b_2}\!
\lambda=\mu_1-\frac{a-C_1} {1-2 b_2}\!
x-\lambda\! sigue una \operatorname{InverseGamma}(\frac{1}{b_2}-1,\frac{a-C_1}{b_2})\!

La distribución de Pearson tipo V es una distribución gamma inversa.

Distribución de Pearson tipo VI[editar]

\lambda=\mu_1 + (a_2-a_1) \frac{m_2+1}{m_2+m_1+2} - a_2\!
\frac{x-\lambda-a_2}{a_2-a_1}\! sigue una :\beta^{\prime}(m_2+1,-m_2-m_1-1)\!

La distribución de Pearson tipo VI es una distribución beta prima o una Distribución F.

Relación con otras Distribuciones[editar]

La familia Pearson subsume las siguientes distribuciones, entre los otros:

Aplicaciones[editar]

Estos modelos son utilizados en los mercados financieros, dado su habilidad para ser parametrizadas de un modo que tiene significado intuitivo para comerciantes de mercado. Un número de modelos está en uso actual que la captura la naturaleza estocástica de la volatilidad de tasas, acciones etcétera. Esta familia de distribuciones puede resultar ser una de lo más importantes. En los Estados Unidos, el Log Pearson III es la distribución predeterminada para el análisis de frecuencias de la inundación.

Notas[editar]

  1. Miller, Jeff; et al. (09-07-2006). «Beta distribution». Earliest Known Uses of Some of the Words of Mathematics. Consultado el 9 de diciembre de 2006.
  2. Miller, Jeff; et al. (07-12-2006). «Gamma distribution». Earliest Known Uses of Some of the Words of Mathematics. Consultado el 9 de diciembre de 2006.
  3. Ord J.K. (1972) p2
  4. Ramsey, Philip H. (01-09-1989). «Critical Values for Spearman's Rank Order Correlation». Consultado el 22 de agosto de 2007.

Fuentes[editar]

Fuentes Primarias[editar]

  • Pearson, Karl (1893). «Contributions to the mathematical theory of evolution [abstract]». Proceedings of the Royal Society 54 (326–330):  pp. 329–333. doi:10.1098/rspl.1893.0079. 
  • Rhind, A. (July/October 1909). «Tables to facilitate the computation of the probable errors of the chief constants of skew frequency distributions». Biometrika 7 (1/2):  pp. 127–147. 

Fuentes Secundarias[editar]

Referencias[editar]

  • Elderton, Sir W.P, Johnson, N.L. (1969) Systems of Frequency Curves. Cambridge University Press.
  • Ord J.K. (1972) Families of Frequency Distributions. Griffin, London.

Plantilla:Distribuciones de Probabilidades