SEIR-Modell

aus Wikipedia, der freien Enzyklopädie
Zur Navigation springen Zur Suche springen
Epidemie-Verlauf im SEIR-Modell für eine Latenzzeit von 1 bis 4 Tagen. Im Vergleich dazu gestrichelt der Verlauf im SIR-Modell. Die Abweichung vom SIR-Modell vergrößert sich mit steigender Latenzzeit. Die Kurven zeigen die Neuinfektionen pro Tag in Abhängigkeit der Zeit in Tagen, relativ zur Populationsgröße. Parameter: Basisreproduktionszahl 2,0, infektiöse Zeit von 4 Tagen, 100 Infektiöse pro Mio. zum Startzeitpunkt.

Als SEIR-Modell bezeichnet man in der mathematischen Epidemiologie einen Ansatz zur Beschreibung der Ausbreitung von ansteckenden Krankheiten. Die Beschreibung ist näher am realen Verlauf als die des SIR-Modells, da hier berücksichtigt wird, dass ein Individuum nach seiner Infektion nicht sofort selbst infektiös ist. Im Gegensatz zu einem Individuum-basierten Modell ist die Beschreibung makroskopisch, d. h. die Population wird als Gesamtheit betrachtet.

Klassisches SEIR-Modell

[Bearbeiten | Quelltext bearbeiten]

Modelldefinition

[Bearbeiten | Quelltext bearbeiten]

Die Population von N Individuen sei zerlegt in die vier Kompartimente S, E, I, R, so dass

Jedes Individuum kann die Prozedur

Susceptible (S) → Exposed (E) → Infectious (I) → Recovered (R)

durchlaufen. Die Ausbreitungsdynamik der Krankheit wird beschrieben durch das System

wobei

die Ableitungen nach der Zeit sind. Betrachtet man die bezüglich der konstanten Population relativen Anteile

nimmt das System die Gestalt

an. Hierbei handelt es sich um ein autonomes nichtlineares System von gewöhnlichen Differentialgleichungen, weshalb das Modell ein dynamisches System beschreibt. Jeder Zustand ohne Exponierte und Infektiöse stellt eine Ruhelage des dynamischen Systems dar, die krankheitsfreies Gleichgewicht genannt wird, zuweilen abgekürzt als DFE für disease free equilibrium. Sollte die effektive Reproduktionszahl den kritischen Wert 1 übersteigen, handelt es sich um eine instabile Ruhelage, aus deren Umfeld sich eine Epidemie entwickelt, woraufhin mit dem Durchlaufen der Epidemie eines der stabilen krankheitsfreien Gleichgewichte angestrebt wird.

Größe Einheit Erklärung
S(t) indv Anteil der Anfälligen, engl. susceptible. Noch nicht infiziert.
E(t) indv Anteil der Exponierten, engl. exposed. Infiziert, aber noch nicht infektiös.
I(t) indv Anteil der Infektiösen, engl. infectious.
R(t) indv Anteil der Erholten, engl. recovered oder resistant. Bzw. verstorben oder nach Symptomen in Quarantäne.
t d Zeit in Tagen, engl. time.
β 1/d Transmissionsrate. Der Kehrwert ist die mittlere Zeit zwischen Kontakten mit Übertragung.
γ 1/d Erholungsrate. Der Kehrwert ist die mittlere infektiöse Zeit.
α 1/d Übergangsrate. Der Kehrwert ist die mittlere Latenzzeit.

Die mittlere Latenzzeit ist die durchschnittliche Zeit, die ein Individuum in der Gruppe E der Exponierten verbringt; diese ist zu unterscheiden von der mittleren Inkubationszeit, denn der Beginn der Infektiosität muss nicht mit dem Beginn der Symptome übereinstimmen.

Für die Transmissionsrate (Übertragungsrate) ist auch die Bezeichnung Kontaktrate geläufig. Eine genauere Überlegung zerlegt diese in ein Produkt , wobei die Transmissionswahrscheinlichkeit und die eigentliche Kontaktrate ist.

