Hartă coropletă

Acum că înțelegi cum se desenează poligoanele, hai să le afișăm pe o hartă. Ține minte: înlocuiești apelul ggplot() cu ggmap(), iar apelul original ggplot() se mută în argumentul base_layer(), după care adaugi stratul de poligoane ca de obicei:

ggmap(corvallis_map_bw,
      base_layer = ggplot(ward_sales,
                          aes(lon, lat))) +
  geom_polygon(aes(group = group, fill = ward))

Încearcă în consolă acum!

Hmm, ceva nu arată bine. Sectoarele 1, 3 și 8 par distorsionate și incorecte. Ce s-a întâmplat? O parte din limitele sectoarelor depășesc granițele hărții. Din cauza setărilor implicite ale ggmap(), datele care ies în afara hărții sunt eliminate înainte de afișare, astfel că unele limite ale poligoanelor dispar, iar punctele rămase se unesc greșit, rezultând forme incorecte.

Nu-ți face griji, există o soluție: ggmap() oferă câteva argumente pentru a controla acest comportament. Argumentele extent = "normal" și maprange = FALSE forțează graficul să folosească intervalul datelor, nu intervalul hărții, pentru a defini limitele de afișare.