Bonjour,
Pour simuler une loi binomiale tu peux déjà utiliser la fonction rbinom() sous R, par exemple si X~bino(15,0.15) tu fais :
x <- rbinom(100, 15, 0.15)
x est un échantillon de 100 simulations de la lois binomiale.
Tu peux ensuite voir sa distribution en faisant :
Tu peux aussi simuler directement la densité en faisant de même avec la fonction dbinom().
Note : pour avoir l'aide de ces fonctions tu peux faire : ?rbinom ou ?dbinom
Sinon tu peux simuler une loi binomiale à l'aide d'une loi uniforme. Avec une loi uniforme il est relativement simple de simuler une bernouilli, et une binomiale est une somme de bernouilli. Je pense que c'est cela qu'ils entendent dans ton énoncé par Monte-Carlo.
Tu trouveras une explication ici : http://math.univ-lille1.fr/~suquet/Polys/simul.pdf si tu ne t'y connais pas.
Bon courage ! Et R est relativement simple à appréhender, tu as une forte communauté en ligne à disposition. Je te conseille d'utiliser RStudio comme IDE pour à la fois éditer ton code et l'exécuter.
Partager