library(lattice)
library(latticeExtra)
library(readr)
data = read_csv2("C:/Users/Tika/Downloads/dkikepadatankelurahan2013.csv")
## ℹ Using "','" as decimal and "'.'" as grouping mark. Use `read_delim()` for more control.
## Rows: 267 Columns: 7
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ";"
## chr (5): NAMA PROVINSI, NAMA KABUPATEN/KOTA, NAMA KECAMATAN, NAMA KELURAHAN,...
## dbl (2): TAHUN, KEPADATAN (JIWA/KM2)
## 
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
colnames(data) = c("Tahun",
                    "Prov",
                    "Kab",
                    "Kec",
                    "Kel",
                    "Luas",
                    "Kepadatan")

xyplot(Kepadatan ~ Luas | Kab,
       data = data,
       layout = c(3,2),
       xlab = "Luas Wilayah",
       ylab = "Kepadatan",
       main = "Kepadatan vs Luas Wilayah")

xyplot(Kepadatan ~ Luas | Kab,
       data = data,
       layout = c(3,2),
       xlab = "Luas Wilayah (KM2)",
       ylab = "Kepadatan (Jiwa/KM2)",
       main = "Kepadatan vs Luas Wilayah") +
  
  layer(panel.abline(h = 20000, lty = "dashed")) +
  
  layer(panel.text(max(x), 21000,
                   "kepadatan tinggi",
                   adj = c(1,0),
                   cex = 0.7))