Tél : (+237) 657 07 98 07
(+ 237) 651 36 96 96
B.P : 749 Yaoundé, Cameroun
Email : info@[Link]
[Link]
Master 1
MAD4037 : Probabilités
Projet de Travail Pratique :
Responsable : NYOBE Samuel
Année académique 2021-2022; Semestre I
L’objectif de ce projet est d’évaluer les compétences de chaque étudiant à écrire des codes de
simulation de variables aléatoires quelconque et leur intégration dans des modèles divers pour
analyse. Il a été vu plusieurs méthodes de simulation de variables aléatoires parmi lesquelles le
changement de variable et la méthode d’inverse de la fonction de répartition.
On fait l’étude sur la modélisation de la croissance en diamètre (en cm) d’un arbre au sein
d’une forêt. Xk (k correspond à une année) est le diamètre (en cm) à l’instant k, la suite
X0 , X1 , · · · , Xk (encore notée (Xk )k≥0 ) est un processus (de Markov) qui modélise l’évolution
annuelle du diamètre d’un arbre où l’état initial X0 suit une loi normale X0 ∼ N (10., 1.2 ) car
les arbres sont généralement inventoriés à partir de 10 cm. Le diamètre à l’instant k dépend
de celui à l’instant k − 1, ainsi, la loi de transition de Xk−1 vers Xk , ∀ k ≥ 1 est inspirée d’un
modèle de croissance celui de Gompertz [Gom25] et peut s’écrire comme suit :
−r
exp −r
Xk = D1−exp Xk−1 k−1 pour k ≥ 1,
où r et D représentent respectivement le taux de croissance et le diamètre maximal d’un arbre.
Dans la suite, prendre r = 0.3 et D = 150 cm. k est le bruit du modèle et est modélisé par la
loi de densité :
(
sin(x) si 0 < x < π2
f (x) =
0 sinon
Consignes :
Vous devez produire des réalisations de ce modèle sur 50 années soit sur 50 pas de temps mis à
part X0 L’année initiale est prise en 2022.
1. En utilisant la méthode de l’inverse de la fonction de répartition, écrire un simulateur de la
loi du bruit k décrite plus haut. En plus de ce simulateur, fournir un code pour le calcul
de la densité, de la probabilité et des quantiles de cette loi. Votre rapport doit inclure les
étapes d’obtention de l’inverse de la fonction de répartition.
2. Testez vos fonctions et représenter des histogrammes d’un jeu tiré de vos simulations et
comparez à la fonction de densité donnée plus haut pour vous assurer du bon fonction-
nement de votre simulateur. Bien vouloir simuler un grand nombre de réalisations.
3. Utilisez ensuite le simulateur du bruit pour simuler des réalisations du modèle de croissance
d’un arbre décrit plus haut. Chaque étudiant doit impérativement re-initialiser le
générateur (seed ) avec son numéro sur la liste officielle des étudiants.
1
4. D’après votre simulation sur les 50 années, en quelle année, l’abre atteindra la plus grande
valeur? Quelle est cette valeur?
5. Simuler enfin 10 réalisations de ce modèle sur les 50 années et les représenter toutes sur
le même graphique. Dans des conditions normales, un arbre doit grandir chaque année
même si cette croissance être très ralentie, d’après vous, ce modèle reflète t-il cette réalité
attendue? Sinon, est-il acceptable néanmoins?
NB: Tous les devoirs doivent remis au délégué et chaque étudiant doit inclure dans un dossier
portant son nom en entier le rapport et les codes R associés à son travail. Bien devoir commentez
vos codes. Ensuite le délégué devra compiler l’ensemble des devoirs dans un fichier compressé
et me l’envoyer. Ceux qui n’auront pas remis leurs devoirs se verront attribuer la note de 00/30.
Délai : 26 Février 2022 à 15h délai de rigueur et à envoyer (uniquement le délégué) à
l’adresse samuellenyobe@[Link].
References
[Gom25] Benjamin Gompertz. On the nature of the function expressive of the law of human
mortality, and on a new mode of determining the value of life contingencies. Philo-
sophical transactions of the Royal Society of London, 115:513–583, 1825.