Skip to content

📈 Previsão do Índice S&P 500 Utilizando LSTM e Mecanismos de Atenção

License

Notifications You must be signed in to change notification settings

gvmossato/alstm-stock-market

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

38 Commits
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Repository files navigation


alstm-logo

python-badge tensoflow-badge python-badge


$${\color[rgb]{1,0.75,0}\text{\Large Previsão do Índice S\&P 500 Utilizando LSTM e Mecanismos de Atenção}}$$

Projeto desenvolvido como Trabalho de Conclusão de Curso durante o último ano de graduação em Engenharia Mecatrônica na Escola Politécnica da Universidade de São Paulo (EP-USP)



View in English (TODO)

🔎 Sumário

Neste repositório:

📜 Sobre

History never repeats itself, but it does often rhyme
Mark Twain

Esse repositório contém o código do estudo para previsão do índice S&P 500, para o qual foi desenvolvido um modelo utilizando células de memória de longo-curto prazo (LSTM) combinadas com mecanismos de atenção. O modelo passou por otimizações com técnicas de grid search e bayesian seach, demonstrando um desempenho promissor na previsão dos preços de fechamento.

O trabalho foi majoritariamente inspirado no artigo Forecasting stock prices with long-short term memory neural network based on attention mechanism (2020) de Jiayu Qiu, Bin Wang e Changjun Zhou, recebendo honras como um dos cinco melhores projetos de 2023 no curso de Engenharia Mecatrônica.

No estudo também foram exploradas aplicações mais práticas do modelo através de técnicas de gestão de banca, avaliando a rentabilidade das previsões em um ambiente controlado e com resultados igualmente promissores.

⁉️ Problemática & Motivação

O aprofundamento financeiro do Brasil, marcado pela liberalização, aumento de crédito e bancarização da população, catalisou um aumento expressivo no número de participantes do mercado financeiro, incluindo também pequenas gestoras e investidores individuais. Entretanto, esse crescimento acompanha também uma ascensão da desinformação, exacerbada pela difusão de conteúdo não verificado em mídias sociais e demais meios.

decision-making-process         fraud

Ao nível de confiança de 95%, o grupo vítimas possuiu proporção (...) significativamente maior apenas em criptomoedas e na opção não tenho investimentos.

Essa dinâmica acaba revelando um recorte populacional que é muitas vezes inexperiente com investimentos e adentra esse universo diretamente através de mercados muito sofisticados (como o de criptomoedas) ou então que ficam expostos a práticas predatórias justamente em virtude da baixa familiaridade com o ambiente financeiro.

Delineado esse contorno, entende-se para a concepção desse projeto que é necessário empoderar o processo de tomada de decisão de investimento dos indivíduos, em especial, daqueles recém chegados ao mercado. Tal cenário então abre espaço para que modelos preditivos possam emergir como auxiliares de investimentos, utilizando informações críticas e atualizadas capazes de influenciar diretamente as tendências do mercado. Nesse contexto, é feito uso do S&P 500: esse é, para além de um índice de ações estadunidenses, um indicador da economia global devido à vasta operação internacional das empresas nesse listadas.

A US investor gets a lot of international exposure by either investing in S&P 500 index funds or actively investing with the S&P 500 as the benchmark. (De, 2013)

A motivação desse trabalho reside então na necessidade de fortalecer o processo decisório no contexto complexo e volátil do mercado de ações. Sob a qual se propõe um modelo que busca ser uma ferramenta auxiliar, não substituindo, mas sim potencializando, o raciocínio estratégico de investidores. Para tal, empresta-se dos recentes desenvolvimentos em LLMs uma simplificação dos mecanismos de atenção para aqui serem aplicados a séries temporais financeiras, provendo especialmente a investidores com recursos limitados e acesso tardio a informações, uma fonte auxiliar para suas decisões.

⚙️ Dados & Pré-Processamento

Fluxograma simplificado do pré-processamento dos dados até serem consumidos pelo modelo:

dados

