Modèle épidémique SIR

Modélisez la dynamique épidémique : voyez comment le taux de transmission, le taux de guérison et R0 façonnent l'épidémie.

Chargement de la simulation interactive...

Leçon

La théorie — Modèle épidémique SIR

Une épidémie a deux fins différentes, et cette page affiche les deux. Le seuil d’immunité collective est le point où chaque infection cesse de produire plus d’un successeur — où la courbe bascule. Le taux d’attaque est la part de la population infectée une fois que tout est réellement terminé. Ce ne sont pas les mêmes nombres, ce ne sont jamais les mêmes nombres, et c’est de leur écart qu’il s’agit ici.

Ce que signifie chaque symbole

R₀
le nombre de reproduction de base, β/γ — combien de personnes une infection contaminerait dans une population encore entièrement réceptive. Un rapport de deux taux : il décrit donc un agent dans un contexte, jamais un agent seul.
s
la fraction réceptive S/N. Presque tout ce qui suit porte sur cette quantité et non sur le temps, ce qui explique que les résultats ne dépendent pas de la vitesse de l’épidémie.
H
le seuil d’immunité collective, 1 − 1/R₀ : la fraction immunisée à laquelle la croissance s’arrête. La page l’affiche.
A
le taux d’attaque — la fraction jamais infectée, une fois I retombé à zéro. La page l’affiche aussi, sous « total infecté ».

D’où vient la formule

  1. Demandons quand l’épidémie croît. Le compartiment infectieux évolue selon dI/dt = βIS/N − γI, qui se factorise en I(βs − γ). Comme I est positif pendant une épidémie, le signe est décidé par la seule parenthèse : l’infection croît tant que βs > γ, c’est-à-dire tant que s > γ/β = 1/R₀. Remarquez qu’il s’agit d’une condition sur la fraction réceptive, où le temps n’apparaît nulle part.
  2. Le point de retournement est donc exactement s = 1/R₀, et la fraction qui n’est plus réceptive à cet instant vaut 1 − 1/R₀. C’est le seuil d’immunité collective, et c’est le nombre qu’affiche la page — 33.3% pour le préréglage grippe, où R₀ = 1.5. C’est aussi là que culmine la courbe des infectieux, ce que le graphique permet de vérifier : le pic n’est pas un hasard de forme, c’est cette condition atteinte.
  3. Vient l’étape que l’on saute d’ordinaire. Au seuil, dI/dt = 0 — mais ce n’est pas I = 0. Une population infectieuse considérable est toujours là, et chacun continue de contaminer au taux βs. Moins d’un successeur chacun, donc les nombres baissent ; mais redescendre d’un pic prend du temps, et chacune de ces infections touche quelqu’un. L’épidémie ne s’arrête pas au seuil. Elle le traverse sur son élan.
  4. Pour trouver où elle s’arrête vraiment, il faut éliminer le temps. Divisons l’équation des réceptifs par celle des guéris : dS/dR = −(β/γ)(S/N) = −R₀s, qui s’intègre en s = s₀·e^(−R₀·r), où r est la fraction guérie. À la fin plus personne n’est infectieux, donc r∞ = 1 − s∞, et la substitution donne la relation de taille finale : s∞ = s₀·e^(−R₀(1 − s∞)). Le taux d’attaque vaut A = 1 − s∞. Elle est transcendante — s∞ figure dedans et dehors de l’exponentielle — donc sans forme close, et c’est pourquoi cette page intègre au lieu d’évaluer une formule.
  5. Mettez les deux nombres côte à côte : l’écart est le dépassement. À R₀ = 1.5 la page affiche un seuil de 33.3% et un total de 58.5% : l’épidémie file 25 points au-delà du point où elle a cessé de croître. Faites varier R₀ et le dépassement n’est pas monotone — il monte à environ 30 points près de R₀ ≈ 2 puis redescend, s’annulant quand R₀ → 1 parce que rien ne se propage, et aux grands R₀ parce que seuil et taux d’attaque se pressent tous deux contre 100%. À R₀ = 5 c’est 80.0% contre 99.3% ; à R₀ = 15, digne de la rougeole, 93.3% contre 100.0%. Le pire dépassement se situe exactement dans la plage de la grippe pandémique et du SARS-CoV-2 initial.

Comment lire ce que vous voyez

Trois courbes et trois nombres. Les courbes sont les compartiments ; celle à surveiller est la bosse des infectieux, dont le sommet est l’étape 2 en train de se produire. Parmi les nombres, « pic d’infectés » et « total infecté » répondent à des questions entièrement différentes : le premier dit combien sont malades en même temps, ce qui intéresse les hôpitaux, le second combien l’ont été à un moment, ce que prédit la relation de taille finale. Pour le préréglage grippe, 6.4% et 58.5%, et le premier ne dit rien du second. Le seuil d’immunité collective au-dessus n’est ni l’un ni l’autre : c’est là où le premier bascule.

