Projekcije širjenja COVID-19 v Sloveniji
Napoved 20.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,15 (podvojitveni čas dnevnih pozitivnih testov 6,3 dni ob 4 uvoženih/dan) in je skoraj enako kot včeraj, kar pomeni, da ukrepi z 9.10., katerih vpliv bi se moral že močno poznati, očitno ne bodo uspeli spustiti R, kar je slaba novica. K sreči smo z ukrepi 16.10. že postavili novo obrambno linijo in z današnjimi ukrepi 20.10. postavljamo še eno. Upam, da bosta zdržali, kajti pritisk virusa je enormen, podvojitveni časi pa kratki, zato se bolnice vse hitreje polnijo. Vpliv ukrepov 16.10. se naj bi začel poznati konec tedna, vpliv ukrepov 20.10. pa naslednji teden. ÄŒe bo današnje dnevno število pozitivnih testov (število bo objavljeno jutri) manjše od 1415, 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. 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 znatna. (slika 1, slika 3, slika 4, slika 5)
Ocenjeno sedemdnevno povprečje efektivnega reprodukcijskega števila je R_ef_7dni=2,17 (včeraj 2,17) s podvojitvenim časom dnevnih pozitivnih testov 6,3 dni (včeraj 6,3 dni). Ocenjeno sedemdnevno povprečje reprodukcijskega števila je R_7dni=2,15 (včeraj 2,16). Sedemdnevno povprečje uvoženih pozitivnih testov je približno 4 uvožene/dan (prejšnji teden isto). (slika 3, slika 4)
Ocenjen delež okuženih in kužnih je sledeč:
V prvem valu se je okužilo približno 1% prebivalstva, kar pomeni približno vsak 100. prebivalec. Na vrhuncu je bil delež kužnih približno 0,35% prebivalstva, par pomeni približno vsak 300. prebivalec. V drugem valu se je do zdaj okužilo približno 2,5% prebivalstva, kar pomeni približno vsak 40. prebivalec. Trenutno je delež kužnih približno 1,5% prebivalstva, kar pomeni približno vsak 70. prebivalec, in hitro raste. Do zdaj se je v obeh valovih skupaj okužilo približno 3,5% prebivalstva, kar pomeni približno vsak 30. prebivalec.
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:
- 980 dnevnih pozitivnih testov (z dnevno fluktuacijo do 1800) 22.10., potem počasen upad (včeraj 980, 22.10.) (slika 1)
- 684 hospitaliziranih 4.11., potem počasen upad (včeraj 699, 4.11.) (slika 2)
- 181 na intenzivni negi 7.11., potem počasen upad (včeraj 187, 7.11) (slika 2)
Kriteriji za 3. paket rdečega seznama ukrepov v 2. valu bodo predvidoma preseženi ob sledečih datumih:
- >360 hospitaliziranih: 22.10. (slika 2)
- >60 na intenzivni: 21.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,16 (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 300 do 1800 okoli trenutne povprečne vrednosti 959, 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,16), 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)