#COMPONENTE IVAE

Descomposición de series temporales (Enfoque Tradicional)

#CARGA DE DATA

library(readxl)
## Warning: package 'readxl' was built under R version 4.1.3
library(tidyr)
## Warning: package 'tidyr' was built under R version 4.1.3
library(dplyr)
## Warning: package 'dplyr' was built under R version 4.1.3
## 
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
library(forecast)
## Warning: package 'forecast' was built under R version 4.1.3
## Registered S3 method overwritten by 'quantmod':
##   method            from
##   as.zoo.data.frame zoo
IVAE_COMPLETO <- read_excel("C:/Users/Administrador/Desktop/Clase Eco.Practica/IVAE_COMPLETO.xlsx", 
    col_names = FALSE, skip = 6, n_max = 10)
## New names:
## * `` -> `...1`
## * `` -> `...2`
## * `` -> `...3`
## * `` -> `...4`
## * `` -> `...5`
## * `` -> `...6`
## * `` -> `...7`
## * `` -> `...8`
## * `` -> `...9`
## * `` -> `...10`
## * `` -> `...11`
## * `` -> `...12`
## * `` -> `...13`
## * `` -> `...14`
## * `` -> `...15`
## * `` -> `...16`
## * `` -> `...17`
## * `` -> `...18`
## * `` -> `...19`
## * `` -> `...20`
## * `` -> `...21`
## * `` -> `...22`
## * `` -> `...23`
## * `` -> `...24`
## * `` -> `...25`
## * `` -> `...26`
## * `` -> `...27`
## * `` -> `...28`
## * `` -> `...29`
## * `` -> `...30`
## * `` -> `...31`
## * `` -> `...32`
## * `` -> `...33`
## * `` -> `...34`
## * `` -> `...35`
## * `` -> `...36`
## * `` -> `...37`
## * `` -> `...38`
## * `` -> `...39`
## * `` -> `...40`
## * `` -> `...41`
## * `` -> `...42`
## * `` -> `...43`
## * `` -> `...44`
## * `` -> `...45`
## * `` -> `...46`
## * `` -> `...47`
## * `` -> `...48`
## * `` -> `...49`
## * `` -> `...50`
## * `` -> `...51`
## * `` -> `...52`
## * `` -> `...53`
## * `` -> `...54`
## * `` -> `...55`
## * `` -> `...56`
## * `` -> `...57`
## * `` -> `...58`
## * `` -> `...59`
## * `` -> `...60`
## * `` -> `...61`
## * `` -> `...62`
## * `` -> `...63`
## * `` -> `...64`
## * `` -> `...65`
## * `` -> `...66`
## * `` -> `...67`
## * `` -> `...68`
## * `` -> `...69`
## * `` -> `...70`
## * `` -> `...71`
## * `` -> `...72`
## * `` -> `...73`
## * `` -> `...74`
## * `` -> `...75`
## * `` -> `...76`
## * `` -> `...77`
## * `` -> `...78`
## * `` -> `...79`
## * `` -> `...80`
## * `` -> `...81`
## * `` -> `...82`
## * `` -> `...83`
## * `` -> `...84`
## * `` -> `...85`
## * `` -> `...86`
## * `` -> `...87`
## * `` -> `...88`
## * `` -> `...89`
## * `` -> `...90`
## * `` -> `...91`
## * `` -> `...92`
## * `` -> `...93`
## * `` -> `...94`
## * `` -> `...95`
## * `` -> `...96`
## * `` -> `...97`
## * `` -> `...98`
## * `` -> `...99`
## * `` -> `...100`
## * `` -> `...101`
## * `` -> `...102`
## * `` -> `...103`
## * `` -> `...104`
## * `` -> `...105`
## * `` -> `...106`
## * `` -> `...107`
## * `` -> `...108`
## * `` -> `...109`
## * `` -> `...110`
## * `` -> `...111`
## * `` -> `...112`
## * `` -> `...113`
## * `` -> `...114`
## * `` -> `...115`
## * `` -> `...116`
## * `` -> `...117`
## * `` -> `...118`
## * `` -> `...119`
## * `` -> `...120`
## * `` -> `...121`
## * `` -> `...122`
## * `` -> `...123`
## * `` -> `...124`
## * `` -> `...125`
## * `` -> `...126`
## * `` -> `...127`
## * `` -> `...128`
## * `` -> `...129`
## * `` -> `...130`
## * `` -> `...131`
## * `` -> `...132`
## * `` -> `...133`
## * `` -> `...134`
## * `` -> `...135`
## * `` -> `...136`
## * `` -> `...137`
## * `` -> `...138`
## * `` -> `...139`
## * `` -> `...140`
## * `` -> `...141`
## * `` -> `...142`
## * `` -> `...143`
## * `` -> `...144`
## * `` -> `...145`
## * `` -> `...146`
## * `` -> `...147`
## * `` -> `...148`
## * `` -> `...149`
## * `` -> `...150`
## * `` -> `...151`
## * `` -> `...152`
## * `` -> `...153`
## * `` -> `...154`
## * `` -> `...155`
## * `` -> `...156`
## * `` -> `...157`
## * `` -> `...158`
## * `` -> `...159`
## * `` -> `...160`
#Generamos serie temporal
data.ivae<-pivot_longer(data = IVAE_COMPLETO[1,],names_to = "vars",cols = 2:160,values_to = "indice") %>% select("indice")
data.ivae.ts<- data.ivae %>% ts(start = c(2009,1),frequency = 12)

