Vidas salvas no Brasil pelo isolamento social (v. 2.1)

Desenvolvido por Paulo J. S. Silva e Claudia Sagastizábal após uma discussão com colegas do CeMEAI

Agradecimentos: Tiago Pereira e Alexandre Delben

O objetivo desta página é apresentar uma estimativa do número de vidas salvas no país pelo isolamento social. Para isso nós fazemos ajustes do parâmetro $R_0$ do modelo SEIR, que representa a taxa de replicação do vírus SARS-CoV-2 (o corona vírus que causa a Covid-19), tentando descobrir se ele varia no tempo. A ideia é buscar identificar tendências na evolução da taxa de propagação do vírus e consequente aceleração ou desaceleração da epidemia depois do início dos protocolos de distanciamento social que foram implementados a partir de 24 de março em vários lugares do país. Quando os dados anteriores ao dia 24 de março não são suficientes para fazer a análise, nós usamos os dados da primeira semana em que conseguimos estimar o comportamento do vírus para representar o período anterior ao isolamento social.

Apresentamos inicialmente os dados para o país todo e depois especializamos os resultados para os estados de São Paulo, Rio de Janeiro e Maranhão (estado de origem do primeiro autor) e para todas as grandes regiões do país. Vamos ficar em casa, vamos salvar vidas.

Obs: estimativas feitas com base nos dados oficiais, veja mais detalhes abaixo, em particular detalhamos a mudanças feitas com relação à versão anterior no final desse documento.

Estimativa de vidas salvas pelo distanciamento social nos próximos 14 dias

Dados atualizados até o dia 30/05/2020.

Os dados a seguir mostram que se mantivermos o isolamento feito na última
semana pelas próximas duas semanas pouparemos uma vida a cada 0.6 minutos.

Já Por região, o tempo em minutos necessário para salvar uma vida será de
1.4 (Norte), 1.9 (Nordeste), 125.2 (Centro_Oeste), 1.1 (Sudeste), 77.5 (Sul).

Total de vidas salvas por dia

Out[4]:
31/05 01/06 02/06 03/06 04/06 05/06 06/06 07/06 08/06 09/06 10/06 11/06 12/06 13/06
Brasil 457 606 778 978 1207 1473 1777 2126 2527 2982 3503 4098 4773 5542
SP 149 198 255 323 400 492 598 719 861 1025 1215 1434 1686 1977
RJ 80 106 138 173 215 263 319 383 456 540 636 746 872 1015
MA 25 33 44 57 72 91 113 139 170 207 252 303 363 434
Norte 161 218 285 367 463 577 711 870 1056 1274 1529 1827 2172 2574
Nordeste 148 196 251 316 390 476 574 685 814 961 1128 1318 1534 1780
Centro_Oeste 2 4 4 5 6 7 9 11 12 15 17 20 23 26
Sudeste 244 322 415 522 644 787 951 1140 1356 1605 1890 2215 2588 3014
Sul 4 5 6 8 10 12 15 18 20 24 27 32 37 42

Obs1: Os resultados estão baseados nos dados oficiais baixados a partir do site Observatório Covid 19. Esses dados sofrem de clara subnotificação e assim as nossas estimativas serão também sub-estimadas. Porém, acreditamos que mesmo assim é possível ter uma ideia da evolução da epidemia e assim ser útil.

Obs2: Devido a essa dificuldade com os dados oficiais, preferimos fazer uma estimativa para alguns dias no futuro, partindo dos dados atuais, a tentar manter informação sobre o total de vidas salvas até hoje. Isso porque estaríamos acumulando informação com imprecisão o que só iria amplificar as incertezas. Pelo mesmo motivo apresentamos abaixo os gráficos de evolução futura por mais 30 dias apenas.

Comentários

O distanciamento social parece ter sido efetivo quando consideramos o Brasil inteiro e, após perder um pouco de força há umas semanas, parece ter melhorado levemente recentemente. Cabe destacar que ele tem perdido força no Centro-Oeste e Sul pnde perderam boa parte do ganho obtido inicialmente. O Norte e Nordeste, que já possuem regiões onde o sistema de UTIs está acima da capacidade ou muito próximo disso, parecem ter entendido a dimensão do problema e passaram a adotar um distanciamento mais efetivo. Cabe, porém, destacar que como estamos usando dados oficiais. Deste modo a desaceleração aparente da transmissão do vírus pode refletir em parte uma dificuldade para testar os casos e levantar dados.

Por fim, vemos que as curvas de uma forma geral foram achatadas. Mas os números de doentes ainda cresce muito, mesmo que mais lentamente. Isso sugere que é imperativo que os governos busquem alternativas de controle da epidemia para não enfretarmos colapsos nos sistemas de saúde em breve.

Resultados Detalhados

Brasil

