Tugas Week 11 Bintel

Visualisasi Trellis Graphics

Menggunakan Package Lattice

library(lattice)

library(readr)
data <- read_delim("C:/Users/nashwa/Downloads/dkikepadatankelurahan2013 (1).csv", 
    delim = ";", escape_double = FALSE, trim_ws = TRUE)

## Rows: 267 Columns: 7
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ";"
## chr (4): NAMA PROVINSI, NAMA KABUPATEN/KOTA, NAMA KECAMATAN, NAMA KELURAHAN
## dbl (3): TAHUN, LUAS WILAYAH (KM2), 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.

data

## # A tibble: 267 × 7
##    TAHUN `NAMA PROVINSI` `NAMA KABUPATEN/KOTA` `NAMA KECAMATAN` `NAMA KELURAHAN`
##    <dbl> <chr>           <chr>                 <chr>            <chr>           
##  1  2013 PROVINSI DKI J… KAB.ADM.KEP.SERIBU    KEP. SERIBU UTR  P. PANGGANG     
##  2  2013 PROVINSI DKI J… KAB.ADM.KEP.SERIBU    KEP. SERIBU UTR  P. KELAPA       
##  3  2013 PROVINSI DKI J… KAB.ADM.KEP.SERIBU    KEP. SERIBU UTR  P. HARAPAN      
##  4  2013 PROVINSI DKI J… KAB.ADM.KEP.SERIBU    KEP. SERIBU SLT  P. UNTUNG JAWA  
##  5  2013 PROVINSI DKI J… KAB.ADM.KEP.SERIBU    KEP. SERIBU SLT  P. TIDUNG       
##  6  2013 PROVINSI DKI J… KAB.ADM.KEP.SERIBU    KEP. SERIBU SLT  P. PARI         
##  7  2013 PROVINSI DKI J… JAKARTA PUSAT         GAMBIR           GAMBIR          
##  8  2013 PROVINSI DKI J… JAKARTA PUSAT         GAMBIR           CIDENG          
##  9  2013 PROVINSI DKI J… JAKARTA PUSAT         GAMBIR           PETOJO UTARA    
## 10  2013 PROVINSI DKI J… JAKARTA PUSAT         GAMBIR           PETOJO SELATAN  
## # ℹ 257 more rows
## # ℹ 2 more variables: `LUAS WILAYAH (KM2)` <dbl>, `KEPADATAN (JIWA/KM2)` <dbl>

colnames(data) = c("Tahun", "Prov", "Kab", "Kec", "Kel", "Luas", "Kepadatan")
data

## # A tibble: 267 × 7
##    Tahun Prov                 Kab                Kec       Kel    Luas Kepadatan
##    <dbl> <chr>                <chr>              <chr>     <chr> <dbl>     <dbl>
##  1  2013 PROVINSI DKI JAKARTA KAB.ADM.KEP.SERIBU KEP. SER… P. P…  0.91      6779
##  2  2013 PROVINSI DKI JAKARTA KAB.ADM.KEP.SERIBU KEP. SER… P. K…  3.76      1705
##  3  2013 PROVINSI DKI JAKARTA KAB.ADM.KEP.SERIBU KEP. SER… P. H…  3.59       628
##  4  2013 PROVINSI DKI JAKARTA KAB.ADM.KEP.SERIBU KEP. SER… P. U…  0.59      3625
##  5  2013 PROVINSI DKI JAKARTA KAB.ADM.KEP.SERIBU KEP. SER… P. T…  1.57      3084
##  6  2013 PROVINSI DKI JAKARTA KAB.ADM.KEP.SERIBU KEP. SER… P. P…  1.39      1968
##  7  2013 PROVINSI DKI JAKARTA JAKARTA PUSAT      GAMBIR    GAMB…  2.58      1350
##  8  2013 PROVINSI DKI JAKARTA JAKARTA PUSAT      GAMBIR    CIDE…  1.26     14584
##  9  2013 PROVINSI DKI JAKARTA JAKARTA PUSAT      GAMBIR    PETO…  1.12     18987
## 10  2013 PROVINSI DKI JAKARTA JAKARTA PUSAT      GAMBIR    PETO…  1.14     14465
## # ℹ 257 more rows

barchart(Kepadatan ~ Kab,
         data=data,
         col="pink",
         main="Kepadatan Penduduk tiap Kabupaten",
         xlab="Kabupaten",
         ylab="Kepadatan Penduduk")

bwplot(Kab ~ Kepadatan,
       data=data,
       col="lavender",
       main="Distribusi Kepadatan Penduduk per Kabupaten",
       xlab="Kepadatan",
       ylab="Kabupaten/Kota")

xyplot(Kepadatan ~ Luas | factor(Kab),
       data=data,
       layout=c(6,1),
       aspect=1,
       col="magenta",
       ylab="Kepadatan Penduduk",
       xlab="Luas Wilayah")

xyplot(Luas ~ Kepadatan, data=data, group=Kab, auto.key=list(space="right"))

Menggunakan Package Lattice Extra

library(latticeExtra)

## Warning: package 'latticeExtra' was built under R version 4.5.3

densityplot(~Luas | Kab, data=data, main= "Densitas Luas",
            col = "red",
            xlab = "Luas",
            ylab = "Densitas")

qqmath(~Kepadatan | Kab, data = data,
       main = "Q-Q Plot Kepadatan Penduduk",
       xlab = "Luas",
       ylab = "Densitas",
       col = "green")

marginal.plot(Kepadatan ~ Luas, data = data, groups = Kab,
              auto.key = TRUE, pch = 24, main = "Marginal Plot Kepadatan Penduduk")