Open rpslimak opened 4 years ago
model Delayed Differencial SIR to vidí (po uvoľnení opatrení) takto ...
@op183 Musím povedať, že @StardustDragon1 má pravdu - je to dosť ťažko čitateľné. (Sorry, k Tvojim emailom som sa zatiaľ nedostal, mal som málo času na tú hobbyepidemiológiu.).
Ak to dobre vidím, tak ku dňu publikácie modelu (4.25) odhaduje vývoj niečoho (infikovaných?) správne, a potom predikuje výbuch infekcií. Samozrejme teraz je každý múdrejši, ale ako vyzerá Tvoj model teraz?
Ja mám podľa SEIR asi takéto niečo: a ten deň som mal toto
(Dátumy ukončenia vlny nepoužívajú ten istý kritérium, preto sú odlišné).
Takto vyzerá dnes (je to aj s legendami). Model predikuje veľmi dobre dnešný stav "pozitívne testovaných". Parametre modelu sú menené vždy pri zmene politiky (zatiaľ tri krát ... 1)rúška a obchody, 2) karanténa, 3) prvé uvoľnenie. Okrem bielej, čo sú reálne dáta, je ostatné model.
ešte to isté, ale bez modelovania "early detected", (odpovedá "delay" modelu IZP)
to je pre porovnanie R0 tak ako ho definuje klasický SEIR model (po uvolnení cca 1.16)
@op183 Ďakujem za nový model aj za vysvetlenie, už je to jasnejšie. Vidím že Tvoj model to nevidí ružovo budúcnosť, aj to že to sedí dodnes pekne.
áno, vždy fixujem model na reálne dáta pokiaľ nepríde k ďalšej "intervencii". Je potrebné si uvedomiť, že modelujeme nie vírus, ale reakciu spoločnosti na jeho správanie. Preto aj akékoľvek predikcie sú len "z brucha", pokiaľ sa nedajú odhadnúť zmeny v parametroch modelu. Posledná zmena == uvoľňovanie, je momentálne len môj "kvalifikovaný odhad". Tie zmeny v parametroch modelu vidíš práve v dolnej časti.
Dolná časť ukazuje "infection rate", to je infekčných dnes / infekčných včera. Farby zodpovedajú vrchným grafom. Sigma9 je disperzia od reálnych dát za posledných 9 dní. (podľa posledných štatistík má výskyt infekcie weibull distribúciu pravdepodobnosti so strednou dobou 8.62 dní, ja od začiatku modelujem latenciu 9 dní, ale exponenciálnu distribúciu. V x-ovej osi sú dni od 6-teho Marca, grid je v týždňoch, sú tam značky prvého v mesiaci. https://github.com/op183/corona-SK-tools . Parametre modelu počítam z reálnych dát predbežne v spreadsheete a potom ručne doladím na minimum sigma9 tými hejblátkami pre lambda a kappa. kappa predstavuje "early detected", teda politiku ktorá sa dotýka len potenciálne infekčných jedincov (napríklad aj povinná karanténa) a lambda, to je politika ktorá sa dotýka každého v populácii (napríklad zavreté obchody alebo rúška).
Ja nie som optimista, stačí si uvedomiť, že to čo vidíme je tá biela čiara (reálne dáta) a tá červená je výpočet nakazených (vrátane nezachytených) jedincov. Mne je jasné, že kým pritvrdením sme šírenie zastavili (vrchol bol okolo 13-teho Apríla, tak opätovné uvoľnenie opatrení nutne znamená opätovný nárast (od 1-vého Mája), ale uvidíme ho reálne na náraste pozitívnych prípadov až od dnes, zajtra, pozajtra ??? ... Ani ten nárast nie je ale žiadna dráma! Naopak, je takmer optimálny. Ak sa ale spustí nákupná turistika do Avionu, alebo pivárne, tak budeme na začiatku a další lockdown sa stane pre kapacitu nemocníc opäť realitou ...
Dnes je v móde sa tešiť z vyliečených a tým podporovať rýchle rušenie obmedzení, ale stačí si uvedomiť, že vyliečení + mŕtvi dnes == nakazení pred mesiacom. Ak vtedy čísla rástli "dramaticky", rovnako "dramaticky" rastú prípady vyliečených dnes. Ďalší mýtus je, že stačí mať pod kontrolou DSS, osady a karanténne centrá, lenže reálne dáta ukazujú, že dnes to už neplatí, vírus už máme v populácii pomerne homogénne.
Zdravim,
tam ako plnite tabulku nakazy_sk, prve 3 su Bratislava v nasledovnom poradi s hodnotami:
529346 - Bratislava - Nové Mesto - hodnota 4 529346 - Bratislava - Nové Mesto - hodnota 3 (duplicita?chyba?) 529320 - Bratislava - Ružinov - hodnota 25
Je to urcite dobre? Zdrojove data ste nedali k dispozicii, je mozne ich zverejnit?
Potom napr. ten Hviezdoslavov - to by mal byt kod obce 501638, ale vy to mate pod Dunajskou Stredou, t.j. 501433...Giraltovce detto.
Košice su zastupene iba 598186 - Košice - Staré Mesto, co je mozno spravne.
Partizanske 505315 je tam tiez duplicitne.
Zverejnite zdrojove tabulky pls.