U svojoj kolumni na Bug.hr dr. Igor Berecki, pedijatar iz Osijeka, posljednjih se tjedana bavi isključivo "koronavirusnim vremenima", kako to naziva, u kojima smo bombardirani vijestima, informacijama, podacima i gomilom brojki.
Većina tih numeričkih podataka kojima smo zatrpani orijentirana je na faktografsko nabrajanje onoga što se već dogodilo: koliko je do sada pozitivnih na SARS-CoV-2, koliko oboljelih od wuhanske groznice, koliko je pacijenata na respiratorima, koliko oporavljenih… A svi mi (sami za sebe ili u raspravama po društvenim mrežama) iz tih brojki pokušavamo donositi zaključke neposredno vezane uz ono što nas zaista zanima: kada ćemo moći u kafić, kada na šišanje, koliko još dana do odlaska na posao, a koliko do godišnjeg odmora na moru… "Ukratko - kada će sve ovo završiti?", pita se dr. Berecki u kolumni koju uz dopuštenje uredništva Buga prenosimo u cijelosti.
Formula umjesto podnaslova: s(t) + i (t) + r (t) = 1
Radujte se narodi jer je Društvo matematičara Amerike (MAA, Mathematical Association of America) izradio statistički model kojim je moguće s prilično velikom pouzdanošću (97% - 99%) procijeniti dinamiku epidemioloških događaja tijekom širenja koronavirusa unutar jedne države – od brzine širenja, preko vrhunca broja oboljelih, sve do jenjavanja epidemije i aproksimativnog datuma njezina završetka.
Svijet će se - prema matematičkom modelu - riješiti 99 % COVID-19 pandemije do 16. lipnja, u Njemačkoj će to biti 15. svibnja, a SAD može očekivati 99 % rješavanje svoje epidemije do 28. svibnja
Model je izrađen za predikciju razvoja COVID-19, kao pomoćni matematički alat za praćenje potvrđenih slučajeva infekcije. Nazvan je SIR, prema tri osnovne ulazne (ovisne) varijable: broju osjetljivih pojedinaca, broju zaraženih i broju oporavljenih (Susceptible, Infected, Recovered). Neovisna varijabla je vrijeme, mjereno u danima. Osnovna dinamika kretanja infekcije aproksimirana je iz podataka o nastanku, širenju i trajanju pandemije hongkonške gripe 1968./69., ponajviše zbog brojnih sličnosti u biologiji uzročnika, načina prenošenja i potencijala zaraznosti između SARS-CoV-2 (uzročnika COVIDA-19) i H3N2 soja A (uzročnika hongkonške gripe).
U modelu su korištena dva uzajamno povezana skupa ovisnih varijabli. Prvi skup su apsolutne vrijednosti (broj ljudi u svakoj grupi, kao funkcija vremena): S = S(t) je broj osjetljivih, I = I(t) je broj zaraženih, R = R(t) je broj oporavljenih pojedinaca. Drugi skup ovisnih varijabli je udio ukupne populacije u svakoj od tri promatrane kategorije; dakle, ako je ukupna populacija = N, tada imamo s(t) = S(t)/N kao osjetljiv udio populacije, i(t) = I(t)/N je zaraženi udio populacije te r(t) = R(t)/N kao oporavljeni udio populacije. U svakom trenutku (t), s(t) + i (t) + r (t) = 1.
Formula umjesto podnaslova: d(i) / d(t) = b*s(t)*i(t) - k*i(t)
Dakle, SIR-model svakoga pojedinca smatra: a) kandidatom za bolest ili b) aktivno bolesnom, zaraznom osobom bez obzira na to ima li ili nema vidljive simptome ili c) osobom koja više nije kandidat za obolijevanje niti je aktivno zarazna (bolesna). Pažljiviji čitatelj će zamijetiti pomalo morbidan detalj: model ne smatra umrle kao zasebnu skupinu. Sa statističko-infektološke točke gledišta, osoba koja je umrla od infekcije ubraja se u grupu oporavljenih (Recovered) jer više ne može oboljeti niti može proširiti bolest. Jedini način na koji pojedinac napušta ''osjetljivu'' skupinu jest ako se zarazi, a jedina dva načina kako prestaje biti ''zaražen'' (i zarazan!) jesu oporavak ili smrt.
U MAA-ovu modelu, vremenska stopa promjene broja osjetljivih S(t) ovisi o: broju pojedinaca koji su trenutno osjetljivi, broju pojedinaca koji su već zaraženi i količini kontakta između osjetljivih i zaraženih (ako pretpostavimo da svaki zaraženi pojedinac ima fiksni broj b kontakata dnevno koji su dovoljni za širenje bolesti). Naravno, nisu svi ti kontakti s osjetljivim pojedincima: ako pretpostavimo homogeno miješanje stanovništva, udio kontakata koji su osjetljivi je u određenoj vremenskoj točki = s(t).
Tako u prosjeku svaki zaraženi pojedinac dnevno stvori određeni broj novih zaraženih jedinki, što je nova vremenski ovisna varijabla b(t), brojka koju se potom matematički korigira faktorom infekcioznosti koronavirusa (ako ste imalo pratili izvješća Stožera i medijske informacije – a znam da jeste – već znate da je to broj R0, koji predstavlja broj osoba koje jedan infektivni pojedinac zarazi u skupini osjetljivih).
Uz logičnu pretpostavku da će se fiksni udio inficirane skupine s vremenom oporaviti (ili umrijeti, ali nećemo sada cjepidlačiti, bitno je da se infekcija ne širi dalje od tog pojedinca), u računu se pojavljuje i varijabla k(t), to jest broj oporavljenih iz inficirane skupine koji se generira u jednom danu. Vrijednosti parametara b i k nije moguće odrediti unaprijed, prije nego što se „s terena“ ne prikupe i ne obrade epidemiološka izvješća iz kojih se izračunava broj novozaraženih tijekom vremena – b(t) i broj oporavljenih u tom vremenu – k(t). No, već nakon desetak dana trajanja epidemije na nekom području te su brojke prilično jasne i stabilne pa ih je moguće uključiti u model i napraviti izračun.
Formula umjesto podnaslova: i(0) = 1.27 x 10- 6
Za usporedbu sa susjedima: Slovenija, Hrvatska i Srbija - pojavnost, broj oboljelih, oblik krivulje i predviđeno trajanje epidemije COVIDA-19
Matematičko-statističke zaljubljenike i one koji žele znati više upućujemo na stranice Mathematical Association of America na kojima je model detaljno razrađen, a nama koji više volimo da nam se sve nacrta preostaje pogledati grafove i krivulje koje je na temeljima MAA-ova modela izradio SUTD (Singapore University of Technology and Design). Već na prvi pogled, podaci koji svakodnevno pristižu iz epidemioloških službi brojnih država svijeta vrlo se uvjerljivo uklapaju u predikcijske krivulje izrađene SIR-modelom.
Kraj je bliži nego što smo se nadali?
A ako im je zaista vjerovati – kraj je bliži nego što smo se nadali! Prema izračunu koji su objavili, u Hrvatskoj će 97 % epidemije proći do 4. svibnja 2020., statistički vrlo uvjerljivih 99 % epidemije će nestati do 15. svibnja 2020., a stopostotni "oporavak" od COVID-19 možemo očekivati do 19. lipnja 2020.
O drugom valu pandemije, onom posebno teškom i nerijetko smrtonosnijem, kakav se redovito zamjećuje iza virusnih pošasti poput španjolske gripe i sličnih, a kakav predviđaju mnogi epidemiolozi i u slučaju COVID-19, za sada nema matematičkih modela. Razlog je jednostavan, a u statističko-matematičkom smislu vrlo složen: do trenutka prognozirane pojave mogućeg drugog vala očekuje se pojava i klinička primjena prvih učinkovitih lijekova i - što je još važnije - cjepiva protiv SARS-CoV-2 virusa. "Jer s dovoljnim postotkom procijepljenosti populacije, sve daljnje matematičke priče o širenju infekcije postaju suvišne", kaže na kraju svog teksta u Bug.hrdr. Igor Berecki.
