Simulace šíření spalniček (Vensim)

From Simulace.info
Revision as of 17:56, 9 June 2019 by Povj01 (talk | contribs) (Definice problému)
Jump to: navigation, search

Zadání

Název simulace: Simulace šíření viru spalniček

Předmět: 4IT495 Simulace systémů

Semestr: LS 2018/2019

Autor: Bc. Jurij Povoroznyk

Typ modelu: Systémová dynamika

Modelovací nástroj: Vensim

Definice problému

V České republice propukla epidemie spalniček. Tato nemoc se k nám dostal od cestovatele z Indie. Celkově bylo nakaženo 50 lidí a další lidé rychle přibývají. Onemocnění se projevuje horečkou, rýmou, kašlem, slzícíma očima a na bukální sliznici jsou bělavé tečky se zarudlým okolím. Virus spalniček se přenáší kapénkovou infekcí. Inkubační doba spalniček je 6–19 dní, průměrně 13 dní. Infikovaní lidé jsou nakažliví ještě 4 až 5 dní před propuknutím této nemoci. Úmrtnost je velmi malá, 3 smrti z 1 000 případů. Cílem tohoto modelu je určit a sledovat průběh této epidemie. Zároveň pomocí tohoto modelu lze určit důležitost karantény a počet vakcín obdržených před vypuknutím epidemie. Údaje získané z této simulace by měly přesvědčit rodiče k očkování svých děti ihned jak to bude možné.

Metoda

Pro zpracování této simulace byl vybrán nástroj Vensim pro simulaci systémové dynamiky. Zadáním specifických hodnot proměnných lze nalézt hraniční hodnoty, optimální hodnoty a sledovat chování ekonomického modelu. U ostatních softwarových nástrojů by došlo k zanedbání prostorového aspektu vnějšího prostředí. Zároveň na takto malé ploše by došlo k zkreslení výsledků zanesením pouze omezeného počtu agentů. Nás především zajimají výsledné faktické hodnoty, které jsou v tomto nástroji daleko více zřejmé a prezentovatelné. Model sleduje vývoj v horizontu několika měsíců.

Model

Diagram stavů a toků

Stock and flow diagram spalnicky.PNG

Diagram kauzálních smyček

Causal diagram spalnicky.PNG

Entity modelu


Míra porodnosti - Pro účel této simulace byla přidaná míra porodnosti, která vyobrazuje celkovou míru porodnosti na jednu osobu. Počáteční hodnota byla nastavena tak, aby byly získány relevantní data. Tento údaj se může změnit pomocí boxu umístěných nad modelem. Minimální hodnota této konstanty je 0, tudíž nedochází k žádnému početí během jednoho měsíce. Maximální možnou hodnotou může být 1, kdy na jednu osobu připadá jedno nově narozené dítě.

Hodnota: míra porodnosti na jednu osobu [0.35]


Narození - Vyjadřuje počet narozených dětí z celkové populace.

Vzorec: celková populace * míra porodnosti


Celková populace - Součet všech lidí vyskytujících se v této konkrétní simulaci. Počet je redukován o mrtvé jedince, kterých je ovšem minimum, jelikož spalničky mají takřka nulovou úmrtnost.

Vzorec: susceptible + infected + recovered + quarantined


Susceptible (náchylný jedinci k onemocnění) - Vyobrazuje počet jedinců, který se můžou nakazit virem spalniček. Tyto údaje jsou v jednotkách - individual, což znázorňuje počet lidí.

Vzorec: INTEG (narození - očkovaný - onemocněl)
Počáteční hodnota: počáteční počet zdravých


Počáteční počet zdravých - Tato konstanta inicializuje do modelu počet jedinců, kteří jsou zdravý a můžou se nakazit tímto virem. Pomocí posuvníku se hodnota může zvýšit či snížit. Slouží pro nastavení simulace na konkrétní situaci z reality.

Hodnota: počáteční počet zdravých [2000]


Denní počet - Nastavuje denní počet setkání jedinců s nakaženými lidmi.

Hodnota: denní počet [4]


