Vizualizace jednorozměrného Gaussova směsného modelu
Teď, když jsi model nafitoval/a do fit_mix_example a extrahoval/a parametry do comp_1, comp_2 a comp_3 (spolu s proportions), je čas vykreslit odpovídající clustery s hustotovým histogramem.
Aby to šlo snáz, je v prostředí definovaná funkce fun_prop(). Tato funkce vrací hodnoty hustoty pro Gaussovo rozdělení podobně jako dnorm, ale navíc přijímá i proporce.
Toto cvičení je součástí kurzu
Mixture Models v R
Pokyny k cvičení
- Vykresli hustotový histogram spolu s hustotou každého clusteru. Nezapomeň, že datový rámec se jmenuje
mix_example. - Použij funkci
stat_function()s argumentemfunnastaveným nafun_proppro zobrazení hustotního rozdělení každého clusteru.
Interaktivní cvičení na vyzkoušení si v praxi
Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.
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]))