Una prueba de la estacionariedad de un AR (2)


17

Considere un proceso AR (2) centrado en la media

Xt=ϕ1Xt1+ϕ2Xt2+ϵt
donde ϵt es el proceso estándar de ruido blanco. Sólo por razones de simplicidad déjame llamar ϕ1=b y ϕ2=a . Centrándome en las raíces de la ecuación de características obtuve
z1,2=b±b2+4a2a
Las condiciones clásicas en los libros de texto son las siguientes: Traté de resolver manualmente (con la ayuda de Mathematica) las desigualdades en las raíces, es decir, el sistemaobteniendo solo¿Sepuede recuperarla tercera condición () agregando las dos soluciones anteriores entre sí obteniendoque a través de algunas consideraciones de signos se convierte en? ¿O me estoy perdiendo una solución?
{|a|<1a±b<1
{|bb2+4a2a|>1|b+b2+4a2a|>1
a±b<1
|a|<1a+b+ab<2a<1|a|<1

Respuestas:


18

Mi conjetura es que la ecuación característica de la que estás saliendo es diferente de la mía. Permítanme proceder en un par de pasos para ver si estamos de acuerdo.

Considere la ecuación

λ2ϕ1λϕ2=0

Si z es una raíz de la ecuación característica "estándar" 1ϕ1zϕ2z2=0 y configurando z1=λ , la pantalla se obtiene reescribiendo la estándar de la siguiente manera:

1ϕ1zϕ2z2=0z2ϕ1z1ϕ2=0λ2ϕ1λϕ2=0
Por lo tanto, una condición alternativa para la estabilidad de unAR(2)es que todas las raíces de la primera pantalla estándentrodel círculo unitario,|z|>1|λ|=|z1|<1.

Utilizamos esta representación para derivar el triángulo de estacionariedad de un proceso AR(2) , es decir, un AR(2) es estable si se cumplen las siguientes tres condiciones:

  1. ϕ2<1+ϕ1
  2. ϕ2<1ϕ1
  3. ϕ2>1

Recuerde que puede escribir las raíces de la primera pantalla (si es real) como

λ1,2=ϕ1±ϕ12+4ϕ22
para encontrar las dos primeras condiciones.

Entonces, el AR(2) es estacionario iff |λ|<1 , por lo tanto (si el λi es real):

1<ϕ1±ϕ12+4ϕ22<12<ϕ1±ϕ12+4ϕ2<2
La mayor de las dosλiestá limitada porϕ1+ϕ12+4ϕ2<2, o:
ϕ1+ϕ12+4ϕ2<2ϕ12+4ϕ2<2ϕ1ϕ12+4ϕ2<(2ϕ1)2ϕ12+4ϕ2<44ϕ1+ϕ12ϕ2<1ϕ1
Analogously, we find that ϕ2<1+ϕ1.

If λi is complex, then ϕ12<4ϕ2 and so

λ1,2=ϕ1/2±i(ϕ12+4ϕ2)/2.
The squared modulus of a complex number is the square of the real plus the square of the imaginary part. Hence,
λ2=(ϕ1/2)2+((ϕ12+4ϕ2)/2)2=ϕ12/4(ϕ12+4ϕ2)/4=ϕ2.
This is stable if |λ|<1, hence if ϕ2<1 or ϕ2>1, as was to be shown. (The restriction ϕ2<1 resulting from ϕ22<1 is redundant in view of ϕ2<1+ϕ1 and ϕ2<1ϕ1.)

Plotting the stationarity triangle, also indicating the line that separates complex from real roots, we get

enter image description here

Produced in R using

phi1 <- seq(from = -2.5, to = 2.5, length = 51) 
plot(phi1,1+phi1,lty="dashed",type="l",xlab="",ylab="",cex.axis=.8,ylim=c(-1.5,1.5))
abline(a = -1, b = 0, lty="dashed")
abline(a = 1, b = -1, lty="dashed")
title(ylab=expression(phi[2]),xlab=expression(phi[1]),cex.lab=.8)
polygon(x = phi1[6:46], y = 1-abs(phi1[6:46]), col="gray")
lines(phi1,-phi1^2/4)
text(0,-.5,expression(phi[2]<phi[1]^2/4),cex=.7)
text(1.2,.5,expression(phi[2]>1-phi[1]),cex=.7)
text(-1.75,.5,expression(phi[2]>1+phi[1]),cex=.7)

this is a very detailed explanation.
Marco

@Christoph: Is there a typo in the answer? Look at equation for λ2. Also, what do you mean by square of a complex number? If z=a+bi then z2=a2b2+2iab. How do you say square of a complex number is "square of the real plus the square of the imaginary part"
shani

1
Thanks, quite right! I was referring to the sqaured modulus, see the edit.
Christoph Hanck

@ChristophHanck, what is your take on Aksakal's answers in these two threads: 1 and 2? Are they in conflict with your answer, and if so, what is the correct answer?
Richard Hardy

I think he is quite right when defining weak stationarity as constancy of the first two moments. Often, and also in the present thread, "stationarity" and "existence of a causal representation", i.e., a summable MA() representation without dependence on the future, are conflated. What my answer therefore more precisely shows is conditions for the existence of the latter.
Christoph Hanck
Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.