Archives pour l'étiquette covid19

Le jeu de la vie et celui des épidémies

Le jeu de la vie, inventé en 1970 par John Conway, n’est pas vraiment un jeu. Ce terme est cependant moins rébarbatif que celui d’automate cellulaire, qui est pourtant plus exact. Il trouve ses origines dans des travaux conduits par John von Neumann dans les années 1940. Nous garderons la métaphore du jeu pour en parler, même si certains trouveront le terme mal adapté quand il s’agit de maladies potentiellement mortelles. L’essentiel est d’aider la compréhension. Voyons quelles en sont les règles.

Les règles du jeu de la vie

Pour jouer, prenez un damier et des pions. Les cases sont considérées comme des cellules ; elles peuvent être mortes ou vivantes. On utilise les pions pour matérialiser les cellules vivantes. Au début du jeu, on place des pions sur n’importe quelle case. On joue ensuite par étapes selon les règles suivantes :

— une cellule morte entourée de trois cellules vivantes ressuscite, sinon elle reste morte ;

— une cellule vivante reste en vie si elle a deux ou trois voisines vivantes, sinon elle meurt.

Bien que l’évolution du jeu soit complètement déterminée par la disposition initiale des cellules, on n’en assiste pas moins à quelques situations qui peuvent paraître surprenantes. Ainsi, en alignant tout simplement trois cellules vivantes les unes à côté des autres, on obtient une situation où les trois cellules se reproduisent, alignées horizontalement puis verticalement et ainsi de suite.

Lorsque trois cellules vivantes sont contiguës, on assiste à une oscillation entre trois cellules en ligne et trois en colonne.

Le jeu des épidémies

 

Ce jeu est loin d’être un simple amusement : il s’agit d’un exemple de ce que l’on nomme « automate cellulaire », particulièrement utile pour modéliser les processus d’expansion des épidémies comme des épizooties. En préalable à ce type d’application, il est nécessaire d’étendre le damier à l’infini. Au départ, toutes les cellules sont saines. On place une cellule infectée puis on « joue » avec la règle probabiliste suivante :

— les cellules voisines de la cellule infectée sont infectées au coup suivant avec la probabilité ;

— la cellule meurt ou est immunisée le coup suivant.

Comment les cellules infectées (en rouge) se multiplient-elles au détriment des cellules saines (en vert) ? Dans cet exemple, la probabilité qu’une cellule voisine d’une cellule infectée soit infectée à son tour est de 25 %. Les cellules mortes ou immunisées sont représentées en bleu.

La question qui intéresse autant les épidémiologistes que le grand public est donc : « Pour quelles valeurs de p, la maladie se propage-t-elle au monde entier ? »

Un modèle probabiliste

Le modèle est ici « probabiliste », et donc on ne peut prédire à l’avance ce qui va se produire dans un cas particulier. Pour avoir une idée rapide de l’évolution moyenne du système, le mieux est de procéder à une simulation. Pour cela, on « joue » selon les règles énoncées ci-dessus en utilisant un générateur de nombres pseudo-aléatoires et on comptabilise le nombre de cellules infectées. En jouant cent fois de suite et en faisant la moyenne des résultats, on obtient une mesure de l’expansion moyenne de l’épidémie.

Taux critique

En dessous d’un certain taux de contamination p, l’épidémie ne s’étend pas. En revanche, au dessus de ce taux, elle envahit le monde entier. Dans le cadre de notre modèle simplifié, le taux critique se situe entre 30 % et 40 %. Une maladie ne devient épidémique que si ce taux est dépassé. Comment ce modèle peut-il être adapté pour bien modéliser différents types d’épidémies ou d’épizooties ? Tout d’abord, on peut modifier le voisinage de chaque cellule, composé ici de huit cellules — les spécialistes parlent de voisinage de Moore, du nom d’Edward Moore, l’un des fondateurs de la théorie des automates. On utilise souvent un voisinage plus simple, dit de von Neumann, constitué des quatre cellules partageant un côté avec la cellule considérée. Avec ce nouveau modèle, le taux critique pour lequel une maladie devient épidémique se situe aux alentours de 60 %. On peut également améliorer le modèle en tenant compte du temps pendant lequel une cellule infectée est contagieuse puis du taux de mortalité et d’immunité ainsi que du temps d’immunité. On arrive ainsi à retrouver la façon dont se sont propagées des épidémies comme la peste dans l’Europe médiévale. Une première vague a tué le tiers de la population en se propageant à partir d’un épicentre situé dans un port, suivie de plusieurs répliques plus faibles, toutes partant du même point. Ces répliques correspondent à la fin de certaines immunités.

