Learning objectives:
Explore sp/sf data abstraction
Practice spatial modelling for aerial data (SIR and BYM)
Display result in tables, plots (ggplot2) and interactive map (leaflet)