Séries Temporais: Gráficos para Sazonalidade

Além da decomposição apresentada no post Séries Temporais: Introdução e Decomposição dos Componentes em R, há outras formas de visualizar a sazonalidade da série. A séries utilizada aqui é a a10 do pacote fpp, e que contém as informações de venda de antibióticos na Austrália no período de 1991 a 2008.


O primeiro passo é instalar o pacote fpp e ggplot2:

library('ggplot2'); # apenas para gerar o autoplot
library('fpp'); # pacote de forecast

Começamos com um gráfico simples da série:

autoplot(a10) +
 ggtitle("Graf. Autoplot para Vendas de Antibioticos na Australia") +
 ylab("$ milhoes") + xlab("Ano")

graf_sazonalidade1

Em seguida, traçamos um gráfico de sazonalidade com a função seasonplot():

# Grafico de sazonalidade
seasonplot(a10, year.labels=TRUE, main="Graf. Seasonplot para Vendas de Antibioticos na Australia");

graf_sazonalidade2Como você pode ver, temos as linhas ano a ano confrontadas. É um gráfico simples de séries, porém com a sazaonalidade confrontada. Podemos observar, por exemplo, que a venda de antibióticos cresce de fevereiro a março, e que janeiro é o pico de vendas. Esse comportamento é padrão para todos os anos, com exceção de 2008. Pode ser o caso dos dados não terem sido completamente coletados (visto que o ano está incompleto).

É possível deixar o gráfico um pouco melhor, com alguns detalhes para facilitar a visualização. Em outras palavras, deixar mais bonito:

# Grafico de sazonalidade com mais detalhes
ggseasonplot(a10, year.labels=TRUE, year.labels.left=TRUE) +
 ylab("$ milhoes") +xlab("Mes")+ ggtitle("Graf. GGSeasonplot para Vendas de Antibioticos na Australia");

graf_sazonalidade3

Outra opção bacana no ggseasonplot() é fazer um gráfico polar:

ggseasonplot(a10, polar=TRUE) +
 ylab("$ milhoes") + xlab("Mes") + ggtitle("Graf. Polar");

graf_sazonalidade4

Por fim, você pode comparar as diferentes estações por mês. No gráfico abaixo, temos as séries divididas pelos meses com o traço azul indicando a média de cada um:

ggsubseriesplot(a10) + ylab("$ milhoes") + xlab("Mes")+
 ggtitle("Graf. GGSubseriesplot")

graf_sazonalidade5

Anúncios

3 comentários

    1. Olá Carlos. Obrigado pelo elogio.
      Não são não, são meros gráficos temporais, porém plotando as “season” umas contra as outras. Para o correlograma, recomendo usar a função ggAcf().
      Abraços

      Curtir

Deixe um comentário

Preencha os seus dados abaixo ou clique em um ícone para log in:

Logotipo do WordPress.com

Você está comentando utilizando sua conta WordPress.com. Sair / Alterar )

Imagem do Twitter

Você está comentando utilizando sua conta Twitter. Sair / Alterar )

Foto do Facebook

Você está comentando utilizando sua conta Facebook. Sair / Alterar )

Foto do Google+

Você está comentando utilizando sua conta Google+. Sair / Alterar )

Conectando a %s