La confrontation avec les données épidémiologiques a permis de montrer que ce type de modèles a une certaine pertinence pour toutes les maladies qui se propagent par contact direct : grippe, tuberculose, coronavirus ou même sida. En revanche, il ne fonctionne plus lorsque la maladie se propage via un agent infectieux, comme dans le cas du paludisme ou du chikungunya.

Géométrie des contagions

Comment considérer maintenant la notion de « cellule voisine » dès que l’on évoque les réseaux de transports aériens, maritimes ou terrestres ? Dans le cas d’une épidémie de grippe humaine, l’aéroport de Paris est voisin de celui de Hong-Kong. Dans le cas d’une épizootie de grippe aviaire, deux élevages fréquentant le même marché aux bestiaux sont voisins. On doit de plus tenir compte des migrations naturelles des oiseaux sauvages. Dans tous ces cas, on retrouve la notion de réseaux.

En modifiant le modèle du jeu, on peut passer du cas où chaque cellule représente un individu à celui où elle représente un domaine où les individus sont en relation constante : un élevage de volaille dans le cas de la grippe aviaire, une ville dans le cas de la tuberculose, du sida ou de la grippe humaine. Ces domaines sont reliés entre eux pour former un réseau. Dans chaque cellule, la modélisation suit une autre logique, celle du modèle « SIR » dû à William Kermack et Anderson Mac Kendrick en 1927 (voir l’article correspondant sur ce blog). Ce modèle compartimente la population en trois classes : S, la classe des individus susceptibles d’attraper la maladie, I, celle de ceux qui en sont infectés (et contagieuses) et R, ceux qui en sont guéris (et immunisés) ou décédés.

Seuil de propagation

On considère l’évolution de ces trois classes dans le temps en fonction de deux taux mesurables expérimentalement. Le premier (a) est le taux de contagion de la maladie pour un infecté, c’est-à-dire la probabilité pour qu’un individu susceptible attrape la maladie après contact avec un individu infecté. Le second taux (b) mesure le passage de l’état I à l’état R.

Après un laps de temps t, on compte a I S t infectés supplémentaires et R augmente de b t. La variation du nombre d’infectés est donc égale à a S – b multiplié par I t. La condition pour que la maladie se propage (et donc donne lieu à une épidémie) est que le nombre de malades infectés augmente, c’est-à-dire que : a S – b > 0. Le quotient b / a a donc valeur de seuil. Si le nombre de sujets susceptibles est strictement inférieur à ce seuil, la maladie ne s’étend pas. Sinon, elle donne lieu à une épidémie (ou à une épizootie).

D’une façon qui peut paraître paradoxale, l’apparition d’une épidémie ne dépend donc pas du nombre de personnes infectées mais du nombre de personnes susceptibles d’attraper la maladie ! Cette remarque justifie à elle seule les politiques de vaccination, même avec un vaccin peu efficace.

Les mathématiques des épidémies

Le passé a connu des épidémies terribles comme les pestes du Moyen-Âge. Avant qu’on comprenne leur mode de transmission, les hommes étaient incapables de s’en prémunir. Quand on le comprit, on put opérer des mises en quarantaine. Enfin, la solution vint avec les vaccinations, qui permettent de réduire le nombre de gens susceptibles de contacter une maladie. Cela suffit pour éviter une épidémie, comme l’explique le modèle SIR. En l’absence de vaccin, on revient aux quarantaines mais avec un peu plus de souplesse.

Le modèle SIR

William Kermack (1898 – 1970) et Anderson Mac Kendrick (1876 – 1943).

William Kermack et Anderson Mac Kendrick ont modélisé les épidémies en 1927. Leur modèle compartimente la population en trois classes : S, la classe des individus susceptibles d’attraper la maladie, I, celle de ceux qui en sont infectés (et contagieux) et R, ceux qui en sont revenus ou morts. Dans les deux cas, ces derniers sont immunisés et ne contamineront plus personne. Le modèle SIR considère l’évolution de ces trois classes dans le temps en fonction de deux taux mesurables expérimentalement. Le premier (a) est le taux de contagion de la maladie pour un infecté, c’est-à-dire la probabilité pour qu’un individu susceptible attrape la maladie après contact avec un individu infecté. Le second taux (b) mesure le passage de l’état I à l’état R. Après un laps de temps Δt, on compte a I S Δt infectés supplémentaires et R augmente de b I Δt. La variation du nombre d’infectés est donc égale à a S – b multiplié par I Δt. La condition pour que la maladie se propage (et donc donne lieu à une épidémie) est que le nombre de malades infectés augmente, c’est-à-dire que a S – b > 0. Le quotient b / a a donc valeur de seuil. Si le nombre de sujets susceptibles est inférieur à ce seuil, la maladie ne s’étend pas. Sinon, elle donne lieu à une épidémie (ou à une épizootie).

