Simuleer een mengsel van twee Gaussische verdelingen
Mixture models kunnen lastig zijn om te begrijpen. Daarom beginnen we met het simuleren van een eenvoudig mengselmodel en bouwen we vandaaruit verder. In deze oefening maak je een Gaussian Mixture Model met twee componenten.
Deze oefening maakt deel uit van de cursus
Mixture Models in R
Oefeninstructies
- Simuleer 500 waarden van een binaire variabele met
sample()met een kans van 0,8 op de waarde 1 en sla het resultaat op in het objectcoin. - Maak daarna een vector
mixturedie, afhankelijk van de waarden incoin, uit twee verschillende Gaussische verdelingen sampelt:- Als de waarden in
coingelijk zijn aan 1, sample dan uit een normale verdeling met een gemiddelde van 5 en sd van 2. - Anders sample je uit een standaardnormale verdeling.
- Als de waarden in
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
# 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)