Visualizando o Modelo de Mistura Gaussiana Univariado
Como você ajustou o modelo em fit_mix_example e extraiu os parâmetros em comp_1, comp_2 e comp_3 (assim como as proportions), vamos agora plotar os clusters correspondentes com o histograma de densidade.
Para facilitar esse objetivo, a função fun_prop() foi definida no ambiente. Essa função fornece os valores de densidade para uma distribuição Gaussiana, como dnorm, mas foi estendida para também aceitar as proporções.
Este exercício faz parte do curso
Modelos de Mistura em R
Instruções do exercício
- Plote o histograma de densidade junto com a densidade de cada cluster. Lembre-se de que o data frame se chama
mix_example. - Use a função
stat_function()com o argumentofunigual afun_proppara desenhar a distribuição de densidade de cada cluster.
Exercício interativo prático
Experimente este exercício completando este código de exemplo.
ggplot(___) + ___(aes(x = x, y = ..density..)) +
stat_function(geom = "line", fun = fun_prop,
args = list(mean = ___[1], sd = ___[2],
proportion = proportions[1])) +
stat_function(geom = "line", fun = fun_prop,
args = list(mean = comp_2[1], sd = comp_2[2],
proportion = ___[2]))+
stat_function(geom = "line", fun = ___,
args = list(mean = comp_3[1], sd = comp_3[2],
proportion = proportions[3]))