Voir la source de Modèle de Gompertz

Le '''modèle de Gompertz''' (ou '''loi de mortalité de Gompertz-Makeham''') établit que le [[Table de mortalité|taux de mortalité]] est la somme de termes indépendants de l'âge (termes de Makeham) et de termes dépendants de l'âge ([[Loi de Gompertz|fonction de Gompertz]]).

Ce modèle suggère également la diminution exponentielle du nombre d'organismes vivants proportionnellement à l'augmentation linéaire de l'âge.

== Histoire ==

La mathématisation de la science de la population progresse au {{s-|XIX|e}}, notamment grâce à la loi de la mortalité établie par [[Benjamin Gompertz]], mais également grâce à la [[loi logistique]] de [[Pierre François Verhulst]], selon laquelle la croissance de la population ralentit.

En effet en 1825, [[Benjamin Gompertz]] propose que la [[force de mortalité]] augmente de façon exponentielle avec l'âge :

:<math> \mu_x = Bc^{x} </math>
où {{mvar|B}} et {{mvar|c}} sont des constantes.

Cependant, se pose le problème de prendre en compte les causes de décès qui ne seraient pas directement liées à l'âge.

C'est ainsi que William Makeham propose d'extrapoler la [[force de mortalité]] aux grands âges sur la base d'une loi de Gompertz modifiée et tenant compte des causes de décès indépendantes de l'âge :
:<math> \mu_x = A + Bc^{x} </math>
où {{mvar|A}} est le risque de mourir pour l'ensemble des causes indépendantes de l'âge.

Dans des conditions où les causes externes de décès sont rares, les termes qui ne dépendent pas de l'âge sont souvent négligeables. On parle alors de la loi de Gompertz, due à Benjamin Gompertz en 1825.

Le modèle est largement utilisé en démographie et [[gérontologie]] pour des prévisions adéquates du taux de mortalité chez certaines espèces (non humaines) et pour comparer les taux de vieillissement actuariels entre et parmi différentes espèces.

Il représente bien la croissance de certaines variables morphologiques (taille, masse corporelle…) d'organismes supérieurs (voir [[#Modélisation de la croissance pour le rat musqué|exemple d'application sur la croissance en masse corporelle du rat musqué en fonction de son âge]]).

== Modèle de Gompertz et dynamique des populations ==

La loi de Gompertz-Makeham décrit la dynamique de la mortalité, qui appartient à la [[dynamique des populations]].

Le modèle a surtout été utilisé pour représenter la croissance de certains organismes.

Il permet de rendre compte de la relation d'[[allométrie]] entre deux variables, en plus de la bonne représentation qu'il offre pour une variable.

Lorsque l'on compare le modèle de Gompertz au [[modèle de Verhulst]], on observe un comportement similaire ([[croissance exponentielle]] de la population) néanmoins dans le cas du modèle de Gompertz, elle est plus rapide.
Ces deux modèles sont concurrents pour la modélisation de la croissance des organismes.

=== Le modèle mathématique de Gompertz ===
Le modèle de Gompertz permet de modéliser la croissance d'une population régulée.

On l'exprime ainsi sous forme d'équation différentielle :
<center><math>\frac{\mathrm{d}x}{\mathrm{d}t}=a x \ln\left(\frac{K}{x}\right)</math></center>

Il est également exprimé sous sa forme intégrée :
<center> <math> x = K \mathrm{e}^{b \mathrm{e}^{-at}} , \ \textrm{ avec } \ b = \ln\left(\frac{x_0}{K}\right) </math> </center>

La forme intégrée du modèle de Gompertz est souvent utilisée pour le calcul numérique, tandis que la [[forme différentielle]] se prête mieux à l'interprétation.

Dans les équations, les paramètres :
* {{mvar|x}} représente la biomasse (taille, masse corporelle…) ;
* {{mvar|t}} le temps ;
* {{mvar|K}} la capacité limite du milieu ;
* {{mvar|a}} une constante.

=== Point d'équilibre du modèle ===