Das Modell macht eine Reihe von Annahmen, welche die Wirklichkeit stark vereinfachen. Im Einzelnen:

  • Es findet eine gute Durchmischung der Bevölkerung statt, so dass jedes mit jedem anderen Individuum in Kontakt kommen kann. Andernfalls müsste man Kontaktnetzwerke oder räumliche Dynamiken beschreiben.
  • Es gibt hinreichend viele Infektiöse, so dass der Verlauf der Epidemie näherungsweise deterministisch stattfindet. Andernfalls müsste man einen stochastischen Prozess beschreiben.
  • Die Werte sind kontinuierlich. Das Modell beschreibt keine diskreten Anzahlen.
  • Sowohl die Dauer der latenten Phase als auch der infektiösen Phase ist exponentialverteilt (siehe den Abschnitt Infektionsalter-Modell unten). Die Transmissionsrate in der infektiösen Phase ist hierbei konstant.
  • Die Krankheit führt entweder zur vollständigen und lebenslangen Immunität oder zum Versterben.
  • Jedes anfällige Individuum ist gleich anfällig für die Krankheit.
  • Ein Nachwuchs neuer Anfälliger und natürliches Versterben tritt nicht auf. Das Modell enthält keine Beschreibung der demografischen Dynamik. Insbesondere ist die Größe der Bevölkerung konstant.

Beziehung zur Basisreproduktionszahl

[Bearbeiten | Quelltext bearbeiten]

Zwischen den Parametern des Modells und der Basisreproduktionszahl lässt sich eine Beziehung herstellen. Damit sich die Krankheit nicht weiter ausbreitet, muss und sein. Einsetzung dieser Bedingungen in die Differentialgleichungen führt zu und , und somit . Diese Gleichgewichtsbedingungen entsprechen der effektiven Reproduktionszahl . Weil aber per Definition gilt, ergibt sich[1]

Gelegentlich wird alternativ so definiert, dass gilt, mit der Bewandtnis, bei Benutzung der Großbuchstaben den Parameter nicht in den Gleichungen mitführen zu müssen.[2][3] Unterscheidungen dieser Art sind durch Normalisierung auf entfernbar, denn das Modell kodiert per se keine diskreten Zahlen.

Man definiert , das ist der immune Anteil der Population. Im Gleichgewicht

ist demnach

Dieser Anteil, oberhalb dessen Herdenimmunität besteht, heißt kritische Immunisierungsschwelle. Diese Schwelle ist nur von der Basisreproduktionszahl abhängig, nicht aber vom gewählten Modell.

Der Parameter , der den Übergang in das Zwischenstadium beschreibt, in der die Person infiziert, aber noch nicht ansteckend ist, spielt bei der Basisreproduktionszahl keine Rolle, solange man die Möglichkeit des Versterbens in der Gruppe der Exponierten vernachlässigt.[2] Falls man dies nicht tut, erhält man eine Formel für wie unten im Abschnitt Einbeziehung demografischer Dynamik angegeben.

Die Beziehung der Parameter zur Basisreproduktionszahl ist also im Allgemeinen für jede Modifikation des Modells neu zu ermitteln. Eine systematische Methode der Bestimmung wurde im Jahr 1989 von Odo Diekmann, Hans Heesterbeek und Hans Metz gefunden,[4] und später von Pauline van den Driessche und James Watmough weiter ausgearbeitet.[5] Bei dieser Herangehensweise wird die Jacobi-Matrix des Systems am krankheitsfreien Gleichgewicht dergestalt zu einer Differenz zerlegt, dass der Minuend die Neuinfektionen charakterisiert und der Subtrahend die sonstigen Übergänge zwischen den Kompartimenten. Das Produkt aus dem Minuend und der Inverse des Subtrahenden bildet eine neue Matrix, von ihnen Next generation matrix genannt. Ihr Spektralradius ist die Basisreproduktionszahl.

Maximal möglicher Anteil an Infizierten

[Bearbeiten | Quelltext bearbeiten]

