SA19008_DescompSerieTemp_3

Descomposición de Series Temporales.

#Carga de datos
library(forecast)
library(readxl)
serie_ivae <- read_excel("C:/Users/Abigail/Desktop/econometria/ejercicios_de_clase/DescompocionSeriesTemp/Eduardo Isaias Sandoval Arias - IVAE_SLV_C.xlsx",
                  col_types = c("skip", "numeric"),
                  skip = 5)

serie_ivae_TS <- ts(data = serie_ivae,
                    start = c(2009, 1),
                    frequency = 12)
serie_ivae_TS %>% 
  autoplot(main = "IVAE, El Salvador 2009-2021[marzo]",
                           xlab = "Años/Meses",
                           ylab = "Indice",colour = 'red')

Modelo Aditivo

Calculo manual

#Calculo de componente de tendencia
library(forecast)
ma2_12 <- ma(serie_ivae_TS, 12, centre = TRUE)
autoplot(serie_ivae_TS,main = "IVAE, El Salvador 2009-2021[marzo]",
           xlab = "Años/Meses",
           ylab = "Indice")+
  autolayer(ma2_12,series = "Tt")

#Calculo de factores estacionales
library(magrittr)
library(forecast)
Yt <- serie_ivae_TS 
Tt <- ma2_12 
SI <- Yt - Tt  
St <- tapply(SI, cycle(SI), mean, na.rm = TRUE) 
St <- St - sum(St) / 12 

#Generar la serie de factores para cada valor de la serie original

St <- rep(St, len = length(Yt)) %>% ts(start = c(2009, 1), frequency = 12) 
autoplot(St,
         main = "Factores Estacionales",
         xlab = "Años/Meses",
         ylab = "Factor Estacional",colour="blue") 

Cálculo del Componente Irregular

It<-Yt-Tt-St
autoplot(It,
         main = "Componente Irregular",
         xlab = "Años/Meses",
         ylab = "It",
         colour="purple")

Descomposición Aditiva (uso de librerías)

#Descomposición Aditiva (usando la libreria stats)

descomposicion_aditiva<-decompose(serie_ivae_TS,type = "additive")
autoplot(descomposicion_aditiva,main="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 = "Descomposición Clásica Aditiva, IVAE")+xlab("Años/Meses")

Modelo Multiplicativo

Calculo manual

#Componente Tendencia Ciclo
Tt<- ma(serie_ivae_TS, 12, centre = TRUE)
autoplot(Tt,main = "Componente Tendencia [Ciclo]", xlab = "Años/Meses",ylab = "Tt",colour="blue")

#Cálculo de Factores Estacionales
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 = "Factores Estacionales",
         xlab = "Años/Meses",
         ylab = "Factor Estacional",colour="red") 

#Cálculo del Componente Irregular
It<-Yt/(Tt*St)
autoplot(It,
         main = "Componente Irregular",
         xlab = "Años/Meses",
         ylab = "It",colour="orange")

Descomposición Multiplicativa (Uso de librerías)

#Descomposición Multiplicativa (libreria stats)
library(stats)
descomposicion_multiplicatica<-decompose(serie_ivae_TS,type = "multiplicative")
autoplot(descomposicion_multiplicatica,main="Descomposición Multiplicativa",xlab="Años/Meses")

#Descomposición Multiplicativa (libreria feasts)
library(tsibble)
library(feasts)
library(ggplot2)
Yt %>% as_tsibble() %>%
  model(classical_decomposition(value, type = "multiplicative")) %>%
  components() %>%
  autoplot() +
  labs(title = "Descomposición Clásica Multiplicativa, IVAE") + xlab("Años/Meses")