Uvod
Uporabljen je kalibriran model SEIR s petimi skupinami.
Shema razširjenega SEIR modela s tremi podskupinami populacije, kjer je možno mešanje med skupinami v oddelkih I. Model je opremljen z možnostjo izločanja cepljenih v oddelkih V. Model je možno razširiti na tak način do 5-ih skupin. Uteži in parametri modelov na shemi niso prikazani.
Nastavitve modela
Simulacija je izvedena na SEIR modelu s 5 starostnimi skupinami. Model je kalibriran na današnje podatke.
Starostne skupine
Starostne skupine so razdeljene glede na tabelo spodaj.
| skupina | število | delež |
|---|---|---|
| 0-24 | 504441 | 24.41% |
| 25-44 | 569430 | 27.56% |
| 45-64 | 596661 | 28.88% |
| 65-74 | 211160 | 10.22% |
| 75> | 184469 | 8.93% |
Bolnišnična obravnava je rezervirana za starostne skupine 45+ v razmerju (0.1, 1.0, 1.5).
Deleži smrti so predpostavljeni samo v skupini 4 in 5. To sta starostni skupini od 65+.
Ref: This represents an increase in deaths from 2.5 to 4.1 per 1000 detected cases. Risk of mortality in patients infected with SARS-CoV-2 variant of concern 202012/1: matched cohort study Robert Challen, Ellen Brooks-Pollock, Jonathan M Read, Louise Dyson, Krasimira Tsaneva-Atanasova, Leon Danon
Dodatno so predpostavljeni večji deleži hospitalizacij v tretji skupini in sicer iz uteži 0.1 na 0.25 ter deleži obravnave na icu iz uteži 0.0 na 0.8. To je zaradi angleške različice virusa, kjer se predvideva hujše poteke bolezni pri nižjih starostnih skupinah (od 45 naprej).
Matrike mešanja
Predpostavljeno je, da bodo sprostitveni ukrepi povečali stopnjo prenosa v starostni skupini otrok in mladostnikov ter delovne populacije. V skupini starostnikov od 65 let naprej se ne bodo spremenile stopnje prenosa.
Povečanje se zgodi od 09.02.2021 naprej.
Matriki stopenj prenosa sta prikazani v nadaljevanju.
| 1 | 1.00 | 1.00 | 1 | 1 |
| 1 | 1.15 | 1.15 | 1 | 1 |
| 1 | 1.15 | 1.15 | 1 | 1 |
| 1 | 1.00 | 1.00 | 1 | 1 |
| 1 | 1.00 | 1.00 | 1 | 1 |
| 1.50 | 1.00 | 1.00 | 1 | 1 |
| 1.25 | 1.15 | 1.15 | 1 | 1 |
| 1.20 | 1.15 | 1.15 | 1 | 1 |
| 1.05 | 1.00 | 1.00 | 1 | 1 |
| 1.05 | 1.00 | 1.00 | 1 | 1 |
Modelski R
Tukaj ostane vse enako kot pri prejšnji simulaciji, le da se osnovno reprodukcijsko število dvigne za 0.5. Predpostavimo, da se reprodukcijsko število dvigne postopoma, s korakom 10 dni za 0.1.
Zaradi različice novega seva je predpostavljeno, da se v 50-ih dnevih poveča iz trenutnega R do R+0.5. Začetek je postavljen na 15.02.2021. Ta hitrost je arbitrarno določena, saj ne poznamo trenutne razširjenosti virusa v populaciji.
S parametri Bw določamo reprodukcijska števila ob različnih časih.
Bw1 = c(1.15,0.5,0.47,0.2,0.1,0.85,0.32,0.33, 0.55, 0.54, 0.55, 0.38, 0.34, 0.35, 0.34, 0.36, 0.42, 0.38, 0.35, 0.33)
casB1 = c(15,7,8,35,35,30,30, 23, 25, 25, 15, 5, 10, 15, 12,10, 12, 10, 10, 10)
Bw2 = c(1.15,0.5,0.47,0.2,0.1,0.85,0.32,0.33, 0.55, 0.54, 0.55, 0.38, 0.34, 0.35, 0.34, 0.36, 0.42, 0.38, 0.35, 0.33, 0.33+0.1/2.9, 0.33+0.2/2.9, 0.33+0.3/2.9, 0.33+0.4/2.9, 0.33+0.5/2.9)
casB2 = c(15,7,8,35,35,30,30, 23, 25, 25, 15, 5, 10, 15, 12,10, 12, 10, 10, 17, 10, 10, 10, 10, 10)
# aprilski ukrepi (ocena)
Bw3 = c(1.15,0.5,0.47,0.2,0.1,0.85,0.32,0.33, 0.55, 0.54, 0.55, 0.38, 0.34, 0.35, 0.34, 0.36, 0.42, 0.38, 0.35, 0.33, 0.33+0.1/2.9, 0.33+0.2/2.9, 0.33+0.3/2.9, 0.33+0.4/2.9, 0.33+0.5/2.9, 1.1/2.9)
casB3 = c(15,7,8,35,35,30,30, 23, 25, 25, 15, 5, 10, 15, 12,10, 12, 10, 10, 17, 10, 10, 10, 10, 15, 10) Iz teh dveh podatkov izračunamo potek R:
Bfun1 = Bt_rect_time(duration_time, Bw1, casB1)
Bfun2 = Bt_rect_time(duration_time, Bw2, casB2)
Bfun3 = Bt_rect_time(duration_time, Bw3, casB3)Izris R
Simulacija
Izračun modelov brez seva in s sevom s cepljenjem
Izris krivulij cepljenja po skupinah
Do 13. 2. 2021 je bilo po podatkih Sledilnika z dvema dozama cepljenih približno 47.000 ljudi. Cepljenje z drugim odmerkom se je začelo 18.1.2021. Tako lahko zaključimo, da je hitrost cepljenja približno 50.000 ljudi na 25 dni. Če to ekstrapoliramo do 31.3.2021, to pomeni približno še 50 dni, ali 150.000 dvojno precepljenih ljudi do konca marca, ali približno 7% precepljene populacije.
Pri tem predpostavimo, da cepimo starostne skupine starejše od 75 let, ki predstavljajo približno 185.000 ljudi. Predpostavimo, da cepimo največ 35% skupine, kar je približno 65.000 ljudi, ostalih 120.000 cepljenih potem predpostavimo v populaciji 65-74, itn.
Možne so seveda še druge strategije.
Izris projekcij modela brez seva in z novim sevom s cepljenjem ter v primeru aprilskih ukrepov
Izris posameznih krivulij modela brez seva in z novim sevom s cepljenjem ter v primeru aprilskih ukrepov