Die Basisreproduktionszahl bestimmt wesentlich den maximal möglichen Anteil von Infizierten an einer gegebenen Population. Bezeichnen wir den Gesamtanteil aller Infizierten mit , also , so folgt aus dem obigen Differentialgleichungssystem

,

oder mit der Basisreproduktionszahl einfach

Mit Umstellung dieser Gleichung nach liefert eine Integration durch Trennung der Variablen

für alle t, wobei der natürliche Logarithmus ist. Insbesondere gilt damit

oder äquivalent

Die Hilfsfunktion hat wegen ein Maximum bei . Der maximal mögliche Anteil der Infizierten wird also für erreicht, d. h. hängt nur von der Basisreproduktionszahl und den Anfangswerten von und ab:

Für eine neu auftretende Erkrankung wie eine Epidemie durch ein unbekanntes Virus gilt und , d. h. der maximal mögliche Anteil der Infizierten an der Population hängt in diesem Fall wie folgt von der Basisreproduktionszahl ab:

Diese Gleichung entspricht der unter SIR-Modell #Maximale Zahl der Infizierten angegebenen Gleichung für das SIR-Modell.

Ausmaß der Epidemie nach ihrem Abklingen

[Bearbeiten | Quelltext bearbeiten]

Im Zusammenhang mit der Basisreproduktionszahl steht auch, welcher Anteil der Population insgesamt infiziert wird, unter Annahme, die Epidemie würde ohne jegliche Quarantäne durchlaufen. Unter Heranziehung der Differentialgleichungen findet man

Für den Anfangswert ist demnach

Bei ist nun , und daher . Daraus resultiert die Gleichung

Dieser Zusammenhang, Final size equation genannt, besitzt eine über das SEIR-Modell hinausgehende Bedeutung. Algebraische Umformung führt zur Gleichung

Mit und lautet die letzte Gleichung , die mit der lambertschen W-Funktion nach umgestellt werden kann, d. h. , was wieder zurückersetzt und nach umgestellt

ergibt.

Für die praktische Berechnung des hier relevanten Teils der W-Funktion betrachtet man das quadratische Taylorpolynom der Funktion an der Stelle −1 und bestimmt von diesem die Nullstelle, auf welche noch einmalig die Fixpunktiteration angewendet wird. Das Resultat ist

Eine gute Näherung der W-Funktion erhält man nun als

wobei die n-te Iteration des Newtonverfahrens

ist. Für alle praktischen Bedürfnisse ist völlig ausreichend.

Exponentielle Anfangsphase

[Bearbeiten | Quelltext bearbeiten]

Am Anfang der Epidemie verläuft die Ausbreitung der Krankheit in guter Näherung exponentiell. Eine Untersuchung dieser anfänglichen Phase gelingt mit der folgenden systematischen Methode, die auch auf andere Modelle anwendbar ist. Der Vektor fasst die vier Anteile zusammen, so dass das autonome System in der abstrakten Form dargestellt werden kann. An einem Zustand nähert die Linearisierung

das System an, wobei mit die Jacobi-Matrix gemeint ist. Der Anfang der Epidemie liegt am krankheitsfreien Gleichgewicht , das als Ruhelage die Gleichung erfüllt, also eine Nullstelle des Vektorfeldes ist. Die Jacobi-Matrix besitzt in dieser Ruhelage den Wert

Die vierte Zeile des Systems darf entfallen, da sie von den übrigen entkoppelt vorliegt, womit sich die Jacobi-Matrix auf drei Zeilen und Spalten verkleinert. Im linearisierten System ist zudem die erste Zeile entkoppelt, womit eigentlich die Betrachtung des kleinen Systems

genügt. Die allgemeine Lösung des linearisierten Systems erhält man mit Standardmethoden. Man kann nachrechnen, dass die Jacobi-Matrix für immer diagonalisierbar ist. Die Lösung ist demzufolge von der Form , wobei der Eigenvektor zum Eigenwert ist und die Konstanten aus der Anfangsbedingung zu ermitteln sind.

