Descomposicion de series temporales

Carga de Datos, IVAE

library(readxl)
library(tidyr)
library(dplyr)
library(forecast)
IVAE_practica <- read_excel("C:/Users/MINEDUCYT/Desktop/Guias de econometria/IVAE practica.xlsx", col_names = FALSE, skip = 6, n_max = 10)

data.ivae<-pivot_longer(data = IVAE_practica[1,],names_to = "vars",cols = 2:160,values_to = "indice") %>% select("indice")
data.ivae.ts<- data.ivae %>% ts(start = c(2009,1),frequency = 12)
data.ivae.ts %>% 
  autoplot(main ="IVAE ENE 2009- MAR 2022",
           xlab="Años/Meses",
           ylab="Indice")

Descomposicion de IVAE [Modelo Aditivo]

library(forecast)
#Descomposición aditiva
descomposicion_aditiva<-decompose(data.ivae.ts,type = "additive")
#gráfico de la descompocición
autoplot(descomposicion_aditiva, main="Descomposicion Aditiva, IVAE ENE 2009-MAR 2022")

#gráfico de la serie original y el componente tendencia ciclo
Yt<-descomposicion_aditiva$x
TC<-descomposicion_aditiva$trend
autoplot(Yt,
         main = "Descomposicion Aditiva, IVAE ENE 2009-MAR 2022")+
  autolayer(TC)

Componente 1 : Agricultura, Ganadería, Silvicultura y Pesca

Carga de datos

library(forecast)
library(readxl)
library(tidyr)
library(dplyr)
componente_1<-pivot_longer(data = IVAE_practica[2,],names_to = "vars",cols = 2:160,values_to = "indice") %>% select("indice")
componente_1.ts<- componente_1 %>% ts(start = c(2009,1),frequency = 12)

Descomposicion de Agricultura, Ganadería, Silvicultura y Pesca [Modelo aditivo]

library(forecast)
#Descomposición aditiva
descomposicion_1<-decompose(componente_1.ts,type = "additive")
#gráfico de la descompocición
autoplot(descomposicion_1, main="Descomposicion Aditiva. Agricultura, Ganadería, Silvicultura y Pesca ENE 2009-MAR 2022")

#gráfico de la serie original y el componente tendencia ciclo
Yt1<-descomposicion_1$x
TC1<-descomposicion_1$trend
autoplot(Yt1,
         main = "Descomposicion Aditiva.Agricultura, Ganadería, Silvicultura y Pesca ENE 2009-MAR 2022")+
  autolayer(TC1)

Componente 2: Índice de Producción Industrial

Carga de datos

library(forecast)
componente_2<-pivot_longer(data = IVAE_practica[3,],names_to = "vars",cols = 2:160,values_to = "indice") %>% select("indice")
componente_2.ts<- componente_2 %>% ts(start = c(2009,1),frequency = 12)

Descomposicion de Índice de Producción Industrial [modelo aditivo]

library(forecast)
#Descomposición aditiva
descomposicion_2<-decompose(componente_2.ts,type = "additive")
#gráfico de la descompocición
autoplot(descomposicion_2, main="Descomposicion Aditiva.Indice de Produccion Industrial ENE 2009-MAR 2022")

#gráfico de la serie original y el componente tendencia ciclo
Yt2<-descomposicion_2$x
TC2<-descomposicion_2$trend
autoplot(Yt2,
         main = "Descomposicion Aditiva.Índice de Producción Industrial ENE 2009-MAR 2022")+
  autolayer(TC2)

Componente 3: Construccion

Carga de datos

library(forecast)
componente_3<-pivot_longer(data = IVAE_practica[4,],names_to = "vars",cols = 2:160,values_to = "indice") %>% select("indice")
componente_3.ts<- componente_3 %>% ts(start = c(2009,1),frequency = 12)

Descomposicion de Construccion [modelo aditivo]

library(forecast)
#Descomposición aditiva
descomposicion_3<-decompose(componente_3.ts,type = "additive")
#gráfico de la descompocición
autoplot(descomposicion_3, main="Descomposicion Aditiva.Construccion ENE 2009-MAR 2022")

#gráfico de la serie original y el componente tendencia ciclo
Yt3<-descomposicion_3$x
TC3<-descomposicion_3$trend
autoplot(Yt3,
         main = "Descomposicion Aditiva.Construccion ENE 2009-MAR 2022")+
  autolayer(TC3)

Componente 4: Comercio, Transporte y Almacenamiento, Actv. alojamiento y servicio de comida

Carga de datos

library(forecast)
componente_4<-pivot_longer(data = IVAE_practica[5,],names_to = "vars",cols = 2:160,values_to = "indice") %>% select("indice")
componente_4.ts<- componente_4 %>% ts(start = c(2009,1),frequency = 12)

Descomposicion de Comercio, Transporte y Almacenamiento, Actv. alojamiento y servicio de comida [modelo aditivo]

library(forecast)
#Descomposición aditiva
descomposicion_4<-decompose(componente_4.ts,type = "additive")
#gráfico de la descompocición
autoplot(descomposicion_4, main="Descomposicion Aditiva.Comercio, Transporte y Almacenamiento, Actv. alojamiento y servicio de comida ENE 2009-MAR 2022")

