Visa resultaten
Det sista – och utan tvekan viktigaste – steget när du skapar en modell är att dela resultaten.
I den här övningen extraherar du skattningarna på länsnivå och visualiserar dem med ggplot2. De slumpmässiga lutningseffekterna på länsnivå behöver adderas till de fasta lutningseffekterna för att få lutningsskattningarna för varje län.
Förutom den här additionen innehåller koden även en sortering av länen efter brottsnivå (lutningsskattningarna) för att göra det lättare att tolka visualiseringen.
Modellen du tidigare anpassade, glmer_out, har laddats in åt dig.
Den här övningen är en del av kursen
Hierarkiska modeller och modeller med blandade effekter i R
Interaktiv övning med praktiskt arbete
Testa den här övningen genom att slutföra den här exempelkoden.
# 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)])