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
Oefeninstructies
- Plot het dichtheidshistogram samen met de dichtheid van elke cluster. Onthoud dat de data frame
mix_exampleheet. - Gebruik de functie
stat_function()met het argumentfungelijk aanfun_propom 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]))