#gráfico de la serie original y el componente tendencia ciclo
Yt4<-descomposicion_4$x
TC4<-descomposicion_4$trend
autoplot(Yt4,
         main = "Descomposicion Aditiva.Comercio, Transporte y Almacenamiento, Actv. alojamiento y servicio de comida ENE 2009-MAR 2022")+
  autolayer(TC4)

Componente 5: Información y Comunicaciones

Carga de datos

library(forecast)
componente_5<-pivot_longer(data = IVAE_practica[6,],names_to = "vars",cols = 2:160,values_to = "indice") %>% select("indice")
componente_5.ts<- componente_5 %>% ts(start = c(2009,1),frequency = 12)

Descomposicion de Informacion y Comunicaciones [modelo aditivo]

library(forecast)
#Descomposición aditiva
descomposicion_5<-decompose(componente_5.ts,type = "additive")
#gráfico de la descompocición
autoplot(descomposicion_5, main="Descomposicion aditiva. Informacion y comunicaciones ENE 2009-MAR 2022")

#gráfico de la serie original y el componente tendencia ciclo
Yt5<-descomposicion_5$x
TC5<-descomposicion_5$trend
autoplot(Yt5,
         main = "Descomposicion Aditiva. Informacion y Comunicaciones ENE 2009-MAR 2022")+
  autolayer(TC5)

Componente 6: Actividades Financieras y de Seguros

Carga de datos

library(forecast)
componente_6<-pivot_longer(data = IVAE_practica[7,],names_to = "vars",cols = 2:160,values_to = "indice") %>% select("indice")
componente_6.ts<- componente_6 %>% ts(start = c(2009,1),frequency = 12)

Descomposicion de Actividades financieras y de Seguros [modelo aditivo]

library(forecast)
#Descomposición aditiva
descomposicion_6<-decompose(componente_6.ts,type = "additive")
#gráfico de la descompocición
autoplot(descomposicion_6, main="Descomposicion aditiva.Actividades financieras  y de Seguros ENE 2009-MAR 2022")

#gráfico de la serie original y el componente tendencia ciclo
Yt6<-descomposicion_6$x
TC6<-descomposicion_6$trend
autoplot(Yt6,
         main = "Descomposicion Aditiva. Actividades financieras y de Seguros ENE 2009-MAR 2022")+
  autolayer(TC6)

Componente 7: Actividades Inmobiliarias

Carga de datos

library(forecast)
componente_7<-pivot_longer(data = IVAE_practica[8,],names_to = "vars",cols = 2:160,values_to = "indice") %>% select("indice")
componente_7.ts<- componente_7 %>% ts(start = c(2009,1),frequency = 12)

Descomposicion de Actividades Inmobilarias [modelo aditivo]

library(forecast)
#Descomposición aditiva
descomposicion_7<-decompose(componente_7.ts,type = "additive")
#gráfico de la descompocición
autoplot(descomposicion_7, main="Descomposicion aditiva.Actividades Inmobilarias ENE 2009-MAR 2022")

#gráfico de la serie original y el componente tendencia ciclo
Yt7<-descomposicion_7$x
TC7<-descomposicion_7$trend
autoplot(Yt7,
         main = "Descomposicion Aditiva. Actividades Inmobilarias ENE 2009-MAR 2022")+
  autolayer(TC7)

Componente 8: Actividades Profesionales, Científicas, Técnicas y otros.

Carga de datos

library(forecast)
componente_8<-pivot_longer(data = IVAE_practica[9,],names_to = "vars",cols = 2:160,values_to = "indice") %>% select("indice")
componente_8.ts<- componente_8 %>% ts(start = c(2009,1),frequency = 12)

Descomposicion de Actividades Profesionales, Científicas, Técnicas y otros.[modelo aditivo]

library(forecast)
#Descomposición aditiva
descomposicion_8<-decompose(componente_8.ts,type = "additive")
#gráfico de la descompocición
autoplot(descomposicion_8, main="Descomposicion aditiva.Actividades Profesionales, Científicas, Técnicas y otros. ENE 2009-MAR 2022")

#gráfico de la serie original y el componente tendencia ciclo
Yt8<-descomposicion_8$x
TC8<-descomposicion_8$trend
autoplot(Yt8,
         main = "Descomposicion Aditiva. Actividades Profesionales, Científicas, Técnicas y otros. ENE 2009-MAR 2022")+
  autolayer(TC8)

Componente 9: Actividades de Administración Pública y Defensa,etc.

Carga de datos

library(forecast)
componente_9<-pivot_longer(data = IVAE_practica[10,],names_to = "vars",cols = 2:160,values_to = "indice") %>% select("indice")
componente_9.ts<- componente_9 %>% ts(start = c(2009,1),frequency = 12)

Descomposicion de Actividades de Administración Pública y Defensa,etc.[modelo aditivo]

library(forecast)
#Descomposición aditiva
descomposicion_9<-decompose(componente_9.ts,type = "additive")
#gráfico de la descompocición
autoplot(descomposicion_9, main="Descomposicion aditiva. Actividades de Administración Pública y Defensa,etc. ENE 2009-MAR 2022")

#gráfico de la serie original y el componente tendencia ciclo
Yt9<-descomposicion_9$x
TC9<-descomposicion_9$trend
autoplot(Yt9,
         main = "Descomposicion Aditiva. Actividades de Administración Pública y Defensa,etc. ENE 2009-MAR 2022")+
  autolayer(TC9)