Miscela di tre distribuzioni gaussiane
Cosa cambia se aggiungiamo un’altra distribuzione alla nostra simulazione? Vedrai che aumentando il numero di componenti la densità di massa si distribuisce per includere l’extra distribuzione, ma la logica resta la stessa dell’esercizio precedente.
Questo esercizio fa parte del corso
Modelli di Mixture in R
Istruzioni dell'esercizio
- Crea
assignments, che assume i valori 0, 1 e 2 con probabilità rispettivamente 0,3, 0,4 e 0,3. - Il data frame
mixturecampiona da una gaussiana conmeanpari a 5 esdpari a 2 quandoassignmentsè 1. Seassignmentsè 2,meanè 10 esdè 1. Altrimenti, è una normale standard. - Traccia l’istogramma con 50 bin.
Esercizio pratico interattivo
Prova a risolvere questo esercizio completando il codice di esempio.
number_observations <- 1000
# Create the assignment object
assignments <- sample(
c(0,1,2), size = number_observations, replace = TRUE, prob = c(0.3, ___, 0.3)
)
# Simulate the GMM with 3 distributions
mixture <- data.frame(
x = ifelse(___ == 1, rnorm(n = number_observations, mean = ___, sd = ___), ifelse(assignments == 2, rnorm(n = number_observations, mean = ___, sd = ___), rnorm(n = ___)))
)
# Plot the mixture
mixture %>%
ggplot() + ___(aes(x = x, y = ..density..), ___ = ___)