Serie Temporal Uca
Serie Temporal Uca
Serie Temporal Uca
SERIES DE TIEMPO
Aspectos básicos y fundamentales
• El operador de rezagos
2
Series Temporales: breve historia
• A inicios del siglo XX, Slutzky y Yule mostraron que series de tiempo
con propiedades similares a series económicas pueden generarse
como sumas o restas (simples o ponderadas) de procesos puramente
aleatorios.
3
Series Temporales: breve historia
4
Componentes de una serie temporal
5
Representaciones gráficas y transformaciones
de una serie temporal
• Primera diferencia: Elimina la tendencia, mantiene los ciclos
estacionales.
• Tasa de crecimiento %: Elimina la tendencia, mantiene los ciclos
estacionales. (Tiende a generar mucha asimetría en los datos o alta
variabilidad)
• Tasa contínua de crecimiento %: Idem caso anterior con la diferencia
que homogeniza más los datos (gracias a los logaritmos).
• Diferencia interanual: Elimina la tendencia y todo factor estacional
• Tasa de crecimiento interanual %: Idem caso anterior, con la
salvedad que tenemos incrementos porcentuales.
• Serie suavizada por media móvil: Elimina la estacionalidad, pero
mantiene la tendencia.
6
Serie Original: PIB trimestral paraguayo en
millones de Gs. costantes
60000000
40000000
pib/t
20000000
0
7
Serie diferenciada: ∆yt = yt − yt−1
2000000 4000000
0
dPIB/t
-6000000-4000000-2000000
8
Serie con tasa de crecimiento ∆%y = yt– yt-1/yt-1 *100
20
10
tasa_pib
0
-10
-20
9
Serie con tasa de crecimiento contínuo
20
10
(diferencia de logaritmos) ∆%yt ≈ (ln yt − ln yt−1) × 100
tasa_pib
0
-10
-20
10
Serie con diferencia interanual o diferencia
estacional ∆4yt ≡ yt − yt−4
4000000
2000000
pib, S4
0
-2000000
11
Serie con tasa de crecimiento interanual
(trimestre actual con mismo trimestre del
periodo anterior) ∆4%yt ≈ (ln yt − ln yt−4) × 100
15
10 5
tasa_4_pib
0 -5
-10
12
pib/ma: x(t)= PIB: window(3 1 0)
20000000 30000000 40000000 50000000 60000000
1995q1
2000q1
2005q1
t
2010q1
Serie suavizada con media móvil
2015q1
13
2020q1
La transformación de datos en STATA para el
PIB trimestral de Paraguay
Transformación Comando en STATA
Serie original PIB
14
Las series de tiempo y
el uso
de
los rezagos
15
El papel del “tiempo” o del “rezago” en los
modelos económicos
• Con mucha frecuencia una serie temporal Y, responde ante un
cambio en X de manera lenta y no instantánea, ese lapso se
denomina rezago.
• Por ejemplo en un modelo de gastos de consumo si se produce un
aumento permanente de 2000 US$, la persona no gastará los 2000
US$ de forma rápida, el impacto del ingreso puede transladarse
varios periodos.
• Yt = constante + 0.4Xt + 0.3Xt−1 + 0.2Xt−2 + ut
• El impacto comporáneo ante un aumento del ingreso de 2000 es un
aumento del gasto de consumo en 800 US$, transcurrido un año, el
aumento inicial de 2000, tiene un impacto de 600 US$ en el gasto de
consumo, y luego de dos años 400 US$, en consecuencia el efecto de
largo plazo es 1800 US$.
• El impacto del “shock” se propaga durante dos años , en
consecuencia el ”shock” tiene un impacto tanto instantáneo y
rezagado.
16
El operador de rezagos
• Un operador de una serie temporal es un “proceso” que transforma una
o más series en nuevas series de tiempo.
• El operador de rezago se denota por L y se define como
• Lxt ≡ xt−1
• L2xt = xt−2
• L3xt = xt−3
• L-1xt = xt+1
17
Ejemplos de aplicación del operador de
rezagos
18
Inverso de un polinomio de grado 1
19
Inverso de un polinomio de grado p
20
El operador de rezagos, ecuaciones en
diferencias y estabilidad de un proceso
21
Operador de rezagos y estabilidad de un
proceso de ecuaciones en diferencias de
orden p
23
24
En los softwares estadísticos…..
Roots of the companion matrix
1
.5
Imaginary
0-.5
-1
-1 -.5 0 .5 1
Real
GRETL
25
Aspectos Generales
de Series
Temporales
26
Series temporales
yt = b0 + b1xt1 + . . .+ bkxtk + ut
27
Características de las series temporales
28
Ejemplos de modelos de series temporales
• Modelo estático de variables con relación contemporánea:
yt = b0 + b1zt + ut
29
Modelo de Rezagos Distribuidos Finitos (RDF)
30
Ejemplo: un modelo dinámico de RDF
• intt = 1.6 + 0.48inft - 0.15inft-1 + 0.32inft-2 +ut
Donde int es la tasa de interés e inf es la tasa de inflación donde
ambos son datos anuales.
31
Series temporales: supuestos para un
estimador insesgado
32
Varianza del estimador BMCO
4. Var(ut|X) = Var(ut) = s2
La varianza del término de error es constante en el tiempo y no
depende de las x.
5. Ausencia de correlación serial:
Corr(ut,us| X)=0 for t ¹ s
Los errores en cualquier momento del tiempo no deben estar
correlacioados
33
Tendencia en una serie temporal
34
Tendencia en una serie temporal
• Tendencia lineal
yt = a0 + a1t + et, t = 1, 2, …
• Tendencia exponencial
log(yt) = a0 + a1t + et, t = 1, 2, …
• Tendencia cuadrática
yt = a0 + a1t + a2t2 + et, t = 1, 2, …
35
Uso de una variable de tendencia en un
modelo de regresión
• Importante: tener en cuenta factores no observados con tendencia
que afectan a la variable endógena pueden estar correlacionados
con nuestras variables explicativas.
36
Ejemplo de una regresion con tendencia
Inversión en viviendas Índice de precios de vivienda
1.05
80000 100000 120000 140000 160000
1
housing price index; 1982 = 1
real housing inv, millions $
.85 .9 .95
60000
.8
1950 1960 1970 1980 1990 1950 1960 1970 1980 1990
1947-1988 1947-1988
37
Ejemplo de una regresion con tendencia
linv Coef. Std. Err. t P>|t| [95% Conf. Interval]
• Para poder analizar la evolución real de una serie se debe ajustar por el factor
estacional. Es decir, se debe trabajar con una serie desestacionalizada.
39
Estacionalidad ejemplos
40
Capturando el efecto estacional en un
modelo de regresión
Ejemplo de cómo capturar alguna influencia estacional
• La ecuación del modelo es suponer que las ventas están en función a los
gastos en publidad y variables dummies estacionales para los trimestres 2,3,y
4, tomando el primer trimestre como referencia.
41
Recordar….
Componente de las series temporales
Yt = Tt + St + It
• Valor observado = Tendencia + Estacionalidad + Irregular
42
INTRODUCCION
PROCESOS
ESTOCÁSTICOS
43
Proceso estocástico estacionario
44
Proceso estocástico estacionario en
Covarianza
• Un proceso estocástico estacionario en covarianza es si E(xt) es constante, Var(xt)
es constante y para todo t, h ≥ 1, Cov(xt, xt+h) depende solo de h y no de t.
• E(xt)=x0
• VAR(xt)=
• Cov(xt, xt+h)=t-h (rezago)
• Un proceso estocástico estacionario en covarianza significa que los segundos
momentos son finitos (constantes) y el valor de la covarianza entre dos periodos
depende solamente de la distancia o rezago entre estos dos periodos de tiempo,
y no del tiempo en el cual se ha calculado la covarianza.
• Al proceso estocástico en covarianza se le conoce como debilmente estacionario,
estacionario de Segundo orden o proceso estocástico estacionario en sentido
amplio.
45
Importancia de la estacionariedad en el
análisis econométrico
• Si una serie no es estacionaria, no se puede analizar su comportamiento
(la media y la varianza serían infinitas).
• Las series económicas que presentan una marcada tendencia creciente o
alta variabilidad no son estacionarias.
• La estacionariedad implica que la variable retorne a una senda estable
en el tiempo ante cualquier shock (reversión a la media).
• Para fines de pronóstico es importante trabajar con series estacionarias
de manera a predecir su comportamiento.
• La estacionariedad require que la varianza de la variable sea constante.
Implica incertidumbre acotada y proceso no explosivo.
46
Modelos Autoregresivos de Media Móvil
(ARMA)
• Recordar que una serie de tiempo {yt}Tt=1 es simplemente la
realización de un proceso estocástico.
• Supongamos que tenemos una serie de tiempo en donde tenemos
las realizaciones del proceso (un valor por período).
y1, y2,.....yt que son los datos observados y queremos proyectar valores
futuros yt+1, yt+2,.....yt+k
47
Tipos de modelos para
especificaciones estacionarias
• Todas las series de tiempo que son estacionarias se pueden modelizar a través de una
estructura tipo ARMA.
• Proceso Autorregresivo (AR): proceso que se explica donde la serie depende sus valores
pasados.
• Proceso de Medias Móviles (MA): proceso que se explica donde la serie depende de sus
innovaciones pasadas (shocks aleatorios = ruido blanco)
• Combinación tanto de procesos autorregresivos y medias móviles (ARMA): proceso que
se explican tanto porque dependen de los valores pasados propios de la serie y de las
innovaciones pasadas.
48
Proceso MA(1) de medias móviles
• Un proceso de medias móviles de orden 1 [MA(1)] se caracteriza por la
siguiente relación:
yt = ϵt + a1ϵ t-1. , t = 1, 2, …
49
Proceso MA(1) de medias móviles
50
Proceso MA(1) de medias móviles
51
Proceso MA(1) de medias móviles
52
Proceso MA(1): Resumen de propiedades del
proceso
Ninguno de los momentos muestrales dependen del tiempo (t), por lo que un
proceso MA(1) siempre será estacionario.
53
Gráfico de autocorrelación de un proceso
MA(1)
Correlograma de un MA(1)
54
Invertibilidad de un proceso MA(1)
• Una condición que hay que verificar es la invertibilidad de un proceso
MA(1).
• Esta condición permitirá que un MA pueda representarse como un
proceso autoregresivo infinito y convergente.
• Convergencia significa que los coeficientes del AR decrecen a cero
conforme avanza el tiempo.
• La idea de la invertibilidad de un MA es que las innovaciones/ruido
blanco pueden invertirse a representaciones de observaciones
pasadas.
• Esto es necesario si queremos pronosticar o modelar una serie
temporal, ya que si el proceso no es invertible las innovaciones del
pasado no podrán estimados, en consecuencia la serie no se podrá
utilizar para fines de pronóstico
• Que un MA pueda representarse como un AR significa que la serie
podrá estimarse
55
Invertibilidad de un proceso MA(1)
56
Proceso MA(q) de orden superior
57
Proceso MA(q) de orden superior
58
Invertibilidad de un Proceso MA(q)
60
Ejemplo: proceso MA(4)
61
Proceso Autorregresivo de Orden 1 AR(1)
62
Proceso Autorregresivo de Orden 1 AR(1)
63
Momentos muestrales de un AR(1)
64
Momentos muestrales de un AR(1)
65
66
0.80 0.60
Autocorrelations of y
0.00 0.20 0.40
-0.20
0 10 20 30
Lag
Bartlett's formula for MA(q) 95% confidence bands
68
Proceso AR(p)
69
Ejemplo: AR(2)
70
¿Es estacionario?
71
Procesos ARMA (Autorregresivo de Media
Móvil)
72
Procesos ARMA (Autorregresivo de Media
Móvil)
73
Identificación de un proceso ARMA
74
Proceso AR(1) en el correlograma
0.50
0.50
0.40
0.40
Partial autocorrelations of X
Autocorrelations of X
0.30
0.30
0.20
0.20
0.10
0.10
0.00
0.00
0 10 20 30 40 0 10 20 30 40
Lag Lag
Bartlett's formula for MA(q) 95% confidence bands 95% Confidence bands [se = 1/sqrt(n)]
75
Proceso MA(1) en el correlograma
76
Ejemplo de un proceso ARMA
77
Autocorrelations of arma11
-0.20 -0.10 0.00 0.10 0.20 0.30
0
5
10
Lag
15
0
5
10
Lag
95% Confidence bands [se = 1/sqrt(n)]
15
20
78
Proceso estacionario con tendencia
• Una serie con tendencia puede no ser estacionaria ya que a medida
que avanza el tiempo la media cambia y puede no ser constante.
• Sin embargo, puede ser estacionaria alrededor de su tendencia en el
tiempo y ser débilmente dependiente (es decir xt y xt-1 se vuelven
independientes cuando h tiende al infinito, Cov [yt , yt+h] → 0,
cuando h → ∞.
• Una serie con tendencia y débilmente dependiente se lo llama serie
estacionaria en tendencia.
79
PROCESOS INTEGRADOS
NO ESTACIONARIOS
MODELO AUTORREGRESIVO
INTEGRADO DE MEDIA MÓVIL
(ARIMA)
80
Procesos Estocásticos No estacionarios
Autorregresivo Integrado de Media Móvil
(ARIMA)
• Hasta ahora hemos visto que podemos modelizar una serie de
tiempo siempre que sea estacionaria.
• Lastimosamente, en el mundo real, pocas series son estacionarias. La
mayoría tienen tendencias estocásticas.
81
Procesos estocásticos no estacionarios:
Caminata aleatoria
82
Procesos estocásticos no estacionarios:
Caminata aleatoria
83
Procesos estocásticos no estacionarios:
Caminata Aleatoria
84
Procesos estocásticos no estacionarios:
Caminata aleatoria con deriva
+ Y0
85
Transformación de series temporales
altamente persistentes
86
Transformación de series temporales
altamente persistentes
87
Transformación de series temporales
altamente persistentes
88
89
190
91
Pronósticos con ARIMA y la metodología Box-
Jenkins
92
Prueba de Raíces Unitarias
Regresiones Espúreas
y
Cointegración
93
Prueba de Raiz Unitaria: Test de Dicky Fuller
94
Prueba de Raíz Unitaria: Test de Dicky Fuller
• La prueba de Dicky - Fuller para comprobar si una serie tiene raíz
unitaria sigue una distribución igual al estadístico t pero se
denomina“ tau”.
20
D.r3 Coef. Std. Err. t P>|t| [95% Conf. Interval]
5 100
0 50 100 150
t
96
Prueba de Raiz Unitaria: Test de Dicky –Fuller
Aumentado - ADF
97
Ejemplo 2: Test de Dicky Fuller
Aumentado(ADF) para la inflación de USA
(Wooldrige, pág 643)
inf
LD. .1383615 .1264026 1.09 0.280 -.1163861 .3931091
L1. -.3103252 .1027077 -3.02 0.004 -.517319 -.1033315
• Una regresion espuria también implica que los shocks no se desvanecen (no son ruido
blanco) y persisten en el tiempo.
• Cuando trabajamos con series de tiempo se debe tener cuidado en inferir resultados
cuando las series son I(1) y tengan tendencias estocásticas.
• Para trabajar con los datos de series de tiempo se deben trabajar con series estacionarias.
99
Cointegración
100
Cointegración: Prueba de Engle-
Granger/Prueba de estacionariedad
• El procedimiento de aplicar la prueba de DFA a los residuos de la
regresión se denomina Prueba de Engle-Granger.
• Los valores asintóticos críticos al 1%, 5% y 10% corresponden a la
tabla de Davidson y Mackinnon del estadístico t (tau) .
• H0: Las series no están cointegradas
• H1: Las series están cointegradas
• El rechazar la hipótesis nula implica que las variables analizadas están
cointegradas, comparten una tendencia en común a largo plazo y
nuestra regresión no es espuria.
• La cointegración es una prueba aplicada sólo cuando trabajamos con
datos no estacionarios. Si las series temporales son no estacionarias
se debe chequear si existe una relación estable entre ellas
(cointegración).
101
Requisitos para la existencia de cointegración
Recordar: cointegración implica relación estable entre las variables de la
regresión estimada
Inconvenientes: Una serie puede ser I(0) y otra I(1) o pueden tener
diferentes orden de integración pero al menos para que exista
cointegración y se descarte la posibilidad de una regresión espúrea, el
residuo de la regresión debe ser estacionario.
102
Ejemplo 3: Cointegración entre r6 y r3, tasas
de los certificados del tesoro de USA a 6 y 3
meses
Realizado la regresión de r6 sobre r3
e
L1. -.6320678 .1120843 -5.64 0.000 -.8540059 -.4101296
LD. -.1460544 .0914803 -1.60 0.113 -.3271945 .0350858
Interpolated Dickey-Fuller
Test 1% Critical 5% Critical 10% Critical
Statistic Value Value Value