unzip(zipfile = "data/1270055001_sa3_2016_aust_shape.zip", exdir = "data/")
sa3 <- st_read(dsn = "data/", layer = "SA3_2016_AUST", stringsAsFactors = FALSE)
-taxstats.sa3 <- merge(x = taxstats, y = sa3, by.x = "Postcode", by.y = "POA_CODE16", all.y = TRUE)
+# This doesn't work
+#taxstats.sa3 <- merge(x = taxstats, y = sa3, by.x = "Postcode", by.y = "POA_CODE16", all.y = TRUE)
# Create a matrix of intersecting postal codes and SA3's