Wyświetlanie wyników
Ostatnim – i zdecydowanie najważniejszym – etapem tworzenia modelu jest prezentacja wyników.
W tym ćwiczeniu wyodrębnisz estymaty na poziomie hrabstw i zwizualizujesz je za pomocą ggplot2. Aby uzyskać estymaty nachylenia dla każdego hrabstwa, należy dodać losowe efekty nachylenia na poziomie hrabstw do efektów stałych.
Poza tym dodawaniem kod zawiera również sortowanie hrabstw według współczynnika przestępczości (estymaty nachylenia), co ułatwia przejrzyste przedstawienie danych.
Model glmer_out, który wcześniej dopasowałeś(-aś), został już wczytany.
To ćwiczenie jest częścią kursu
Hierarchiczne i mieszane modele efektów w R
Interaktywne ćwiczenie praktyczne
Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.
# 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)])