Diese allgemeine Lösung gilt auch für sonderbare Anfangswerte. Einer der Summanden fällt jedoch rasch auf null ab, weshalb man ein wenig später für die Neuinfektionen eigentlich nur eine exponentielle Kurve beobachtet, womit sich die Lösungen ebenfalls zu Exponentialfunktionen vereinfachen. Das Verhältnis von Exponierten zu Infizierten strebt mit dem Abfallen des Summanden einen speziellen Wert an. Weil infolge aufgrund gilt, finden sich die natürlichen Anfangswerte

Weiterhin ergibt sich nun, dass die Wachstumskonstante als Eigenwert der Jacobi-Matrix die quadratische Gleichung

erfüllen muss, womit eine Beziehung zwischen der Wachstumskonstante und den Parametern des Modells gefunden ist.[6] Wie bei jedem exponentiellen Wachstumsvorgang steht die Wachstumskonstante in eins-zu-eins-Beziehung zur anschaulicheren Vervielfachungszeit

speziell zur Verdopplungszeit .

Die Linearisierung an der Ruhelage ist überdies für die Untersuchung der Stabilität von Bedeutung. Besitzt die Jacobi-Matrix keinen Eigenwert mit Realteil null, spricht man von einer hyperbolischen Ruhelage. Laut dem Linearisierungssatz, einem Grundresultat der Stabilitätstheorie dynamischer Systeme, besitzt das linearisierte System in der Nähe der hyperbolischen Ruhelage dasselbe qualitative Verhalten wie das ursprüngliche. Somit ist eine hyperbolische Ruhelage genau dann stabil, wenn sich sämtliche Eigenwerte links der imaginären Achse befinden. Weil die vorliegende Jacobi-Matrix allerdings den Eigenwert null enthält, müsste man weitergehende Mittel zur Untersuchung der Stabilität einsetzen. Eine Umgehung dieses Problems bietet das später beschriebene Modell mit demografischer Dynamik im speziellen Fall, wo die Geburtenrate mit der Sterberate übereinstimmt. Dann ist N konstant bezüglich der Zeit, womit man das äquivalente System definieren darf, in dem N eine explizite Konstante ist. Der Eigenwert null tritt nun nicht mehr auf. Damit einhergehend verbleibt aber aufgrund des Versterbens der Resistenten nur noch die vollständig anfällige Population als einziges krankheitsfreies Gleichgewicht. Ferner tritt im instabilen Fall das stabile endemische Gleichgewicht in Erscheinung, das im einfachen Modell nicht existiert.

Die Untersuchung zur Stabilität wurde von van den Driessche in allgemeiner Weise durchgeführt. Sie kam zum Resultat, dass das krankheitsfreie Gleichgewicht bei einem deterministischen Kompartiment-Modell unter gewissen natürlichen Voraussetzungen an die Modelldefinition erwartungsgemäß immer dann instabil ist, wenn die Basisreproduktionszahl über dem kritischen Wert 1 liegt, und zumindest lokal asymptotisch stabil, wenn sie darunter liegt. Ebenso findet man bei Erläuterungen dieser Untersuchung aufgrund der besagten Problematik zunächst die Beschränkung auf solche Modelle vor, wo ein einziges krankheitsfreies Gleichgewicht existiert.

Beispielrechnung

[Bearbeiten | Quelltext bearbeiten]

Es folgt eine Beispielrechnung zu einer Parameterbelegung, wie sie im März 2020 für die COVID-19-Pandemie in Deutschland abgeschätzt wurde.[7] Die Ausbreitung verläuft bezüglich einer Basisreproduktionszahl von 2,4, was der Unterlassung wesentlicher Quarantäne entspricht. Zur numerischen Lösung des Anfangswertproblems genügt das Euler-Verfahren.

from numpy import array as vector

# Explizites Euler-Verfahren
def euler_method(f, t0, x0, t1, h):
    t = t0; x = x0
    a = [[t, x]]
    for k in range(0, 1 + int((t1 - t0)/h)):
        t = t0 + k*h
        x = x + h*f(t, x)
        a.append([t, x])
    return a