Conforme a imagem:

  1. Fonte: o conjunto de dados utilizado nesse projeto é um recorte histórico abrangendo de 03 de janeiro de 1983 a 01 de setembro de 2023, retirado do Yahoo! Finance.

  2. Redução de Ruído: transformada wavelet com a família de funções Coiflets até terceira ordem, a qual é particularmente eficaz na redução de ruído em sinais não estacionários como os de preços de ações.

  3. Normalização: para assegurar que todas as variáveis tenham o mesmo peso durante o treinamento da rede, aplicamos a normalização Z-Score. Isso coloca todas as variáveis na mesma escala, neutralizando o efeito de disparidades nas magnitudes de preços e volumes transacionais.

  4. Repartição: os dados são divididos em três segmentos: treino, validação e teste, na proporção de 95/2.5/2.5. Escolhemos essa divisão para permitir um ajuste fino dos hiperparâmetros (usando o conjunto de validação) e uma avaliação honesta do desempenho do modelo (usando o conjunto de teste). Ao manter a ordem temporal (sem embaralhamento), respeitamos a sequência natural dos eventos no mercado de ações.

  5. Janelamento: implementamos uma técnica de janela deslizante com tamanho de 20 dias (aproximadamente um mês em dias úteis), que constitui o passo temporal do nosso modelo. Dentro dessa janela, seis séries temporais distintas — abertura, fechamento, máxima, mínima, fechamento ajustado e volume (em quantidades) — são fornecidas ao modelo. Com esses dados ele prediz o valor de fechamento no vigésimo primeiro dia.

🦾 Otimização

Como citado, realizamos a tunagem de hiperparâmetros em duas fases distintas: grid search e bayesian search. No primeiro, buscamos explorar deterministicamente as redondezas do modelo apresentado no artigo base utilizado. Assim, com um total de 48 combinações únicas de hiperparâmetros avaliadas utilizando a metodologia de validação cruzada com k-dobras obtivemos 144 execuções distintas (k=3).

Na sequência, o bayesian search foi implementado com base nas combinações de hiperparâmetros mais promissoras do grid search, em uma tentativa de refinar a rede. Nesse ponto, decidimos fixar o tamanho do estado oculto em 20, baseando-nos nas descobertas do grid search e na sintonia com o tamanho de entrada de dados (20 dias úteis). Ao longo de 100 configurações de redes, avaliamos e validamos diferentes modelos, totalizando 300 execuções.

Um resumo dos testes encontra-se na tabela abaixo:

Hiperparâmetro Grid Search Bayesian Search
Tamanho do Estado Oculto 10, 20, 50, 100 Não testado
Taxa de Aprendizado 0,001; 0,01; 0,1 0,0001 a 0,01 (0,00018)
Tamanho do Lote 64, 128, 256, 512 64, 128, 256, 512, 1024
Taxa de Dropout Não testado 0% a 30% (12,241%)

🧠 Rede & Treinamento

A arquitetura proposta para o modelo de previsão do índice S&P 500 incorpora um total de 20 células LSTM, isto é, para uma janela móvel de 20 dias são lidos OHLC, fechamento ajustado e volume (em quantidades), para então se predizer o fechamento do 21º dia.

modelo

Como mostra a imagem anterior, após o processamento pelas células LSTM, as saídas são submetidas ao mecanismo de soft attention. Esse mecanismo avalia as contribuições de cada célula LSTM e pondera sua influência, permitindo que a importância de momentos distintos no tempo seja diferenciada, ao invés de focar apenas na informação mais recente. Isso se baseia na premissa de que eventos passados dentro da janela de tempo podem ter relevância semelhante ou até maior do que os mais recentes.

A "camada de atenção" então agrega as saídas das células LSTM ponderadas em um vetor de contexto que concentra as informações relevantes detectadas pela rede. Esse vetor de contexto é passado por uma camada de dropout, com uma taxa de desativação de 12,241%, antes de ser apresentado à última camada da rede.

A etapa final do modelo é composta por uma camada densa com um único neurônio, cuja função de ativação linear é adequada para tarefas de regressão como a previsão de índices de ações. Esse neurônio processa o vetor e produz o output final da rede: a previsão do valor de fechamento do S&P 500.

Por fim, o treinamento da rede ocorre ao longo de 2000 epochs com um batch size de 128. O algoritmo ADAM é o escolhido para otimização, operando com uma taxa de aprendizado de 0,00018 - esses parâmetros foram selecionados com base nos resultados da busca em grid, da busca bayesiana e da análise da curva de aprendizado.

📈 Resultados

Note

A seção de resultados se prolonga por algumas dezenas de páginas da monografia, então, não sendo pertinente trazer todos os resultados, abordamos aqui um recorte conveniente do que alcançamos com o modelo.

