Математические модели могут прогнозировать, как инфекционные заболевания прогрессируют, чтобы показать вероятный исход эпидемии ( в том числе среди растений ) и помочь информировать общественное здравоохранение и вмешательства в здоровье растений. Модели используют основные предположения или собранные статистические данные вместе с математикой, чтобы найти параметры для различных инфекционных заболеваний и использовать эти параметры для расчета эффектов различных вмешательств, таких как программы массовой вакцинации . Моделирование может помочь решить, какие вмешательства следует избегать, а какие следует опробовать, или может предсказать будущие модели роста и т. д.
Моделирование инфекционных заболеваний — это инструмент, который использовался для изучения механизмов распространения заболеваний, прогнозирования будущего течения вспышки и оценки стратегий борьбы с эпидемией. [1]
Первым ученым, который систематически пытался количественно оценить причины смерти, был Джон Граунт в своей книге « Естественные и политические наблюдения, сделанные над счетами смертности » в 1662 году. Счета, которые он изучал, представляли собой списки цифр и причин смерти, публиковавшиеся еженедельно. Анализ Граунтом причин смерти считается началом «теории конкурирующих рисков», которая, по словам Дейли и Гани [1] , является «теорией, которая в настоящее время прочно укоренилась среди современных эпидемиологов».
Самый ранний отчет о математическом моделировании распространения болезней был выполнен в 1760 году Даниилом Бернулли . Получив медицинское образование, Бернулли создал математическую модель для защиты практики прививок от оспы . [2] Расчеты по этой модели показали, что всеобщая прививка от оспы увеличит продолжительность жизни с 26 лет 7 месяцев до 29 лет 9 месяцев. [3] Работа Даниила Бернулли предшествовала современному пониманию теории микробов . [4]
В начале 20 века Уильям Хамер [5] и Рональд Росс [6] применили закон действия масс для объяснения поведения эпидемий.
В 1920-х годах появились компартментальные модели. Модель эпидемии Кермака-Маккендрика (1927) и модель эпидемии Рида-Фроста (1928) описывают взаимосвязь между восприимчивыми , инфицированными и иммунными лицами в популяции. Модель эпидемии Кермака-Маккендрика успешно предсказывала поведение вспышек, очень похожее на то, что наблюдалось во многих зарегистрированных эпидемиях. [7]
В последнее время агентные модели (ABM) используются вместо более простых компартментальных моделей . [8] Например, эпидемиологические ABM используются для информирования о вмешательствах в здравоохранение (нефармацевтических) против распространения SARS-CoV-2 . [9] Эпидемиологические ABM, несмотря на их сложность и необходимость высокой вычислительной мощности, подвергаются критике за упрощение и нереалистичные предположения. [10] [11] Тем не менее, они могут быть полезны для информирования о решениях относительно мер смягчения и подавления в случаях, когда ABM точно откалиброваны. [12]
Модели хороши ровно настолько, насколько хороши предположения, на которых они основаны. Если модель делает прогнозы, которые не соответствуют наблюдаемым результатам, а математика верна, то первоначальные предположения должны быть изменены, чтобы сделать модель полезной. [13]
«Стохастический» означает быть или иметь случайную величину. Стохастическая модель — это инструмент для оценки вероятностных распределений потенциальных результатов, допускающий случайные вариации в одном или нескольких входных данных с течением времени. Стохастические модели зависят от случайных вариаций в риске воздействия, заболевании и другой динамике болезни. Статистическое распространение болезни на уровне агента в малых или больших популяциях может быть определено стохастическими методами. [14] [15] [16]
При работе с большими группами населения, как в случае туберкулеза, часто используются детерминированные или компартментальные математические модели. В детерминированной модели индивидуумы в популяции распределяются по разным подгруппам или компартментам, каждый из которых представляет определенную стадию эпидемии. [17]
Скорость перехода из одного класса в другой математически выражается как производные, поэтому модель формулируется с использованием дифференциальных уравнений. При построении таких моделей необходимо предполагать, что численность популяции в отсеке дифференцируема по времени и что эпидемический процесс является детерминированным. Другими словами, изменения в популяции отсека можно рассчитать, используя только историю, которая использовалась для разработки модели. [7]
Распространенное объяснение роста эпидемий заключается в том, что 1 человек заражает 2, те 2 заражают 4 и так далее и так далее, причем число инфицированных удваивается с каждым поколением. Это аналогично игре в пятнашки , где 1 человек помечает 2, те 2 помечают 4 других, которые никогда не были помечены, и так далее. По мере развития этой игры она становится все более безумной, поскольку помеченные пробегают мимо ранее помеченных, чтобы выследить тех, кто никогда не был помечен. Таким образом, эта модель эпидемии приводит к кривой, которая растет экспоненциально, пока не рухнет до нуля, поскольку все население инфицировано. т.е. никакого коллективного иммунитета и никакого пика и постепенного спада, как это наблюдается в реальности. [18]
Эпидемии можно моделировать как заболевания, распространяющиеся по сетям контактов между людьми. Такая сеть может быть представлена математически с помощью графа и называется контактной сетью. [19] Каждый узел в контактной сети является представлением индивидуума, а каждая связь (ребро) между парой узлов представляет контакт между ними. Связи в контактных сетях могут использоваться для передачи заболевания между индивидуумами, и каждое заболевание имеет свою собственную динамику поверх своей контактной сети. Сочетание динамики заболевания под влиянием вмешательств, если таковые имеются, на контактной сети может быть смоделировано с помощью другой сети, известной как сеть передачи. В сети передачи все связи отвечают за передачу заболевания. Если такая сеть является локально древовидной сетью, то есть любое локальное соседство в такой сети принимает форму дерева , то базовое воспроизведение может быть записано в терминах средней степени избыточности сети передачи, такой что:
где — средняя степень (средняя степень) сети, а — второй момент распределения степеней сети передачи . Однако не всегда просто найти сеть передачи из контактной сети и динамики заболевания. [20] Например, если контактную сеть можно аппроксимировать графом Эрдеша–Реньи с пуассоновским распределением степеней , а параметры распространения заболевания определены в приведенном выше примере, например, — скорость передачи на человека, а средний инфекционный период заболевания составляет , то базовое число воспроизводства равно [21] [22], поскольку для распределения Пуассона.
Базовое репродуктивное число (обозначаемое как R 0 ) является мерой того, насколько передается болезнь. Это среднее число людей, которых один заразный человек заразит в течение своей инфекции. Это количество определяет, будет ли инфекция увеличиваться субэкспоненциально , вымрет или останется постоянной: если R 0 > 1, то каждый человек в среднем заражает более одного другого человека, поэтому болезнь будет распространяться; если R 0 < 1, то каждый человек в среднем заражает менее одного человека, поэтому болезнь вымрет; и если R 0 = 1, то каждый человек в среднем заразит ровно одного другого человека, поэтому болезнь станет эндемичной: она будет распространяться по всей популяции, но не увеличится и не уменьшится. [23]
Инфекционное заболевание считается эндемичным , когда оно может поддерживаться в популяции без необходимости внешнего воздействия. Это означает, что в среднем каждый инфицированный человек заражает ровно одного другого человека (если больше, то число инфицированных людей будет расти субэкспоненциально , и начнется эпидемия , если меньше, то болезнь исчезнет). В математических терминах это:
Базовое репродуктивное число ( R 0 ) заболевания, предполагая, что все восприимчивы, умноженное на долю населения, которое фактически восприимчиво ( S ), должно быть равно единице (поскольку те, кто не восприимчив, не фигурируют в наших расчетах, поскольку они не могут заразиться болезнью). Обратите внимание, что это соотношение означает, что для того, чтобы заболевание находилось в эндемическом устойчивом состоянии , чем выше базовое репродуктивное число, тем ниже должна быть доля восприимчивой популяции, и наоборот. Это выражение имеет ограничения, касающиеся доли восприимчивости, например, R 0 равно 0,5, подразумевает, что S должно быть равно 2, однако эта доля превышает размер популяции. [ необходима цитата ]
Предположим, что распределение возраста прямоугольное, а возраст инфицирования имеет такое же распределение для каждого года рождения. Пусть средний возраст инфицирования будет A , например, когда индивидуумы моложе A восприимчивы, а те, кто старше A, иммунны (или заразны). Тогда можно показать простым рассуждением, что доля популяции, которая восприимчива, определяется как:
Мы повторяем, что L — это возраст, в котором в этой модели предполагается, что каждый индивидуум умрет. Но математическое определение эндемичного устойчивого состояния можно переформулировать так:
Поэтому, в силу транзитивного свойства :
Это обеспечивает простой способ оценки параметра R0 с использованием легкодоступных данных.
Для популяции с экспоненциальным распределением возраста ,
Это позволяет определить базовое репродуктивное число заболевания с учетом A и L в любом типе распределения популяции.
Компартментальные модели формулируются как цепи Маркова . [24] Классическая компартментальная модель в эпидемиологии — это модель SIR, которая может использоваться как простая модель для моделирования эпидемий. Также используются несколько других типов компартментальных моделей.
В 1927 году WO Kermack и AG McKendrick создали модель, в которой они рассматривали фиксированную популяцию только с тремя отсеками: восприимчивый, ; инфицированный, ; и выздоровевший, . Отсеки, используемые для этой модели, состоят из трех классов: [25]
Существует множество модификаций модели SIR, в том числе те, которые включают в себя рождение и смерть, где после выздоровления иммунитет отсутствует (модель SIS), где иммунитет сохраняется только в течение короткого периода времени (SIRS), где существует латентный период болезни, когда человек не заразен ( SEIS и SEIR ), и где младенцы могут рождаться с иммунитетом (MSIR). [ необходима ссылка ]
Математические модели должны интегрировать растущий объем данных , генерируемых по взаимодействиям хозяина и патогена . Многие теоретические исследования динамики популяции , структуры и эволюции инфекционных заболеваний растений и животных, включая людей, посвящены этой проблеме. [26]
Темы исследований включают:
Если доля населения, обладающего иммунитетом, превышает уровень коллективного иммунитета к заболеванию, то заболевание больше не может сохраняться в популяции, и его передача прекращается. [27] Таким образом, заболевание может быть устранено из популяции, если достаточное количество людей обладает иммунитетом либо благодаря вакцинации, либо благодаря выздоровлению от предшествующего воздействия заболевания. Например, искоренение оспы , последний случай дикой природы был зарегистрирован в 1977 году, и сертификация искоренения местной передачи 2 из 3 типов дикого полиовируса (тип 2 в 2015 году, после последнего зарегистрированного случая в 1999 году, и тип 3 в 2019 году, после последнего зарегистрированного случая в 2012 году). [28]
Уровень коллективного иммунитета будет обозначаться q . Напомним, что для стабильного состояния: [ необходима цитата ]
По очереди,
что приблизительно составляет: [ необходима цитата ]
S будет (1 − q ), поскольку q — это доля населения, которая имеет иммунитет, а q + S должно быть равно единице (поскольку в этой упрощенной модели каждый либо восприимчив, либо имеет иммунитет). Тогда:
Помните, что это пороговый уровень. Вымирание от передачи произойдет только в том случае, если доля иммунных лиц превысит этот уровень из-за программы массовой вакцинации.
Мы только что рассчитали критический порог иммунизации (обозначаемый q c ). Это минимальная доля населения, которая должна быть иммунизирована при рождении (или близко к рождению), чтобы инфекция в популяции исчезла.
Потому что доля окончательной численности популяции p , которая никогда не заражается, может быть определена как:
Следовательно,
Решая для , получаем:
Если используемая вакцина недостаточно эффективна или требуемый охват не может быть достигнут, программа может не превысить q c . Такая программа защитит вакцинированных людей от болезни, но может изменить динамику передачи. [ необходима цитата ]
Предположим, что часть населения q (где q < q c ) привита при рождении от инфекции с R 0 > 1. Программа вакцинации изменяет R 0 на R q , где
Это изменение происходит просто потому, что в популяции теперь меньше восприимчивых людей, которые могут быть инфицированы. R q — это просто R 0 минус те, кто в обычных условиях был бы инфицирован, но не может быть инфицирован сейчас, поскольку они обладают иммунитетом.
Вследствие этого более низкого базового числа воспроизводства средний возраст инфицирования A также изменится на некоторое новое значение A q у тех, кто остался невакцинированным.
Вспомним соотношение, связывающее R 0 , A и L . Предполагая, что продолжительность жизни не изменилась, теперь: [ необходима цитата ]
Но R 0 = L / A, поэтому:
Таким образом, программа вакцинации может повысить средний возраст инфицирования, а невакцинированные лица будут испытывать меньшую силу инфицирования из-за присутствия вакцинированной группы. Для заболевания, которое приводит к большей клинической тяжести у пожилых людей, невакцинированная часть населения может испытать болезнь относительно позже в жизни, чем это произошло бы при отсутствии вакцины.
Если программа вакцинации приводит к тому, что доля иммунных лиц в популяции превышает критический порог в течение значительного периода времени, передача инфекционного заболевания в этой популяции прекращается. Если ликвидация происходит везде в одно и то же время, то это может привести к искоренению . [ необходима цитата ]
Модели имеют преимущество в том, что они одновременно изучают несколько результатов, а не делают один прогноз. Модели показали высокую степень надежности в прошлых пандемиях, таких как SARS , SARS-CoV-2 , [29] свиной грипп , MERS и Эбола . [30]
{{cite journal}}
: Цитировать журнал требует |journal=
( помощь ){{cite journal}}
: CS1 maint: DOI inactive as of November 2024 (link){{cite book}}
: CS1 maint: location (link)