Generalized SEIR model and a new reproduction number reflecting the real epidemic dynamics

Technische Samenvatting: Gegeneraliseerd SEIR-model en een nieuw reproductienummer dat de werkelijke epidemische dynamiek weerspiegelt

Probleemstelling
Bestaande wiskundige modellen voor epidemische dynamiek, zoals het klassieke SIR-model, houden vaak geen rekening met kritieke factoren zoals de incubatieperiode (de blootgestelde compartiment), verborgen (asymptomatische) infecties, herinfecties en de effecten van vaccinaties en pasgeborenen. Eerdere door de auteur voorgestelde modellen adresseerden verborgen dynamiek en herinfecties, maar verwaarloosden het blootgestelde compartiment ($E$). Bovendien berusten standaarddefinities van het effectieve reproductienummer ($R_t$) uitsluitend op zichtbare (symptomatische) gevallen, wat de werkelijke snelheid van infectie verspreiding mogelijk niet nauwkeurig weergeeft wanneer een aanzienlijk deel van de populatie asymptomatisch is of niet geregistreerd staat. Dit artikel adresseert de behoefte aan een uitgebreider model dat het blootgestelde compartiment integreert met verborgen dynamiek en stelt een nieuwe metriek voor die de werkelijke epidemische trends beter weerspiegelt.

Methodologie
De auteur stelt een Gegeneraliseerd SEIR-model (Susceptible-Exposed-Infectious-Removed) voor, bestaande uit zes gekoppelde differentiaalvergelijkingen. Dit model breidt eerder werk uit door een compartiment voor blootgestelde personen ($E$) te introduceren en zowel de infectieuze ($I$) als de verwijderde ($R$) populaties op te delen in zichtbare (geregistreerde) en verborgen (niet-geregistreerde/asymptomatische) subcompartimenten.

Belangrijke kenmerken van het wiskundige kader zijn onder meer:

• Zes Differentiaalvergelijkingen: Die de dynamiek van Susceptibele ($S$), Blootgestelde ($E$), Zichtbaar Infectieuze ($I^v$), Verborgen Infectieuze ($I^h$), Zichtbaar Verwijderde ($R^v$) en Verborgen Verwijderde ($R^h$) personen beheersen.
• Parameters: Het model incorporeert 14 onbekende parameters per epidemische golf, waaronder infectiesnelheden ($\alpha$), verwijderingssnelheden ($\rho$), afname van immuniteit ($\delta$), de zichtbaarheidscoëfficiënt ($\beta$) en een coëfficiënt ($\mu$) die de netto veranderingssnelheid van de susceptibele populatie vertegenwoordigt door geboorten, vaccinaties en herinfecties.
• Parameteridentificatie: De auteur stelt een methode voor om deze 14 parameters te identificeren met behulp van de methode van kleinste kwadraten, waarbij het model wordt aangepast aan geobserveerde dagelijkse, wekelijkse of maandelijkse aantallen geregistreerde gevallen ($D_j$). Het model relateert de afgeleide van de zichtbare verwijderde gevallen aan de geobserveerde gegevens.
• Nieuw Reproductienummer ($R_\tau^E$): Een nieuw reproductienummer wordt gedefinieerd op basis van de ratio van blootgestelde personen op tijdstip $t$ versus $t-\tau$. Cruciaal is dat het artikel aantoont dat deze snelheid kan worden geschat met enkel de ratio van dagelijkse (of wekelijkse) geregistreerde gevallen, waardoor de zichtbare gegevens effectief worden gekoppeld aan de dynamiek van de verborgen blootgestelde populatie.

Belangrijkste Bijdragen

1. Gegeneraliseerde Modelstructuur: De introductie van een zes-vergelijkingen SEIR-model dat gelijktijdig rekening houdt met de incubatieperiode, verborgen infecties, herinfecties en demografische veranderingen (geboorten/vaccinaties).
2. Parameterschattingsmethode: Een specifieke procedure voor het schatten van de 14 parameters van het model met behulp van standaard epidemiologische gegevens (gecumuleerde of dagelijkse casus-aantallen) zonder directe observatie van verborgen gevallen te vereisen.
3. Herzien Reproductienummer: Het voorstel van een reproductienummer ($R_\tau^E$) dat de toenemende of afnemende snelheid van blootgestelde personen weerspiegelt. Het artikel stelt vast dat de ratio van opeenvolgende dagelijkse (of wekelijkse) geregistreerde casus-aantallen dient als een geldige estimator voor deze snelheid, wat een praktisch instrument biedt voor epidemische beheersingsbeslissingen.
4. Evenwichtsanalyse: Het artikel biedt een analyse van evenwichtspunten, waarbij wordt opgemerkt dat een niet-nul evenwicht (endemische staat) alleen mogelijk is wanneer de netto groei van de susceptibele populatie ($\mu$) niet nul is.

Resultaten
Het model werd toegepast om de kink (pertussis) epidemie in Engeland van 2023 tot 2025 te simuleren.

• Simulatiescenario's:
  • $\mu = 0$: Simulaties die uitgaan van geen netto groei in de susceptibele populatie (bijv. perfecte vaccinatie die geboorten balanceert) resulteerden in een monotone afname van gevallen na een piek, wat wijst op een beperkte epidemische duur.
  • $\mu > 0$: Simulaties die een positieve netto groeisnelheid incorporeerden (weerspiegelend dat geboortecijfers de vaccinatie-effectiviteit overtreffen) produceerden oscillerende epidemische golven met een periode van ongeveer 2500 dagen (6,85 jaar).
• Validatie: Het model met $\mu = 7,48 \text{ dag}^{-1}$ vertoonde een sterke overeenstemming met historische gegevens en recente observaties (2025) van kinkgevallen in Engeland.
• Voorspelling: Op basis van het gekalibreerde model voorspelt het artikel dat het aantal nieuwe kinkgevallen in Engeland opnieuw zal beginnen toe te nemen in juli–september 2027. Het model geeft aan dat het aantal blootgestelde personen onmiddellijk zal stijgen na het minimum in de dagelijkse geregistreerde gevallen.

Betekenis en Claims
Het artikel beweert dat het voorgestelde Gegeneraliseerde SEIR-model een realistischerere weergave biedt van epidemische dynamiek door de integratie van de incubatieperiode en verborgen transmissie. De primaire betekenis ligt in de mogelijkheid om:

• Toekomstige epidemische golven (specifiek de kink-resurgentie in Engeland in 2027) met hoge nauwkeurigheid te voorspellen met enkel geregistreerde casusgegevens.
• Een praktische, op data gebaseerde definitie van het reproductienummer te bieden die rekening houdt met verborgen dynamiek, waardoor epidemiologen de werkelijke snelheid van infectie verspreiding kunnen monitoren met behulp van standaard surveillancegegevens.
• Aan te tonen dat hoewel intensieve vaccinatie theoretisch de duur van een epidemie kan verkorten, het mogelijk onvoldoende is om langdurige oscillaties te voorkomen als geboortecijfers en afname van immuniteit niet adequaat worden gecompenseerd.

De auteur merkt op dat het grote aantal parameters (14) de optimalisatie bemoeilijkt, wat suggereert dat toekomstig onderzoek gebruik moet maken van High Performance Computing (HPC) en AI-methoden voor grootschalige parameteridentificatie. De studie concludeert dat het model er succesvol in is geslaagd de kloof tussen theoretische compartimentale dynamiek en observeerbare real-world data te overbruggen.