Zobrazení výsledků
Posledním – a pravděpodobně nejdůležitějším – krokem při tvorbě modelu je sdílení výsledků.
V tomto cvičení extrahuješ odhady na úrovni okresů a vykreslíš je pomocí ggplot2. Aby ses dostal/a k odhadům sklonu pro každý okres, je potřeba přičíst náhodné efekty sklonu k fixním efektům sklonu.
Kód kromě tohoto součtu také řadí okresy podle míry kriminality (odhadů sklonu), aby bylo snazší data přehledně vizualizovat.
Model glmer_out, který jsi dříve natrénoval/a, je již načten.
Toto cvičení je součástí kurzu
Hierarchické a smíšené modely s náhodnými efekty v R
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 the fixed-effect slope for Year2
Year2_slope <- fixef(___)['Year2']
# Extract out the random-effect slopes for county
county_slope <- ranef(___)$County
# Create a new column for the slope
county_slope$slope <- county_slope$Year2 + Year2_slope
# Use the row names to create a county name column
county_slope$county <- rownames(county_slope)
# Create an ordered county-level factor based upon slope values
county_slope$county_plot <- factor(county_slope$county,
levels = county_slope$county[order(county_slope$slope)])