Visualisering av ett univariat gaussiskt blandningsmodell
Nu när du anpassat modellen till fit_mix_example och extraherat parametrarna till comp_1, comp_2 och comp_3 (samt proportions), är det dags att plotta motsvarande kluster tillsammans med densitetshistogrammet.
För att underlätta detta har funktionen fun_prop() definierats i miljön. Den här funktionen returnerar densitetsvärden för en gaussisk fördelning, liknande dnorm, men är utökad till att även ta emot proportioner.
Den här övningen är en del av kursen
Blandningsmodeller i R
Övningsinstruktioner
- Rita densitetshistogrammet tillsammans med densiteten för varje kluster. Kom ihåg att dataramen heter
mix_example. - Använd funktionen
stat_function()med argumentetfunsatt tillfun_propför att rita densitetsfördelningen för varje kluster.
Interaktiv övning med praktiskt arbete
Testa den här övningen genom att slutföra den här exempelkoden.
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]))