Essa é uma estimativa do número de vidas que serão salvas no país por termos mantido o isolamento social na última semana e se continuarmos firmes nos próximos 14 dias. Veja que o número cresce muito rapidamente à medida que os dias passam, enfatizando a necessidade de medidas de mitigação.

Vidas que serão salvas nos próximos 14 dias = 457 (31/05), 606 (01/06), 778 (02/06), 978 (03/06), 1207 (04/06), 1473 (05/06), 1777 (06/06), 2126 (07/06), 2527 (08/06), 2982 (09/06), 3503 (10/06), 4098 (11/06), 4773 (12/06), 5542 (13/06).

Veja a evolução da taxa de reprodução do vírus e como ela caiu depois do início do isolamento social. Veja também como já achatamos a curva. Já há um grande ganho. Assim, ainda há muito o que fazer e a janela de oportunidade fica cada vez mais estreita.

Out[6]:
R0 antes do início da quarentena              = 2.27
R0 na última semana                           = 1.49
Estimativa de óbitos (com taxa de reprodução da última semana) =
    [1186, 1247, 1312, 1379, 1451, 1525, 1604, 1687, 1773, 1865, 1961, 2061, 2167, 2278]
Estimativa de óbitos (com taxa de reprodução de antes da quarentena) =
    [1643, 1853, 2090, 2357, 2658, 2998, 3381, 3813, 4300, 4847, 5464, 6159, 6940, 7820]

São Paulo

Essa é uma estimativa do número de vidas que serão salvas no estado de São Paulo por termos mantido o isolamento social na última semana e se continuarmos firmes nos próximos 14 dias. Veja que o número cresce rapidamente à medida que os dias passam.

Vidas que serão salvas nos próximos 14 dias = 149 (31/05), 198 (01/06), 255 (02/06), 323 (03/06), 400 (04/06), 492 (05/06), 598 (06/06), 719 (07/06), 861 (08/06), 1025 (09/06), 1215 (10/06), 1434 (11/06), 1686 (12/06), 1977 (13/06).

Veja a evolução da taxa de reprodução do vírus e como ela caiu depois do início do isolamento social. Veja também como já achatamos a curva. Já há um grande ganho seguido de um período de leve piora que foi depois revertido. Mas o número de caso ainda não diminui e portanto ainda há muito o que fazer.

Out[9]:
R0 antes do início da quarentena              = 2.51
R0 na última semana                           = 1.29
Estimativa de óbitos (com taxa de reprodução da última semana) =
    [221, 228, 235, 242, 250, 257, 265, 274, 282, 291, 300, 309, 319, 329]
Estimativa de óbitos (com taxa de reprodução de antes da quarentena) =
    [370, 426, 490, 565, 650, 749, 863, 993, 1143, 1316, 1515, 1743, 2005, 2306]

Rio de Janeiro

Essa é uma estimativa do número de vidas que serão salvas no estado do Rio de Janeiro por termos mantido o isolamento social na última semana e se continuarmos firmes nos próximos 14 dias. Veja que o número cresce rapidamente à medida que os dias passam.

Vidas que serão salvas nos próximos 14 dias = 80 (31/05), 106 (01/06), 138 (02/06), 173 (03/06), 215 (04/06), 263 (05/06), 319 (06/06), 383 (07/06), 456 (08/06), 540 (09/06), 636 (10/06), 746 (11/06), 872 (12/06), 1015 (13/06).

Veja a evolução da taxa de reprodução do vírus e como ela caiu depois do início do isolamento social. Veja também como já achatamos a curva. Já há um grande ganho, mas houve uma pequena piora nos últimos dias e ainda há muito o que fazer.

Out[12]:
R0 antes do início da quarentena              = 2.26
R0 na última semana                           = 1.59
Estimativa de óbitos (com taxa de reprodução da última semana) =
    [252, 268, 284, 302, 321, 341, 362, 385, 409, 435, 462, 491, 521, 553]
Estimativa de óbitos (com taxa de reprodução de antes da quarentena) =
    [332, 374, 422, 475, 536, 604, 681, 768, 865, 975, 1098, 1237, 1393, 1568]

Maranhão

Essa é uma estimativa do número de vidas que serão salvas no estado do Maranhão por termos mantido o isolamento social na última semana e se continuarmos firmes nos próximos 14 dias. Veja que o número cresce rapidamente à medida que os dias passam.

Vidas que serão salvas nos próximos 14 dias = 25 (31/05), 33 (01/06), 44 (02/06), 57 (03/06), 72 (04/06), 91 (05/06), 113 (06/06), 139 (07/06), 170 (08/06), 207 (09/06), 252 (10/06), 303 (11/06), 363 (12/06), 434 (13/06).

