Mistura de três distribuições Gaussianas
O que muda se incluirmos outra distribuição na nossa simulação? Você vai ver que aumentar o número de componentes espalha a densidade de massa para incluir a distribuição extra, mas a lógica continua a mesma do exercício anterior.
Este exercício faz parte do curso
Modelos de Mistura em R
Instruções do exercício
- Crie
assignments, que assume os valores 0, 1 e 2 com probabilidades de 0, 0.4 e 0.3, respectivamente. - O data frame
mixtureamostra de uma Gaussiana commean5 esd2 quandoassignmentsé 1. Seassignmentsé 2, entãomeané 10 esdé 1. Caso contrário, é uma distribuição normal padrão. - Plote o histograma com 50 bins.
Exercício interativo prático
Experimente este exercício completando este código de exemplo.
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..), ___ = ___)