De resultaten weergeven
De laatste, en misschien wel belangrijkste stap bij het bouwen van een model, is het delen van de resultaten.
In deze oefening haal je de schattingen op county-niveau op en visualiseer je ze met ggplot2. De random-effect hellingen op county-niveau moeten worden opgeteld bij de fixed-effect hellingen om de hellingschatting per county te krijgen.
Naast deze optelling sorteert de code de counties op misdaadcijfer (de hellingschattingen) om de visualisatie duidelijker te maken.
Het model dat je eerder hebt gefit, glmer_out, is voor je ingeladen.
Deze oefening maakt deel uit van de cursus
Hiërarchische en Mixed-Effects-modellen in R
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
# 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)])