Suppose
Une population fermée et bien mélangée — pas de naissances, pas de décès, pas de voyages, et chaque réceptif également susceptible de croiser chaque infectieux. Des β et γ constants : aucun changement de comportement, aucune saison, aucune intervention en cours de route. Une immunité définitive, personne ne retourne au réservoir réceptif. Et un modèle déterministe, raison pour laquelle une poignée d’infectieux déclenche toujours une épidémie ici, alors qu’en réalité les petits nombres s’éteignent souvent par hasard même quand R₀ > 1.
Ne tient plus quand
L’expression « seuil d’immunité collective » invite à une lecture que les mathématiques ne soutiennent pas : qu’une épidémie livrée à elle-même s’y arrête. Elle ne s’y arrête pas — elle s’arrête au taux d’attaque, que l’étape 5 montre pouvoir être 25 points plus haut. Le seuil est une cible vaccinale, un niveau atteignable sans épidémie ; atteint par l’infection, il coûte le seuil et le dépassement en plus. Cette distinction n’avait rien d’académique en 2020. Deux mises en garde de plus pour qui prend ces nombres au sérieux : le bon mélange est l’hypothèse la moins réaliste ici, et l’assouplir abaisse généralement seuil et taux d’attaque, parce que les individus très connectés sont infectés tôt et cessent de faire pont ; et cette page démarre avec dix infectieux plutôt qu’un germe infinitésimal, d’où ses 58.5% juste au-dessus des 58.3% de la limite classique s₀ → 1.

R₀ est un rapport de deux taux, donc jamais une propriété du seul pathogène 🖖

Le modèle SIR divise la population en Susceptibles, Infectés et Rétablis. Sa dynamique est régie par des équations différentielles non linéaires : dS/dt = -βIS/N et dI/dt = βIS/N - γI. R₀ = β/γ est le taux de reproduction de base. Si R₀ > 1, la maladie se propage de façon exponentielle avant d'atteindre l'immunité collective.

Pourquoi une épidémie s'arrête avant d'infecter tout le monde 🖖

La courbe des infections monte, atteint un pic, puis redescend — non pas parce que l'agent pathogène disparaît, mais parce que tant de personnes ont déjà été infectées que le virus rencontre sans cesse des individus déjà immunisés. C'est pourquoi une épidémie s'éteint bien avant d'atteindre tout le monde : avec R₀ = 1.5, seuls environ 58% de la population sont infectés, et non 100%.

Le théorème du seuil est né d'un chimiste aveugle 🖖

Le modèle SIR a été publié en 1927 par William Kermack et Anderson McKendrick. Trois ans plus tôt, une explosion de laboratoire avait rendu Kermack totalement aveugle à l'âge de 26 ans ; il a pourtant contribué à démontrer le théorème du seuil épidémique, résolvant une grande partie des calculs de tête, dit-on. Leur article est resté largement ignoré pendant des décennies avant de devenir un pilier de l'épidémiologie moderne.

Problème entièrement résolu

  1. Une infection ayant R₀ = 3 se propageant au-delà de la fraction immunisée 5 étapes

    Une infection a un R₀ = 3. Trouvez la fraction qui doit être immunisée pour stopper sa propagation — et montrez que l'épidémie continue de progresser au-delà de ce point.

    1. R₀ est le nombre moyen de nouveaux cas qu'un cas génère dans une population entièrement réceptive : le taux de transmission divisé par le taux de guérison. C'est le rapport de deux taux, c'est pourquoi il n'a ni unité ni échelle de temps.

    2. Une épidémie progresse tant que chaque cas se remplace. Une fois qu'une partie de la population est immunisée, seule la part de sujets réceptifs compte, de sorte que la croissance se poursuit exactement tant que cette part dépasse 1/R₀.

    3. Inversez-le. Le simulateur ci-dessus affiche ce seuil, et c'est tout le fondement de la politique de vaccination — il n'est jamais nécessaire d'atteindre tout le monde.

    4. Mais le seuil est l'endroit où la croissance s'arrête, pas là où l'épidémie s'arrête. Les personnes infectées sont toujours contagieuses à ce moment-là, de sorte que les cas continuent de s'accumuler au-delà.

    5. Ce dépassement explique pourquoi le pic et le seuil sont des chiffres différents, et pourquoi une épidémie non atténuée se termine avec une proportion d'infectés bien supérieure au seuil.

    Réponse

    66.7%. Le pic survient avec environ 30% d'infectés en même temps, et — le point le plus facilement ignoré — l'épidémie se poursuit au-delà du seuil au lieu de s'y arrêter, de sorte que la taille finale dépasse largement ce niveau. L'immunité collective est une affirmation sur le moment où la courbe s'inverse, et non sur la fin de l'épidémie. Et comme le seuil est 1 − 1/R₀, il est extrêmement sensible au sommet : la rougeole avec un R₀ ≈ 15 nécessite une couverture de 93.3%, c'est pourquoi la rougeole réapparaît en premier dès que la couverture vaccinale baisse.

Parcours

Croissance limitée

Références (3)

Exemples de problèmes

  • Rougeole (R₀≈15) - Rougeole : extrêmement contagieuse (R₀ ≈ 15), entraînant des flambées rapides.
  • Grippe (R₀≈1,5) - Grippe saisonnière : propagation modérée (R₀ ≈ 1.5), avec un pic durant les mois d'hiver.
  • COVID Delta - COVID-19 Delta : modélisé avec une phase d'incubation SEIR (R₀ ≈ 5.0).
  • Subcritique - Cas sous-critique : R₀ < 1.0, la maladie s'éteint naturellement sans provoquer d'épidémie.