Simuler un mélange de deux distributions gaussiennes
Les modèles de mélange peuvent être difficiles à saisir. Commençons donc par en simuler un simple, puis nous irons plus loin. Dans cet exercice, vous allez créer un modèle de mélange gaussien à deux composantes.
Cette activité fait partie du cours
Modèles de mélange en R
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 decoin, de la façon suivante :- Lorsque la valeur de
coinest 1, échantillonnez une loi normale de moyenne 5 et d'écart-type 2. - Sinon, échantillonnez une loi normale standard.
- Lorsque la valeur 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)