Počet setkání - Tento výpočet určí možnost setkání zdravých jedinců s nakaženými. Výpočet byl omezen pouze na lidi, kteří nejsou v karanténě. Vyléčený lidí nemůžou být nakažení, ovšem můžou sloužit jako přenášeči, proto jsou v tomto vzorci využiti.

Hodnota: (denní počet * Susceptible) * (Infected / (celková populace - Quarantined))


Míra nakažlivosti - Procenta znázorňují možnost nákazy touto nemocí. Hodnota je nastavená především tak, že někteří lidí mají silnější imunitu, proto šance na nakažení není 100%, ale pouze 60%.

Hodnota: míra nakažlivosti [0.6]


Onemocněl - Znázorňuje počet lidí za jeden měsíc, kteří onemocněli. Výpočet je omezen pouze na infikovaný subjekty. Jedinci nacházející se v karanténě nemůžou nakazit zdravé lidi.

Vzorec: (míra nakažlivosti * počet setkání) + (Infected - Quarantined)


Infected (infikovaný jedinci) - Vyobrazuje počet jedinců, který jsou nakažení virem spalniček. Zde je počáteční hodnota počet nakažených subjektů, kteří nejsou v karanténě a neví, že mají tento virus v těle.

Vzorec: INTEG (onemocněl - karanténa - mrtvý - vyléčen)
Počáteční hodnota: počáteční počet nakažených


Počáteční počet nakažených - Počet nakažených jedinců virem spalniček při spuštění simulace.

Počáteční hodnota: počáteční počet nakažených [50]


Přístupnost vakcíny - Udává procentuální dostupnost vakcíny proti spalničkám. Ne všichni mají možnost tuto vakcínu získat. Ať už kvůli sociální situaci nebo v dané oblasti není dostatečný počet vakcín.

Počáteční hodnota: přístupnost vakcíny [0.5]


Množství vakcíny - Výpočet sloužící pro určení indexu lidí z celkového počtu zdravých subjektů, kteří mají zájem o očkování proti spalničkám.

Vzorec: (onemocněl / Susceptible) * přístupnost vakcíny


Nebezpečí výskytu - Konstanta udává nebezpečí výskytu spalniček v dané simulaci. Minimální hodnota je 0, kde nebezpečí je nulové. Maximální je 1, což označuje 100% možnost výskytu spalniček. Tato hodnota slouží primárně jako motivace lidí nechat se očkovat proti tomuto viru.

Vzorec: (onemocněl / Susceptible) * přístupnost vakcíny


Síla očkovaných 1, 2 nebo žádnou vakcínou - Tyto konstanty slouží pro udání síly vakcíny, která se aplikuje na subjekt. Síla vakcíny se liší podle počtu aplikovaných dávek a efektivnosti daného vzorku. Údaje byly získány ze statistického úřadu. [1] Pokud subjekt nemá žádnou vakcínu, je šance na onemocnění 90 %. Proto byla nastavená hodnota 0.1. [2]

Hodnoty: 1 vakcína [0.87], 2 vakcíny [0.95], žádná vakcína [0.1]


Účinnost - Zde se propočítá, jaká je síla vakcíny a kterou možnost si lidé nejčastěji volí. Konkrétní možnost si zvolí na základě strachu z výskytu této nemoci v jejich oblasti.

Vzorec: IF THEN ELSE(nebezpečí výskytu > 0.6, síla očkovaných 1 vakcínou , 0) + IF THEN ELSE(nebezpečí výskytu > 0.8, síla očkovaných 2 vakcínama, 0 ) + IF THEN ELSE(nebezpečí výskytu <= 0.6, síla neočkováno, 0)


Míra úmrnosti - Údává procentuální možnost úmrtí na tuto nemoc. Úmrtnost u tohoto konkrétního viru je velmi nízká, proto byla nastavená odpovídající hodnota. [3]

Hodnota: míra úmrtnosti [0.001]


Mrtvý - Počet mrtvých osob z infikovaných jedinců.

Vzorec: Infected * míra úmrtnosti


Míra zotavení - Udává jaká je šance na zotavení. Spalničky nejsou příliš smrtelné, tzn. že všichni co projdou 13 denní nemocí, tak se uzdraví. V této simulace byly vyloučený zdravotní potíže, které můžou být i po absolvování tohoto viru.

