Комментарии 22
Я спрашиваю потому, что все модели, что описывались в популярных статьях это какая-то фикция. Чем могут помочь графики без привязки к географии? На мой делитантский взглят тут больше подходят модели типа погодных, движения воздуха/жидкости или симуляции ядерных взрывов. Другой перспективный (опять-же на мой взгляд) подход будет симуляция систем частиц. Все такие модели схожи требованием огромных вычислительных мощьностей.
Проблема не в самой модели. Модель содержит максимально абстрактные вероятностные переменные. Проблема в том, для каких групп величина этих переменных считается одной и той же, и из статистики какого качества эта величина берётся. Естественно, при ленивом моделировании, оказывается, что выборка — заведомо велика, и получается, почти в прямом смысле, средняя температура по больнице, включая морг.
На практике же инфекция может достаточно свободно распространяться внутри некоторых кластеров и крайне неохотно передается между кластерами, что существенно замедляет ее распространение.
У нас на работе два отдела (абонентский и производственно-технический) переболели осенью практически в полном составе и почти одновременно. А других отделов это не коснулось. Во всяком случае — в тот раз.
Здесь на Хабре один автор показал, что в условиях ограниченного числа знакомств скорость распространение инфекции быстро замедляется и из экспонентной сначала становится линейной, а потом — константной.
Собственно говоря, по статистике большинства стран это прекрасно видно — сначала экспотенциальный рост сменяется линейным, а потом после прохождения пика несколько успокаивается и долго тянется «толстый хвост». И это не затухание эпидемии, а именно ее тлеющее распространение.
Для этого мы добавляем соответствующие компартменты на каждый отдел офиса. Допустим два отдела. Тогда изначально будет две группы S (S1 и S2) и две группы I (I1 и I2), а R оставим одну (она не на кого не влияет). Из S1 индивиды попадают в I1, а потом в общую R. Из S2 в I2, а потом в общую R. На интенсивность перехода из S1 в I1 будет сильно влиять сама I1 и совсем немного I2 (так как контакт между отделами всё таки возможен). И соответственно наоборот: на интенсивность перехода из S2 в I2 сильно влияет I2 и слегка I1.
Так получается описываемая Вами кластеризация. Вот, что у меня быстро получилось по такой схеме. Эпидемия возникает только в одном отделе. В другом отделе она начинается после (просачивается таки).
![image](https://habrastorage.org/getpro/habr/comment_images/c1a/2f5/801/c1a2f58011de815880de1598ab446dd0.png)
Про тлеющее распространение могу сказать только, что замечал подобное явление в моделях, в которых каким-либо образом восстанавливается пул восприимчивых. Либо когда иммунитет не образуется вовсе (SIS), либо утрачивается со временем (SIRS), либо модели с рождаемостью (новорожденные по умолчанию восприимчивы).
Если отдела два, то между ними возникает канал передачи инфекции. если три, то таких каналов уже три. Если четыре — то уже шесть.
А ведь кроме отделов одного предприятия существуют другие предприятия, с которыми данное предприятие может контактировать, а может и нет. И семьи работников, которые работают на этих других предприятиях и ходят в общественные места.
Т.е. простое усложнение структуры — тупиковая ветвь. Потому что по мере усложнения структуры модель усложняется нелинейно.
Значит, модель должна описывать сложную иерархическую структуру, но при этом быть сама простой. Вроде как фрактальные объекты — описываются просто, но имеют бесконечно сложную структуру.
Так и тут — система представляет собой фрактал — люди взаимодействуют между собой. Но каждый человек за время инфицирования может инфицировать ограниченное число контактов. И то только с некоторой вероятностью.
Но он входит в кластер, который взаимодействует с другими кластерами. И эти кластеры тоже заражают друг-друга. Причем вероятность заражения нелинейна во времени — чем больше иефицированных в кластере, темм больше вероятность передачи заразы в другие кластеры. Но тоже не во все — ведь число контактов каждого кластера тоже ограничено.
И т.д. от уровня к уровню.
Но мне кажется, что сначала стоит попробовать создать модель для ограниченного числа контактов между людьми. И посмотреть, как это скажется. Только надо учесть, что пространство общения псевдоэвклидово. Т.е. число общих контактов у контактирующих друг с другом выше, чем у неконтактирующих.
Ваш вопрос натолкнул меня на хорошую идею, благодарю. Конкретно описанные в статье модели не могут учитывать географию. И за это как раз часто и ругают SIR модель. Однако я хотел сказать, что компартментальные модели — это неисчерпаемый источник новых моделей, которые могут учитывать различные, необходимые в конкретной ситуации факторы.
Например, географию. Можно условно разделить моделируемый географический участок на такие же группы. Допустим на страны/регионы/города или ещё мельче. Для каждого региона будет свой компартмент S и соответственно I и R, которые имеют свои вероятности переходов и коэффициенты воздействия на другие географические участки.
С поведением также. Объединяем индивидов популяции в группы по схожему поведению. Например, носящие маски и игнорирующие маски. Каждая группа, конечно, имеет свой коэффициент заразности и восприимчивости.
Графики же самый простой способ визуализации. Действительно, если будет географическое разделение, можно визуализировать распространение цветами на карте.
![](https://habrastorage.org/webt/qf/nh/fn/qfnhfnuli12lu1zun0ryvhmrqyc.gif)
https://commons.wikimedia.org/wiki/File:COVID-19_Ukraine_epidemic_dynamics_model.gif
А откуда берутся вторая волна, третья?.. Хотя бы на качественном уровне. Типа начался спад, све расслабились… Или как?
habr.com/ru/post/551682/#comment_22908330
Учитывается ли группа не переболевших, но по не выясненной ещё причине не восприимчивые?
Но я думаю, это не та причина, которая обуславливает волны COVID.
![image](https://habrastorage.org/getpro/habr/comment_images/23f/b88/1c5/23fb881c517a3d936f53398d951241de.png)
Если модель учитывает рождаемость и смертность, то волны происходят из-за появления новорожденных, у которых нет иммунитета. Но это тоже пока не про COVID — слишком мало времени ещё прошло, чтобы говорить о таких волнах.
Классической причиной возникновения волн гриппа являются сезонные колебания температуры, которые изменяют возможности распространения вируса. Также возможно влияет введение карантинных мер на определённый период. Их эффективность это отдельный разговор и повод для споров. Сейчас, конечно, не об этом. В любом случае меняется эффективность заражения.
Либо используют модель, учитывающую общую рождаемость и смертность (то есть некоторые добавления/вычитания индивидов из нужных компартментов), и делают смертность на стадии болезни больше чем на других.
И где хоть одна модель, которая объясняет "волны"?
Модель не может объяснить никакие волны. Объяснить волны должна гипотеза, по которой подбирается эмпирический закон изменения параметров модели.
Модель не в состоянии заменить понимание механики процессов, которые она моделирует. Она служит только задаче получения одних возможных данных из других, имеющихся.
Скажем, эмпирическая модель колебания пружинного маятника в виде затухающей синусоидальной функции координаты от времени, построенная на основании только наблюдения за маятником, никак не объясняет, под воздействием каких сил эти колебания происходят и почему затухают. Она может только показать, сколько времени понадобится для прекращения колебаний конкретного маятника, если оттянуть его в некое положение. А вот предсказать поведение любого маятника мы сможем только если будем знать, почему он колеблется (соответственно, как действующие на груз силы зависят от свойств пружины и положения маятника в пространстве относительно планеты).
habr.com/ru/post/551682/#comment_22908330
Такие модели, могут продемонстрировать волны заболеваемости.
Чисто интуитивно это выглядит так, как будто после высокого пика эпидемия должна взять передышку. Видимо должны постепенно изменится круги общения людей, чтобы появились новые связи
О факторах реально влияющих на волны ковида вопрос отдельный. Я указал возможность влияния карантинных мер. Всё таки я не исключаю сезонность температуры и действительно интересна Ваша идея про круг общения. Стоит подумать как реализовать это в рамках компартментов.
Думаю, здесь работает совокупность факторов, поэтому получается примерно одно и то же время спада эпидемии по странам.
Конструирование эпидемиологических моделей