Eesti geotahkdiagramming (geofacets)

Mõtlesin teel Visualising Data Hackathonile, et võiks siin Eesti geo-facet’id ära teha. Algus venis veidi, ja see osutus nii kole lihtsaks, et sain valmis enne kui üritus pihta algas (NB! ma olen graafiku välja venitanud siin, muidu ei ole ta kuidagi Eesti moodi - maakonnapiirid tuleks uuesti ümber teha, et nad paremini sellisesse gridi mahuks).

library(tidyverse)
library(geofacet)

  estgrid <- data.frame(
    name = c("Harju maakond", "Lääne-Viru maakond", "Ida-Viru maakond", "Jõgeva maakond", "Järva maakond", "Hiiu maakond", "Lääne maakond", "Rapla maakond", "Tartu maakond", "Viljandi maakond", "Saare maakond", "Pärnu maakond", "Valga maakond", "Põlva maakond", "Võru maakond"),
    code = c("Harju", "Lääne-Viru", "Ida-Viru", "Jõgeva", "Järva", "Hiiu", "Lääne", "Rapla", "Tartu", "Viljandi", "Saare", "Pärnu", "Valga", "Põlva", "Võru"),
    row = c(1, 1, 1, 2, 2, 2, 2, 2, 3, 3, 3, 3, 4, 4, 5),
    col = c(4, 5, 6, 6, 5, 1, 3, 4, 6, 5, 1, 4, 5, 6, 6),
    stringsAsFactors = FALSE
  )

tqqtus <- read_csv("http://www.ut.ee/~iseppo/tqqtus_pikk.csv")


ggplot(tqqtus, aes(x=aasta, y=töötus))+
  geom_line(color="red")+
  theme_bw()+
  facet_geo(~Maakond, grid=estgrid, label="Maakond")