+ scale_fill_viridis(name = "",
+ limits = c(min(sa3$incomeperearningcapita, na.rm = TRUE),max(sa3$incomeperearningcapita, na.rm = TRUE))) +
+ scale_color_viridis(name = "",
+ limits = c(min(sa3$incomeperearningcapita, na.rm = TRUE),max(sa3$incomeperearningcapita, na.rm = TRUE))) +
+ coord_sf(datum = NA) + # Work around https://github.com/tidyverse/ggplot2/issues/2071 to remove gridlines
+ labs(title = paste0(x," \nincome distribution"),
+ subtitle = "2014/15, in 1000s AUD",
+ caption = "\nSource: Australian Taxation Office") +
+ theme_economist() +
+ theme(legend.position = "bottom",
+ legend.text = element_text(angle = 45, hjust = 1, size = 8),
+ axis.text = element_blank(),
+ axis.ticks = element_blank())
+ print(plot)
+})
+}
+
+saveGIF(plots(),movie.name = "AUCitiesIncomeDistribution.gif", interval = 2, loop = 2)