def SEIR_model(alpha, beta, gamma):
    def f(t, x):
        s, e, i, r = x
        return vector([
            -beta*s*i,
            beta*s*i - alpha*e,
            alpha*e - gamma*i,
            gamma*i
        ])
    return f

def SEIR_simulation(alpha, beta, gamma, e0, i0, days, step=0.1):
    x0 = vector([1.0 - e0 - i0, e0, i0, 0.0])
    f = SEIR_model(alpha, beta, gamma)
    return euler_method(f, 0, x0, days, step)

def diagram(simulation):
    import matplotlib.pyplot as plot
    plot.style.use('fivethirtyeight')
    figure,axes = plot.subplots()
    figure.subplots_adjust(bottom = 0.15)
    axes.grid(linestyle = ':', linewidth = 2.0, color = "#808080")
    t,x = zip(*simulation())
    s, e, i, r = zip(*x)
    axes.plot(t, s, color = "#0000cc")
    axes.plot(t, e, color = "#ffb000", linestyle = '--')
    axes.plot(t, i, color = "#a00060")
    axes.plot(t, r, color = "#008000", linestyle = '--')
    plot.show()

def simulation1():
    N = 83200000 # Einwohnerzahl von Deutschland 2019/2020
    R0 = 2.4; gamma = 1/3.0
    return SEIR_simulation(
        alpha = 1/5.5, beta = R0*gamma, gamma = gamma,
        e0 = 40000.0/N, i0 = 10000.0/N, days = 140)

diagram(simulation1)
Dynamics of an infectious disease, according to the basic SEIR model

Die vier Anteile, jeweils abhängig von der Zeit in Tagen.
s in Blau, e in Gelb gestrichelt, i in Magenta, r in Grün gestrichelt.

Ersichtlich ist an diesem Beispiel, dass die Epidemie aufgrund der Aufheizung durch die Infektiösen auch noch nach dem Erreichen der kritischen Immunisierungsschwelle

weiterläuft. Insgesamt würden sich 88 % der Bevölkerung mit der Krankheit anstecken. Die Epidemie ließe sich allerdings spätestens nach Erreichen der kritischen Immunisierungsschwelle durch eine ca. einmonatige strenge Quarantäne stoppen.

Modellerweiterungen

[Bearbeiten | Quelltext bearbeiten]

Einbeziehung demografischer Dynamik

[Bearbeiten | Quelltext bearbeiten]
Flussdiagramm

Unter der Annahme einer konstanten Geburtenrate und konstanten Sterberate wurde das erweiterte Modell

aufgestellt. Gegenüber dem einfachen SEIR-Modell beschreibt dieses Modell auch einen langfristigen endemischen Verlauf, bei dem es zu einer Schwingung des Anteils der Anfälligen kommen kann, bis sich dieser ins endemische Gleichgewicht eingependelt hat. Für verschwindende Latenzzeit geht das Modell ins korrespondierende SIR-Modell über, wo dieses Phänomen bereits auftritt.

Sofern sich die Geburten- und Sterberate voneinander unterscheiden, ist die Populationsgröße

keine Konstante mehr, sondern verläuft gemäß

Die Beziehung zur Basisreproduktionszahl ist

Ist das Wachstum der Bevölkerung oder ihr Schwund über den betrachteten Zeitraum verschwindend gering, kann man die Geburten- und Sterberate zur Vereinfachung gleich setzen. In diesem Fall pendelt sich der Anteil der Anfälligen genau in den Gegenwert der kritischen Immunisierungsschwelle ein.

Zeitabhängige Transmissionsrate

[Bearbeiten | Quelltext bearbeiten]

Vorgänge wie Verhaltensänderungen, Quarantäne und Saisonalität bewirken eine Veränderung der Transmissionsrate. Diese Umstände finden ihre Berücksichtigung in der Modellierung der Transmissionsrate als zeitabhängige Funktion , wobei das übrige Modell identisch beibehalten wird.[8] Die einfachsten Ansätze für die Saisonalität nehmen die Transmissionsrate zum Beispiel als Sinusschwingung an, mit Berg in den kälteren und Tal in den wärmeren Monaten.[9]