#Graficamos
data.ivae.ts %>% 
  autoplot(main ="IVAE ENE 2009- MAR 2022",
           xlab="Años/Meses",
           ylab="Indice")

##DESCOMPOSICIÓN ADITIVA IVAE

library(forecast)
#Descomposición aditiva
descompocicion_aditiva<-decompose(data.ivae.ts,type = "additive")
#gráfico de la descompocición
autoplot(descompocicion_aditiva, main="Descompisicion")

#gráfico de la serie original y el componente tendencia ciclo
Yt<-descompocicion_aditiva$x
TC<-descompocicion_aditiva$trend
autoplot(Yt,
         main = "Componente IVAE")+
  autolayer(TC)
## Warning: Removed 12 row(s) containing missing values (geom_path).

##CARGA DE DATA AGRICULTURA, GANADERIA, SILVICULTURA Y PESCA

data_agri<-pivot_longer(data = IVAE_COMPLETO[2,],names_to = "vars", cols = 2:160,values_to = "indice") %>% select("indice")
data_agri_ts<-data_agri %>% ts(start = c(2009,1),frequency = 12)

##DESCOMPOSICIÓN ADITIVA AGRICULTURA, GANADERIA, SILVICULTURA Y PESCA

descomposicion_aditiva2<-decompose(data_agri_ts,type = "additive")
Yt2<-descomposicion_aditiva2$x
Tc2<-descomposicion_aditiva2$trend

autoplot(Yt2,main="Agricultura, Ganaderia, Silvicultura y Pesca",xlab="Años/Meses",ylab = "indice") + autolayer(Tc2)
## Warning: Removed 12 row(s) containing missing values (geom_path).

##CARGA DE DATA INDICE DE PRODUCCOÓN INDUSTRIAL: INDUSTRIAS MANUFACTURERAS, EXPLOTACIÓN DE MINAS Y CANTERAS Y OTRAS ACTIVIDADES INDUSTRIALES

data_IPI<-pivot_longer(data = IVAE_COMPLETO[3,],names_to = "vars", cols = 2:160,values_to = "indice") %>% select("indice")
data_IPI_ts<-data_IPI %>% ts(start = c(2009,1),frequency = 12)

##DESCOMPOSICIÓN ADITIVA DE PRODUCCOÓN INDUSTRIAL

descomposicion_aditiva3<-decompose(data_IPI_ts,type = "additive")
Yt3<-descomposicion_aditiva3$x
Tc3<-descomposicion_aditiva3$trend

autoplot(Yt3,main=" Indice de Producción Industrial: Industrias Manufactureras, Explotacion de Minas y Canteras y Otras Actividades Industriales",xlab="Años/Meses",ylab = "indice") + autolayer(Tc3)
## Warning: Removed 12 row(s) containing missing values (geom_path).

##CARGA DE DATA DE CONSTRUCCIÓN

data_constru<-pivot_longer(data = IVAE_COMPLETO[4,],names_to = "vars", cols = 2:160,values_to = "indice") %>% select("indice")
data_constru_ts<-data_constru %>% ts(start = c(2009,1),frequency = 12)

##DESCOMPOSICIÓN ADITIVA DE CONSTRUCCIÓN

descomposicion_aditiva4<-decompose(data_constru_ts,type = "additive")
Yt4<-descomposicion_aditiva4$x
Tc4<-descomposicion_aditiva4$trend

autoplot(Yt4,main=" Construccion",xlab="Años/Meses",ylab = "indice") + autolayer(Tc4)
## Warning: Removed 12 row(s) containing missing values (geom_path).

##CARGA DE DATA DE COMERCIO, TRANSPORTE Y ALMACENAMIENTO, ACTIVIDADES DE ALOJAMIENTO Y DE SERVICIO DE COMIDAS

data_comer<-pivot_longer(data = IVAE_COMPLETO[5,],names_to = "vars", cols = 2:160,values_to = "indice") %>% select("indice")
data_comer_ts<-data_comer %>% ts(start = c(2009,1),frequency = 12)

##DESCOMPOSICIÓN ADITIVA DE COMERCIO, TRANSPORTE Y ALMACENAMIENTO, ACTIVIDADES DE ALOJAMIENTO Y DE SERVICIO DE COMIDAS

descomposicion_aditiva5<-decompose(data_comer_ts,type = "additive")
Yt5<-descomposicion_aditiva5$x
Tc5<-descomposicion_aditiva5$trend