Qualitativamente, os resultados no conjunto de teste e o retorno acumulado ao longo do período são bastante fidedignos ao que se observou no mercado para à época:

test-results         return-results

Optamos ainda por fazer um estudo comparativo do mecanismo de atenção, testando variações desse: uma rede sem atenção, uma rede com a atenção como proposta no artigo de referência (benchmark) e a "atenção clássica", proposta no artigo inaugural Attention Is All You Need:

attention-results

É observável, portanto, a influência e consequente melhora do desempenho do modelo com o uso do mecanismo clássico. Quando passamos, todavia, a avaliar quantitativamente o modelo frente ao benchmark, encontramos um problema inicial:

Benchmark Modelo
Data mínima (teste) 2019-05-17 2022-09-23
Data máxima (teste) 2019-07-01 2023-09-01
Observações 31 237
Preço Mínimo 2751,53 3577,03
Preço Máximo 2971,41 4588,96
Amplitude de Preços 219,88 1011,93
Preço Médio 2870,27 4095,33
Volatilidade Anualizada 11,28% 17,58%

Como sintetiza a tabela anterior, o conjunto de teste do nosso modelo encontrava-se em um mercado muito mais complexo: mais observações com maiores amplitudes de preço e uma volatilidade notadamente superior. Assim, a comparação direta não poderia ser considerada justa, optamos então por normalizar a comparação das métricas de erro por dois métodos, preço médio e amplitude:

Benchmark Modelo
Métricas
RMSE 0,3475 19,5238
MAE 0,1935 13,9011
$R^2$ 0,8783 0,9940
Normalização pelo Preço Médio
RMSE 0,00012107 0,00476734
MAE 0,00006742 0,00339437
Normalização pela Amplitude de Preços
RMSE 0,00158041 0,01929363
MAE 0,00088003 0,01373717

Todavia, mesmo após a normalização, com exceção do $R^2$, não superamos o benchmark. Não obstante, esses resultados forneceram um insight valioso: embora não estejamos acertando adequadamente o preço de fechamento, estamos fazendo uma leitura muito satisfatória da tendência. Ora, vamos então aprofundar essa análise:

trend-results         confusion-matrix

Constatamos que o modelo de fato parece seguir muito bem as oscilações, prevendo com consistência quando o mercado irá subir ou cair:

  • Dado que o índice subiu, acertamos 84,72% das vezes.
  • Dado que o índice caiu, acertamos 80,43% das vezes.

