Articulo de referencia

Correlograma

Un gráfico que muestra 100 números aleatorios con una función sinusoidal "oculta" y una autocorrelación (correlograma) de la serie en la parte inferior. En el análisis de datos,...

Un gráfico que muestra 100 números aleatorios con una función sinusoidal "oculta" y una autocorrelación (correlograma) de la serie en la parte inferior.

En el análisis de datos, un correlograma es un gráfico de estadísticas de correlación . Por ejemplo, en el análisis de series temporales , un gráfico de las autocorrelaciones de la muestra.rh{\displaystyle r_{h}}versush{\displaystyle h\,}(los desfases temporales) es un autocorrelograma . Si se grafica la correlación cruzada , el resultado se llama un correlograma cruzado .

El correlograma es una herramienta comúnmente utilizada para comprobar la aleatoriedad de un conjunto de datos . Si los datos son aleatorios, las autocorrelaciones deberían ser cercanas a cero para cualquier intervalo de tiempo. Si no son aleatorios, una o más autocorrelaciones serán significativamente distintas de cero.

Además, los correlogramas se utilizan en la etapa de identificación del modelo para los modelos de series temporales autorregresivas de media móvil de Box-Jenkins . Las autocorrelaciones deben ser cercanas a cero para indicar aleatoriedad; si el analista no verifica la aleatoriedad, la validez de muchas de las conclusiones estadísticas se vuelve cuestionable. El correlograma es una excelente manera de verificar dicha aleatoriedad.

En el análisis multivariado , las matrices de correlación mostradas como imágenes con mapas de color también pueden denominarse "correlogramas" o "corrgramas". [ 1 ] [ 2 ] [ 3 ]

Aplicaciones

El correlograma puede ayudar a proporcionar respuestas a las siguientes preguntas: [ 4 ]

  • ¿Son aleatorios los datos?
  • ¿Está una observación relacionada con una observación adyacente?
  • ¿Una observación está relacionada con una observación dos veces eliminada? (etc.)
  • ¿La serie temporal observada es ruido blanco ?
  • ¿La serie temporal observada es sinusoidal?
  • ¿La serie temporal observada es autorregresiva?
  • ¿Cuál es el modelo apropiado para la serie temporal observada?
  • ¿Es el modelo?
Y=constante+error{\displaystyle Y={\text{constante}}+{\text{error}}}
¿Válido y suficiente?
  • ¿Es la fórmula?sY¯=s/norte{\displaystyle s_{\bar {Y}}=s/{\sqrt {N}}}¿válido?

Importancia

La aleatoriedad (junto con el modelo fijo, la variación fija y la distribución fija) es uno de los cuatro supuestos que suelen subyacer a todos los procesos de medición. El supuesto de aleatoriedad es de vital importancia por las siguientes tres razones:

  • La mayoría de las pruebas estadísticas estándar dependen del azar. La validez de las conclusiones de la prueba está directamente relacionada con la validez del supuesto de azar.
  • Muchas fórmulas estadísticas de uso común dependen del supuesto de aleatoriedad, siendo la más común la fórmula para determinar el error estándar de la media muestral:
sY¯=s/norte{\displaystyle s_{\bar {Y}}=s/{\sqrt {N}}}

donde s es la desviación estándar de los datos. Si bien se usa con frecuencia, los resultados obtenidos con esta fórmula no tienen valor a menos que se cumpla el supuesto de aleatoriedad.

  • Para datos univariados, el modelo predeterminado es
Y=constante+error{\displaystyle Y={\text{constante}}+{\text{error}}}

Si los datos no son aleatorios, este modelo es incorrecto e inválido, y las estimaciones de los parámetros (como la constante) se vuelven absurdas e inválidas.

Estimación de autocorrelaciones

El coeficiente de autocorrelación en el retardo h viene dado por

rh=doh/do0{\displaystyle r_{h}=c_{h}/c_{0}\,}

donde c h es la función de autocovarianza

doh=1nortet=1norteh(YtY¯)(Yt+hY¯){\displaystyle c_{h}={\frac {1}{N}}\sum _{t=1}^{Nh}\left(Y_{t}-{\bar {Y}}\right)\left(Y_{t+h}-{\bar {Y}}\right)}

y c 0 es la función de varianza

do0=1nortet=1norte(YtY¯)2{\displaystyle c_{0}={\frac {1}{N}}\sum _{t=1}^{N}\left(Y_{t}-{\bar {Y}}\right)^{2}}

El valor resultante de r h estará comprendido entre −1 y  +1.

Estimación alternativa

Algunas fuentes pueden utilizar la siguiente fórmula para la función de autocovarianza:

doh=1norteht=1norteh(YtY¯)(Yt+hY¯){\displaystyle c_{h}={\frac {1}{Nh}}\sum _{t=1}^{Nh}\left(Y_{t}-{\bar {Y}}\right)\left(Y_{t+h}-{\bar {Y}}\right)}

Si bien esta definición presenta menos sesgos , la formulación (1/ N ) posee algunas propiedades estadísticas deseables y es la forma más utilizada en la literatura estadística. Para más detalles, consulte las páginas 20 y 49-50 de Chatfield.

A diferencia de la definición anterior, esta definición nos permite calculardoh{\displaystyle c_{h}}de una manera un poco más intuitiva. Consideremos el ejemploY1,,Ynorte{\ Displaystyle Y_ {1}, \ puntos, Y_ {N}}, dóndeYiRnorte{\displaystyle Y_{i}\in \mathbb {R} ^{n}}parai=1,,norte{\displaystyle i=1,\dots ,N}. Entonces, deja

