1. Učit se
    2. /
  2. Kurzy
    3. /
  3. Vizualizace geoprostorových dat v R
Connected

cvičení

Příklad intervalové škály

Vrátíme se k tvému grafu podílu obyvatelstva ve věku 18 až 24 let:

tm_shape(prop_by_age) +
  tm_raster("age_18_24", palette = vir) +
  tm_legend(position = c("right", "bottom"))

Graf byl problematický, protože většina podílů spadala do nejnižší barevné úrovně, a v důsledku toho nebylo v grafu příliš vidět. Jedním ze způsobů, jak tento problém vyřešit, je rozdělit rozsah proměnné nikoli do intervalů stejné délky, ale do smysluplnějších kategorií.

Začneme replikací výchozích přihrádek tmap: pět kategorií rozdělených pomocí "pretty" breaks. Pak můžeš vyzkoušet několik dalších metod dělení proměnné do intervalů. Přímé použití funkce classIntervals() ti rychle ukáže, kde budou hranice, ale nejlepší způsob, jak sadu hranic otestovat, je vytvořit graf.

(Mimochodem, dalším způsobem, jak takový problém řešit, je najít transformaci proměnné tak, aby byly intervaly stejné délky na transformované škále užitečnější.)

Pokyny

100 XP
  • Zavolej classIntervals() na values(prop_by_age[["age_18_24"]]) s argumenty n = 5 a style = "pretty". Vidíš problém? 130 770 buněk mřížky skončí v první přihrádce.
  • Teď zavolej classIntervals() stejně jako výše, ale s style = "quantile".
  • Použij stejně velké přihrádky tak, že předáš argumenty n a style do vrstvy tm_raster() ve svém grafu.
  • Vytvoř histogram z values(prop_by_age[["age_18_24"]]). Kde bys hranice umístil/a?
  • Definuj vlastní hranice v tm_raster() zadáním breaks = c(0.025, 0.05, 0.1, 0.2, 0.25, 0.3, 1).
  • Ulož finální graf jako leaflet pomocí tmap_save() s filename "prop_18-24.html".