Frente a esses resultados, optamos então por tentar validar o modelo em uma abordagem um pouco mais prática, donde surge a iniciativa de aplicar estratégias de gestão de banca para operar no mercado considerando as previsões. Foram exploradas diversas estratégias (Martingale, Paroli, D'Alembert, etc.) em um ambiente simulado simplificado, cujas hipóteses adotadas foram:

  1. Livre de custos
  2. Liquidez e volume suficientes no mercado
  3. Operações ao preço de fechamento
  4. Sem alavancagem
  5. Compra e short selling são igualmente complexos

Destarte, mediante a previsão do modelo e a estratégia de gestão escolhida, o investidor entra comprado ou vendido no ativo, ganhando ou perdendo consoante a variação do índice no período. Os resultados para uma das estratégias mais rentáveis — Paroli — pode ser visto no gráfico abaixo:

paroli-results

Finalmente, o quadro de resultados por estratégia de gestão de banca fica expresso por:

bet-results

Embora os resultados sejam promissores, em especial o de estratégias mais agressivas quanto ao retorno potencial, elas também revelam uma exposição exagerada desse operador, que acabou se beneficiado de um período de mercado particularmente favorável, mas que possivelmente não resistiria a momentos de crise.

Diante disso, concluímos que estratégias de gestão de risco moderadas, como a Propocional de 25%, são mais factíveis, tanto em retornos quanto em exposição, equilibrando melhor o potencial de lucro com a minimização de riscos.

🌎 Aplicação

Como um todo, a implementação do projeto pode ser segmentada entre duas grandes frentes:

  1. Modelo, que compreende basicamente a tudo que fora exposto até aqui, como a rede, o treinamento, as validações técnicas e práticas, etc.

  2. Aplicação, uma plataforma web, hospedada em um outro repositório dedicado, para proporcionar uma interface intuitiva e acessível o suficiente a fim de permitir que usuários comuns pudessem usufruir das predições do modelo sem conhecimento técnico em programação.

front-repo-banner

A aplicação permaneceu operante até meados de abril de 2024, sendo incorporados os dados mais recentes disponíveis à época a cada seis meses, em treinamentos incrementais automáticos com 15 epochs e integração ao IBM Cloudant, provedor do banco de dados NoSQL utilizado.

👨‍💻 Uso & Código

Para executar o modelo carregando pesos de sessões de treinamento anteriores (se necessário) ou definir testes para ajustes de hiperparâmetros, utilize:

poetry run model [-t {grid,bayes}] [-w]

Parâmetros opcionais:

  • -t, --tuning: especifica o método de otimização a ser executado. Se não especificado, o ajuste de parâmetros não será realizado. Aceita:

    • grid: utiliza o grid search para otimizar os parâmetros.

    • bayes: utiliza o bayesian search para otimizar os parâmetros.

  • -w, --load-weights: carrega os pesos salvos da sessão de treinamento mais recente.

    • Padrão: False (não carrega os pesos automaticamente).

Warning

As configurações para cada tipo de ajuste devem ser definidas diretamente no código, em ./alstm_stock_market/run.py


Para executar a aplicação para realizar previsões com os dados mais recentes disponíveis, sincronização com a nuvem e treinamentos incrementais automáticos (se necessário), utilize:

poetry run app

Note que para integração com o banco de dados será necessário especificar as varáveis de ambiente requeridas pelo serviço de nuvem: DATABASE, CLOUDANT_USERNAME, CLOUDANT_PASSWORD e CLOUDANT_HOST.


Já com relação ao código em si, a árvore de arquivos do projeto está organizada como:

📦alstm_stock_market
 ┣ 📂images
 ┣ 📂logs
 ┣ 📂src
 ┃ ┣ 📂app
 ┃ ┃ ┗ 📜app.py
 ┃ ┣ 📂data
 ┃ ┃ ┗ 📜preprocessor.py
 ┃ ┣ 📂helpers
 ┃ ┃ ┣ 📂calendars
 ┃ ┃ ┃ ┗ 📜us.cal
 ┃ ┃ ┣ 📜plotter.py
 ┃ ┃ ┗ 📜utils.py
 ┃ ┣ 📂manager
 ┃ ┃ ┣ 📜manager.py
 ┃ ┃ ┗ 📜strategies.py
 ┃ ┣ 📂model
 ┃ ┃ ┣ 📂weights
 ┃ ┃ ┣ 📜evaluator.py
 ┃ ┃ ┣ 📜model.py
 ┃ ┃ ┗ 📜params.py
 ┗ 📜run.py

Sua subdivisão em módulos concentra as distintas operações do código em diretórios dedicados:

  • 📂src/model/: contém os arquivos referentes ao modelo em si, como arquitetura, hiperparâmetros da rede e métricas de avaliação. Na subpasta weights encontram-se os arquivos .h5 com os pesos do modelo após treinamentos.

  • 📂src/app/: contém a lógica que permite ao modelo ser executado em produção, conforme trabalhado na seção 🌎 Aplicação. As rotinas de treinamentos incrementais e comunicação com a nuvem — IBM Cloudant — também se encontram aqui.

  • 📂src/manager/: contém os arquivos referentes à gestão de banca. Subdivisão do código implementada para avaliar o desempenho do modelo em um cenário ainda controlado, mas mais próximo da prática, operando com distintas estratégias frente às previsões.

  • 📂src/data/: contém os arquivos referentes a todo o pipeline de dados exposto, capaz de lidar com cada um dos casos de uso (treinamento inicial, treinamentos adicionais, uso em produção, etc.).

  • 📂src/helpers/: contém os arquivos gerais e de uso compartilhado entre os demais módulos, funções e métodos auxiliares.


Por fim, a pasta images é utilizada para salvar os plots em formatados vetorizados, se desejável, enquanto logs armazena registros de execução do módulo de aplicação.

🤝 Colaboradores

Este projeto foi desenvolvido por Gabriel Mossato em colaboração com Paulino Fonseca, ambos à época graduandos sob orientação do Prof. Dr. Oswaldo Luiz do Valle Costa, pertencente ao Departamento de Engenharia Elétrica da Escola Politécnica da Universidade de São Paulo (EP-USP).


⬆️ Voltar ao Início

About

📈 Previsão do Índice S&P 500 Utilizando LSTM e Mecanismos de Atenção

Topics

Resources

License

Stars

Watchers

Forks

Languages