SEIR model s petimi skupinami

Uvod

Izvedene so kalibracija modela SEIR s petimi skupinami okuženih na dejanske podatke.

Simulacija 1

Starostne skupine

## [1] "<1 "    "1-4"    "5-9 "   "10-14 " "15-19 " "20-24 "
## [1] "25-29 " "30-34 " "35-39 " "40-44 "
## [1] "45-49"  "50-54 " "55-59 " "60-64 "
## [1] "65-69 " "70-74 "
## [1] "75-79 " "80-84 " "85+ "
## [1] "Deleži:"
## [1] 0.24414409 0.27559808 0.28877759 0.10219920 0.08928104

Dvignemo mešanje pri starostnih skupinah 25-65 za 15%. Možno je mešanje vseh skupin med sabo.

Matrika mešanja
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

Bolnišnična obravnava je rezervirana za starostne skupine 45+ v razmerju (0.1, 1.0, 1.5).

Deleži smrti so samo v skupini 4 in 5. To sta starostni skupini od 65+.

Določitev R modela

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.33)

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

#pdat = izracun_modela_V05_01(Bfun, duration_time, w=w, param)
pdat = izracun_modela_dejanski_okuzeni_V5_01(Bfun, sidat, duration_time, w, param)

Izris projekcij modela vseh skupin posamično

Izris projekcij modela vseh skupin skupaj