Descomposición de series temporales

library(forecast)
library(readxl)
serie_IVAE <- read_excel("C:/Users/gusta_000/Desktop/Econometria/IVAE_SLV_C.xlsx", col_types = c("skip", "numeric"),skip = 5)
library(magrittr)
serie_IVAE.ts <- ts(data = serie_IVAE,start = c(2009, 1),frequency = 12)

serie_IVAE.ts %>% autoplot(main = "Gráfico 1: IVAE, El Salvador 2009-2021[marzo]",xlab = "Años/Meses",ylab = "Indice")

Modelo Aditivo

componente de tendencia

ma2_12 <- ma(serie_IVAE.ts, 12, centre = TRUE)
autoplot(serie_IVAE.ts,main = "Gráfico 2: IVAE, El Salvador 2009-2021[marzo]",xlab = "Años/Meses", ylab = "Indice")+autolayer(ma2_12,series = "Tt")

Cálculo de los Factores Estacionales [Componente St]

library(magrittr)
YT<- serie_IVAE.ts 
TT<- ma2_12 
SI<- YT - TT 
St<- tapply(SI, cycle(SI), mean, na.rm = TRUE)
St<- St - sum(St) / 12 
St<-rep(St, len = length(YT)) %>% ts(start = c(2009, 1), frequency = 12) 
autoplot(St,main = "Gráfico 3: Factores Estacionales IVAE 2009 - 2021 ", xlab = "Años/Meses", ylab = "Factor Estacional")

Cálculo del Componente Irregular It

It<-YT-TT-St
autoplot(It,
         main = "Gráfico 4: Componente Irregular, IVAE 2009 - 2021",
         xlab = "Años/Meses",
         ylab = "It")

Descomposición Aditiva

aditiva_descomposicion<-decompose(serie_IVAE.ts,type = "additive")
autoplot(aditiva_descomposicion,main="Gráfico 5: Descomposición Aditiva",xlab="Años/Meses")

Descomposición Aditiva usando libreria feasts

library(tsibble)
library(feasts)
library(ggplot2)
YT %>% as_tsibble() %>%
  model(
    classical_decomposition(value, type = "additive")
  ) %>%
  components() %>%
  autoplot() +
  labs(title = "Gráfico 6: Descomposición Clásica Aditiva, IVAE")+xlab("Años/Meses")

MODELO MULTIPLICATIVO

Componente Tendencia Ciclo [Tt=TCt]

TT<- ma(serie_IVAE.ts, 12, centre = TRUE)
autoplot(TT,main = "Gráfico 7: Componente Tendencia [Ciclo]", xlab = "Años/Meses",ylab = "TT")

Cálculo de Factores Estacionales [St]

SI<-YT/TT 
St <- tapply(SI, cycle(SI), mean, na.rm = TRUE) 
St <- St*12/sum(St) 
St <-
  rep(St, len = length(YT)) %>% ts(start = c(2009, 1), frequency = 12) 
autoplot(St,
         main = "Gráfico 8: Factores Estacionales",
         xlab = "Años/Meses",
         ylab = "Factor Estacional") 

Cálculo del Componente Irregular [It]

It<-YT/(TT*St)
autoplot(It,
         main = "Gráfico 9: Componente Irregular",
         xlab = "Años/Meses",
         ylab = "It")

Descomposición Multiplicativa (usando la librería stats)

descomposicion_multiplicativa<-decompose(serie_IVAE.ts,type = "multiplicative")
autoplot(descomposicion_multiplicativa,main="Gráfico 10: Descomposición Multiplicativa",xlab="Años/Meses")

Descomposición Multiplicativa usando libreria feasts

library(tsibble)
library(feasts)
library(ggplot2)

YT %>% as_tsibble() %>%
  model(
    classical_decomposition(value, type = "multiplicative")
  ) %>%
  components() %>%
  autoplot() +
  labs(title = "Gráfico 11: Descomposición Clásica Multiplicativa, IVAE")+xlab("Años/Meses")