library(readxl)
library(tidyr)
library(dplyr)
library(forecast)
IVAE_03_22 <- read_excel("IVAE_03_22.xlsx",
col_names = FALSE, skip = 6, n_max = 10)
data.ivae<-pivot_longer(data = IVAE_03_22[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_aditiva<-decompose(data.ivae.ts,type = "additive")
#Gráfico de la descomposición
autoplot(descomposicion_aditiva, main="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 = "IVAE, ENE 2009 - MAR 2022",
xlab="Años/Meses",
ylab="Indice")+
autolayer(TC)
data.agricultura<-pivot_longer(data = IVAE_03_22[2,],names_to = "vars",cols = 2:160,values_to = "indice") %>% select("indice")
data.agricultura.ts<- data.agricultura %>% ts(start = c(2009,1),frequency = 12)
data.agricultura.ts %>%
autoplot(main ="Agricultura, Ganadería, Silvicultura y Pesca, enero 2009 - marzo 2022",
xlab="Años/Meses",
ylab="Indice")
descomposicion_aditiva2<-decompose(data.agricultura.ts,type = "additive")
#Gráfico de la descomposición
autoplot(descomposicion_aditiva2, main="Agricultura, Ganadería, Silvicultura y Pesca, enero 2009 - marzo 2022")
#Gráfico de la serie original y el componente tendencia ciclo
Yt_agricultura<-descomposicion_aditiva2$x
TC_agricultura<-descomposicion_aditiva2$trend
autoplot(Yt_agricultura,
main = "Agricultura, Ganadería, Silvicultura y Pesca, enero 2009 - marzo 2022",
xlab="Años/Meses",
ylab="Indice")+
autolayer(TC_agricultura)
data.industria<-pivot_longer(data = IVAE_03_22[3,],names_to = "vars",cols = 2:160,values_to = "indice") %>% select("indice")
data.industria.ts<- data.industria %>% ts(start = c(2009,1),frequency = 12)
data.industria.ts %>%
autoplot(main ="Industrias Manufactureras, Explotación de Minas y Canteras y Otras Actividades Industriales, enero 2009 - marzo 2022",
xlab="Años/Meses",
ylab="Indice")
descomposicion_aditiva3<-decompose(data.industria.ts,type = "additive")
#Gráfico de la descomposición
autoplot(descomposicion_aditiva3, main="Industrias Manufactureras, Explotación de Minas y Canteras y Otras Actividades Industriales")
#Gráfico de la serie original y el componente tendencia ciclo
Yt_industria<-descomposicion_aditiva3$x
TC_industria<-descomposicion_aditiva3$trend
autoplot(Yt_industria,
main = "Industrias Manufactureras, Explotación de Minas y Canteras y Otras Actividades Industriales",
xlab="Años/Meses",
ylab="Indice")+
autolayer(TC_industria)
data.construccion<-pivot_longer(data = IVAE_03_22[4,],names_to = "vars",cols = 2:160,values_to = "indice") %>% select("indice")
data.construccion.ts<- data.construccion %>% ts(start = c(2009,1),frequency = 12)
data.construccion.ts %>%
autoplot(main ="Construcción",
xlab="Años/Meses",
ylab="Indice")
descomposicion_aditiva4<-decompose(data.construccion.ts,type = "additive")
#Gráfico de la descomposición
autoplot(descomposicion_aditiva4, main="Construcción")
#Gráfico de la serie original y el componente tendencia ciclo
Yt_construccion<-descomposicion_aditiva4$x
TC_construccion<-descomposicion_aditiva4$trend
autoplot(Yt_construccion,
main = "Construcción",
xlab="Años/Meses",
ylab="Indice")+
autolayer(TC_construccion)
data.comercio<-pivot_longer(data = IVAE_03_22[5,],names_to = "vars",cols = 2:160,values_to = "indice") %>% select("indice")
data.comercio.ts<- data.comercio %>% ts(start = c(2009,1),frequency = 12)
data.comercio.ts %>%
autoplot(main ="Comercio, Transporte y Almacenamiento, Actividades de Alojamiento y de Servicio de Comidas",
xlab="Años/Meses",
ylab="Indice")
descomposicion_aditiva5<-decompose(data.comercio.ts,type = "additive")
#Gráfico de la descomposición
autoplot(descomposicion_aditiva5, main="Comercio, Transporte y Almacenamiento, Actividades de Alojamiento y de Servicio de Comidas")
#Gráfico de la serie original y el componente tendencia ciclo
Yt_comercio<-descomposicion_aditiva5$x
TC_comercio<-descomposicion_aditiva5$trend
autoplot(Yt_comercio,
main = "Comercio, Transporte y Almacenamiento, Actividades de Alojamiento y de Servicio de Comidas",
xlab="Años/Meses",
ylab="Indice")+
autolayer(TC_comercio)
data.comunicaciones<-pivot_longer(data = IVAE_03_22[6,],names_to = "vars",cols = 2:160,values_to = "indice") %>% select("indice")
data.comunicaciones.ts<- data.comunicaciones %>% ts(start = c(2009,1),frequency = 12)
data.comunicaciones.ts %>%
autoplot(main ="Información y Comunicaciones",
xlab="Años/Meses",
ylab="Indice")
descomposicion_aditiva6<-decompose(data.comunicaciones.ts,type = "additive")
#Gráfico de la descomposición
autoplot(descomposicion_aditiva6, main="Información y Comunicaciones")
#Gráfico de la serie original y el componente tendencia ciclo
Yt_comunicaciones<-descomposicion_aditiva6$x
TC_comunicaciones<-descomposicion_aditiva6$trend
autoplot(Yt_comunicaciones,
main = "Información y Comunicaciones",
xlab="Años/Meses",
ylab="Indice")+
autolayer(TC_comunicaciones)
data.financias<-pivot_longer(data = IVAE_03_22[7,],names_to = "vars",cols = 2:160,values_to = "indice") %>% select("indice")
data.financias.ts<- data.financias %>% ts(start = c(2009,1),frequency = 12)
data.financias.ts %>%
autoplot(main ="Actividades Financieras y de Seguros",
xlab="Años/Meses",
ylab="Indice")
descomposicion_aditiva7<-decompose(data.financias.ts,type = "additive")
#Gráfico de la descomposición
autoplot(descomposicion_aditiva7, main="Actividades Financieras y de Seguros")
#Gráfico de la serie original y el componente tendencia ciclo
Yt_financias<-descomposicion_aditiva7$x
TC_financias<-descomposicion_aditiva7$trend
autoplot(Yt_financias,
main = "Actividades Financieras y de Seguros",
xlab="Años/Meses",
ylab="Indice")+
autolayer(TC_financias)
data.inmobiliarias<-pivot_longer(data = IVAE_03_22[8,],names_to = "vars",cols = 2:160,values_to = "indice") %>% select("indice")
data.inmobiliarias.ts<- data.inmobiliarias %>% ts(start = c(2009,1),frequency = 12)
data.inmobiliarias.ts %>%
autoplot(main ="Actividades Inmobiliarias",
xlab="Años/Meses",
ylab="Indice")
descomposicion_aditiva8<-decompose(data.inmobiliarias.ts,type = "additive")
#Gráfico de la descomposición
autoplot(descomposicion_aditiva8, main="Actividades Inmobiliarias")
#Gráfico de la serie original y el componente tendencia ciclo
Yt_inmobiliarias<-descomposicion_aditiva8$x
TC_inmobiliarias<-descomposicion_aditiva8$trend
autoplot(Yt_inmobiliarias,
main = "Actividades Inmobiliarias",
xlab="Años/Meses",
ylab="Indice")+
autolayer(TC_inmobiliarias)
data.profesionales<-pivot_longer(data = IVAE_03_22[9,],names_to = "vars",cols = 2:160,values_to = "indice") %>% select("indice")
data.profesionales.ts<- data.profesionales %>% ts(start = c(2009,1),frequency = 12)
data.profesionales.ts %>%
autoplot(main ="Actividades Profesionales, Científicas, Técnicas, Administrativas, de Apoyo y Otros Servicios",
xlab="Años/Meses",
ylab="Indice")
descomposicion_aditiva9<-decompose(data.profesionales.ts,type = "additive")
#Gráfico de la descomposición
autoplot(descomposicion_aditiva9, main="Actividades Profesionales, Científicas, Técnicas, Administrativas, de Apoyo y Otros Servicios")
#Gráfico de la serie original y el componente tendencia ciclo
Yt_profesionales<-descomposicion_aditiva9$x
TC_profesionales<-descomposicion_aditiva9$trend
autoplot(Yt_profesionales,
main = "Actividades Profesionales, Científicas, Técnicas, Administrativas, de Apoyo y Otros Servicios",
xlab="Años/Meses",
ylab="Indice")+
autolayer(TC_profesionales)