Veja a evolução da taxa de reprodução do vírus e como ela caiu depois do início do isolamento social. Veja também como já achatamos a curva. Já há um grande ganho, mas ainda há muito o que fazer pois a taxa de reprodução do vírus ainda está alta.

Out[15]:
R0 antes do início da quarentena              = 2.72
R0 na última semana                           = 1.67
Estimativa de óbitos (com taxa de reprodução da última semana) =
    [46, 50, 53, 57, 62, 66, 71, 76, 82, 88, 94, 101, 109, 117]
Estimativa de óbitos (com taxa de reprodução de antes da quarentena) =
    [71, 83, 97, 114, 134, 157, 184, 215, 252, 295, 346, 404, 472, 551]

Norte

Houve uma queda significativa na região norte, porém o patamar atual ainda é alto. É importante buscar formas de mitigação mais efetivas para que a situação não se deteriore ainda mais.

Vidas que serão salvas nos próximos 14 dias = 161 (31/05), 218 (01/06), 285 (02/06), 367 (03/06), 463 (04/06), 577 (05/06), 711 (06/06), 870 (07/06), 1056 (08/06), 1274 (09/06), 1529 (10/06), 1827 (11/06), 2172 (12/06), 2574 (13/06).

Veja a evolução da taxa de reprodução do vírus e como ela caiu depois do início do isolamento social. Veja também como já achatamos a curva. Já há um grande ganho, mas com o número crescente de casos e ainda há muito o que fazer.

Out[18]:
R0 antes do início da quarentena              = 2.70
R0 na última semana                           = 1.55
Estimativa de óbitos (com taxa de reprodução da última semana) =
    [273, 289, 306, 323, 342, 362, 383, 405, 429, 454, 480, 507, 537, 567]
Estimativa de óbitos (com taxa de reprodução de antes da quarentena) =
    [434, 507, 591, 690, 805, 939, 1094, 1275, 1485, 1728, 2009, 2334, 2709, 3141]

Nordeste

Essa é uma estimativa do número de vidas que serão salvas na região Nordeste por termos mantido o isolamento social na última semana e se continuarmos firmes nos próximos 14 dias. Veja que o número cresce rapidamente à medida que os dias passam.

Vidas que serão salvas nos próximos 14 dias = 148 (31/05), 196 (01/06), 251 (02/06), 316 (03/06), 390 (04/06), 476 (05/06), 574 (06/06), 685 (07/06), 814 (08/06), 961 (09/06), 1128 (10/06), 1318 (11/06), 1534 (12/06), 1780 (13/06).

Veja a evolução da taxa de reprodução do vírus e como ela caiu depois do início do isolamento social. Veja também como já achatamos a curva. Já há um grande ganho, mas ainda há muito o que fazer.

Out[21]:
R0 antes do início da quarentena              = 2.29
R0 na última semana                           = 1.47
Estimativa de óbitos (com taxa de reprodução da última semana) =
    [361, 378, 397, 416, 436, 457, 479, 503, 527, 552, 579, 607, 637, 667]
Estimativa de óbitos (com taxa de reprodução de antes da quarentena) =
    [509, 574, 648, 732, 826, 933, 1053, 1188, 1341, 1513, 1707, 1925, 2171, 2447]

Centro-Oeste

Essa é uma estimativa do número de vidas que serão salvas na região Centro-Oeste. O número é baixo porque houve um aumento recente da taxa de reprodução do vírus, espero que essa tendência seja revertida em breve.

Vidas que serão salvas nos próximos 14 dias = 2 (31/05), 4 (01/06), 4 (02/06), 5 (03/06), 6 (04/06), 7 (05/06), 9 (06/06), 11 (07/06), 12 (08/06), 15 (09/06), 17 (10/06), 20 (11/06), 23 (12/06), 26 (13/06).

Veja a evolução da taxa de reprodução do vírus e como ela caiu depois do início do isolamento social mas subiu de volta depois. Houve apenas um pequeno achatamento a curva.

Out[24]:
R0 antes do início da quarentena              = 1.93
R0 na última semana                           = 1.68
Estimativa de óbitos (com taxa de reprodução da última semana) =
    [21, 22, 24, 26, 28, 30, 32, 34, 37, 39, 42, 45, 48, 52]
Estimativa de óbitos (com taxa de reprodução de antes da quarentena) =
    [23, 26, 28, 31, 34, 37, 41, 45, 49, 54, 59, 65, 71, 78]

Sudeste

Essa é uma estimativa do número de vidas que serão salvas na região Sudeste por termos mantido o isolamento social na última semana e se continuarmos firmes nos próximos 14 dias. Veja que o número cresce rapidamente à medida que os dias passam.

Vidas que serão salvas nos próximos 14 dias = 244 (31/05), 322 (01/06), 415 (02/06), 522 (03/06), 644 (04/06), 787 (05/06), 951 (06/06), 1140 (07/06), 1356 (08/06), 1605 (09/06), 1890 (10/06), 2215 (11/06), 2588 (12/06), 3014 (13/06).

