Simuler un mélange de deux distributions gaussiennes
Les modèles de mélange peuvent être difficiles à appréhender. Commençons donc par simuler un modèle de mélange simple, puis nous irons plus loin. Dans cet exercice, vous allez créer un modèle de mélange gaussien à deux composantes.
Cet exercice fait partie du cours
<cours>Modèles de mélange en R</cours>Instructions de l’exercice
- Simulez 500 observations d’une variable binaire avec
sample()en utilisant une probabilité de 0,8 pour la valeur 1, et enregistrez le résultat dans l’objetcoin. - Créez ensuite un vecteur
mixturequi échantillonne deux distributions gaussiennes différentes selon les valeurs decoinde la manière suivante :- Lorsque les valeurs de
coinvalent 1, échantillonnez une loi normale de moyenne 5 et d’écart-type 2. - Sinon, échantillonnez une loi normale standard.
- Lorsque les valeurs de
Exercice interactif pratique
Essayez cet exercice en complétant ce code d’exemple.
# Create coin object
coin <- ___(c(0,___), size = ___,
replace = TRUE, prob = c(0.2, ___))
# Sample from two different Gaussian distributions
___ <- ifelse(___ == 1, rnorm(n = ___, mean = ___, sd = ___),
rnorm(n = ___))
# Check the first elements
head(mixture)