incógnita=[Y1Y¯YnorteY¯]Rnorte×norte{\displaystyle X={\begin{bmatrix}Y_{1}-{\bar {Y}}&\cdots &Y_{N}-{\bar {Y}}\end{bmatrix}}\in \mathbb {R} ^{n\times N}}

Luego calculamos la matriz de Gram.Q=incógnitaincógnita{\displaystyle Q=X^{\top }X}. Finalmente,doh{\displaystyle c_{h}}se calcula como la media muestral de lah{\displaystyle h}diagonal deQ{\displaystyle Q}. Por ejemplo, el0{\displaystyle 0}diagonal (la diagonal principal) deQ{\displaystyle Q}tienenorte{\displaystyle N}elementos, y su media muestral corresponde ado0{\displaystyle c_{0}}. El1{\displaystyle 1}primera diagonal (a la derecha de la diagonal principal) deQ{\displaystyle Q}tienenorte1{\displaystyle N-1}elementos, y su media muestral corresponde ado1{\displaystyle c_{1}}, etcétera.

Inferencia estadística con correlogramas

Ejemplo de correlograma a partir de una muestra de 400 puntos de un proceso autorregresivo de primer orden con una correlación de 0,75 entre puntos adyacentes, junto con los intervalos de confianza del 95 % (representados en negro alrededor de las estimaciones de correlación y en rojo alrededor de cero), calculados mediante las ecuaciones de esta sección. La línea azul discontinua muestra la función de autocorrelación real del proceso muestreado.
20 correlogramas a partir de muestras de 400 puntos del mismo proceso aleatorio que en la figura anterior.

En el mismo gráfico se pueden trazar límites superiores e inferiores para la autocorrelación con nivel de significancia.α{\displaystyle \alpha \,}:

B=±z1α/2Smi(rh){\displaystyle B=\pm z_{1-\alpha /2}SE(r_{h})\,}conrh{\displaystyle r_{h}\,}como la autocorrelación estimada en el retardoh{\displaystyle h\,}.

Si la autocorrelación es mayor (menor) que este límite superior (inferior), la hipótesis nula de que no hay autocorrelación en y más allá de un retardo dado se rechaza a un nivel de significancia de α{\displaystyle \alpha \,}Esta prueba es aproximada y supone que la serie temporal es gaussiana .

En lo anterior, z 1− α /2 es el cuantil de la distribución normal ; SE es el error estándar, que se puede calcular mediante la fórmula de Bartlett para procesos MA( ℓ ):

Smi(r1)=1norte{\displaystyle SE(r_{1})={\frac {1}{\sqrt {N}}}}
Smi(rh)=1+2i=1h1ri2norte{\displaystyle SE(r_{h})={\sqrt {\frac {1+2\sum _{i=1}^{h-1}r_{i}^{2}}{N}}}}parah>1.{\displaystyle h>1.\,}

En el ejemplo representado, podemos rechazar la hipótesis nula de que no existe autocorrelación entre puntos temporales separados por desfases de hasta 4. Para la mayoría de los periodos más largos, no se puede rechazar la hipótesis nula de que no existe autocorrelación.

Tenga en cuenta que existen dos fórmulas distintas para generar los intervalos de confianza:

1. Si el correlograma se utiliza para comprobar la aleatoriedad (es decir, que no existe dependencia temporal en los datos), se recomienda la siguiente fórmula:

±z1α/2norte{\displaystyle \pm {\frac {z_{1-\alpha /2}}{\sqrt {N}}}}

donde N es el tamaño de la muestra , z es la función cuantil de la distribución normal estándar y α es el nivel de significancia . En este caso, las bandas de confianza tienen un ancho fijo que depende del tamaño de la muestra.

2. Los correlogramas también se utilizan en la etapa de identificación del modelo para ajustar modelos ARIMA . En este caso, se asume un modelo de media móvil para los datos y se deben generar las siguientes bandas de confianza:

±z1α/21norte(1+2i=1kri2){\displaystyle \pm z_{1-\alpha /2}{\sqrt {{\frac {1}{N}}\left(1+2\sum _{i=1}^{k}r_{i}^{2}\right)}}}

donde k es el retardo. En este caso, los intervalos de confianza aumentan a medida que aumenta el retardo.

Software

Los correlogramas están disponibles en la mayoría de las bibliotecas estadísticas de uso general.

Correlogramas:

Corrgramas:

Referencias

  1. Friendly, Michael (19 de agosto de 2002). "Corrgramas: representaciones exploratorias para matrices de correlación" (PDF) . The American Statistician . 56 (4). Taylor & Francis : 316–324 . doi : 10.1198/000313002533 . Consultado el 19 de enero de 2014 .
  2. 1 2 "CRAN – Paquete corrgram" . cran.r-project.org . 29 de agosto de 2013. Consultado el 19 de enero de 2014 .
  3. 1 2 "Quick-R: Correlogramas" . statmethods.net . Consultado el 19 de enero de 2014 .
  4. "1.3.3.1. Gráfico de autocorrelación" . www.itl.nist.gov . Consultado el 20 de agosto de 2018 .
  5. "Visualización § Gráfico de autocorrelación" .

Lecturas adicionales

  • Hanke, John E.; Reitsch, Arthur G.; Wichern, Dean W. Previsión empresarial (7.ª  ed.). Upper Saddle River, NJ: Prentice Hall.
  • Box, GEP; Jenkins, G. (1976). Análisis de series temporales: pronóstico y control . Holden-Day.
  • Chatfield, C. (1989). El análisis de series temporales: una introducción (cuarta  ed.). Nueva York, NY: Chapman & Hall.
  • Gráfico de autocorrelación

Dominio público Este artículo incorpora material de dominio público del Instituto Nacional de Estándares y Tecnología.