Veja a evolução da taxa de reprodução do vírus e como ela caiu depois do início do isolamento social. Veja também como já achatamos a curva. Já há um grande ganho seguida de uma piora bem leve que foi revertido últimos dias. Porém o número de casos ainda cresce e assim ainda há muito o que fazer.

Out[27]:
R0 antes do início da quarentena              = 2.31
R0 na última semana                           = 1.37
Estimativa de óbitos (com taxa de reprodução da última semana) =
    [488, 509, 530, 552, 576, 600, 625, 651, 679, 707, 737, 768, 800, 833]
Estimativa de óbitos (com taxa de reprodução de antes da quarentena) =
    [732, 831, 945, 1074, 1220, 1387, 1576, 1791, 2035, 2312, 2627, 2983, 3388, 3847]

Sul

Esse é número de vidas perdidas na região Sul com o aumento forte recente dos casos. Espero que possamos reverter isso.

Vidas que serão salvas nos próximos 14 dias = 4 (31/05), 5 (01/06), 6 (02/06), 8 (03/06), 10 (04/06), 12 (05/06), 15 (06/06), 18 (07/06), 20 (08/06), 24 (09/06), 27 (10/06), 32 (11/06), 37 (12/06), 42 (13/06).

Veja a evolução da taxa de reprodução do vírus e como ela caiu depois do início do isolamento social seguido de um período de quase estabilidade. Nesse momento passou-se a impressão que a epidemia estava sendo controlada, mas logo depois há uma subida forte, voltando para uma situação desconfortavelmente próxima à situação original. O número de infectados ainda é baixo e portanto ainda há como reverter de novo a situação, mas isso deve ser feito rapidamente.

Out[30]:
R0 antes do início da quarentena              = 1.96
R0 na última semana                           = 1.49
Estimativa de óbitos (com taxa de reprodução da última semana) =
    [18, 19, 20, 21, 22, 23, 24, 25, 27, 28, 30, 31, 33, 35]
Estimativa de óbitos (com taxa de reprodução de antes da quarentena) =
    [22, 24, 26, 29, 32, 35, 39, 43, 47, 52, 57, 63, 70, 77]

Observações técnicas

Esta página foi escrita com Jupyter e é executada diariamente para atualizar as informações. O código foi escrito em Julia e JuMP. Nós estimamos o $R_0$ ajustando os dados oficiais ao modelo SEIR, mas permitindo que o $R_t$ varie no tempo. Para evitar oscilações bruscas e overfitting, usamos uma regularização de variação total. O problema de otimização não linear é então resolvido usando o solver Ipopt.

Melhorias da versão 2.1

Introduzimos uma janela rolante para suavizar os dados diários, eles tem muita variação, principalmente nos finais de semana. Estamos usando uma janela de 14 dias (tem que ser múltiplo de 7 devido ao comportamento periódico nas semanas). Corrigir um erro que fazia o programa estimar mal o total de casos ativos que infectam (obrigado Tiago por chamar atenção para isso). Isso tinha pouca influência nos resultados porque estava sendo compensado na estimativa da taxa de mortalidade.

Melhorias da versão 2.0

A versão 2.0 incorpora alguma melhorias com respeito à versão original da página. Com o crescimento rápido do número de casos, a versão anterior tinha uma tendência a "esquecer" os dados antigos, pois os valores eram muito pequenos quanto comparados com os novos. Isso estava gerando uma subestimação da taxa de reprodução passada do vírus e, consequentemente, dava a impressão que o isolamento social estava sendo menos efetivo do que de fato tem sido em diminuir as taxas propagação da Covid-19. Para evitar esse problema, trocamos o ajuste baseado em erros absolutos por outro baseado em erros relativos. Vimos também que era necessário desconsiderar os primeiros dias da epidemia, pois a taxa de reprodução do vírus tende variar rápido demais. Isso deve ocorrer porque quando aparecem os primeiros doentes eles aparecem de forma pouco ordenada, aos "soluços", à medida que o sistema de testagem limpa a fila original.

Agradecimento

A pesquisa dos autores é apoiada pela FAPESP, através CEPID-CeMEAI (FAPESP 2013/07375-0) e Temático (2018/24293-0), e pelo CNPq (304301/2019-1 e 306089/2019-0).

Nota de isenção de responsabilidade

Esta página não é uma publicação oficial da Unicamp, seu conteúdo não foi examinado e/ou editado por esta instituição. A responsabilidade por seu conteúdo é exclusivamente do autor.

Disclaimer

This web site is not a official publication of Unicamp. Its content was not edited or supervised by the University. The responsibility for its contents lies entirely with the author.