Indicaciones:

  • Aplique el filtro HP para cada uno de los países de Centroamérica, para los periodos disponibles para cada país

  • represente sus resultados de forma gráfica:

      1. Usando la versión personalizada y
      1. usando la representación rápida.
    • En todos los gráficos indique el país y el periodo al que corresponden los datos.

Version personalizada

library(readxl)
library(mFilter)
library(forecast)

datos <- read_excel("datos_PIB_trim_CA.xlsx", skip = 4)

paises <- colnames(datos)[-1]   

for (pais in paises) {
  

  
  serie <- as.numeric(datos[[pais]])
  
  serie <- na.interp(serie)   
  
  PIB_trim <- ts(serie,
                 start = c(1990,1), 
                 frequency = 4)
  
  hp_result <- hpfilter(PIB_trim, freq = 1600)
  
  cycle <- hp_result$cycle
  
  plot(cycle, type = "l",
       ylab = "Ciclo",
       xlab = "Tiempo")
  
  title(main = paste("Componente Cíclico -", pais), cex.main = 1.3)
  
}

Representacion rapida

plot(hp_result, ask = FALSE)
title(main = paste("Filtro HP -", pais,
                   "\nPeriodo:", start(PIB_trim)[1], "-", end(PIB_trim)[1]))