Modelagem preditiva de séries temporais hierárquicas no Python

Introdução

Frequentemente, séries temporais podem ser desagregadas em vários atributos de interesse. Por exemplo1, o número total de veículos produzidos no país ao longo do tempo pode ser dividido em veículos comerciais leves, caminhões, ônibus, etc. Cada uma dessas categorias pode ainda ser desagregada em outros níveis, por regiões, por fabricante, etc., caracterizando o que pode ser chamado de "séries temporais hierárquicas".

Essa riqueza de informação e dados possibilita (e pode ser de interesse) gerar previsões desagregadas das séries, de modo que os pontos de previsão das séries desagregadas possam ser analisados individualmente e que, quando agregados de alguma forma, sejam coerentes com os valores agregados da série.

Aprenda a coletar, processar e analisar dados na formação de Do Zero à Análise de Dados Econômicos e Financeiros com Python.

Quer aprender a como criar modelos de previsão? Veja nosso curso Modelagem e Previsão usando Python.

Previsão Hierárquica

Existem várias abordagens para gerar previsões hierárquicas com séries temporais, a mais comum e intuitiva é a bottom-up (de baixo para cima), onde primeiro geramos as previsões para cada série no nível inferior e, em seguida, somamos os pontos de previsão para obter previsões da série acima.

Por outro lado, a abordagem top-down (de cima para baixo) envolve primeiro a geração de previsões para a série agregada e, em seguida, desagregá-la na hierarquia utilizando proporções (geralmente baseadas na série histórica).

Ainda existem outras abordagens como a do Minimum Trace e suas variantes, sendo que cada uma tem seus prós e contras. Para se aprofundar no tema veja Forecasting: principles and practice de Hyndman, R.J., & Athanasopoulos, G. (2021), que apresenta detalhes sobre as abordagens disponíveis.

Exemplo: dados de produção da ANFAVEA

Mensalmente a ANFAVEA disponibiliza séries temporais da produção, licenciamento, exportação, etc. de veículos. Os dados são desagregados pelas categorias citadas acima e neste exemplo iremos explorar alguns métodos de gerar modelos de previsão para as séries de produção de veículos.

Os dados são coletados via SGS do BCB usando a biblioteca python-bcb, no período de jan/1993 até jul/2024. A separação de treino/teste ocorre nos doze últimos meses da série.

Bibliotecas

Para a análise, vamos utilizar as seguintes bibliotecas Python:

    • python-bcb: Para acessar os dados do SGS.
    • hierarchicalforecast: Para realizar previsões hierárquicas.
    • statsforecast: Para utilizar modelos de séries temporais.
    • pandas: Para manipulação de dados.
    • plotnine: Para gerar gráficos.

Coleta dados

O primeiro passo é coletar os dados do SGS. Utilizaremos a função sgs.get() para buscar as séries temporais de produção de veículos, desde janeiro de 1993 até julho de 2024.

Tratamento de Dados

Após coletar os dados, devemos transformá-los em um formato adequado para a análise. Isso inclui:

  • Limpar e padronizar os nomes das colunas.
  • Converter a coluna de data para o formato datetime.
  • Criar uma coluna 'top_level' com o valor 'Total' para indicar o nível superior da hierarquia.
  • Transformar o DataFrame para o formato 'long', separando as categorias de veículos em diferentes linhas.

Criação do Modelo

Para realizar as previsões, vamos utilizar o objeto StatsForecast. Esse objeto permite utilizar diferentes modelos de séries temporais, como AutoARIMA, AutoETS e Naive.

Após gerar as previsões com os diferentes modelos, precisamos reconciliá-las para garantir a coerência entre os níveis da hierarquia. Para isso, vamos utilizar o objeto HierarchicalReconciliation.

Após a reconciliação das previsões, precisamos avaliá-las. Para isso, vamos utilizar o objeto HierarchicalEvaluation. O objetivo é comparar as previsões com os dados reais do conjunto de teste e calcular o erro quadrático médio (MSE).

AutoARIMA AutoETS Naive AutoARIMA/BottomUp AutoETS/BottomUp Naive/BottomUp AutoARIMA/TopDown_method-forecast_proportions AutoETS/TopDown_method-forecast_proportions Naive/TopDown_method-forecast_proportions AutoARIMA/MinTrace_method-ols AutoETS/MinTrace_method-ols Naive/MinTrace_method-ols
level metric
Overall mse-scaled 0.785841 0.876167 1.0 0.818236 0.859989 1.0 0.772272 0.888246 1.0 0.783855 0.875968 1.0
top_level mse-scaled 0.738763 0.87002 1.0 0.798764 0.840054 1.0 0.738763 0.87002 1.0 0.751004 0.862252 1.0
top_level/bottom_level mse-scaled 0.841087 0.883381 1.0 0.841087 0.883381 1.0 0.811595 0.909635 1.0 0.822405 0.892064 1.0

 

Podemos também avaliar as previsões visualmente, comparando os dados de previsão com os dados de teste.

 

Conclusões

A modelagem preditiva de séries temporais hierárquicas é uma ferramenta poderosa para analisar e prever dados que se estruturam em diferentes níveis de granularidade. Com a utilização das bibliotecas Python e dos métodos descritos neste artigo, você pode aplicar essa técnica para obter insights e tomar decisões mais precisas em seus projetos.

Quer aprender mais?

Clique aqui para fazer seu cadastro no Boletim AM e baixar o código que produziu este exercício, além de receber novos exercícios com exemplos reais de análise de dados envolvendo as áreas de Data Science, Econometria, Machine Learning, Macroeconomia Aplicada, Finanças Quantitativas e Políticas Públicas diretamente em seu e-mail.

Compartilhe esse artigo

Facebook
Twitter
LinkedIn
WhatsApp
Telegram
Email
Print

Comente o que achou desse artigo

Outros artigos relacionados

Coletando dados para monitoramento climático com Python

As condições climatológicas influenciam desde a safra de grãos até a decisão de um vendedor ambulante levar seu carrinho para a praia ou não. Por sua importância e impactos na economia do país, neste exercício mostramos como coletar e elaborar análises de dados sobre o clima usando o Python.

Coletando dados de secas e queimadas no Brasil com Python

Neste artigo exploramos fontes públicas de dados sobre secas e queimadas no Brasil. Mostramos como acessar, coletar e preparar os dados para elaboração de análises. Usamos a linguagem Python para desenvolver uma rotina automatizada.

Como analisar demonstrações contábeis usando IA

Neste post, vamos explorar como utilizar o modelo de linguagem Gemini do Google para analisar demonstrações contábeis anuais da Eletrobras e extrair informações relevantes para tomada de decisão. Através de um código Python, vamos importar os dados direto da CVM, conectar com o Gemini e gerar resumos sobre as contas das demonstrações e perspectivas futuras sobre as finanças da empresa.

como podemos ajudar?

Preencha os seus dados abaixo e fale conosco no WhatsApp

Boletim AM

Preencha o formulário abaixo para receber nossos boletins semanais diretamente em seu e-mail.