Un effet de seuil

D’une façon qui peut paraître paradoxale, l’apparition d’une épidémie ne dépend donc pas du nombre de personnes infectées, mais du nombre de personnes susceptibles d’attraper la maladie !

La vaccination permet de passer en dessous du seuil qui permet une épidémie.

Cette remarque justifie à elle seule les politiques de vaccination, même avec un vaccin peu efficace. Pour éviter une épidémie, c’est le nombre d’individus susceptibles d’attraper la maladie qu’il faut diminuer ! Ce nombre dépend de chaque maladie. Pour la rougeole, pour éviter l’épidémie, une couverture vaccinale de 95 % est nécessaire, ce qui n’est plus assuré en France du fait de campagnes obscurantistes antivaccins. La France a donc la honte d’exporter la rougeole dans des pays comme la Suisse ou le Costa Rica.

Calcul du seuil

Le calcul du seuil est très délicat. Pour cela l’usage est d’introduire un autre coefficient, appelé le taux de reproduction de base et noté R0 (R zéro) qui est le nombre de cas secondaires produits par un individu infectieux moyen au cours de sa période d’infectiosité, dans une population entièrement constituée de personnes susceptibles d’attraper la maladie. D’après la définition même, il y a épidémie si et seulement si R0 > 1. Un beau résultat … qui ne signifie rien si on ne sait pas évaluer ce R! Ainsi, au début de l’épidémie de Covid-19, il était de 3,3 selon l’institut Pasteur et descendu à 0,5 après le confinement. Ces simples résultats montrent que la méthode ne vaut que pour un territoire uniforme et qu’une autre logique doit s’appliquer quand les cas sont rares. C’est d’ailleurs bien le cas puisqu’on essaye de confiner les foyers de contaminations isolés. En revanche, les modèles utilisés supposent une population homogène, ce qui n’est absolument pas le cas. Les personnes affaiblies par le grand âge ou des comorbidités risquent plus de contracter une forme grave de la maladie tandis que les jeunes sont souvent plus insouciants et, de ce fait, risquent davantage d’être contaminés mais avec des formes bénignes ou même asymptomatiques. Le modèle SIR est-il applicable dans le cas du Covid ? Sans doute non car il suppose un taux de contamination indépendant des conditions climatiques et homogène dans la population ce qui semble faux. Sans doute faudrait-il considérer la saisonnalité et des groupes d’âge, ou de la quantité de virus à laquelle on a été exposé, qui implique la charge virale. D’autre part, les données collectées pour calculer le R0 sont sujettes à caution. En particulier, les chiffres en provenance de Chine au début de l’épidémie laissaient prévoir une létalité très faible. Tout au plus, de 0,04 % ce qui n’est guère inquiétant. Elle fut en fait bien plus importante et très variable selon les pays, d’autant plus que les méthodes de collectes des décès différentes selon les pays biaisent les résultats. Pour conclure, les modèles sont importants d’un point de vue théorique mais peuvent être peu fiables en pratique. S’y fier aveuglément peut être dangereux, il importe de garder son bon sens et ne pas tuer une population atteinte sous prétexte de la sauver. La survie dépend aussi de l’économie ! De même, comment évaluer les traitements quand 98 % des personnes infectées guérissent sans le moindre traitement ? Le problème se pose également pour l’évaluation d’un vaccin.

Les mesures de confinement

En l’absence de vaccin, il reste qu’on peut essayer de diminuer le taux a ce qui augmente le seuil b / a. Pour cela, on peut éviter une quarantaine stricte et se contenter de réduire les interactions sociales. Dans le cas du Covid19, il semble que limiter les interactions sociales à 5 par jour suffit.

Le théorème du moustique

Le modèle SIR justifie également le théorème du moustique découvert par Ronald Ross (1857 – 1932) en 1911, et selon lequel il n’est pas besoin d’éliminer tous les moustiques pour éradiquer le paludisme, il suffit d’en faire passer la population sous un certain seuil. Cette découverte a précédé sa justification théorique au moyen du modèle SIR. Auparavant, le lien entre les marécages et le paludisme était connu, comme le montrent les tentatives de drainage des marais Pontins près de Rome de l’Antiquité jusqu’au XIXe siècle. De même, le paludisme a disparu de France, où il était autrefois endémique dans les régions humides comme la Sologne, le marais Poitevin et même autour de Port-Royal des Champs, au cours du XIXe siècle grâce au drainage et non à la consommation de Quinquina, un apéritif à base de Quinine, un préventif du paludisme, comme cela a été parfois affirmé !

Voir aussi : https://blogs.futura-sciences.com/lehning/2020/03/19/le-jeu-de-la-vie-et-celui-des-epidemies/