Exemplo Séries Temporais: PIB Brasileiro

Licença

This work is licensed under the Creative Commons Attribution-ShareAlike 4.0 International License. To view a copy of this license, visit http://creativecommons.org/licenses/by-sa/4.0/ or send a letter to Creative Commons, PO Box 1866, Mountain View, CA 94042, USA.

License: CC BY-SA 4.0

Citação

Sugestão de citação: FIGUEIREDO, Adriano Marcos Rodrigues. Exemplo Séries Temporais: PIB Brasileiro. Campo Grande-MS, Brasil: RStudio/Rpubs, 2023. Disponível em https://rpubs.com/amrofi/exemplo_pib_brasil.

1 Dados e plot inicial

Os dados vem do software. O pacote ecoseries foi removido do CRAN, e portanto, tentaremos com rbcb ou sidrar. Neste exemplo utiliza-se a série n. 22109 do PIB a preços de mercado (SCN-2010, Trimestral, (1995=100), dados dessazonalizados ), de 1996T01 até 2022T04.

library(rbcb)
bacen<-rbcb::get_series(22109)
bacen<-bacen[,2]
attach(bacen)
pib.ts<-ts(bacen, start = c(1996,1),frequency=4)
plot(pib.ts,main="PIB a preços de mercado 
     SCN-2010 (Trimestral) (1995=100)",
     sub="BCB série 22109",type = "o",col="black",lwd=2,lty=1, 
     ylab = "Índice 1995=100",xlab="trimestre")

2 No pacote lattice

library(lattice)
xyplot.ts(pib.ts,
          xlab = "Datas - Anos", ylab = "Índice (1995=100)",
          main = "PIB a preços de mercado - Série dessazonalizada",
          sub="Fonte: BCB, 2023. Utilizando o pacote rbcb.")

3 No ggplot2:

library(ggplot2)
datas <- seq(as.Date(paste(c(start(pib.ts),1), collapse="/")), 
             by = "quarter", length.out = length(pib.ts))
dados.df <- data.frame(datas, bacen) # dataframe 
colnames(dados.df)<-c("datas","PIB")
## plot
ggplot(dados.df) +
  aes(x = datas, y = PIB) +
  geom_line(colour = "#112446") +
  labs(
    x = "Datas - Anos",
    y = "Índice (1995=100)",
    title = "PIB a preços de mercado - Série dessazonalizada",
    subtitle = "BCB série 22109, SCN-2010 (Trimestral) ",
    caption = "Fonte: BCB, 2023. Utilizando o pacote rbcb."
  ) +
  ggthemes::theme_stata()

4 No dygraphs

library(dygraphs)
dygraph(dados.df, 
        main = "PIB a preços de mercado - Série dessazonalizada BCB série 22109 <br> SCN-2010 (Trimestral)",
        xlab = "Datas - Anos", ylab = "Índice (1995=100)") %>% 
    dyAxis("x", drawGrid = TRUE) %>% 
    dyEvent("2008-07-01", "2008T3", labelLoc = "bottom") %>%
    dyEvent("2009-01-01", "2009T1", labelLoc = "bottom") %>% 
    dyEvent("2014-01-01", "2014T1", labelLoc = "bottom") %>% 
    dyEvent("2016-10-01", "2016T4", labelLoc = "bottom") %>% 
    dyEvent("2020-04-01", "2020T2", labelLoc = "bottom") %>% 
    dyOptions(drawPoints = TRUE, pointSize = 2)

5 Séries em diferenças

# PIB
dpib<-diff(pib.ts,1)
dpib2<-diff(pib.ts,2)
plot(dpib2,
     main="Séries de Diferenças do PIB",
     type = "o",col="black",lwd=2,lty=1, 
     ylab = "Indice",xlab="trimestre")
lines(dpib,type="o",col = "red",lwd=2,lty=2)
legend("bottomleft",c("d(PIB,2)", "d(PIB)"),
       cex=0.7,lwd=2,lty=1:2,col=c(1,2))