Zu bemerken ist, dass das Differentialgleichungssystem mit der expliziten Zeitabhängigkeit kein autonomes mehr ist und damit nicht mehr direkt ein dynamisches System vorliegt. Man kann aus dem System allerdings durch Hinzunahme der Gleichung künstlich ein autonomes gewinnen.

Interagierende Subpopulationen

[Bearbeiten | Quelltext bearbeiten]

Da soziale Kontakte, Schwere und Letalität der Erkrankung sowie die Wirksamkeit prophylaktischer Maßnahmen sich zwischen interagierenden Subpopulationen deutlich unterscheiden können, werden im „Interacting Subpopulation SEIR model“ individuelle SEIR Modelle für jede Subgruppe konstruiert und durch Kontakt-Interaktionen zwischen den Gruppen verknüpft[10]. Solche Modelle werden zum Beispiel zur Modellierung der Covid-19-Pandemie verwendet, um personalisierte, beschleunigte Impfstrategien zu entwickeln[11], welche eine Verkürzung der Pandemie und eine Reduktion von Fallzahlen und Todesfällen voraussagen.

Verallgemeinerungen

[Bearbeiten | Quelltext bearbeiten]

Infektionsalter-Modell

[Bearbeiten | Quelltext bearbeiten]

Die im SEIR-Modell vorausgesetzte Exponentialverteilung der Latenzzeit und der infektiösen Zeit mag recht künstlich erscheinen. Ein allgemeines, von dieser Voraussetzung befreites Modell beschreibt den zeitlichen Verlauf des Anteils der Anfälligen durch die Gleichung[12]

Hierbei ist die Infektionsrate im Infektionsalter . Mit dem Infektionsalter ist nicht das Alter eines Individuums gemeint, sondern die seit der Infektion vergangene Zeit. ist die Zeitableitung von .

Das Modell ist wie folgt interpretierbar. Ursache für die Neuinzidenz sind in der Vergangenheit infizierte Infektiöse. Um in der Zeit zurück kam hinzu, die sich nun im Infektionsalter befinden und daher zu beitragen. Das Integral über die gesamte Vergangenheit ist schließlich noch mit zu multiplizieren, weil zum aktuellen Zeitpunkt nur noch der Anteil anfällig ist.

Per Definition von ergibt sich die Basisreproduktionszahl nun gemäß der Formel

Die normierte Funktion lässt sich als Dichte einer Wahrscheinlichkeitsverteilung auffassen, die man Verteilung der Generationszeit, englisch generation interval distribution, nennt. Ihr Erwartungswert, die Generationszeit, ist der mittlere Abstand zwischen dem Zeitpunkt der Infektion des Spenders und dem Zeitpunkt der Infektion des Empfängers der Krankheit. Im Modell bleibt die Final size equation gültig, womit das Ausmaß der Epidemie allein von der Basisreproduktionszahl abhängt, nicht jedoch von der Dichte.

Zur Anfangsphase der Epidemie ergibt sich aus dem Modell für eine vollständig anfällige Population ( in guter Näherung) mit dem Ansatz zudem die epidemiologische Euler-Lotka-Gleichung

Diese besagt, bei Kenntnis der Dichte und der Wachstumskonstante ist die Basisreproduktionszahl ohne Kenntnis der Infektionsrate ermittelbar.

Um zu verstehen, in welchem Bezug die Modelle stehen, zerlegt man die Infektionsrate in das Produkt , wobei die Transmissionsrate ist, und die Wahrscheinlichkeit infektiös zu sein. Im Folgenden ist die Transmissionsrate konstant. Beginnt die infektiöse Zeit sofort und ist deren Dauer zufällig verteilt gemäß der Verteilungsfunktion , gilt . Speziell für ergibt sich das SIR-Modell, wobei die mittlere infektiöse Dauer nichts anderes als der Erwartungswert der Verteilung ist.

Dauert nun die Latenzzeit zufällig lang gemäß Verteilungsfunktion , und die danach beginnende infektiöse Zeit zufällig lang gemäß Verteilungsfunktion , erhält man