autoplot(Yt5,main=" Comercio, Transporte y Almacenamiento, Actividades de Alojamiento y de Servicio de Comidas",xlab="Años/Meses",ylab = "indice") + autolayer(Tc5)
## Warning: Removed 12 row(s) containing missing values (geom_path).

##CARGA DE DATA DE INFORMATICA Y COMUNICACIONES

data_info<-pivot_longer(data = IVAE_COMPLETO[6,],names_to = "vars", cols = 2:160,values_to = "indice") %>% select("indice")
data_info_ts<-data_info %>% ts(start = c(2009,1),frequency = 12)

##DESCOMPOSICIÓN ADITIVA DE INFORMATICA Y COMUNICACIONES

descomposicion_aditiva6<-decompose(data_info_ts,type = "additive")
Yt6<-descomposicion_aditiva6$x
Tc6<-descomposicion_aditiva6$trend

autoplot(Yt6,main="Información y Comunicaciones",xlab="Años/Meses",ylab = "indice") + autolayer(Tc6)
## Warning: Removed 12 row(s) containing missing values (geom_path).

##CARGA DE DATA DE ACTIVIDADES FINANCIERAS Y DE SEGUROS

data_Finan<-pivot_longer(data = IVAE_COMPLETO[7,],names_to = "vars", cols = 2:160,values_to = "indice") %>% select("indice")
data_Finan_ts<-data_Finan %>% ts(start = c(2009,1),frequency = 12)

##DESCOMPOSICIÓN ADITIVA DE ACTIVIDADES FINANCIERAS Y DE SEGUROS

descomposicion_aditiva7<-decompose(data_Finan_ts,type = "additive")
Yt7<-descomposicion_aditiva7$x
Tc7<-descomposicion_aditiva7$trend

autoplot(Yt7,main="Actividades Financieras y de Seguros",xlab="Años/Meses",ylab = "indice") + autolayer(Tc7)
## Warning: Removed 12 row(s) containing missing values (geom_path).

##CARGA DE DATA DE ACTIVIDADES INMOBILIARIAS

data_inmobi<-pivot_longer(data = IVAE_COMPLETO[8,],names_to = "vars", cols = 2:160,values_to = "indice") %>% select("indice")
data_inmobi_ts<-data_inmobi %>% ts(start = c(2009,1),frequency = 12)

##DESCOMPOSICIÓN ADITIVA DE ACTIVIDADES INMOBILIARIAS

descomposicion_aditiva8<-decompose(data_inmobi_ts,type = "additive")
Yt8<-descomposicion_aditiva8$x
Tc8<-descomposicion_aditiva8$trend

autoplot(Yt8,main="Actividades Inmobiliarias",xlab="Años/Meses",ylab = "indice") + autolayer(Tc8)
## Warning: Removed 12 row(s) containing missing values (geom_path).

##CARGA DE DATA DE ACTIVIDADES PROFESIONALES, CIENTIFICAS, TÉCNICAS, ADMINISTRATIVAS, APOYO Y OTROS SERVICIOS

data_acti<-pivot_longer(data = IVAE_COMPLETO[9,],names_to = "vars", cols = 2:160,values_to = "indice") %>% select("indice")
data_acti_ts<-data_acti %>% ts(start = c(2009,1),frequency = 12)

##DESCOMPOSICIÓN ADITIVA DE ACTIVIDADES PROFESIONALES, CIENTIFICAS, TÉCNICAS, ADMINISTRATIVAS, APOYO Y OTROS SERVICIOS

descomposicion_aditiva9<-decompose(data_acti_ts,type = "additive")
Yt9<-descomposicion_aditiva9$x
Tc9<-descomposicion_aditiva9$trend

autoplot(Yt9,main="Actividades Profesionales, Científicas, Técnicas, Administrativas, de Apoyo y Otros Servicios",xlab="Años/Meses",ylab = "indice") + autolayer(Tc9)
## Warning: Removed 12 row(s) containing missing values (geom_path).

##CARGA DE DATA DE ACTIVIDADES ADMINISTRATIVAS PÚBLICA Y DEFENSA, ENSEÑANZA, SALUD Y ASISTENCIA SOCIAL

data_admon<-pivot_longer(data = IVAE_COMPLETO[10,],names_to = "vars", cols = 2:160,values_to = "indice") %>% select("indice")
data_admon_ts<-data_admon %>% ts(start = c(2009,1),frequency = 12)

##DESCOMPOSICIÓN ADITIVA DE ACTIVIDADES ADMINISTRATIVAS PÚBLICA Y DEFENSA, ENSEÑANZA, SALUD Y ASISTENCIA SOCIAL

descomposicion_aditiva10<-decompose(data_admon_ts,type = "additive")
Yt10<-descomposicion_aditiva10$x
Tc10<-descomposicion_aditiva10$trend

autoplot(Yt10,main="Actividades de Administración Pública y Defensa, Enseñanza, Salud y Asistencia Socials",xlab="Años/Meses",ylab = "indice") + autolayer(Tc10)
## Warning: Removed 12 row(s) containing missing values (geom_path).