stringtranslate.com

Математическое моделирование инфекционных заболеваний

Математические модели могут прогнозировать, как инфекционные заболевания прогрессируют, чтобы показать вероятный исход эпидемии ( в том числе среди растений ) и помочь информировать общественное здравоохранение и вмешательства в здоровье растений. Модели используют основные предположения или собранные статистические данные вместе с математикой, чтобы найти параметры для различных инфекционных заболеваний и использовать эти параметры для расчета эффектов различных вмешательств, таких как программы массовой вакцинации . Моделирование может помочь решить, какие вмешательства следует избегать, а какие следует опробовать, или может предсказать будущие модели роста и т. д.

История

Моделирование инфекционных заболеваний — это инструмент, который использовался для изучения механизмов распространения заболеваний, прогнозирования будущего течения вспышки и оценки стратегий борьбы с эпидемией. [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, которая может использоваться как простая модель для моделирования эпидемий. Также используются несколько других типов компартментальных моделей.

Модель SIR

Диаграмма модели SIR с начальными значениями , а также показателями заражения и выздоровления
Анимация модели 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]

Смотрите также

Ссылки

  1. ^ ab Daley DJ, Gani J (2005). Моделирование эпидемий: Введение . Нью-Йорк: Cambridge University Press.
  2. ^ Hethcote HW (2000). «Математика инфекционных заболеваний». SIAM Review . 42 (4): 599–653. Bibcode : 2000SIAMR..42..599H. doi : 10.1137/S0036144500371907. S2CID  10836889.
  3. ^ Блоуэр С., Бернулли Д. (2004). «Попытка нового анализа смертности, вызванной оспой, и преимуществ вакцинации для ее предотвращения». Обзоры по медицинской вирусологии . 14 (5): 275–88. doi :10.1002/rmv.443. PMID  15334536. S2CID  8169180.
  4. ^ «Теория микробов — обзор | Темы ScienceDirect».
  5. ^ Хамер В. (1928). Эпидемиология старая и новая . Лондон: Kegan Paul.
  6. ^ Росс Р. (1910). Профилактика малярии. Даттон.
  7. ^ ab Brauer F, Castillo-Chávez C (2001). Математические модели в популяционной биологии и эпидемиологии . Нью-Йорк: Springer.
  8. ^ Eisinger D, Thulke HH (апрель 2008 г.). «Формирование пространственной структуры способствует искоренению инфекционных заболеваний». Журнал прикладной экологии . 45 (2): 415–423. Bibcode : 2008JApEc..45..415E. doi : 10.1111/j.1365-2664.2007.01439.x. PMC 2326892. PMID  18784795 . 
  9. ^ Адам Д. (апрель 2020 г.). «Специальный отчет: Моделирование, определяющее реакцию мира на COVID-19». Nature . 580 (7803): 316–318. Bibcode :2020Natur.580..316A. doi : 10.1038/d41586-020-01003-6 . PMID  32242115. S2CID  214771531.
  10. ^ Squazzoni F, Polhill JG, Edmonds B, Ahrweiler P, Antosz P, Scholz G и др. (2020). «Вычислительные модели, имеющие значение во время глобальной вспышки пандемии: призыв к действию». Журнал искусственных обществ и социального моделирования . 23 (2): 10. doi : 10.18564/jasss.4298 . hdl : 10037/19057 . ISSN  1460-7425. S2CID  216426533.
  11. ^ Шридхар Д., Маджумдер М.С. (апрель 2020 г.). «Моделирование пандемии». BMJ . 369 : m1567. doi : 10.1136/bmj.m1567 . PMID  32317328. S2CID  216074714.
  12. ^ Maziarz M, Zach M (октябрь 2020 г.). «Моделирование на основе агентов для прогнозирования эпидемии SARS-CoV-2 и оценки вмешательства: методологическая оценка». Журнал оценки в клинической практике . 26 (5): 1352–1360. doi :10.1111/jep.13459. PMC 7461315. PMID  32820573 . 
  13. ^ Huppert A, Katriel G (2013). «Математическое моделирование и прогнозирование в эпидемиологии инфекционных заболеваний». Клиническая микробиология и инфекция . 19 (11): 999–1005. doi : 10.1111/1469-0691.12308 . PMID  24266045.
  14. ^ Тембин Х. «COVID-19: Перспектива игры типа «усредненное поле», основанная на данных. Игры». Games Journal. doi : 10.3390/g11040051 . hdl : 10419/257469 . {{cite journal}}: Цитировать журнал требует |journal=( помощь )
  15. ^ Nakamura GM, Monteiro AC, Cardoso GC, Martinez AS (февраль 2017 г.). "Эффективный метод комплексного вычисления распространения эпидемий на уровне агентов в сетях". Scientific Reports . 7 (1): 40885. arXiv : 1606.07825 . Bibcode :2017NatSR...740885N. doi :10.1038/srep40885. ISSN  2045-2322. PMC 5247741 . PMID  28106086. 
  16. ^ Nakamura GM, Cardoso GC, Martinez AS (февраль 2020 г.). «Улучшенные уравнения восприимчивость–инфекционность–восприимчивость к эпидемиям на основе неопределенностей и автокорреляционных функций». Royal Society Open Science . 7 (2): 191504. Bibcode :2020RSOS....791504N. doi :10.1098/rsos.191504. ISSN  2054-5703. PMC 7062106 . PMID  32257317. 
  17. ^ Dietz K (1967). «Эпидемии и слухи: обзор». Журнал Королевского статистического общества. Серия A (общая) . 130 (4): 505–528. doi :10.2307/2982521. JSTOR  2982521.
  18. ^ Майер Б.Ф., Брокманн Д. (2020). «Эффективное сдерживание объясняет субэкспоненциальный рост недавних подтвержденных случаев COVID-19 в Китае». Science . 368 (6492): 742–746. Bibcode :2020Sci...368..742M. doi :10.1126/science.abb4557. PMC 7164388 . PMID  32269067. 
  19. ^ Сетевая наука Альберта-Ласло Барабаши.
  20. ^ Kenah E, Robins JM (сентябрь 2007 г.). «Второй взгляд на распространение эпидемий в сетях». Physical Review E. 76 ( 3 Pt 2): 036113. arXiv : q-bio/0610057 . Bibcode : 2007PhRvE..76c6113K. doi : 10.1103/PhysRevE.76.036113. ISSN  1539-3755. PMC 2215389. PMID 17930312  . 
  21. ^ Пастор-Саторрас Р., Кастеллано С., Ван Мигем П., Веспиньяни А. (31.08.2015). «Эпидемические процессы в сложных сетях». Reviews of Modern Physics . 87 (3): 925–979. arXiv : 1408.2701 . Bibcode :2015RvMP...87..925P. doi : 10.1103/RevModPhys.87.925 . S2CID  14306926.
  22. ^ K Rizi A, Faqeeh A, Badie-Modiri A, Kivelä M (2022-04-20). «Распространение эпидемий и отслеживание цифровых контактов: эффекты гетерогенного смешивания и неудачи карантина». Physical Review E. 105 ( 4): 044313. arXiv : 2103.12634 . Bibcode : 2022PhRvE.105d4313R. doi : 10.1103/PhysRevE.105.044313 . PMID  35590624. S2CID  232320251.
  23. ^ «Основное репродуктивное число — обзор | Темы ScienceDirect».
  24. ^ Косма Шализи (15 ноября 2018 г.). «Данные в пространстве и времени; Лекция 21: Модели отсеков» (PDF) . Университет Карнеги — Меллона . Получено 19 сентября 2020 г. .
  25. ^ Kermack WO, McKendrick AG (1991). "Вклад в математическую теорию эпидемий--I. 1927". Bulletin of Mathematical Biology . 53 (1–2): 33–55. Bibcode : 1927RSPSA.115..700K. doi : 10.1007/BF02464423 . JSTOR  94815. PMID  2059741.
  26. ^ Brauer F (2017). «Математическая эпидемиология: прошлое, настоящее и будущее». Моделирование инфекционных заболеваний . 2 (2): 113–127. doi :10.1016/j.idm.2017.02.001. PMC 6001967. PMID  29928732 . 
  27. ^ Бриттон Т., Болл Ф., Трапман П. (2020). «Математическая модель раскрывает влияние гетерогенности популяции на коллективный иммунитет к SARS-CoV-2». Science . 369 (6505): 846–849. Bibcode :2020Sci...369..846B. doi :10.1126/science.abc6810. PMC 7331793 . PMID  32576668. 
  28. ^ Поллард А. Дж., Бийкер Э. М. (2021). «Руководство по вакцинологии: от основных принципов к новым разработкам». Nature Reviews Immunology . 21 (2): 83–100. doi :10.1038/s41577-020-00479-7. PMC 7754704. PMID 33353987  . 
  29. ^ Ренц А., Видерспик Л., Дрэгер А. (2020). «FBA раскрывает гуанилаткиназу как потенциальную цель для противовирусной терапии против SARS-CoV-2». Биоинформатика . 36 (Приложение_2): i813–i821. doi :10.1093/bioinformatics/btaa813. PMC 7773487. PMID  33381848 . 
  30. ^ Costris-Vas C, Schwartz EJ, Smith? RJ (ноябрь 2020 г.). «Прогнозирование COVID-19 с использованием прошлых пандемий в качестве руководства: насколько надежными были математические модели тогда и насколько надежными они будут сейчас?». Mathematical Biosciences and Engineering . 17 (6): 7502–7518. doi : 10.3934/mbe.2020383 (неактивен 24.04.2024). PMID  33378907.{{cite journal}}: CS1 maint: DOI inactive as of April 2024 (link)

Дальнейшее чтение

Внешние ссылки

Программное обеспечение