La Serie de Gram–Charlier A (nombrada en honor a Jørgen Pedersen Gram y Carl Charlier), y la serie de Edgeworth (en honor a Francis Ysidro Edgeworth) son series que aproximan una distribución de probabilidad en términos de sus cumulantes. Las series son la misma, pero el orden de sus términos varían (y por consiguiente su exactitud para truncar la serie).
Serie de Gram–Charlier A
La idea fundamental de estas expansiones es escribir una función característica con función de densidad de probabilidad que se aproxime en términos de la función característica de una función de una distribución con propiedades conocidas y recuperar la función inicial mediante la inversa de la transformada de Fourier.
Seaf la función característica de la distribución cuya función de densidad es F, y κr sus cumulantes. Se expande en términos de una distribución conocida con densidad de probabilidad
, función característica
y cumulantes estándarγr. Es común escoger la distribución normal como
, pero también es posible escoger otras. Por definición de los cumulantes se tiene la siguiente identidad:
![{\displaystyle f(t)=\exp \left[\sum _{r=1}^{\infty }(\kappa _{r}-\gamma _{r}){\frac {(it)^{r}}{r!}}\right]\psi (t)\,.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/024a3d7b6213d4771f53bd6e8f03e4a41c6d2104)
Por las propiedades de la transformada de Fourier, (it)rψ(t) es la transformada de Fourier de (−1)r Dr
(x), donde D es el operador diferencial respecto a x. Por tanto, se obtiene para F la expansión
![{\displaystyle F(x)=\exp \left[\sum _{r=1}^{\infty }(\kappa _{r}-\gamma _{r}){\frac {(-D)^{r}}{r!}}\right]\Psi (x)\,.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1e48162ecb9f83b722186f57e309cad9eb8e0265)
Si se elige
como la normal con media y varianza dadas por F, es decir, media μ = κ1 y varianza σ2 = κ2, entonces la expansión es
![{\displaystyle F(x)=\exp \left[\sum _{r=3}^{\infty }\kappa _{r}{\frac {(-D)^{r}}{r!}}\right]{\frac {1}{{\sqrt {2\pi }}\sigma }}\exp \left[-{\frac {(x-\mu )^{2}}{2\sigma ^{2}}}\right]\,.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/47ddea32a50d31462215ec3e3bc502a9bfff6462)
Al expandir el exponencial y agrupando lo términos según el orden de las derivadas se obtiene la serie de Gram–Charlier A. Al incluir solo primeros dos términos de corrección de la normal se obtiene
![{\displaystyle F(x)={\frac {1}{{\sqrt {2\pi }}\sigma }}\exp \left[-{\frac {(x-\mu )^{2}}{2\sigma ^{2}}}\right]\left[1+{\frac {\kappa _{3}}{3!\sigma ^{3}}}H_{3}\left({\frac {x-\mu }{\sigma }}\right)+{\frac {\kappa _{4}}{4!\sigma ^{4}}}H_{4}\left({\frac {x-\mu }{\sigma }}\right)\right]\,,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ca67001f98a9010d8a23f519c5a2cee534ce6f17)
con H3(x) = x3 − 3x y
H4(x) = x4 − 6x2 + 3 (Polinomios de Hermite).
Nótese que esta expresión no es necesariamente positiva y por lo tanto, no es una distribución de probabilidad válida. La serie de Gram–Charlier A diverge en muchos casos. La serie converge solo si F(x) decrece más rápido que exp(−x2/4) hacia infinito (Cramér 1957). Cuando no converge, la serie no puede ser una serie asintótica porque no es posible estimar el error de la expansión. por esta razón, la serie de Edgeworth se prefiere a la de Gram–Charlier A.
Serie de Edgeworth
Edgeworth desarrolló una expansión similar a partir del teorema del límite central. La ventaja de la serie de Edgeworth es que se controla el error, siendo así una verdadera serie asintótica.
Sea {Xi} una secuencia de variables aleatorias independiente e idénticamente distribuidas de media μ y varianza σ2, y sea Yn su suma estandarizada:
![{\displaystyle Y_{n}={\frac {1}{\sqrt {n}}}\sum _{i=1}^{n}{\frac {X_{i}-\mu }{\sigma }}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/51e318ff27a84e2580667df87d2dee59904cc8f7)
Sea Fn la función de distribución acumulada de las variables Yn. Entonces por el teorema del límite central,
![{\displaystyle \lim _{n\to \infty }F_{n}(x)=\Phi (x)\equiv \int _{-\infty }^{x}{\tfrac {1}{\sqrt {2\pi }}}e^{-{\frac {1}{2}}q^{2}}dq}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c44f99c89b9e3647e3ee0b7961a9372fa2bae38d)
para todo x, siempre que la media y la varianza sean finitas.
Ahora asuma que las variables aleatorias Xi tienen media μ, varianza σ2 y cumulantes de mayor orden κr=σrλr. Expandiendo en términos de la distribución normal estándar
![{\displaystyle \Psi (x)={\frac {1}{\sqrt {2\pi }}}\exp(-{\tfrac {1}{2}}x^{2})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4c60054294291691054da4d6550d0ca5b0afadf1)
Entonces las diferencias entre cumulantes en la expresión de la función característica fn(t) de Fn son
![{\displaystyle \kappa _{1}^{F(n)}-\gamma _{1}=0\,,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/21422be24f4b328f21f210df247acda99fa0ecff)
![{\displaystyle \kappa _{2}^{F(n)}-\gamma _{2}=0\,,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/65f55b9060a9017a5285b2f7e61b9ed834f907bf)
![{\displaystyle \kappa _{r}^{F(n)}-\gamma _{r}={\frac {\kappa _{r}}{\sigma ^{r}n^{r/2-1}}}={\frac {\lambda _{r}}{n^{r/2-1}}};\quad r\geq 3\,.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5f3416521b570c342b7e858cb531c1f80ffd4cf4)
La serie de Edgeworth se desarrolla de manera similar a la de Gram–Charlier A series, solo que ahora los términos se asocian según la potencia de n. Así,
![{\displaystyle f_{n}(t)=\left[1+\sum _{j=1}^{\infty }{\frac {P_{j}(it)}{n^{j/2}}}\right]\exp(-t^{2}/2)\,,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5c08972d7bbbad7d52852c12accece1feafc2d29)
donde Pj(x) es un polinomio de grado 3j. Nuevamente, luego de invertir la transformada de Fourier, la función de distribución Fn es
![{\displaystyle F_{n}(x)=\Phi (x)+\sum _{j=1}^{\infty }{\frac {P_{j}(-D)}{n^{j/2}}}\Phi (x)\,.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5b3c392ffd0e72d09bc79334486d09b1b40e8c3e)
Los primeros cinco términos de la expansión son[1]
![{\displaystyle {\begin{aligned}F_{n}(x)=&\Phi (x)\\&-{\frac {1}{n^{1/2}}}{\bigg (}{\tfrac {1}{6}}\lambda _{3}\,\Phi ^{(3)}(x){\bigg )}\\&+{\frac {1}{n}}{\bigg (}{\tfrac {1}{24}}\lambda _{4}\,\Phi ^{(4)}(x)+{\tfrac {1}{72}}\lambda _{3}^{2}\,\Phi ^{(6)}(x){\bigg )}\\&-{\frac {1}{n^{3/2}}}{\bigg (}{\tfrac {1}{120}}\lambda _{5}\,\Phi ^{(5)}(x)+{\tfrac {1}{144}}\lambda _{3}\lambda _{4}\,\Phi ^{(7)}(x)+{\tfrac {1}{1296}}\lambda _{3}^{3}\,\Phi ^{(9)}(x){\bigg )}\\&+{\frac {1}{n^{2}}}{\bigg (}{\tfrac {1}{720}}\lambda _{6}\,\Phi ^{(6)}(x)+{\big (}{\tfrac {1}{1152}}\lambda _{4}^{2}+{\tfrac {1}{720}}\lambda _{3}\lambda _{5}{\big )}\Phi ^{(8)}(x)\\&\qquad \quad +{\tfrac {1}{1728}}\lambda _{3}^{2}\lambda _{4}\,\Phi ^{(10)}(x)+{\tfrac {1}{31104}}\lambda _{3}^{4}\,\Phi ^{(12)}(x){\bigg )}\\&+O(n^{-5/2})\,.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/36203131c39444e6fd040c36a5447a02f869d653)
Acá Φ(j)(x) es la j-ésima derivada de Φ(·) en x. Blinnikov y Moessner (1998) dieron un algoritmo sencillo para calcular los términos de mayor orden de la expansión.
Notas
Bibliografía
- H. Cramér. (1957). Mathematical Methods of Statistics. Princeton University Press, Princeton.
- D. L. Wallace. (1958). "Asymptotic approximations to distributions". Annals of Mathematical Statistics, 29: 635–654.
- M. Kendall & A. Stuart. (1977), The advanced theory of statistics, Vol 1: Distribution theory, 4th Edition, Macmillan, New York.
- P. McCullagh (1987). Tensor Methods in Statistics. Chapman and Hall, London.
- D. R. Cox and O. E. Barndorff-Nielsen (1989). Asymptotic Techniques for Use in Statistics. Chapman and Hall, London.
- P. Hall (1992). The Bootstrap and Edgeworth Expansion. Springer, New York.
- S. Blinnikov and R. Moessner (1998). Expansions for nearly Gaussian distributions. Astronomy and astrophysics Supplement series, 130: 193–205.
- J. E. Kolassa (2006). Series Approximation Methods in Statistics (3rd ed.). (Lecture Notes in Statistics #88). Springer, New York.