Peta per Kelurahan

Package

library(tidyverse) #ggplot2
library(sf)
library(viridis)
library(patchwork)

Data kelurahan

shp_smg <- st_read("datasets/KELURAHAN_AR/KELURAHAN_AR.shp")

data_smg_3 <- read_csv("datasets/data_smg_3.csv")

data_smg_3_peta <- left_join(shp_smg,
                             data_smg_3)

Peta kelurahan

data_smg_3_peta |> 
  ggplot(aes(fill = stunting_22)) +
  geom_sf() +
  scale_fill_viridis(direction = -1)

ggsave("ggsave/peta_kelurahan_1.png", width = 7, height = 7)

data_smg_3_peta |> 
  ggplot(aes(fill = prev_stu_22)) +
  geom_sf() +
  scale_fill_viridis(direction = -1) +
  geom_sf_text(aes(label = WADMKD),
               size = 1)

ggsave("ggsave/peta_kelurahan_2.png", width = 7, height = 7)