He visto varias veces que las personas rechazan el nulo en una prueba aumentada de Dickey-Fuller , y luego afirman que muestra que su serie es estacionaria (desafortunadamente, no puedo mostrar las fuentes de estas afirmaciones, pero imagino que existen afirmaciones similares aquí y allá en uno u otro diario).
Sostengo que es un malentendido (que el rechazo del nulo de una raíz unitaria no es necesariamente lo mismo que tener una serie estacionaria, especialmente porque las formas alternativas de no estacionariedad rara vez se investigan o incluso se consideran cuando se realizan tales pruebas).
Lo que busco es:
a) un buen contraejemplo claro al reclamo (puedo imaginar un par en este momento, pero apuesto a que alguien que no sea yo tendrá algo mejor de lo que tengo en mente). Podría ser una descripción de una situación específica, tal vez con datos (simulados o reales; ambos tienen sus ventajas); o
b) un argumento convincente de por qué el rechazo en un Dickey-Fuller aumentado debería considerarse como un factor estacionario
(o incluso ambos (a) y (b) si te sientes inteligente)
fuente
Respuestas:
Aquí hay un ejemplo de una serie no estacionaria que ni siquiera una prueba de ruido blanco puede detectar (y mucho menos una prueba de tipo Dickey-Fuller):
Sí, esto puede ser sorprendente, pero esto no es ruido blanco .
La mayoría de los ejemplos de contadores no estacionarios se basan en una violación de las dos primeras condiciones estacionarias: tendencias deterministas (media no constante) o series temporales de raíz unitaria / heteroscedastica (varianza no constante). Sin embargo, también puede tener procesos no estacionarios que tienen media constante y varianza, pero violan la tercera condición: la función de autocovarianza (ACVF) debe ser constante en el tiempo y una función desolamente.| s - t |c o v ( xs, xt) El | s-t |
La serie de tiempo anterior es un ejemplo de una serie de este tipo, que tiene media cero, varianza unitaria, pero el ACVF depende del tiempo. Más precisamente, el proceso anterior es un proceso MA (1) localmente estacionario con parámetros tales que se convierte en ruido blanco espurio (ver Referencias a continuación): el parámetro del proceso MA cambia con el tiempoXt= εt+ θ1εt - 1
donde es el tiempo normalizado. La razón por la que esto parece ruido blanco (aunque, por definición matemática, claramente no lo es), es que el ACVF variable en el tiempo se integra a cero con el tiempo. Dado que el ACVF de muestra converge al ACVF promedio, esto significa que la autocovarianza de muestra (y la autocorrelación (ACF)) convergerán en una función que parece ruido blanco. Entonces, incluso una prueba de Ljung-Box no podrá detectar esta no estacionariedad. El documento (descargo de responsabilidad: soy el autor) sobre Pruebas de ruido blanco contra alternativas localmente estacionarias propone una extensión de las pruebas de Box para hacer frente a dichos procesos localmente estacionarios.u = t / T
Para obtener más código R y más detalles, consulte también esta publicación de blog .
Actualización después del comentario de mpiktas :
Es cierto que esto podría parecer un caso teóricamente interesante que no se ve en la práctica. Estoy de acuerdo en que es poco probable que vea un ruido blanco espurio en un conjunto de datos del mundo real directamente, pero lo verá en casi cualquier residuo de un ajuste de modelo estacionario. Sin entrar en demasiados detalles teóricos, imagine un modelo general variable en el tiempo con una función de covarianza variable en el tiempo . Si ajusta un modelo constante , esta estimación estará cerca del promedio de tiempo del modelo verdadero ; y, naturalmente, los residuos ahora estarán cerca de , que por construcción deθ ( u ) γθ( k , u ) θˆ θ ( u ) θ ( u ) - θˆ θˆ se integrará a cero (aproximadamente). Ver Goerg (2012) para más detalles.
Veamos un ejemplo
Entonces ruido fraccional con el parámetro (dado que creemos que todo está bien y tenemos un modelo estacionario). Verifiquemos los residuos: d <0,5reˆ= 0.23 reˆ< 0.5
Se ve bien, ¿verdad? Bueno, el problema es que los residuos son ruido blanco espurio . ¿Cómo puedo saber? Primero puedo probarlo
y segundo, sabemos por la literatura que los datos del anillo de los árboles son, de hecho, ruido fraccionado localmente estacionario: ver Goerg (2012) y Ferreira, Olea y Palma (2013) .
Esto muestra que mi, ciertamente, un ejemplo teórico, en realidad está ocurriendo en la mayoría de los ejemplos del mundo real.
fuente
Ejemplo 1
Se sabe que los procesos de raíz unitaria con un componente MA negativo fuerte conducen a pruebas de ADF con un tamaño empírico muy superior al nominal (por ejemplo, Schwert, JBES 1989 ).
Ejemplo 2
Dependiendo del tipo de cambio de varianza, la prueba ADF aún se rechazará con frecuencia. En mi ejemplo a continuación, tenemos una ruptura de varianza hacia abajo, lo que hace que la prueba "crea" que la serie converge, lo que lleva a un rechazo del nulo de una raíz unitaria.
(Como comentario aparte, la prueba ADF "pierde" su distribución nula asintótica fundamental en presencia de heterocedasticidad incondicional).
fuente
La prueba de raíz unitaria es notoriamente difícil. Por lo general, usar una prueba no es suficiente y debe tener mucho cuidado con los supuestos exactos que usa la prueba.
La forma en que se construye el ADF lo hace vulnerable a una serie que son tendencias simples no lineales con ruido blanco agregado. Aquí hay un ejemplo:
Aquí tenemos la tendencia exponencial y vemos que el ADF funciona bastante mal. Acepta el nulo de raíz unitaria el 30% del tiempo y lo rechaza el 70% del tiempo.
Por lo general, el resultado de cualquier análisis no es afirmar que la serie es estacionaria o no. Si los métodos utilizados en el análisis requieren estacionariedad, la suposición errónea de que la serie es estacionaria cuando en realidad no lo es, generalmente se manifiesta de una forma u otra. Así que personalmente miro todo el análisis, no solo la parte de prueba de raíz unitaria. Por ejemplo, OLS y NLS funcionan bien para datos no estacionarios, donde la no estacionariedad está en la media, es decir, la tendencia. Entonces, si alguien afirma erróneamente que la serie es estacionaria y aplica OLS / NLS, esta afirmación podría no ser relevante.
fuente
sapply(oo, "[[","p.value")
?