Les mathématiques des épidémies

Jusqu'où la courbe des infectés et des décès va-t-elle grimper? Quelle est l'efficacité des mesures de confinement décidées dans la plupart des pays? Peut-on atteindre l'immunité collective? Comment mener le déconfinement? Y aura-t-il une seconde vague? Qu'est-ce que le taux de reproduction, annoncé lundi soir à 0,6 par le gouvernement, alors qu'il valait 4 avant le confinement?
Les modèles mathématiques des épidémies donnent des éléments pour répondre à ce type de questions. Je voudrais expliquer ce que les modèles les plus simples nous disons. Des versions plus élaborées permettent de faire des prévisions, ont convaincu nos dirigeants politiques de confiner la population, malgré la crise économique qui va en résulter. Sur verra en particulier le rôle du nombre de reproductions de la base R0, qui est le nombre moyen de personnes à qui un infecté transmet la maladie au début de l'épidémie. Sa valeur exacte fait encore débat pour le Covid-19. Mais surtout elle dépend non seulement de la façon dont la maladie se transmet, mais aussi du nombre de contacts que chaque individu à chaque jour. C'est pour faire baisser ce nombre en dessous de 1 (on verra pourquoi ci-dessous) que l'on nous confine. Une vidéo complète cet article pour les plus curieux. Certains des passages soulignés renvoient à une version plus complète (avec plus de mathématiques) de cet article, disponible en ligne.
Un modèle simple pour cerner la situation
Le modèle de Reed-Frost est un des plus anciens modèles mathématiques des épidémies, il date de 1929. Il est très simpliste, mais il permet d'introduire des notions essentielles et d'obtenir une formule mathématique importante. Les individus sont de trois types: S comme sensible (être infecté), I comme infecté et infectieux (capable d'infecter un sensitive), R comme remis ou retiré (soit guéri, soit mort). Dans ce modèle, dit «SIR», le temps est discret (sur progresse semaine par semaine, par exemple); la population est supposée grande, sa taille est n. Au début, il y a n-1 individus de type S, 1 individu de type I et 0 de type R. Un individu qui est infecté une semaine donnée infectera chaque S avec les réponses p la semaine suivante, puis guérira (devenant R). L'épidémie se poursuit tant qu'il y a des infectés (que je n'est pas nul), puis elle s'arrête.
Deux remarques: on néglige pour simplifier la phase d'incubation; on suppose qu'un individu de type R, s'il n'est pas mort (ce qui heureusement est le cas de l'immense majorité des R), est immunisé - en réalité, à ce stade, on ne sait pas encore beaucoup de choses sur l'immunité au Covid-19.
Le nombre de reproductions de base, noté R0, est le nombre moyen de susceptibles qu'un infecté infecté au début de l'épidémie, c'est-à-dire dire «quand presque toute la population est sensible».
Dans le modèle de Reed-Frost, combien vaut R0? Donnons des valeurs à nos deux paramètres. Supposons que n = 1 000, et que p = 0,0025 (soit 0,25%). Le premier infecté a autour de lui n-1 ≃ n = 1 000 individus sensibles. Puisqu'il infecte chacun d'eux avec la réponse p, ce nombre moyen vaut ici R0 = n × p = 2,5.
Il est assez clair que si R0 <1, il n'y a pas d'épidémie majeure avec un petit nombre d'infectés initiaux, de même si R0 = 1. Par contre si R0> 1, un seul infecté initial peut déclencher une épidémie majeure, qui touche une fraction importante de la population.
Quelle fraction de la population totale sera touchée par l'épidémie?
Une question importante, que vous voyez beaucoup de gens au vu de la croissance exponentielle du nombre d'individus touchés, et qui nourrit les discussions entre épidémiologistes modélisateurs et responsables politiques, est «quelle fraction de la population totale sera touchée par l'épidémie» ? Si l'on admet qu'une personne guérie est immunisée, il est clair que l'épidémie va s'arrêter tôt ou tard. En outre, tout le monde n'est pas touché. Lorsqu'une fraction de la population est immunisée, dans le modèle de Reed-Frost, pour savoir combien de susceptibles un infecté infecté en moyenne, il faut multiplier p par le nombre S de susceptibles restant, et bien avant que ce nombre ne s ' annule, le produit p × S passe en dessous de 1, et alors l'épidémie s'arrête.
Voici ce que donnent 10000 simulations dans les cas R0 = 0,95, et R0 = 2,5. Lorsque R0 <1, aucune épidémie majeure n'a lieu. Dans le cas R0 = 2,5, une certaine fraction des simulations (qui ne dépend pas de la taille de la population) n'aboutissent pas à une épidémie majeure, tandis que la proportion d'individus infectés dans le cas d'une épidémie majeure se concentre quand la taille de la population augmente autour d'une certaine valeur (qui augmente avec R0, comme le sont les simulations pour plusieurs valeurs de R0
La hauteur de chaque barre indique le nombre de simulations produisant la proportion de la population totale qui est infectée, indiquée horizontalement. n est la taille totale de la population. Étienne Pardoux, auteur fourni
Quelle est la fraction de la population touchée par une épidémie majeure? Supposons que R0> 1. Dans la réalité, il y a plusieurs infectés initiaux, au moins en Europe de l'Ouest et aux États-Unis, suffisamment pour que l'épidémie majeure soit inévitable. On peut montrer que l'épidémie touche au total une fraction τ de la population, qui pour une grande taille de population, est proche de la solution de l'équation 1-τ = e-R0 × τ. Cette équation a toujours la solution τ = 0, et si R0> 1, il y a une autre solution, qui correspond à une épidémie majeure. Sur les valeurs numériques suivantes:
pour R0 = 2, τ = 79%; pour R0 = 3, τ = 94%; pour R0 = 4, τ = 98%; et pour R0 = 5, τ = 99%.
La preuve de ce résultat, dans le cas du modèle de Reed-Frost est assez élémentaire En fait sur peut être établi dans un cadre plus général, mais avec quand même les trois hypothèses suivantes:
Pas d'immunité au départ (naturelle ou par vaccination). Vrai pour le SARS-CoV-2, pas pour la grippe saisonnière. Si la moitié de la population était immunisée, la taille totale serait très réduite!
Une communauté homogène, ce qui n'est vrai ni pour le SRAS-CoV-2 ni pour aucune maladie. Quelles hétérogénéités sont à prendre en compte dépend de chaque maladie et de son mode de transmission. On considère qu'il faut réduire le nombre d'infectés I d'environ 10% à 20% par rapport aux prédictions du modèle homogène.
Sur une supposé que le comportement des individus ne change pas au cours de l'épidémie. Bien sûr les mesures prises par les autorités (fermeture des écoles et des lieux publics, confinement) changent la donne.
On peut tracer la courbe de τ en fonction de R0.
Taille totale de l'épidémie en fonction de R0. Notre modèle est en bleu. Si on prend en compte l'hétérogénéité, en rouge pointillé; et enfin si la moitié de la population était immunisée, en jaune. Étienne Pardoux, auteur fourni
Comment réduire R0, le nombre de reproductions de base? On peut factoriser R0: R0 = p × c × l, où est la réponse de transmission à chaque «contact», c est le nombre de «contacts» par jour et la durée de la période d'infection. Alors, les mesures de prévention visent à réduire R0 en obtenu:
la réponse de transmission à chaque «contact» p: masques, lavage des mains;
le nombre de «contacts» par jour c: confinement, éviter les transports publics et les grands rassemblements;
la durée de la période d'infection l: diagnostic rapide, isolation des infectés, traitement efficace (si disponible).
Il n'y aura pas d'épidémie majeure si le nombre de reproductions de base peut être réduit à une valeur inférieure à 1 grâce à diverses mesures. Conclusion: il faudrait des masques pour tout le monde (réduction de p) et tester massivement (réduction de l, comme le montre l'exemple de la Corée).
Si une proportion v de la population est vaccinée, alors R0 est remplacé par Rv = (1-v) R0 (chaque tentative d'infecteur ”ne réussit pas à obtenir les réponses 1-v): Rv est inférieur à 1 si v> 1-1 / R0. Deux conclusions sur la vaccination (qui n'est malheureusement pas encore disponible pour le Covid-19): d'une part, si on connaît R0, on connaît la fraction de la population qu'il faut vacciner pour être sûr qu'il n 'y aura pas d'épidémie majeure - par exemple, si R0 vaut 2,5, il faut vacciner 60% de la population, 75% si R0 vaut 4, 80% si R0 vaut 5. D'autre part, en se vaccinant , sur se protège soi-même, mais surtout sur protège les autres.
Observateur la courbe de propagation de l'épidémie ne suffit pas pour estimer R0 car le taux de croissance de la courbe dépend de deux facteurs: R0 et le «temps de génération» G, c'est-à-dire la durée entre le moment où l'on est infecté et celui où l'on infecte.
Le taux de croissance croît avec R0, décroît avec G.Il existe des formules qui relient r, R0 et la loi des réponses de G. Par une méthode de «suivi des contacts» (recherche de qui a infecté qui et quand; ceci a été pratiqué avec succès aux Contamines en février, pour juguler ce foyer d'épidémie), sur peut obtenir des informations sur la loi éprouvée de G, d'où, en combinant avec reconnu de r, sur peut estimer R0!
La dynamique de l'épidémie
Jusqu'ici sur un surtout discuté la fraction de la population touchée par l'épidémie. Pour décrire l'évolution de l'épidémie, dans le cas d'une grande population, et partant d'une situation où une toute petite fraction de la population est touchée, on peut utiliser un modèle déterministe SIR, qui s'écrit sous la forme d'un système d'équations différentielles, dont je présente ci-dessous graphiquement quelques solutions représentatives. On voit sur les quatre premières courbes que plus R0 est grand, plus la vague des infectés arrivent tôt et monte haut. En outre, la réduction de la hauteur du pic est plus marquée que la réduction de la «taille totale de l'épidémie», l'aire sous la courbe), lorsque l'on réduit R0, voiture la courbe s'étale.
Si on part avec R0 = 2,5, puis que l'on applique un confinement au bout de la troisième semaine, l'effet de réduction de la vague est très importante, même si le taux de reproduction après confinement reste supérieur à 1. Dans la réalité, sur n'observe pas une modification si rapide de la courbe, à la fois parce que le confinement ne se rencontre pas bien en place instantanément, et surtout parce que notre modèle ne prend pas en compte les effets de retard inhérents à l'épidémie (temps de génération G, période d'incubation).
L'évolution du nombre d'infectés. À gauche, sans confinement et pour R0 = 3 (en noir), R0 = 2,5 (en bleu), R0 = 2 (en rouge), R0 = 1,5 (en vert). À droite, R0 = 2,5 en noir puis confinement après la 3ᵉ semaine, soit un changement de R0: R0 = 1,5 (en bleu), R0 = 1 (en rouge), R0 = 0,85 (en vert) . Étienne Pardoux, auteur fourni
L'épidémie actuelle du Covid-19
Concernant le Covid-19 en France, il y a encore beaucoup d'incertitudes quant à la valeur de R0. Certaines estimations tournent autour de 2,5, d'autres donnent des valeurs plus importantes, entre 3 et 5! Il est probable que R0 ne prenne pas la même valeur dans les zones peu peuplées, et en région parisienne par exemple où beaucoup de gens se pressent dans les transports en commun. Si le gouvernement n'avait pas pris de mesures fortes, on peut penser qu'entre 60 et 80% de la population aurait été touchée. Et surtout la vague serait arrivée très vite et serait montée très haut, submergeant le système de santé et augmentant la mortalité.
Ceci ne règle pas le problème de comment sortir du confinement. Une des difficultés est qu'à ce jour sur n'a pas une bonne estimation du nombre de personnes en France qui ont été touchées par le virus, c'est-à-dire I (t) + R (t), si t est la date d'aujourd'hui. Sur sait combien ont eu des symptômes assez sérieux pour se rendre à l'hôpital, mais une des particularités de cette maladie est qu'elle est bénigne et même non décélée dans la majorité des cas Et sur ne connaît pas avec précision la proportion de ces cas. Donc on ne sait pas aujourd'hui quelle fraction de la population est immunisée, informations qui seront importantes pour la sortie du confinement.
Dès que l'on dispose de tests sérologiques fiables de détection des anticorps, sur pourra tester par sondage un certain nombre de Français, pour avoir cette information.
Que se passe-t-il dans deux mois, trois mois, six mois? Aurons-nous un retour à la normale relativement rapide, comme cela à l'air d'être le cas en Chine? Au niveau européen? Au niveau de la planète? Y aura-t-il des rebonds? Chacun peut formuler des conjectures, mais je pense que personne, à ce stade, n'a de certitude.
Les épidémies et la recherche mathématique
Les épidémies constituant un champ d'application important des mathématiques. Les événements que nous vivons attirent l'attention du public sur cette branche de la modélisation mathématique, et sont aussi l'occasion pour les spécialistes d'accélérer leurs travaux, pour tenter de répondre à des questions qui conditionnent la survie de dizaines, voire de plusieurs milliers de gens. Bien sûr, le plus important est la recherche de médicaments efficaces et d'un vaccin, mais les modèles mathématiques restent indispensables pour les décideurs. Voici quelques-uns des problèmes sur lesquels travaillent les mathématiciens et statisticiens qui recherchent les modèles des épidémies:
L'hétérogénéité de la population (mais vues et lieux de travail, répartition spatiale, les voyages…) et ses effets.
Faut-il privilégier les modèles les plus réalistes, ou les plus simples (dont on peut mieux tirer des conclusions)? Leurs réponses sont-elles très différentes?
Quel est l'effet des diverses mesures préventives?
Que peut-on estimer à partir des données disponibles, et quelles données supplémentaires sont utiles utiles?