Projekcije širjenja COVID-19 v Sloveniji

Napoved 18.10.2020

Misel meseca: Slovenci smo športni fenomen, bodimo to tudi pri Covidu in pokažimo svetu, da znamo stopiti skupaj ter virus premagati in biti tudi pri tem najboljši.

Uvod

Uvodna pojasnila so podana tukaj. Za radovedne in željne znanja smo pripravili Malo šolo epidemije: epizoda 1, epizoda 2, epizoda 3.

Podatke pridobivamo od COVID-19 Sledilnika in se jim najlepše zahvaljujemo. Brez teh obsežnih kvalitetnih podatkov, ki jih vsakodnevno pretočimo v svoj model, se ne bi dalo napovedati nič pametnega, tako pa so naše napovedi izredno točne (z ozirom na področje napovedovanja razvoja epidemije).

Analiza stanja

Ocenjeno reprodukcijsko število R je 2,10 (podvojitveni čas dnevnih pozitivnih testov 6,5 dni ob 4 uvoženih/dan). Prvič po dveh tednih je R nekoliko padel in posledično podvojitveni čas nekoliko narastel, kar nakazuje, da bodo ukrepi z 9.10., katerih vpliv naj bi v prihodnjih dneh naraščal, morebiti uspeli kompenzirati naraščajoče posledice izgube sledenja stikov zaradi presežene kapacitete epidemiološkega nadzora. ÄŒe bi uspeli ukrepi z 9.10. pomembno zmanjšati R, bi bila to zelo dobra novica, kajti potem bi imeli ukrepi z 16.10. lažje delo pri preostalem spuščanju R pod 0,95 (ob 4 uvoženih/dan), potrebnem za obvladovanje epidemije. Je pa današnji delež pozitivnih testov izredno visok, kar 19,3%, in bi to lahko pomenilo, da je dnevno število pozitivnih testov v primerjavi s prejšnjimi dnevi podcenjeno in da posledično še vedno drvimo v smeri naraščajočega R in krajšanja podvojitvenega časa. Kaj se dogaja z R in pri kateri vrednosti se bo ustalil, se bo izkristaliziralo v naslednjih dneh. Dokler se R ne bo ustalil, bo negotovost podanih napovedi velika. ÄŒe bo današnje dnevno število pozitivnih testov (število bo objavljeno jutri) manjše od 325, bo R nekoliko padel, podvojitveni čas pa nekoliko narastel, sicer bo ravno obratno.

ÄŒe družba kot celota ne bomo uspeli bistveno zmanjšati širjenja Covida z odgovornim ravnanjem vseh, ne vidimo razloga, zakaj razvoj epidemije ne bi sledil predstavljenim neugodnim krivuljam. ÄŒe epidemije ne bomo uspeli sami pravočasno ustaviti, nas bo ustavila ona, kajti naraščajoče eksponentne funkcije se ne da premagati. Upam, da bomo tokrat dovolj zreli in pametni, da se bomo znali sami spraviti dol z eksponentno naraščajoče krivulje, drugače nam bo morala spet pomagati država z ukrepi. Zato dosledno izvajajmo samozaščitne ukrepe, vsi in vedno. Epidemija bo začela pojenjati, če bo R manjši od 0,95 (ob 4 uvoženih/dan). ÄŒuvajte se, kajti verjetnost, da srečate kužno osebo, je že znatna. (slika 1, slika 3, slika 4, slika 5)

Ocenjeno sedemdnevno povprečje efektivnega reprodukcijskega števila je R_ef_7dni=2,11 (včeraj 2,14) s podvojitvenim časom dnevnih pozitivnih testov 6,5 dni (včeraj 6,4 dni). Ocenjeno sedemdnevno povprečje reprodukcijskega števila je R_7dni=2,10 (včeraj 2,12). Sedemdnevno povprečje uvoženih pozitivnih testov je približno 4 uvožene/dan (prejšnji teden isto). (slika 3, slika 4)

Napoved razvoja epidemije

Ob (optimistični) predpostavki, da ukrepi 16.10. polno primejo in zaustavijo širjenje epidemije (reprodukcijsko število pod 0,9), se nakazujejo naslednje maksimalne številke ob navedenih datumih:

• 950 dnevnih pozitivnih testov (z dnevno fluktuacijo do 1800) 22.10., potem počasen upad (včeraj 980, 22.10.) (slika 1)
• 693 hospitaliziranih 4.11., potem počasen upad (včeraj 729, 4.11.) (slika 2)
• 190 na intenzivni negi 7.11., potem počasen upad (včeraj 195, 7.11) (slika 2)

Bolj kot bodo ukrepi z 9.10. uspeli spustiti R (če ga bodo sploh uspeli spustiti), nižje bodo te maksimalne številke.

Kriteriji za 2. paket rdečega seznama ukrepov v 2. valu bodo predvidoma preseženi ob sledečih datumih:

• >170 14/100.000: že preseženo 12.10. (slika 1)
• >300 hospitaliziranih: 20.10. (slika 2)
• >50 na intenzivni: 19.10. (slika 2)

V napovedi za pozitivne teste (slika 1) in bolnice (slika 2) smo predpostavili, da ukrepi (z "ukrepi" so mišljeni ukrepi in njihovo izvajanje) z 9.10. ne bodo uspeli zmanjšati R, da pa bodo ukrepi z 16.10. dovolj učinkoviti in bodo uspeli R spustiti na 0,9, do takrat pa bo epidemija sledila trenutni dinamiki z R=2,12 (pri tej vrednosti dobimo najboljše ujemanje s podatki) in 4 uvoženimi/dan. V tem primeru bo število dnevnih pozitivnih testov v prihodnjem obdobju fluktuiralo v območju od 200 do 1500 okoli trenutne povprečne vrednosti 771, ki bo sprva še naraščala do 22.10. (na sliki je zaradi sedemdnevnega povprečenja vrh že nekoliko prej) in nato začela padati. V prognozi za bolnice so predstavljeni tudi scenarij za malo verjeten primer, da bodo ukrepi z 9.10. prijeli (R=0,9), scenarij za primer, da tudi novi ukrepi ne bodo prijeli (po 16.10. še naprej R=2,12), ter predvidoma konservativen scenarij z R=2,5 po 9.10. Namen predstavitve različnih scenarijev je, da se seznanimo z negotovostmi in karakterističnimi časi.

Še prognoza na osnovi epidemioloških podatkov (stanje 11.10.2020), ki jih tedensko objavlja NIJZ. S slike 6 je razvidno, da naraščajo vse kategorije razen "povezani z importiranimi". Vse kategorije razen "importirani" in "povezani z importiranimi" naraščajo še vedno vse hitreje (krivulje se krivijo v smeri navzgor), kar nakazuje, da imamo eksponentno rast, kar je slabo.

Iz tabele 1 (stanje 15.10.2020) je razvidno, da je smrtnost po starostnih kategorijah v drugem valu epidemije glede na pozitivne teste približno 50% manjša kot v prvem valu (upoštevan je tipičen dvotedenski časovni zamik od pozitivnega testa do smrti), predvidoma zaradi boljšega zdravljenja.

Slika 1: Z "ukrepi" so mišljeni ukrepi in njihovo izvajanje

Slika 2: Z "ukrepi" so mišljeni ukrepi in njihovo izvajanje

 Slika 3:

 Slika 4: (Stanje 14.10.2020, zaradi preseženih kapacitet epidemiološke službe novejši podatki niso na razpolago)

 Slika 5:

Slika 6: (Stanje 11.10.2020)

Tabela 1: (Stanje 15.10.2020)