Zobrazení výsledků pro chlamydie
V předchozím cvičení jsi na illinoiská data o chlamydiích aplikoval/a GLMER. V tomto cvičení si ukážeme několik způsobů, jak výsledky zobrazit. Tyto metody můžeš využít k přípravě přehledů modelu pro klienta nebo dokument, ve kterém výsledky popisuješ. Doporučuji ti ale také samostatně prozkoumávat výstupy modelu a vytvářet vlastní způsoby jejich prezentace – originální přístup tě jako datového analytika může příjemně odlišit od ostatních!
Postup bude následující:
- Prozkoumej odhady modelu.
- Vykresli data a pro každou věkovou skupinu přizpůsob
glm. I když výsledky nejsou úplně totožné s výstupyglmer(), tato aproximace pomáhá výsledky přehledně vizualizovat.
Toto cvičení je součástí kurzu
Hierarchické a smíšené modely s náhodnými efekty v R
Pokyny k cvičení
- Extrahuj odhady pevných efektů z
model_outpomocífixef(). - Extrahuj odhady náhodných efektů z
model_outpomocíranef(). - Spusť kód pro vykreslení dat pomocí metod z balíčku
ggplot2.
Interaktivní cvičení na vyzkoušení si v praxi
Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.
# Extract out fixed effects
___
# Extract out random effects
___
# Run code to see one method for plotting the data
ggplot(data = il_data_2,
aes(x = year, y = count, group = county)) +
geom_line() +
facet_grid(age ~ . ) +
stat_smooth(method = "glm",
method.args = list(family = "poisson"),
se = FALSE,
alpha = 0.5) +
theme_minimal()