Speziell für und ergibt sich das SEIR-Modell. Die Erwartungswerte sind die mittlere Latenzzeit und die mittlere infektiöse Zeit.

Für einen nicht-zufälligen infektiösen Zeitraum mit Beginn im Infektionsalter und Ende in , also , vereinfacht sich das Modell stattdessen zu

eine retardierte Differentialgleichung. Die Beziehung zur Basisreproduktionszahl ist dann .

Stochastische Varianten

[Bearbeiten | Quelltext bearbeiten]
Waagerecht die Dichte der Verteilung des Zeitpunktes, an dem die Epidemie zu einer Wahl von Parametern im CTMC-Modell erlischt. Mit wachsender Basisreproduktionszahl trennt sich die Dichte in einen kleinen und einen großen Berg auf. Von einer Population von tausend Individuen sind zu Beginn vier infektiös.

Zur Beschreibung von Epidemie-Verläufen im Bereich weniger Infizierter wurde zu den klassischen Modellen jeweils eine direkte Entsprechung als zeitkontinuierlicher Markow-Prozess aufgestellt (CTMC: Continuous time Markov chain). In diesen stochastischen Modellen kann eine Epidemie vorläufig zum Erliegen kommen, obwohl die Reproduktionszahl größer als eins ist.

Für einen Markow-Prozess sei wie üblich

die Übergangswahrscheinlichkeit vom Zustand in den Zustand . Für das SEIR-Modell ist und[13]

Ein Verfahren zur Simulation dieser Markow-Prozesse ist der Gillespie-Algorithmus. Gemäß den Gesetzmäßigkeiten von Markow-Prozessen ist die Verweildauer im Zustand exponentialverteilt mit Verteilungsfunktion , wobei die Ereignisrate die Summe der Einzelraten ist, das heißt

Mittels der Inversionsmethode ist für eine Standardzufallszahl die Verweildauer bzw. der Zeitschritt demnach

Eine zweite Standardzufallszahl wählt aus, welches der paarweise disjunkten Ereignisse stattfinden soll:

  • Eine Infektion geschieht mit Wahrscheinlichkeit ,
  • ein Übergang ins Infektiöse mit Wahrscheinlichkeit ,
  • eine Erholung mit der restlichen Wahrscheinlichkeit .
Epidemie-Verläufe im DTMC-Modell. Aufgetragen ist jeweils die Anzahl der Infektiösen. Zu Beginn sind es zehn von tausend. Die schwarze Kurve zeigt den Erwartungswert, die gestrichelte den Wert im klassischen Modell.

Für einen kleinen Zeitschritt entsteht aus dem CTMC-Modell eine Näherung in Form eines zeitdiskreten Markow-Prozesses (DTMC: Discrete time Markov chain). Ist der Zeitschritt hinreichend klein, so dass mit ihm höchstens ein Ereignis stattfindet, darf man den Term ignorieren und erhält die Übergangswahrscheinlichkeiten analog zum CTMC-Modell.[14]

Die Simulation gestaltet sich besonders einfach. Eine Standardzufallszahl wählt aus, welches der paarweise disjunkten Ereignisse im Zeitschritt stattfinden soll:

  • Eine Infektion geschieht mit Wahrscheinlichkeit ,
  • ein Übergang ins Infektiöse mit Wahrscheinlichkeit ,
  • eine Erholung mit Wahrscheinlichkeit ,
  • Verharren im Zustand mit der restlichen Wahrscheinlichkeit.
  • Modellrechner für den Verlauf einer Epidemie nach dem SEIR-Modell mit änderbaren Input-Variablen (englisch)

Einzelnachweise

