Simulacija vpliva odpiranja šol

Uvod

Uporabljen je kalibriran model SEIR s petimi starostnimi skupinami za simulacijo vpliva odpiranja šol s 25-01-2021.

Simulacija

Simulacija je izvedena na SEIR modelu s 5 starostnimi skupinami. Model je kalibriran na današnje podatke. Cepljenje še ni upoštevano.

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+.

Matrike mešanja

Predpostavljeno je, da bo odprtje šol povečalo stopnjo prenosa v starostni skupini otrok in mladostnikov za največ 25%. Vendar se bo hkratno povečala stopnja prenosa med skupino 1 in ostalimi skupinami v razmerju: 15%, 10%, 5% in 5% za skupine 2, 3, 4 in 5.

Povečanje se bo zgodilo od 25-01-2021 naprej.

Matriki stopenj prenosa sta prikazani v nadaljevanju.

Matrika mešanja1
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
Matrika mešanja2
1.25 1.00 1.00 1 1
1.15 1.15 1.15 1 1
1.10 1.15 1.15 1 1
1.05 1.00 1.00 1 1
1.05 1.00 1.00 1 1

Modelski R

S parametri Bw določamo reprodukcijska števila ob različnih časih.

#Bw = c(1.15,0.5,0.47,0.2,0.1,0.85,0.32,0.33,0.54,0.5,0.65,0.34,0.34,0.35, 0.30, 0.40, 0.33)
Bw = 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.33,  0.38, 0.42, 0.35)

casB = c(15,7,8,35,35,30,30, 23, 25, 25, 15, 5, 10, 15, 10, 15, 14, 10)

Iz teh dveh podatkov izračunamo potek R:

Bfun = Bt_rect_time(duration_time, Bw, casB)

Izris R

Izračun modela skupaj z dejanskimi podatki

# brez spremembe
pdat1 = izracun_modela_dejanski_okuzeni_V5_02(Bfun, sidat, duration_time, w, param1)

# s spremembo
pdat2 = izracun_modela_dejanski_okuzeni_V5_02(Bfun, sidat, duration_time, w, param2)

pdat = pdat2

Razlike med situacijo brez šol in s šolami

Situacija brez šol je označena z 1, situacija s šolami pa kot situacija 2.

Hospitalizirani

Izris relativne razlike med krivuljami

Obravnave na intenzivni negi

Izris relativne razlike med krivuljami

Smrti

Izris relativne razlike med krivuljami