Aan de slagGa gratis aan de slag

Een univariate Gaussian Mixture Model visualiseren

Je hebt het model gefit in fit_mix_example en de parameters geëxtraheerd in comp_1, comp_2 en comp_3 (en ook de proportions). Laten we nu de bijbehorende clusters samen met het dichtheidshistogram plotten.

Om dit makkelijker te maken is de functie fun_prop() in de omgeving gedefinieerd. Deze functie geeft de dichtheidswaarden voor een Gaussiaanse verdeling, zoals dnorm, maar is uitgebreid zodat ook de verhoudingen worden geaccepteerd.

Deze oefening maakt deel uit van de cursus

Mixture Models in R

Cursus bekijken

Oefeninstructies

  • Plot het dichtheidshistogram samen met de dichtheid van elke cluster. Onthoud dat de data frame mix_example heet.
  • Gebruik de functie stat_function() met het argument fun gelijk aan fun_prop om de dichtheidsverdeling voor elke cluster te tekenen.

Praktische interactieve oefening

Probeer deze oefening eens door deze voorbeeldcode in te vullen.

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]))
Code bewerken en uitvoeren