[Bearbeiten | Quelltext bearbeiten]
  1. Institute for Disease Modeling: SEIR Model (Memento des Originals vom 7. Mai 2020 im Internet Archive)  Info: Der Archivlink wurde automatisch eingesetzt und noch nicht geprüft. Bitte prüfe Original- und Archivlink gemäß Anleitung und entferne dann diesen Hinweis.@1@2Vorlage:Webachiv/IABot/www.idmod.org (abgerufen am 12. April 2020)
  2. a b Pauline van den Driessche: Reproduction numbers of infectious disease models. In: Infectious Disease Modelling, Band 2, KeAi Publishing, August 2017, S. 288–303. doi:10.1016/j.idm.2017.06.002. PMC 6002118 (freier Volltext).
  3. Odo Diekmann, Hans Heesterbeek, Tom Britton: Mathematical tools for understanding infectious disease dynamics. Princeton University Press, 2013, S. 35.
  4. Odo Diekmann, J. A. P. Heesterbeek, J. A. J. Metz: On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations. In: Journal of Mathematical Biology, Band 28, 1990, S. 365–382. doi:10.1007/BF00178324.
  5. Pauline van den Driessche, James Watmough: Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. In: Mathematical Biosciences, Band 180, Nr. 1–2, November–December 2002, S. 29–48. doi:10.1016/s0025-5564(02)00108-6.
  6. Junling Ma: Estimating epidemic exponential growth rate and basic reproduction number. In: Infectious Disease Modelling, Band 5, 2020, S. 129–141. doi:10.1016/j.idm.2019.12.009.
  7. Stellungnahme der Deutschen Gesellschaft für Epidemiologie (DGEpi) zur Verbreitung des neuen Coronavirus (SARS-CoV-2). (PDF) Deutsche Gesellschaft für Epidemiologie, 18. März 2020, abgerufen am 26. März 2020. Die der Beispielrechnung zugrunde liegende Version vom 19. März ist nicht mehr auffindbar. Hier ersatzweise die Version vom Vortag. Die aktuelle Version siehe hier (Memento des Originals vom 29. Mai 2020 im Internet Archive)  Info: Der Archivlink wurde automatisch eingesetzt und noch nicht geprüft. Bitte prüfe Original- und Archivlink gemäß Anleitung und entferne dann diesen Hinweis.@1@2Vorlage:Webachiv/IABot/www.dgepi.de.
  8. Gerardo Chowell, Cécile Viboud, Lone Simonsen, Seyed M. Moghadas: Characterizing the reproduction number of epidemics with early subexponential growth dynamics. In: Journal of the Royal Society Interface, Band 13, Nr. 123, 17. August 2016. doi:10.1098/rsif.2016.0659.
  9. M. Keeling, P. Rohani: Modeling Infectious Diseases in Humans and Animals. Abschnitt 5.2. (S. 159): Modeling forcing in childhood infectious diseases: Measles.
  10. Patrick Hunziker: Personalized-dose Covid-19 vaccination in a wave of virus Variants of Concern: Trading individual efficacy for societal benefit. In: Precision Nanomedicine. 24. Juli 2021, ISSN 2639-9431, doi:10.33218/001c.26101 (precisionnanomedicine.com [abgerufen am 18. August 2021]).
  11. Patrick Hunziker: Vaccination strategies for minimizing loss of life in Covid-19 in a Europe lacking vaccines. Infectious Diseases (except HIV/AIDS), 1. Februar 2021, doi:10.1101/2021.01.29.21250747 (medrxiv.org [abgerufen am 18. August 2021]).
  12. Fred Brauer: Mathematical epidemiology: Past, present, and future. In: Infectious Disease Modelling, Band 2, Nr. 2, 2017, S. 113–127. doi:10.1016/j.idm.2017.02.001. PMC 6001967 (freier Volltext).
  13. William Tritch, Linda J. S. Allen: Duration of a minor epidemic. In: Infectious Disease Modelling, Band 3, 2018, S. 60–73. doi:10.1016/j.idm.2018.03.002. PMC 6326226 (freier Volltext).
  14. Linda J. S. Allen: An Introduction to Stochastic Epidemic Models. In: Fred Brauer, Pauline van den Driessche, Jianhong Wu: Mathematical Epidemiology. In: Lecture Notes in Mathematics, Band 1945. Springer, Berlin/Heidelberg 2008, S. 81–130. doi:10.1007/978-3-540-78911-6 3.