LoslegenKostenlos starten

Univariates Gaußsches Mischmodell visualisieren

Da du das Modell in fit_mix_example geschätzt und die Parameter in comp_1, comp_2 und comp_3 (sowie die proportions) extrahiert hast, lass uns nun die zugehörigen Cluster zusammen mit dem Dichtehistogramm plotten.

Um dir das zu erleichtern, wurde die Funktion fun_prop() in der Umgebung definiert. Diese Funktion liefert die Dichtewerte für eine gaußsche Verteilung, ähnlich wie dnorm, ist aber erweitert, sodass sie auch die Anteile akzeptiert.

Diese Übung ist Teil des Kurses

<Kurs>Mischungsmodelle in R</Kurs>
Kurs ansehen

Übungsanweisungen

  • Zeichne das Dichtehistogramm zusammen mit der Dichte jedes Clusters. Denk daran, dass der Data Frame mix_example heißt.
  • Verwende die Funktion stat_function() mit dem Argument fun gleich fun_prop, um die Dichteverteilung für jeden Cluster zu zeichnen.

Interaktive praktische Übung

Versuche dich an dieser Übung, indem du diesen Beispielcode vervollständigst.

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 bearbeiten und ausführen