ComeçarComece de graça

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

Ver curso

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 mixture amostra de uma Gaussiana com mean 5 e sd 2 quando assignments é 1. Se assignments é 2, então mean é 10 e sd é 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..), ___ = ___)
Editar e executar o código