Un état d'équilibre de la population est observé quand la population n'évolue pas. Les points d'équilibre sont les valeurs {{math|''x''{{exp|*}}}} pour lesquelles <math> \dot{x} = \tfrac{\mathrm{d}x}{\mathrm{d}t} = 0 </math>. On trouve deux points d'équilibre : {{math|''x''{{indexp|1|*}} {{=}} 0}} et {{math|''x''{{indexp|2|*}} {{=}} ''K''}}.

{{Boîte déroulante|titre= Calcul des points d'équilibre|contenu= 

Un point d'équilibre est atteint quand <math> \dot{x} = \tfrac{\mathrm{d}x}{\mathrm{d}t} = 0 </math>, soit :
:<math>\frac{\mathrm{d}x}{\mathrm{d}t}= 0 \longleftrightarrow a x \ln\left(\frac{K}{x}\right)=0 </math>

Deux solutions possibles : soit {{math|''ax'' {{=}} 0}}, soit {{math|ln({{sfrac|''K''|''x''}}) {{=}} 0}}.

* '''Pour le premier point''' :

:<math> ax = 0  
\text{, }
a\ne 0 \Longleftrightarrow x=0 </math>, on notera ce point d'équilibre {{math|''x''{{indexp|1|*}}}}.
* '''Pour le deuxième point : '''

:<math>\ln\left(\frac{K}{x}\right)=0 \Longleftrightarrow  \frac{K}{x}=1 \Longleftrightarrow x = K </math>,

on notera ce point d'équilibre {{math|''x''{{indexp|2|*}}}}.}}

=== Stabilité locale ===
On étudie la stabilité au point d'équilibre  {{math|''x''{{indexp|1|*}}}} et {{math|''x''{{indexp|2|*}}}}.
Ce qui veut dire que l'on va déterminer pour tout {{mvar|x}} proche de {{math|''x''{{exp|*}}}} , si l'on '''se rapproche''' ou si l'on '''s'éloigne''' de {{math|''x''{{exp|*}}}}.

On observe le signe de {{sfrac|d''f''|d''x''}} :
* si {{math|{{sfrac|d''f''|d''x''}} > 0}} alors {{math|''x''{{exp|*}}}} est instable ;
* si {{math|{{sfrac|d''f''|d''x''}} < 0}} alors {{math|''x''{{exp|*}}}} est stable.

Pour cela, on calcule la [[dérivée partielle]] en ces points :
<center> <math>\frac{\mathrm{d}f}{\mathrm{d}x}=a \times \ln(K)-a \times \ln(x)-a</math>. </center>

Il apparait que :
* {{math|''x''{{indexp|1|*}}}} est instable : pour des {{mvar|x}} proches de {{math|''x''{{indexp|1|*}}}}, on '''s'éloigne''' du point d'équilibre {{math|''x''{{indexp|1|*}}}}.
* {{math|''x''{{indexp|2|*}}}} est stable : pour des {{mvar|x}} proches de {{math|''x''{{indexp|2|*}}}}, on '''se rapproche''' du point d'équilibre {{math|''x''{{indexp|2|*}}}}.

{{Démonstration| contenu = 
Au point d'équilibre {{math|''x''{{indexp|1|*}} {{=}} 0}}, on a {{sfrac|d''f''|d''x''}} non défini (car {{math|ln(0)}} n’existe pas) mais <math>\lim_{x \to 0} \frac{\mathrm{d}f}{\mathrm{d}x}= +\infty >0 \Longrightarrow x_1^*</math> est instable.

Au point d'équilibre {{math|''x''{{indexp|2|*}} {{=}} ''K''}}, on a {{math|{{sfrac|d''f''|d''x''}} {{=}} –''a'' < 0}}, donc {{math|''x''{{indexp|2|*}}}} est stable.}}

=== Modélisation de la croissance pour le rat musqué ===

Le modèle de Gompertz représente bien la croissance en [[masse corporelle]] du [[rat musqué]] ou sa taille en fonction de son âge{{Référence souhaitée|date=6 août 2023}}.

''Ajustement du modèle de Gompertz aux données sur la masse corporelle avec <math> x_0=16,\text{ }a=0,036 \text{ et } K=760 </math>.''

L'équation s'écrit donc <math> x(t)=760\times \exp\left(\ln\left(\frac{16}{760}\right) \mathrm{e}^{-0,036t}\right)</math>

[[Fichier:Exemple modele gompertz.jpg|center|thumb|Exemple de représentation de données et ajustement au modèle de Gompertz.]]


=== Modélisation de la croissance des tumeurs ===

La courbe de Gompertz est utilisée pour ajuster les données de la croissance des tumeurs. En fait, les tumeurs sont des populations de cellules dans un espace confiné où la disponibilité des éléments nutritifs est limitée. En notant la taille de la tumeur {{math|''X''(''t'')}}, il est utile d'écrire la courbe de Gompertz comme suit :
:<math>X(t) = K \times \mathrm{e}^{(b \times \mathrm{e}^{-at})}, \ \textrm{ avec } \ b = \ln\left(\frac{x_0}{K}\right) </math>
où {{math|''X''(0)}} est la taille de la [[tumeur]] au moment de l'observation de départ, et {{mvar|K}} est la taille maximale qui peut être atteinte avec les éléments nutritifs disponibles, c’est-à-dire la population asymptotique de cellules tumorales.

== Quelques exemples d'applications ==
* Applications en [[microbiologie]] et alimentation pour décrire les phases de croissance microbienne (voir : {{ISBN|978-2-84444-558-2}})
* Applications dans des modèles de prédictions de croissance

== Voir également ==

=== Liens internes ===
* [[Allométrie]]
* [[Croissance exponentielle]]
* [[Dynamique des populations]]
* [[Modèle de Verhulst]]

=== Liens externes ===
* {{Lien web|url= http://hermandenis.perso.sfr.fr/219.%20Ann.%20XIX.%20Mesure%20de%20la%20longevite.pdf <!-- Paramètre obligatoire --> |titre= MESURE DE LA LONGÉVITÉ DES ANIMAUX ET DES ÊTRES HUMAINS<!-- Paramètre obligatoire --> |auteur= Herman Denis}}
* {{Bases}}

=== Bibliographie ===
* {{Ouvrage|langue=fr|auteur1=Alain Pavé|titre=Modélisation des systèmes vivants|sous-titre=de la cellule à l'écosystème|lieu=Paris|éditeur=[[éditions Lavoisier]]|date=octobre 2012|pages totales=634|isbn=978-2-7462-3911-1}}
* {{Ouvrage|langue=en|prénom1=James Dickson|nom1=Murray|titre=Mathematical Biology<!-- Paramètre obligatoire -->|sous-titre=Spatial Models and Biomedical Applications|tome=II|lieu=New York etc.|éditeur=Springer|année=2011|numéro d'édition=3|pages totales=811|isbn=978-0-387-95228-4|isbn2=0-387-95228-4|lire en ligne=https://books.google.com/books?id=Dk8pFkP2cKYC&printsec=frontcover}}
* {{Ouvrage|langue=fr|prénom1=Graziella|nom1=Caselli|titre=Démographie : Analyse et Synthèse <!-- Paramètre obligatoire -->|sous-titre=Population et société|volume=VI|éditeur=Institut National d’Études Démographiques|année=|pages totales=585|isbn=}}
* {{Ouvrage|langue=fr|prénom1=Graziella|nom1=Caselli|titre=Démographie : Analyse et Synthèse <!-- Paramètre obligatoire -->|sous-titre=Les Déterminants de la mortalité|volume=III|lieu=Paris|éditeur=Institut national d’études démographiques|année=2002|pages totales=483|isbn=2-7332-2013-6}}
* {{Ouvrage|langue=fr|prénom1=Alain|nom1=Branger|titre=Microbiochimie et Alimentation <!-- Paramètre obligatoire -->|lieu=Dijon|éditeur=Educagri|année=2012|pages totales=344|passage=113-114|isbn=978-2-84444-558-2|lire en ligne=https://books.google.com/books?id=tvxQg3cEJyoC&printsec=frontcover}}

{{Portail|statistiques|Démographie}}

{{DEFAULTSORT:Modele de Gompertz}}
[[Catégorie:Statistiques]]
[[Catégorie:Démographie]]
[[Catégorie:Sciences actuarielles]]