Hodnota: míra zotavení [1]


Vyléčen - Výpočet vyléčených osob.

Vzorec: míra zotavování * Infected


Recovered (vyléčení a imunní jedinci) - Počet subjektů, kteří jsou vyléčení ze spalniček, nebo jsou imunní kvůli absolvování nebo vakcíně.

Vzorec: INTEG (očkovaný + vyléčen + zotavování v karanténě)
Počáteční hodnota: počet vyléčených [0]


Míra karantény - Procentuální počet infikovaných lidí jde do karantény. Číslo zobrazující povinnost karantény, ať už kvůli vyhlášce nebo individuální preferencí. Toto číslo jde regulovat posuvníkem umístěným nad modelem.

Hodnota: míra karantény [0.6]


Karanténa - Výpočet infikovaných lidí, kteří jdou do karantény.

Vzorec: míra karantény * Infected


Quarantined (jedinci nacházející se v karanténě) - Počet subjektů, kteří se nachází v karanténě. Tito jedinci nemůžou nakazit zdravý jedince.

Vzorec: INTEG (Quarantined + karanténa - zotavování v karanténě)
Počáteční hodnota: počet jedinců v karanténě [0]


Zotavování v karanténě - Průběh léčby je identický jako u lidí, kteří nejsou v karanténě.

Vzorec: Quarantined * míra zotavování


Posuvníky určené k regulaci konstant

Posuvníky níže slouží pro úpravu primárních konstant, které jsou důležité pro získání výsledných dat z této simulace. Primárně jde o míru karantény a nebezpečí výskytu, jenž mění preferenci lidí, zda se jít naočkovat proti spalničkám či nikoliv. Boxy usnadňují regulaci a umožňují měnit hodnoty, které vstupují do jednotlivých vzorců.

Posuvniky spalnicky.PNG

Ovlivňující faktory

Z následujících dvou obrázku je viditelné, jaké faktory ovlivňují celkový počet infikovaných a vyléčených (imunních) jedinců.

Ovlivnujici faktory recovered.PNG Ovlivnujici faktory infected.PNG

Výsledek

Výsledek simulace bude porovnán mezi třemi různými situacemi. Níže jsou vyobrazeny hodnoty operativních proměnných, který se mění podle zkoumaných situací, které jsou popsány v zadání.

Normální situace

  • Počáteční počet zdravých - 2 000
  • Počáteční počet nakažených - 50
  • Míra karantény - 0.6
  • Nebezpečí výskytu - 0.2
  • Míra porodnosti - 0.35


Normal situation spalnicky table.PNG
Normal situation spalnicky graph.PNG

Větší nebezpečí výskytu spalniček

  • Počáteční počet zdravých - 2 000
  • Počáteční počet nakažených - 50
  • Míra karantény - 0.6
  • Nebezpečí výskytu - 0.9
  • Míra porodnosti - 0.35


Danger spalnicky table.PNG
Danger spalnicky graph.PNG

Větší důraz na karanténu nemocných lidí

  • Počáteční počet zdravých - 2 000
  • Počáteční počet nakažených - 50
  • Míra karantény - 0.9
  • Nebezpečí výskytu - 0.2
  • Míra porodnosti - 0.35


Karantena spalnicky table.PNG
Karantena spalnicky graph.PNG

Závěr

Kód

Diagram stavů a toků
Diagram kauzálních smyček

Zdroje

  1. LIMBERKOVÁ, Radomíra. (2019) Kdy jste imunní vůči onemocnění spalničkami? Nerozumíte laboratornímu nálezu? [online]. [cit. 2019-09-06] Dostupné z: http://www.szu.cz/tema/prevence/kdy-jste-imunni-vuci-onemocneni-spalnickami-nerozumite
  2. SHERMAN, Alter. (2015) 90 percent chance of getting measles, if you are exposed and unvaccinated [online]. [cit. 2019-09-06] Dostupné z: https://www.childrensdayton.org/the-hub/90-percent-chance-getting-measles-if-you-are-exposed-and-unvaccinated
  3. ROZALIO, s. (2018) Spalničky [online]. [cit. 2019-09-06] Dostupné z: http://www.rozalio.cz/o-nemocech/spalnicky