Social Vulnerability Analysis to COVID-19 in DKI Jakarta Using Factor and Cluster Analysis

Sun, 26 Feb 2023 00:00:00 +0000

Materi ini diproduksi oleh tim dari Algoritma untuk Analisis Kerentanan Sosial terhadap COVID-19 di DKI Jakarta Menggunakan Analisis Faktor dan Klaster. Materi berikut hanya ditujukan untuk kalangan terbatas, meliputi individu/personal yang menerima materi ini secara langsung dari lembaga pelatihan. Materi ini dilarang untuk direproduksi, didistribusikan, diterjemahkan, atau diadaptasikan dalam bentuk apapun di luar izin dari individu dan organisasi yang berkepentingan.

Algoritma adalah pusat pendidikan Data Science di Jakarta. Kami mengadakan workshop dan program pelatihan untuk membantu para profesional dan pelajar untuk mendapatkan keahlian dalam berbagai bidang dalam ruang lingkup Data Science: data visualization, machine learning, data modeling, statistical inference, dan lain-lainnya.

Sebelum masuk ke dalam materi dan menjalankan kode-kode di dalam materi ini, silakan anda melihat bagian Library and Setup untuk melihat dan memastikan semua persyaratan dasar untuk mengikuti materi ini sudah terpenuhi termasuk package-package yang diperlukan. Pada bagian Tujuan Pembelajaran anda dapat melihat secara umum apa saja yang akan dipelajari dalam modul materi ini. Kami harap materi ini akan bermanfaat bagi karir ataupun menambah keahlian peserta.

Preface

Pendahuluan

Pandemi COVID-19 merupakan bentuk bahaya yang memiliki potensi mengancam segala aspek kehidupan masyarakat seperti sosial, kesehatan, dan psikologis. Dampak pada masyarakat Indonesia tidak berbeda jauh dengan negara lain yang juga mengalami COVID 19. Kondisi kerentanan sosial menjadi realitas nyata yang terjadi pada masyarakat dalam menghadapi pandemi COVID-19.

Di era pemulihan ekonomi ini, kerentanan dapat dipandang dari dua arah. Dimana sisi positif kerentanan ini dapat kita pandang sebagai suatu masalah yang dapat diselesaikan dalam ranah pihak CSR. Di sinilah peran industri untuk pro aktif dalam mencari peluang dengan memanfaatkan kerentanan yang ditimbulkan oleh COVID-19. Sementara sisi lainnya, kerentanan ini dapat kita pandang sebagai permasalahan yang disebabkan oleh berbagai aspek yang diakibatkan oleh COVID-19, sehingga dapat ditanggulangi dengan mengintegrasikan perencanaan penanggulangan resiko dengan program pengembangan lainnya.

Materi ini bertujuan memberikan pemahaman kepada peserta workshop terkait penggunaan unsupervised learning yang dapat digunakan untuk project berbasis social behaviour. Beberapa metode yang dapat diimplementasikan diantaranya analisis faktor dalam membentuk faktor-faktor laten kerentanan sosial serta PCA dalam melakukan ekstraksi variabel menjadi komponen yang akan membentuk index kerentanan sosial masing-masing daerah. Hasil index ini dapat kita manfaatkan untuk melakukan visualisasi demografis dan profiling masing-masing daerah menggunakan fuzzy c-means clustering sehingga insight yang diambil akan lebih bermakna. Adapun setelah mempelajari materi ini peserta diharapkan dapat memahami kegunaan dan potensi pemanfaatan analisis faktor dan fuzzy clustering sesuai dengan proses bisnis yang ada di bidang industri yang dijalani oleh peserta.

Library dan Setup

Untuk dapat mengikuti materi ini, peserta diharapkan sudah menginstall beberapa packages di bawah ini. Apabila package tersebut belum terinstall, silahkan jalankan chunk atau baris kode berikut. Apabila sudah ter-install, lewati chunk berikut dan muat package yang dibutuhkan dengan menjalankan chunk selanjutnya.

`1`	`# install.packages(c("readr","tidyverse", "stringr", "rmarkdown", "GGally", "psych", "GPArotation", "factoextra", "BAMMtools", "ppclust", "gridExtra", "leaflet"))`

#data cleaning
library(readxl) # untuk membaca data excel
library(readr) # untuk membaca data
library(tidyverse) # untuk data manipulation
library(stringr) # untuk persiapan data

#data analysis
library(GGally) # untuk membuat matriks korelasi
library(psych) # untuk melakukan analisis KMO
library(GPArotation) # untuk melakukan factor analysis
library(factoextra) # untuk mengekstrak dan visualisasi Hasil analisis data multivariat
library(BAMMtools) # untuk membuat natural breaks pada pemetaan
library(ppclust) # untuk membuat Probabilistic and Possibilistic Cluster Analysis

#data visualization
library(gridExtra) # untuk membuat mixed plot dalam satu hasil
library(leaflet) # untuk visualisasi interaktif

Tujuan Pembelajaran

Tujuan utama dari workshop ini adalah untuk memberikan pengenalan yang komprehensif mengenai tools dan perangkat lunak yang digunakan untuk mengetahui kelompok yang memiliki resiko lebih tinggi (rentan) dan seberapa cepat pulih terhadap dampak COVID-19, yakni sebuah open-source populer: R. Adapun materi ini akan mencakup:

Introductory Module:

Dasar Bahasa Pemrograman R
- Pengenalan bahasa pemrograman R
- Bekerja dengan RStudio
- Fungsi standar pada R
- Data manipulasi dan persiapan data menggunakan R tidyverse
Introduction to Statistics
- Statistika Deskriptif
  - Ukuran penyebaran data (Measure of spread)
  - Hubungan antar data (Variable Relationship)
- Statistika Inferensial
  - Uji hipotesis

**Main Module: **

Data wrangling dan EDA
- Gambaran umum variabel penyusun indeks kerentanan sosial
- Melihat summary, plot distribusi, serta korelasi antar variabel
Mengecek asumsi dalam analisis faktor
- Uji multikolinearitas pada seluruh variabel
- Kaiser–Meyer–Olkin (KMO) - Measure of Sampling Adequacy (MSA)
- The Bartlett’s Test statistics
Membentuk faktor-faktor penyusun kerentanan sosial
- Factor Analysis dalam membentuk faktor-faktor penyusun kerentanan sosial menggunakan fungsi fa()
- fa.parallel() untuk menemukan jumlah faktor optimal
- Melakukan ekstraksi dan rotating faktor
- Melakukan ekstraksi dengan menghasilkan faktor umum menggunakan fa.diagram()
- Pembobotan dalam agregasi indeks.
Melakukan pembentukan Indeks Kerentanan Sosial
- Melakukan Principle Component Analysis (PCA) untuk mengambil komponen yang tereduksi.
- Membentuk SoVI (Social Vulnerability Index) masing-masing kecamatan
Visualisasi Geospasial Kerentanan Sosial Kecamatan di Provinsi DKI Jakarta
- Indeks kerentanan sosial dikategorikan menjadi tiga yaitu : rendah, sedang dan tinggi.Pengkategorian menggunakan Natural Breaks pada R menggunakan getJenksBreaks()
- Membuat visualisasi geocoding berupa choropleth dengan leaflet()
Fuzzy C-Means Clustering terhadap Kerentanan Sosial
- Pemilihan nilai k optimum dengan elbow method
- FCM menggunakan fcm()
- Melakukan cluster profiling

Ground Truth Story

Sampai saat ini, pandemi COVID-19 belum juga usai. Provinsi DKI Jakarta menjadi provinsi dengan jumlah kasus positif COVID-19 terbanyak di Indonesia. Meski disadari bahwa virus COVID-19 dapat menginfeksi siapa saja, tetapi beberapa kelompok orang memiliki tingkat risiko yang lebih tinggi untuk terpapar. Sebagai contoh, jika dilihat pada data kematian akibat COVID-19 per kelompok umur per tanggal 30 September 2022, penduduk rentan seperti usia lansia (≥ 60 tahun) menyumbang kematian sebesar 47,6% dari total kematian akibat COVID-19.

Selain itu, pada kasus positif COVID-19 dan sedang dirawat atau isolasi, proporsi penduduk berjenis kelamin perempuan lebih banyak dibandingkan laki-laki (Kemenkes, 2021). Tingkat kerentanan sosial masyarakat telah banyak dilakukan melalui indeks kerentanan sosial. Kerentanan sosial wilayah menggambarkan kerapuhan sosial dari suatu wilayah akibat pengaruh dari adanya bahaya. Faktor-faktor kerentanan sosial menjadi salah satu hal yang krusial di tengah pandemi COVID-19. Hal ini dikarenakan gabungan faktor kerentanan sosial tersebut dapat meningkatkan risiko terkena COVID-19 dan risiko pulih dari kondisi setelah terpapar COVID-19.

Tujuan Analisis

Untuk mengetahui kelompok yang memiliki resiko lebih tinggi (rentan) dan seberapa cepat pulih terhadap dampak COVID-19 menjadi penting, mengingat hal ini dapat meningkatkan kesadaran, mencegah terinfeksi COVID-19, dan menurunkan tingkat kematian.

Data Loading

Pada analisis ini, dilakukan observasi pada level kecamatan dengan jumlah observasi sebanyak 44 kecamatan di Provinsi DKI Jakarta. Hal ini diakibatkan observasi level kecamatan memiliki data yang terbaru dan paling lengkap. Adapun variabel yang digunakan berasal dari berbagai sumber dengan tanggal terakhir diambil yaitu 30 April 2021, yaitu:

Website Open Data DKI Jakarta berupa variabel:
- persentase penduduk perempuan,
- persentase kepala keluarga berjenis kelamin perempuan,
- persentase penduduk lansia,
- jumlah pengangguran,
- persentase penduduk wajib KTP,
- jumlah daerah rawan banjir,
- persentase pekerja di sektor informal, dan
- jumlah puskesmas.
Website Open Data COVID-19 Provinsi DKI Jakarta berupa variabel:
- jumlah kasus positif COVID-19,
- jumlah penduduk yang dirawat akibat COVID-19, dan
- jumlah penduduk yang melakukan isolasi mandiri akibat COVID-19.
Publikasi Kabupaten/Kota dalam Angka oleh BPS untuk setiap kota dan kabupaten di Provinsi DKI Jakarta berupa variabel:
- jumlah RS umum

Berikut adalah 14 variabel yang digunakan:

1
2

rentan <- readxl::read_xlsx("dataset/df_kerentanansosial.xlsx")
rentan

#> # A tibble: 44 × 15
#>    Kecamatan perem…¹ KK_pe…² lansia penga…³ wajib…⁴ banjir infor…⁵    RS puske…⁶
#>    <chr>       <dbl>   <dbl>  <dbl>   <dbl>   <dbl>  <dbl>   <dbl> <dbl>   <dbl>
#>  1 Cakung       49.2    18.1   6.22    2909    70.2     31    6.50     1       9
#>  2 Cempaka …    50.1    23.0  10.9      544    75.1      0   13.9      2       4
#>  3 Cengkare…    49.3    18.5   6.59    2839    70.8     55    6.79     1      10
#>  4 Cilandak     50.4    20.8   9.31     584    73.7      2   10.9      3       6
#>  5 Cilincing    49.5    20.5   5.98    2746    74.4     38   11.3      1      10
#>  6 Cipayung     49.6    17.1   7.00    2388    71.2     15   14.3      0      11
#>  7 Ciracas      49.7    18.7   7.75    1818    71.6      8    7.19     1       6
#>  8 Duren Sa…    50.1    19.6   9.27    2071    72.2      2    6.43     2      12
#>  9 Gambir       49.5    25.1  13.2      702    76.5      0   20.0      1       4
#> 10 Grogolpe…    50.3    22.9  11.4     1187    75.1     28   15.1      3      10
#> # … with 34 more rows, 5 more variables: positiv_cov <dbl>, rawat <dbl>,
#> #   isoman <dbl>, Kode_shp <dbl>, Kode <chr>, and abbreviated variable names
#> #   ¹perempuan, ²KK_perempuan, ³pengangguran, ⁴wajib_ktp, ⁵informal, ⁶puskesmas

Berikut adalah deskripsi data untuk rentan:

Kecamatan = Nama kecamatan pada masing-masing kota administratif pada Provinsi DKI Jakarta
perempuan = persentase penduduk perempuan
KK_perempuan = persenyase kepala keluarga berjenis kelamin perempuan
lansia = persentase penduduk lansia
pengangguran = jumlah pengangguran
wajib_ktp = persentase penduduk wajib KTP
banjir = jumlah daerah rawan banjir
informal = persentase pekerja di sektor informal
RS = jumlah RS umum
puskesmas = jumlah puskesmas
positif_cov = jumlah kasus positif akibat COVID-19 (data 30 april 2021)
rawat = jumlah penduduk yang dirawat akibat COVID-19 (data 30 april 2021)
isoman = jumlah penduduk yang melakukan isolasi mandiri akibat COVID-19 (data 30 april 2021)
kode_shp = kode wilayah untuk kebutuhan pemetaan peta.

Data Cleansing

Melakukan inspeksi untuk melihat keseluruhan data.

`1`	`glimpse(rentan)`

#> Rows: 44
#> Columns: 15
#> $ Kecamatan    <chr> "Cakung", "Cempaka Putih", "Cengkareng", "Cilandak", "Cil…
#> $ perempuan    <dbl> 49.2096, 50.0816, 49.3437, 50.4252, 49.4534, 49.6205, 49.…
#> $ KK_perempuan <dbl> 18.05886, 23.00326, 18.54968, 20.79268, 20.51346, 17.1236…
#> $ lansia       <dbl> 6.217495, 10.852327, 6.585898, 9.310706, 5.975218, 7.0006…
#> $ pengangguran <dbl> 2909, 544, 2839, 584, 2746, 2388, 1818, 2071, 702, 1187, …
#> $ wajib_ktp    <dbl> 70.20066, 75.12540, 70.84496, 73.73298, 74.35980, 71.1721…
#> $ banjir       <dbl> 31, 0, 55, 2, 38, 15, 8, 2, 0, 28, 12, 32, 2, 0, 4, 4, 0,…
#> $ informal     <dbl> 6.503218, 13.910844, 6.791649, 10.869546, 11.266610, 14.2…
#> $ RS           <dbl> 1, 2, 1, 3, 1, 0, 1, 2, 1, 3, 3, 2, 0, 1, 3, 2, 4, 2, 1, …
#> $ puskesmas    <dbl> 9, 4, 10, 6, 10, 11, 6, 12, 4, 10, 6, 12, 7, 13, 11, 9, 8…
#> $ positiv_cov  <dbl> 11740, 4998, 13036, 8201, 7899, 10437, 9529, 14329, 3796,…
#> $ rawat        <dbl> 82, 50, 109, 62, 73, 63, 73, 139, 25, 67, 80, 55, 25, 73,…
#> $ isoman       <dbl> 102, 17, 146, 114, 63, 124, 108, 145, 42, 43, 175, 87, 19…
#> $ Kode_shp     <dbl> 38, 5, 15, 28, 12, 42, 41, 39, 1, 16, 31, 35, 8, 20, 29, …
#> $ Kode         <chr> "IDN.7.4.1_1", "IDN.7.2.1_1", "IDN.7.1.1_1", "IDN.7.3.1_1…

# ubah tipe data Kode_Shp sebagai character
rentan <- rentan %>% 
  mutate(Kode_shp = as.character(Kode_shp))
glimpse(rentan)

#> Rows: 44
#> Columns: 15
#> $ Kecamatan    <chr> "Cakung", "Cempaka Putih", "Cengkareng", "Cilandak", "Cil…
#> $ perempuan    <dbl> 49.2096, 50.0816, 49.3437, 50.4252, 49.4534, 49.6205, 49.…
#> $ KK_perempuan <dbl> 18.05886, 23.00326, 18.54968, 20.79268, 20.51346, 17.1236…
#> $ lansia       <dbl> 6.217495, 10.852327, 6.585898, 9.310706, 5.975218, 7.0006…
#> $ pengangguran <dbl> 2909, 544, 2839, 584, 2746, 2388, 1818, 2071, 702, 1187, …
#> $ wajib_ktp    <dbl> 70.20066, 75.12540, 70.84496, 73.73298, 74.35980, 71.1721…
#> $ banjir       <dbl> 31, 0, 55, 2, 38, 15, 8, 2, 0, 28, 12, 32, 2, 0, 4, 4, 0,…
#> $ informal     <dbl> 6.503218, 13.910844, 6.791649, 10.869546, 11.266610, 14.2…
#> $ RS           <dbl> 1, 2, 1, 3, 1, 0, 1, 2, 1, 3, 3, 2, 0, 1, 3, 2, 4, 2, 1, …
#> $ puskesmas    <dbl> 9, 4, 10, 6, 10, 11, 6, 12, 4, 10, 6, 12, 7, 13, 11, 9, 8…
#> $ positiv_cov  <dbl> 11740, 4998, 13036, 8201, 7899, 10437, 9529, 14329, 3796,…
#> $ rawat        <dbl> 82, 50, 109, 62, 73, 63, 73, 139, 25, 67, 80, 55, 25, 73,…
#> $ isoman       <dbl> 102, 17, 146, 114, 63, 124, 108, 145, 42, 43, 175, 87, 19…
#> $ Kode_shp     <chr> "38", "5", "15", "28", "12", "42", "41", "39", "1", "16",…
#> $ Kode         <chr> "IDN.7.4.1_1", "IDN.7.2.1_1", "IDN.7.1.1_1", "IDN.7.3.1_1…

Melakukan pengecekan missing value

Melakukan pengecekan missing value pada data untuk setiap variabel.

`1`	`anyNA(rentan)`

#> [1] FALSE

1
2

# pak albert
is.null(rentan)

#> [1] FALSE

1
2

# pak anthony
colSums(is.na(rentan))

#>    Kecamatan    perempuan KK_perempuan       lansia pengangguran    wajib_ktp 
#>            0            0            0            0            0            0 
#>       banjir     informal           RS    puskesmas  positiv_cov        rawat 
#>            0            0            0            0            0            0 
#>       isoman     Kode_shp         Kode 
#>            0            0            0

Insight: data yang kita punya tidak memiliki missing value

Melakukan pengecekan duplikasi data

Cek terlebih dahulu apakah terdapat duplikasi pada data yang ada.

1
2
3

rentan %>% 
  duplicated() %>% 
  sum()

#> [1] 0

1
2

# pak mula
anyDuplicated(rentan)

#> [1] 0

Insight: data rentan tidak memiliki observasi yang duplikat.

Data rentan sudah bersih dari missing value dan juga duplikasi.

Exploratory Data Analysis

Sebelum dilakukannya analisis lebih lanjut, alangkah lebih baik kita mengeksplorasi data kerentanan sosial dalam bentuk visualisasi.

1. Gambaran umum variabel penyusun indeks kerentanan sosial

Pada step ini dilakukan untuk memahami data terlebih dahulu dan mencoba mengumpulkan sebanyak mungkin wawasan dari data. Hal ini dapat dilakukan menggunakan five number summary

Exploratory Data Analysis refers to the critical process of performing initial investigations on data so as to discover patterns,to spot anomalies,and to check the summary statistics and graphical representations

1
2

# Please run the code down below
summary(rentan)

#>   Kecamatan           perempuan      KK_perempuan       lansia      
#>  Length:44          Min.   :48.27   Min.   :14.78   Min.   : 5.873  
#>  Class :character   1st Qu.:49.43   1st Qu.:18.93   1st Qu.: 7.844  
#>  Mode  :character   Median :49.67   Median :21.11   Median : 8.845  
#>                     Mean   :49.71   Mean   :21.09   Mean   : 8.997  
#>                     3rd Qu.:49.98   3rd Qu.:22.94   3rd Qu.: 9.776  
#>                     Max.   :50.98   Max.   :27.49   Max.   :13.219  
#>   pengangguran      wajib_ktp         banjir         informal     
#>  Min.   : 207.0   Min.   :70.20   Min.   : 0.00   Min.   : 6.148  
#>  1st Qu.: 596.2   1st Qu.:72.06   1st Qu.: 0.00   1st Qu.: 9.543  
#>  Median : 931.0   Median :73.04   Median : 9.00   Median :10.540  
#>  Mean   :1244.9   Mean   :73.57   Mean   :12.86   Mean   :12.412  
#>  3rd Qu.:1666.5   3rd Qu.:74.52   3rd Qu.:18.75   3rd Qu.:16.318  
#>  Max.   :3468.0   Max.   :79.73   Max.   :59.00   Max.   :25.738  
#>        RS          puskesmas      positiv_cov        rawat       
#>  Min.   :0.000   Min.   : 3.00   Min.   :  242   Min.   :  0.00  
#>  1st Qu.:1.000   1st Qu.: 6.00   1st Qu.: 5698   1st Qu.: 49.75  
#>  Median :1.500   Median : 7.50   Median : 7966   Median : 62.50  
#>  Mean   :1.818   Mean   : 7.75   Mean   : 7587   Mean   : 65.61  
#>  3rd Qu.:3.000   3rd Qu.: 9.25   3rd Qu.: 9622   3rd Qu.: 77.00  
#>  Max.   :7.000   Max.   :14.00   Max.   :14329   Max.   :142.00  
#>      isoman         Kode_shp             Kode          
#>  Min.   :  0.00   Length:44          Length:44         
#>  1st Qu.: 41.75   Class :character   Class :character  
#>  Median : 62.50   Mode  :character   Mode  :character  
#>  Mean   : 67.93                                        
#>  3rd Qu.: 88.50                                        
#>  Max.   :175.00

Insight: Berdasarkan hasil summary(rentan):

Rata-rata persentase penduduk perempuan untuk masing-masing kecamatan adalah sebesar 49.71%. Ternyata 21.09% adalah KK/kepala keluarga.
Ada sebanyak 59 titik banjir di suatu kecamatan tertentu sebagai nilai maksimum dan juga rata-rata daerah rawan banjir di DKI Jakarta sebanyak 13 titik.
[Pak Albert] positif covid bnyak tp yg dirawat dan isoman dikit
[Pak Dary] Jumlah kasus positif akibat covid-19 yang terbesar sebanyak 14329 orang sedangkan penduduk yang dirawat sebanyak 142 orang

2. Melihat visualisasi antar beberapa variabel

Pertama dilakukan pemeriksaan nilai korelasi antar variabel dengan melihat scatter plot antar masing-masing nilai variabel berpasangan pada data. Dari library GGally

1
2

#rentan[baris, kolom]
rentan[,c(2:13)]

#> # A tibble: 44 × 12
#>    perempuan KK_pe…¹ lansia penga…² wajib…³ banjir infor…⁴    RS puske…⁵ posit…⁶
#>        <dbl>   <dbl>  <dbl>   <dbl>   <dbl>  <dbl>   <dbl> <dbl>   <dbl>   <dbl>
#>  1      49.2    18.1   6.22    2909    70.2     31    6.50     1       9   11740
#>  2      50.1    23.0  10.9      544    75.1      0   13.9      2       4    4998
#>  3      49.3    18.5   6.59    2839    70.8     55    6.79     1      10   13036
#>  4      50.4    20.8   9.31     584    73.7      2   10.9      3       6    8201
#>  5      49.5    20.5   5.98    2746    74.4     38   11.3      1      10    7899
#>  6      49.6    17.1   7.00    2388    71.2     15   14.3      0      11   10437
#>  7      49.7    18.7   7.75    1818    71.6      8    7.19     1       6    9529
#>  8      50.1    19.6   9.27    2071    72.2      2    6.43     2      12   14329
#>  9      49.5    25.1  13.2      702    76.5      0   20.0      1       4    3796
#> 10      50.3    22.9  11.4     1187    75.1     28   15.1      3      10    8298
#> # … with 34 more rows, 2 more variables: rawat <dbl>, isoman <dbl>, and
#> #   abbreviated variable names ¹KK_perempuan, ²pengangguran, ³wajib_ktp,
#> #   ⁴informal, ⁵puskesmas, ⁶positiv_cov

1
2
3

ggpairs(rentan[,c(2:13)], upper = list(continuous = wrap("cor", size=2)))+
  theme(axis.text = element_text(colour = "black", size = 4),
        strip.text = element_text(size = 5))

Insight:

Nilai korelasi variabel lansia dan kk_perempuan besarnya adalah 0.813. Hal ini mengindikasikan adanya hubungan yang kuat antara 2 variabel tersebut.
KK_perempuan dan juga perempuan memiliki distribusi normal. Variabel yang lain memiliki distribusi yg skewed

Penyusun Kerentanan Sosial

Pada bagian ini dilakukan tiga teknik analisis untuk melihat bagaimana kerentanan sosial yang diakibatkan COVID-19.

1. Digunakan metode Factor Analysis untuk membentuk faktor umum dan interpretasi pada faktor umum dan menentukan variabel dari masing masing komponen yang terbentuk.
1. Digunakan metode PCA untuk meringkas informasi serta akan dipilih beberapa PC untuk dapat merangkum informasi yang dibutuhkan.
1. Kemudian melakukan pengklusteran menggunakan metode non-hierarchical clustering (Fuzzy C-Means Clustering)

Analisis Faktor / Factor Analysis dan Asumsi

Analisis faktor termasuk salah satu teknik statistika yang dapat digunakan untuk memberikan deskripsi yang relatif sederhana melalui reduksi jumlah variabel yang disebut faktor.

Analisis faktor dipergunakan untuk mereduksi data atau meringkas, dari variabel lama yang banyak diubah menjadi sedikit variabel baru yang disebut faktor, dan masih memuat sebagian besar informasi yang terkandung dalam variabel asli.

Analisis faktor dipergunakan dalam kondisi sebagai berikut :

Mengenali atau mengidentifikasi faktor, yang menjelaskan korelasi antara suatu set variabel.
Mengenali atau mengidentifikasi suatu set variabel baru yang independent dan lebih sedikit jumlahnya untuk menggantikan suatu set variabel asli yang saling berkorelasi di dalam analisis multivariat selanjutnya.
Mengidentifikasi suatu set variabel yang penting dari suatu set variabel yang lebih banyak jumlahnya untuk dipergunakan dalam analisis multivariat selanjutnya.

Karakteristik Analisis Faktor

Variabel-variabel dalam satu faktor mempunyai korelasi yang tinggi
Korelasi dengan variabel-variabel pada faktor lain relatif rendah
Tiap-tiap kelompok dari variabel mewakili suatu konstruksi dasar yang disebut faktor.

Melakukan uji multikolinearitas pada seluruh variabel.

Tujuan dari membentuk plot korelasi adalah untuk menguji kolinearitas dari data. Asumsi pada analisis faktor yaitu seluruh variabel harus memiliki korelasi, tetapi tidak terlalu tinggi. Silahkan gunakan ggcorr() untuk memudahkan melihat heat map korelasi antar data.

Kita akan mengecek bahwa variabel-variabel yang kita punya tidak memiliki korelasi yg tingi
Variabel-variabel tsb harus mempunyai nilai korelasi (gaboleh 0)

`1`	`ggcorr(rentan, label = T, label_round = 3, label_size = 2)`

Correlation plot di atas yang digunakan untuk menguji multikolinearitas dari data. Pada analisis ini, digunakan batas (>0,85) sebagai batas korelasi ekstrem terhadap setiap variabel.

Insight: Variabel-variabel yang kita punya semuanya miliki nilai korelasi (tidak ada yang tidak berhubungan) dan juga tidak ada yang melebihi 0.85. Sehingga asumsi multikolinearitas terpenuhi.

Kaiser–Meyer–Olkin (KMO) - Measure of Sampling Adequacy(MSA)

KMO-MSA (proporsi korelasi) adalah uji yang dilakukan untuk menentukan kelayakan (appropriateness) dari suatu analisis faktor yang akan dilakukan.

Skala uji KMO-MSA berkisar antara 0 sampai 1.

Jika nilai KMO-MSA lebih rendah dari 0.5, maka analisis faktor tidak layak dilakukan. Sedangkan jika nilai KMO lebih besar dari 0.5, maka analisis faktor layak dilakukan.

Untuk dapat melakukan pengujian KMO-MSA, pada R disediakan suatu function KMO() dari library psych.

Pertama membentuk matriks korelasi dari variabel numerik pada data

`1`	`rentan_cor <- cor(rentan[,2:13])`

Mencari nilai KMO-MSA dengan KMO() dari library psych menggunakan matriks korelasi yang sudah dibuat sebelumnya

1
2

rentan_KMO <- KMO(rentan_cor)
rentan_KMO

#> Kaiser-Meyer-Olkin factor adequacy
#> Call: KMO(r = rentan_cor)
#> Overall MSA =  0.74
#> MSA for each item = 
#>    perempuan KK_perempuan       lansia pengangguran    wajib_ktp       banjir 
#>         0.59         0.65         0.74         0.68         0.75         0.85 
#>     informal           RS    puskesmas  positiv_cov        rawat       isoman 
#>         0.71         0.54         0.87         0.72         0.83         0.75

Insight:

Untuk nilai overall MSA : 0.74 -> data yang kita miliki sudah cukup untuk dilakukan analisis faktor
Dari hasil output nilai MSA untuk masing-masing variabel: Tidak ada nilai yang besarnya < 0.5, sehingga semua variabel layak untuk masuk ke dalam analisis faktor

The Bartlett’s Test statistics

Bartlett’s test digunakan untuk mengkonfirmasi bahwa data berkorelasi secara berpasangan (Mickey & Sharma, 1997). Dengan kata lain uji bartlett digunakan untuk mengetahui homogenitas variance.

Hipotesis uji
- H0 = Varians di antara masing-masing variabel adalah sama (variance homogen)
- H1 = Setidaknya satu variabel memiliki varians yang tidak sama dengan yang lain (variance tidak homogen)
Pengujian yang diharapkan adalah data memiliki variance yang tidak homogen (tolak H0).

Di R fungsi yang dapat digunakan adalah bartlett.test() dari base. nilai alpha = 0.05

`1`	`bartlett.test(rentan[,c(2:13)])`

#> 
#> 	Bartlett test of homogeneity of variances
#> 
#> data:  rentan[, c(2:13)]
#> Bartlett's K-squared = 4496.9, df = 11, p-value < 0.00000000000000022

Apabila nilai p-value < alpha = tolak H0

Insight: Nilai p-value < 0.05/ nilai alpha. Sehingga data kita memiliki variance yang tidak homogen/ berbeda satu sama lainnya.

Membentuk Analisis Faktor

Untuk membentuk faktor-faktor penyusun kerentanan sosial digunakan metode analisis faktor, fungsi yang tersedia di R adalah fa() dari library psych.

Parameter fungsi fa()

r: Correlation matrix dari data awal

1
2

fa_rentan <- fa(r = rentan_cor)
fa_rentan

#> Factor Analysis using method =  minres
#> Call: fa(r = rentan_cor)
#> Standardized loadings (pattern matrix) based upon correlation matrix
#>                MR1      h2    u2 com
#> perempuan     0.03 0.00083 0.999   1
#> KK_perempuan -0.43 0.18150 0.818   1
#> lansia       -0.53 0.27661 0.723   1
#> pengangguran  0.65 0.41698 0.583   1
#> wajib_ktp    -0.62 0.38288 0.617   1
#> banjir        0.56 0.31535 0.685   1
#> informal     -0.57 0.32381 0.676   1
#> RS            0.19 0.03505 0.965   1
#> puskesmas     0.70 0.48876 0.511   1
#> positiv_cov   0.97 0.93146 0.069   1
#> rawat         0.75 0.56067 0.439   1
#> isoman        0.78 0.60599 0.394   1
#> 
#>                 MR1
#> SS loadings    4.52
#> Proportion Var 0.38
#> 
#> Mean item complexity =  1
#> Test of the hypothesis that 1 factor is sufficient.
#> 
#> The degrees of freedom for the null model are  66  and the objective function was  9.22
#> The degrees of freedom for the model are 54  and the objective function was  4.63 
#> 
#> The root mean square of the residuals (RMSR) is  0.18 
#> The df corrected root mean square of the residuals is  0.2 
#> 
#> Fit based upon off diagonal values = 0.81

Secara default, objek factor analysis yang terbentuk menghasilkan sejumlah 1 faktor. Jumlah ini dapat saja belum optimal, sehingga perlu untuk mencari jumlah faktor yang optimal dari data yang kita miliki.

Menentukan jumlah faktor optimal

Secara teori, kita dapat mengambil faktor sebanyak variabel yang kita miliki. Tapi hal tersebut tentunya berkebalikan dengan tujuan analisis faktor.

Sehingga, akan dicari faktor laten dalam kumpulan variabel. Oleh karena itu, kita hanya akan tertarik pada faktor-faktor yang menjelaskan proporsi variasi yang substansial dalam data.

Langkah di balik penentuan jumlah faktor optimal:

Mencari eigen value: variance/informasi yang dirangkum oleh tiap faktor dengan eigen().

1
2

ev <- eigen(rentan_cor) # get eigenvalues
ev$values

#>  [1] 4.99667306 2.46039233 1.12792938 0.99712076 0.65274905 0.49540729
#>  [7] 0.37580580 0.36473044 0.24043800 0.13934914 0.10192060 0.04748414

Membentuk scree plot: grafik yang menunjukkan relasi antara faktor dengan nilai Eigennya

Di R terdapat sebuah fungsi yaitu fa.parallel(). Fungsi ini mencari jumlah faktor optimal dari nilai eigen yang berurutan, kemudian penurunan tajam dalam plot menunjukkan jumlah komponen atau faktor yang sesuai untuk diekstraksi. Parameter: x = data frame yang berisikan variabel numerik fa = fa untuk factor analysis

1
2

nfactors <- fa.parallel(x = rentan[,c(2:13)],
                        fa = "fa")

#> Parallel analysis suggests that the number of factors =  2  and the number of components =  NA

`1`	`nfactors`

#> Call: fa.parallel(x = rentan[, c(2:13)], fa = "fa")
#> Parallel analysis suggests that the number of factors =  2  and the number of components =  NA 
#> 
#>  Eigen Values of 
#> 
#>  eigen values of factors
#>  [1]  4.52  1.69  0.50  0.24 -0.04 -0.07 -0.11 -0.18 -0.35 -0.42 -0.60 -0.66
#> 
#>  eigen values of simulated factors
#>  [1]  1.28  0.75  0.55  0.38  0.21  0.09 -0.04 -0.16 -0.29 -0.38 -0.49 -0.60
#> 
#>  eigen values of components 
#>  [1] 5.00 2.46 1.13 1.00 0.65 0.50 0.38 0.36 0.24 0.14 0.10 0.05
#> 
#>  eigen values of simulated components
#> [1] NA

Insight: Jumlah faktor yang akan kita gunakan untuk kebutuhan analisis sebanyak 2 buah faktor.

*** END OF DAY 2 *** *** START OF DAY 3 ***

Membentuk Analisis Faktor dengan Jumlah Faktor yang ditentukan

Selanjutnya akan kembali dilakukan analisis faktor menggunakan jumlah faktor optimal yang sudah didapatkan.

Parameter tambahan fungsi fa():

r = Correlation matrix dari data awal
nfactors: banyaknya faktor yang akan dibentuk
rotate : tipe rotasi. Rotasi meminimalkan kompleksitas faktor untuk membuat struktur lebih sederhana untuk ditafsirkan. Beberapa tipe rotasi yang dapat dipilih:

Rotasi ortogonal (varimax, quartimax, equamax) -> Jika ingin memulai dengan asumsi faktor saling tidak berkorelasi
Rotasi oblique (oblimin, promax) -> asumsi faktor saling berkorelasi

# jumlah faktor optimal yang didapat
nfact <- nfactors$nfact

# factor analysis dengan tipe rotasi varimax
fa_rentan_opt <- fa(r = rentan_cor,
                    nfactors = 2, #jumlah faktor optimal yang didapat dari fa_paralel
                    rotate= "varimax")
fa_rentan_opt

#> Factor Analysis using method =  minres
#> Call: fa(r = rentan_cor, nfactors = 2, rotate = "varimax")
#> Standardized loadings (pattern matrix) based upon correlation matrix
#>                MR1   MR2   h2      u2 com
#> perempuan     0.21  0.57 0.37  0.6313 1.3
#> KK_perempuan -0.24  0.80 0.70  0.2972 1.2
#> lansia       -0.33  0.90 0.93  0.0705 1.3
#> pengangguran  0.60 -0.20 0.40  0.5982 1.2
#> wajib_ktp    -0.58  0.18 0.37  0.6297 1.2
#> banjir        0.56 -0.07 0.32  0.6826 1.0
#> informal     -0.43  0.61 0.55  0.4521 1.8
#> RS            0.31  0.32 0.20  0.7994 2.0
#> puskesmas     0.69 -0.11 0.48  0.5152 1.1
#> positiv_cov   1.00 -0.03 1.01 -0.0074 1.0
#> rawat         0.87  0.16 0.78  0.2229 1.1
#> isoman        0.79 -0.06 0.63  0.3660 1.0
#> 
#>                        MR1  MR2
#> SS loadings           4.36 2.38
#> Proportion Var        0.36 0.20
#> Cumulative Var        0.36 0.56
#> Proportion Explained  0.65 0.35
#> Cumulative Proportion 0.65 1.00
#> 
#> Mean item complexity =  1.3
#> Test of the hypothesis that 2 factors are sufficient.
#> 
#> The degrees of freedom for the null model are  66  and the objective function was  9.22
#> The degrees of freedom for the model are 43  and the objective function was  2.13 
#> 
#> The root mean square of the residuals (RMSR) is  0.08 
#> The df corrected root mean square of the residuals is  0.1 
#> 
#> Fit based upon off diagonal values = 0.96

Dari output di atas dengan menggunakan faktor sebanyak 2 didapatkan beberapa informasi sebagai berikut:

SS Loadings: Sum square loadings (jumlah kuadrat masing-masing nilai loading)
Proportion Var: Memberitahu kita berapa banyak varians/informasi keseluruhan yang diperhitungkan faktor dari semua variabel.
Cumulative Var: Nilai akumulasi dari proportion var.
Proportion Explained: Jumlah variansi relatif yang terjelaskan. Merupakan pembagian dari proportion var/ sum(proportion var).
Cumulative Proportion: Jumlah akumulasi dari proportion explained.
MR: Factor Loading -> kekuatan hubungan dari masing-masing variabel terhadap faktor utama.
h2: jumlah varians dalam item/variabel yang dijelaskan oleh faktor-faktor (yang dipertahankan) -> communality
u2: merupakan nilai 1-h2 adalah residual variance.
com: item complexity.

Dengan menggunakan 2 faktor, maka informasi keseluruhan yang diperhitungkan faktor dari semua variabel sebanyak 0.56

Pembentukan Faktor Umum dan Interpretasi

Salah satu tujuan dari factor analysis adalah melakukan ekstraksi dengan menghasilkan faktor utama Faktor utama yang dihasilkan nantinya akan dijadikan sebagai faktor-faktor pembentuk kerentanan sosial pada data. Pembentukan faktor umum dilakukan dengan menggunakan fungsi fa.diagram() dengan parameter:

fa.results = object dari output factor analysis
rsize = layout tampilan (opsional)

1
2

fa.diagram(fa.results = fa_rentan_opt,
           rsize = 2)

Arah dari masing-masing variabel ditentukan berdasarkan pengaruhnya terhadap kerentanan sosial. Arah positif berpengaruh terhadap peningkatan kerentanan sosial, sebaliknya arah negatif menurunkan tingkat kerentanan sosial.

Dari plot di atas dilakukan pengekstrakan terhadap 2 faktor, yaitu faktor 1 dan faktor 2

Interpretasi Faktor:

Faktor kesehatan, kependudukan/ Socio-health
- Faktor ini terdiri dari 7 variabel pembentuk, yaitu positiv_cov, rawat, isoman, puskesmas, pengangguran, wajib_ktp, dan banjir.
- Pada faktor ini menjelaskan 36% varians yang membentuk kerentanan sosial.
- 6 variabel berpengaruh positif dan 1 variabel berpengaruh negatif dalam kerentanan sosial.
- Insight dari variabel positiv_cov, semakin meningkat angka positif terhadap covid, hal tersebut meningkatkan kerentanan suatu daerah tertentu.
- Insight dari variabel wajib_ktp, masyarakat yang memiliki KTP lebih mudah mendapatkan akses terhadap pendataan, sehingga hal tersebut dapat menurunkan kerentanan sosial daerah.
Faktor demografi, perempuan, dan pekerjaan
- Faktor ini terdiri dari 5 variabel pembentuk, yaitu KK_perempuan, lansia, perempuan, informal, RS.
- Pada faktor ini menjelaskan 20% varians yang membentuk kerentanan sosial.
- 5 variabel berpengaruh positif dan 0 variabel berpengaruh negatif dalam kerentanan sosial.
- Insight dari variabel lansia, semakin banyak penduduk lansia hal tersebut mengakibatkan kerentanan daerah semakin tinggi
- Insight dari variabel perempuan, perempuan meningkatkan kerentanan suatu daerah.

Indeks Kerentanan Sosial

Ada banyak metode untuk melakukan pembobotan dalam agregasi indeks. Beberapa di antaranya sudah pernah dilakukan oleh beberapa peneliti:

(Cutter et al., 2003) menggunakan pembobot tertimbang,
(Solangaarachchi et al., 2012) menggunakan pembobot berdasarkan kontribusi terhadap varians terjelaskan,
(Siagian et al., 2014) menggunakan pembobot tidak sama, dan
(Rygel et al., 2006) menggunakan pareto-ranking. Tidak ada kesepakatan yang membatasi penggunaan pembobotan yang berbeda.

Dalam analisis ini, digunakan pembobotan tidak sama untuk melakukan agregasi indeks. Pembobot didapatkan dari pembagian antara varians terjelaskan masing-masing komponen dibagi dengan total varians terjelaskan atau nilai proportion explained dari object fa_rentan_opt. Dengan cara ini komponen yang memiliki varians terjelaskan yang lebih tinggi berpengaruh lebih besar terhadap kerentanan sosial. Perhitungan kerentanan sosial mengikuti persamaan sebagai berikut:

$$SoVI=|(bobot_{faktor_1}*komponen_1)+...+(bobot_{faktor_n}*komponen_n)|$$ $$SoVI=|(0.65*komponen_1)+(0.35*komponen_2)|$$

SoVI merujuk peda Social Vulnarability Index.

Melakukan PCA pada data `rentan`

Principal Component Analysis (PCA)

PCA merupakan suatu algoritma unsupervised learning yang bertujuan untuk dimentionality reduction. Salah satu tujuan dari PCA adalah untuk mengubah variabel asli yang saling berkorelasi satu sama lainnya menjadi satu set variabel baru yang memiliki korelasi lebih kecil atau saling bebas.

Ide dasar dari PCA adalah untuk membuat sumbu (axis) baru yang dapat menangkap informasi sebesar mungkin. Sumbu baru ini adalah yang dinamakan sebagai Principal Component (PC). Untuk melakukan dimensionality reduction, kita akan memilih beberapa PC untuk dapat merangkum informasi yang dibutuhkan. Nilai masing-masing PC dibentuk dari matriks varian-kovarians dari suatu set variabel melalui kombinasi linear dari variabel-variabel tersebut.

$$ PC1 = a_{11}x_1 + a_{12}x_2 + … + a_{1p}x_p $$ di mana:

vektor $a_1p$ adalah eigen vector dari suatu matriks varian kovarians
Nilai $x_i$ adalah nilai suatu data pada baris ke-i

Hal-hal yang menjadi pedoman dalam menentukan jumlah PC, yaitu:

1. Mempertahankan PC yang cukup untuk menjelaskan persentase yang spesifik dari total variansi/informasi
1. Mempertahankan PC yang memiliki eigen value lebih besar dibandingkan rata-rata eigen value dari setiap PC

PCA bertujuan untuk melakukan reduksi data dari sekumpulan variabel (Mickey & Sharma, 1997).

Pada analisis ini, PCA digunakan untuk menghitung nilai SoVI masing-masing kecamatan yang ada pada dataset dengan cara mengambil sejumlah komponen hasil ekstraksi dari PCA. Untuk itu akan dilakukan inspeksi untuk pca$x yaitu nilai hasil proyeksi titik ke PC untuk tiap baris. Digunakan untuk mendapatkan nilai data yang baru.

scale. = TRUE -> meminta function prcomp untuk melakukan scaling pada data.

# membuat object PCA
rentan_pca<-prcomp(rentan[,2:13], scale. = T)
# mendapatkan nilai data baru untuk setiap komponen
as.data.frame(rentan_pca$x)

#>            PC1         PC2           PC3         PC4          PC5         PC6
#> 1   3.48456033  1.56233715 -0.7262031306 -0.64374876 -0.163406655 -0.74142914
#> 2  -2.61111513 -0.60474909  0.4623937304  0.46275127 -0.410487137  0.14726787
#> 3   4.69523879  0.59195276 -1.1135726052 -0.54173232 -1.069298346 -1.26127490
#> 4  -0.07868070 -1.03628501  1.8524291410  0.96318139 -0.221424833  0.04214496
#> 5   1.73955227  1.04632938 -2.1244807456 -0.26018450 -0.575303645  0.28200629
#> 6   2.40101072  1.04152301  0.1490534767 -1.43692714  0.430527948  0.07293507
#> 7   1.47299263  0.98088246  0.8909485733 -0.12320638 -0.112644504 -0.68748910
#> 8   3.65189428 -1.34080735  1.0195621716 -0.18171083  0.833926734  0.68497539
#> 9  -3.85810044 -0.79893795 -0.5297227000 -0.48534633 -0.015113992 -0.74881476
#> 10 -0.32625959 -1.99896649 -0.6896076783  0.23807582 -0.829430129  1.03202752
#> 11  2.32288113 -0.13471501  1.7152986424  0.77615806  0.462729597 -1.56262929
#> 12  1.43693184 -0.50581706 -1.0658360613 -0.73017377 -0.413635373  0.55337595
#> 13 -2.23131123  1.11093067 -0.2564729792 -1.35134454  0.332982277 -0.03771696
#> 14  2.06904285  2.03022952  0.1334168134 -0.83122203  1.131629143  1.60172454
#> 15 -1.00940775 -2.09519893  0.6252760231 -0.06247669  0.749362864  0.99176351
#> 16  0.75549234  0.01366052  1.0553934229  0.22434269  0.180822524  0.52727196
#> 17  1.17120823 -0.57828383  1.0732616691  1.19262510  0.991410531  0.19302265
#> 18 -1.16926742 -2.19515212  0.5938420953  0.92533747 -3.503834527  0.42561428
#> 19  0.71400790 -0.51636508 -0.7618743472 -0.68790299 -0.535352749 -0.37498532
#> 20  0.66977733  0.86489300  1.6526938008 -0.31658631  0.300355568  0.53129242
#> 21 -3.15653638  5.24745844 -0.7101244918  1.51018729  0.115313056 -0.13649751
#> 22 -2.55772946  3.90784192  0.3965305074  1.81267331 -0.678143059  1.13022869
#> 23  0.57826303  0.78898797 -3.1501664312  1.77212274  0.387939358  0.74058395
#> 24  1.91908109 -0.20008249  0.0696986694 -0.17269433 -0.089099525 -0.27334988
#> 25  1.19797823  0.20286494  0.5629103133  0.72892971 -0.374757875 -0.68179857
#> 26 -1.22109471  1.71524920  0.6440751859 -0.71813540 -0.332483201  0.09010989
#> 27 -1.01682652 -0.01570862 -0.2125766600 -1.06744533 -0.658894427  0.74542732
#> 28 -3.77920641 -1.03195936 -0.4335916309  0.46714341  0.419938379 -0.75983914
#> 29 -1.11452399  1.34439781  0.3318911686 -0.17200816 -0.206826726 -1.01389663
#> 30 -0.05873287  0.69593554  0.3356780518 -0.25660730  0.860836924  0.78852424
#> 31 -0.34496790 -0.19586614  0.5680489211 -0.68223887  0.008939626  0.34164596
#> 32  2.19698054 -0.19466645  1.0548372637  0.10057658  0.162867266 -0.14215674
#> 33  0.80652458  1.56043783  0.9096708340 -0.03121866 -0.486329281 -0.57429236
#> 34  1.13940601 -0.71320586 -0.1105783337  0.86872584  0.027719201 -1.12141646
#> 35 -0.35318806  0.97289851  1.3358545106 -0.11538835 -0.900059066  0.27063810
#> 36  1.36190064 -2.44887092  0.1229781253  2.36144592  1.069984258  0.44754357
#> 37 -3.58600353 -0.89269648 -1.3454499528  1.81886283  0.777286784 -0.71508794
#> 38 -2.60005266 -1.60224334 -0.0006642107 -1.81537960  0.067163191 -0.02390341
#> 39 -2.21807215 -1.83516590  0.0311274943  0.54999613  1.352077250 -0.84975348
#> 40 -3.56728777 -1.95501559 -0.5361500483 -1.42888784 -0.155764761  0.07911948
#> 41 -1.22433435  0.37359418 -1.0565102557 -1.61793584  1.326002021 -0.18906513
#> 42 -1.72504403  0.07419043 -0.5079317552 -1.27044157  0.038640929 -0.64927437
#> 43  4.64139040 -1.83810853 -2.4714698905  0.71849580 -0.147790038  0.24376350
#> 44 -0.61837209 -1.39772762  0.2161133027 -0.49068753 -0.148375580  0.58166401
#>             PC7          PC8          PC9         PC10        PC11         PC12
#> 1   0.012117983  0.236421545 -0.732919742 -0.084711871 -0.22610393 -0.095257111
#> 2   0.551614496  0.332221385 -0.490948264  0.392511458 -0.06738762 -0.112767371
#> 3  -0.316734949  0.007986329  0.327179857 -0.064118113  0.02166961  0.033720971
#> 4  -0.520082615  0.630615822 -0.140783257 -0.069422341  0.32765856  0.260842858
#> 5  -0.610990438 -0.053163022 -0.140483498  0.535396676  0.41391011  0.202238880
#> 6  -1.510745747  0.384439805  0.109622319  0.111577894 -0.57123147 -0.047214407
#> 7   0.077645805  0.993953000 -0.474404432  0.064816040 -0.17435874  0.212487724
#> 8   0.632755287  1.548845233  0.871243438  0.060502812 -0.32624830  0.004914267
#> 9  -0.147307939  0.611472391  0.350290517 -0.740857835 -0.41056797 -0.083587931
#> 10 -0.137454436 -0.531248127  0.177108736 -0.183597847 -0.21485374 -0.463019304
#> 11 -0.992746166  0.521392664 -0.008322405 -0.256275091  0.23489421  0.061772413
#> 12 -0.395628343 -0.465278754 -0.055430637 -1.079024373  0.20008366  0.103251900
#> 13  0.027419195 -0.461185254 -0.274399597 -0.125255987 -0.03503136 -0.058923152
#> 14  0.460807460 -0.321820840 -0.485397446  0.150778723 -0.18275363  0.351702172
#> 15 -1.608857391 -0.537192009  0.834675817 -0.140441565  0.06475569  0.198580310
#> 16  0.575321610  0.237319667  0.271420560 -0.283755030  0.28452951 -0.105363439
#> 17  0.801154972 -0.036439018 -0.214001536  0.255376924 -0.42570689 -0.327915319
#> 18 -0.042779684  0.066384639  0.350459435  0.126983737 -0.42661622  0.081365015
#> 19  0.647008708  0.193608026  0.141037884  0.100587585  0.26691401 -0.388404916
#> 20  0.175322894 -0.100857800  0.285232932  0.123908640 -0.27130871 -0.332064862
#> 21 -0.009528516 -0.269888388  1.016844122 -0.274309927 -0.28384145 -0.067435405
#> 22 -0.283355450 -0.186875810  0.113413452  0.179839413 -0.11626426  0.393277672
#> 23 -0.452248165  1.738169013 -0.264065314  0.104902503  0.33217872 -0.166833578
#> 24  0.446709572 -0.313627937 -0.017142734  0.295443779  0.25410103 -0.005488873
#> 25 -0.796035328 -0.613509374 -0.349219966 -0.171680757 -0.17087305 -0.041445735
#> 26  0.130245827 -0.364312420 -0.041753178 -0.312645521  0.67633690 -0.175314371
#> 27  0.336869785  1.163504596 -0.746206127  0.138773322  0.21510217  0.231571356
#> 28 -0.212140262 -0.199723445 -0.634096679 -0.003234272 -0.59795957 -0.063372530
#> 29  0.851470624 -0.534115492 -0.065248404 -0.332718042  0.48311562  0.262127426
#> 30  0.534008829 -0.507265427 -0.188117714 -0.468647511  0.06360698 -0.223664110
#> 31 -0.344784974 -0.762761837  0.454485151  0.762689180  0.28211137  0.091729659
#> 32  0.057104696  0.071026649  0.851614475 -0.407883635  0.53246947 -0.132709203
#> 33 -0.379317854 -0.077420592 -0.480411688  0.180826601 -0.40498217 -0.127484503
#> 34  1.390061776 -0.490018033  0.935895825  0.325176019 -0.18108126  0.435469514
#> 35  0.249374867 -0.287547191 -0.769988023  0.325057864  0.11969088 -0.322515954
#> 36  0.535325263 -0.353248400 -0.719451459 -0.283278722  0.09408529  0.233658315
#> 37  0.121582539  0.428379630  0.075721893 -0.211958889  0.22986636 -0.155627493
#> 38  0.402433863  0.019034508 -0.710962577 -0.726939057 -0.28919976  0.411774538
#> 39 -1.024558973 -0.692401188 -0.455343923  0.801455412  0.22388622  0.005126438
#> 40  0.228133745  0.697606134  0.630507762  0.449615761  0.05341087  0.180955425
#> 41  0.092146590 -0.025762576  0.613804100  0.257215678 -0.34524199  0.006401819
#> 42  0.626166714 -0.137449560  0.149281572  0.466551965  0.29753308 -0.134197368
#> 43  0.395243655 -1.338480975 -0.025294586  0.004848276 -0.34151263  0.064062553
#> 44 -0.572749524 -0.220787565 -0.075446663  0.005920122  0.39121441 -0.196424291

Pertanyaan Pak Salomo : kalau misalnya data ada variable yg MSA nya di Bawah 0.5 dan harus didrop, Ketika pembentukan PC diikutkan atau dibuang juga? Harus di drop langsung dari awal.

Sesuai dengan hasil dari factor analysis dengan objek fa_rentan_opt(), akan kita ambil sebanyak 2 principal component untuk kebutuhan pemetaan SoVI masing-masing kecamatan di Provinsi DKI Jakarta.

1
2

pc_keep <- as.data.frame(rentan_pca$x[,1:2])
pc_keep

#>            PC1         PC2
#> 1   3.48456033  1.56233715
#> 2  -2.61111513 -0.60474909
#> 3   4.69523879  0.59195276
#> 4  -0.07868070 -1.03628501
#> 5   1.73955227  1.04632938
#> 6   2.40101072  1.04152301
#> 7   1.47299263  0.98088246
#> 8   3.65189428 -1.34080735
#> 9  -3.85810044 -0.79893795
#> 10 -0.32625959 -1.99896649
#> 11  2.32288113 -0.13471501
#> 12  1.43693184 -0.50581706
#> 13 -2.23131123  1.11093067
#> 14  2.06904285  2.03022952
#> 15 -1.00940775 -2.09519893
#> 16  0.75549234  0.01366052
#> 17  1.17120823 -0.57828383
#> 18 -1.16926742 -2.19515212
#> 19  0.71400790 -0.51636508
#> 20  0.66977733  0.86489300
#> 21 -3.15653638  5.24745844
#> 22 -2.55772946  3.90784192
#> 23  0.57826303  0.78898797
#> 24  1.91908109 -0.20008249
#> 25  1.19797823  0.20286494
#> 26 -1.22109471  1.71524920
#> 27 -1.01682652 -0.01570862
#> 28 -3.77920641 -1.03195936
#> 29 -1.11452399  1.34439781
#> 30 -0.05873287  0.69593554
#> 31 -0.34496790 -0.19586614
#> 32  2.19698054 -0.19466645
#> 33  0.80652458  1.56043783
#> 34  1.13940601 -0.71320586
#> 35 -0.35318806  0.97289851
#> 36  1.36190064 -2.44887092
#> 37 -3.58600353 -0.89269648
#> 38 -2.60005266 -1.60224334
#> 39 -2.21807215 -1.83516590
#> 40 -3.56728777 -1.95501559
#> 41 -1.22433435  0.37359418
#> 42 -1.72504403  0.07419043
#> 43  4.64139040 -1.83810853
#> 44 -0.61837209 -1.39772762

Selanjutnya, melakukan perhitungan persamaan indeks kerentanan sosial. Sehingga dihasilkan SoVI atau Social Vulnerability Index untuk masing-masing kecamatan.

$$SoVI=|(0.65*komponen_1)+(0.35*komponen_2)|$$

abs(-5)

#> [1] 5

1
2
3

SoVI <- pc_keep %>% # object pc yang kita ambil sebanyak 2
  mutate(SoVI = 0.65*abs(PC1) + 0.35*abs(PC2))
SoVI

#>            PC1         PC2      SoVI
#> 1   3.48456033  1.56233715 2.8117822
#> 2  -2.61111513 -0.60474909 1.9088870
#> 3   4.69523879  0.59195276 3.2590887
#> 4  -0.07868070 -1.03628501 0.4138422
#> 5   1.73955227  1.04632938 1.4969243
#> 6   2.40101072  1.04152301 1.9251900
#> 7   1.47299263  0.98088246 1.3007541
#> 8   3.65189428 -1.34080735 2.8430139
#> 9  -3.85810044 -0.79893795 2.7873936
#> 10 -0.32625959 -1.99896649 0.9117070
#> 11  2.32288113 -0.13471501 1.5570230
#> 12  1.43693184 -0.50581706 1.1110417
#> 13 -2.23131123  1.11093067 1.8391780
#> 14  2.06904285  2.03022952 2.0554582
#> 15 -1.00940775 -2.09519893 1.3894347
#> 16  0.75549234  0.01366052 0.4958512
#> 17  1.17120823 -0.57828383 0.9636847
#> 18 -1.16926742 -2.19515212 1.5283271
#> 19  0.71400790 -0.51636508 0.6448329
#> 20  0.66977733  0.86489300 0.7380678
#> 21 -3.15653638  5.24745844 3.8883591
#> 22 -2.55772946  3.90784192 3.0302688
#> 23  0.57826303  0.78898797 0.6520168
#> 24  1.91908109 -0.20008249 1.3174316
#> 25  1.19797823  0.20286494 0.8496886
#> 26 -1.22109471  1.71524920 1.3940488
#> 27 -1.01682652 -0.01570862 0.6664353
#> 28 -3.77920641 -1.03195936 2.8176699
#> 29 -1.11452399  1.34439781 1.1949798
#> 30 -0.05873287  0.69593554 0.2817538
#> 31 -0.34496790 -0.19586614 0.2927823
#> 32  2.19698054 -0.19466645 1.4961706
#> 33  0.80652458  1.56043783 1.0703942
#> 34  1.13940601 -0.71320586 0.9902360
#> 35 -0.35318806  0.97289851 0.5700867
#> 36  1.36190064 -2.44887092 1.7423402
#> 37 -3.58600353 -0.89269648 2.6433461
#> 38 -2.60005266 -1.60224334 2.2508194
#> 39 -2.21807215 -1.83516590 2.0840550
#> 40 -3.56728777 -1.95501559 3.0029925
#> 41 -1.22433435  0.37359418 0.9265753
#> 42 -1.72504403  0.07419043 1.1472453
#> 43  4.64139040 -1.83810853 3.6602417
#> 44 -0.61837209 -1.39772762 0.8911465

Melakukan penggabungan variabel kecamatan, kode_shp atau kode area, dan juga nilai SoVI yang sudah didapatkan

Selanjutnya, dilakukan pengabungan nilai hasil perhitungan Social Vulnerability Index dengan nama-nama kecamatan dan kode area yang ada pada dataset rentan untuk kebutuhan visualisasi geospasial.

`1`	`glimpse(rentan)`

#> Rows: 44
#> Columns: 15
#> $ Kecamatan    <chr> "Cakung", "Cempaka Putih", "Cengkareng", "Cilandak", "Cil…
#> $ perempuan    <dbl> 49.2096, 50.0816, 49.3437, 50.4252, 49.4534, 49.6205, 49.…
#> $ KK_perempuan <dbl> 18.05886, 23.00326, 18.54968, 20.79268, 20.51346, 17.1236…
#> $ lansia       <dbl> 6.217495, 10.852327, 6.585898, 9.310706, 5.975218, 7.0006…
#> $ pengangguran <dbl> 2909, 544, 2839, 584, 2746, 2388, 1818, 2071, 702, 1187, …
#> $ wajib_ktp    <dbl> 70.20066, 75.12540, 70.84496, 73.73298, 74.35980, 71.1721…
#> $ banjir       <dbl> 31, 0, 55, 2, 38, 15, 8, 2, 0, 28, 12, 32, 2, 0, 4, 4, 0,…
#> $ informal     <dbl> 6.503218, 13.910844, 6.791649, 10.869546, 11.266610, 14.2…
#> $ RS           <dbl> 1, 2, 1, 3, 1, 0, 1, 2, 1, 3, 3, 2, 0, 1, 3, 2, 4, 2, 1, …
#> $ puskesmas    <dbl> 9, 4, 10, 6, 10, 11, 6, 12, 4, 10, 6, 12, 7, 13, 11, 9, 8…
#> $ positiv_cov  <dbl> 11740, 4998, 13036, 8201, 7899, 10437, 9529, 14329, 3796,…
#> $ rawat        <dbl> 82, 50, 109, 62, 73, 63, 73, 139, 25, 67, 80, 55, 25, 73,…
#> $ isoman       <dbl> 102, 17, 146, 114, 63, 124, 108, 145, 42, 43, 175, 87, 19…
#> $ Kode_shp     <chr> "38", "5", "15", "28", "12", "42", "41", "39", "1", "16",…
#> $ Kode         <chr> "IDN.7.4.1_1", "IDN.7.2.1_1", "IDN.7.1.1_1", "IDN.7.3.1_1…

dokumentasi

untuk mengambil keseluruhan nilai ~ (.)
mengambil nilai yang berlainan/berkebalikan menggunakan tanda !
cbind -> column combine

rentan_new <- rentan %>% 
  select_if(~!is.numeric(.)) %>% # ambil kolom Kecamatan dan kode_shp
  cbind(SoVI = SoVI[,3]) # gabungkan dengan kolom SoVI
  
rentan_new

#>                   Kecamatan Kode_shp         Kode      SoVI
#> 1                    Cakung       38  IDN.7.4.1_1 2.8117822
#> 2             Cempaka Putih        5  IDN.7.2.1_1 1.9088870
#> 3                Cengkareng       15  IDN.7.1.1_1 3.2590887
#> 4                  Cilandak       28  IDN.7.3.1_1 0.4138422
#> 5                 Cilincing       12  IDN.7.5.1_1 1.4969243
#> 6                  Cipayung       42  IDN.7.4.2_1 1.9251900
#> 7                   Ciracas       41  IDN.7.4.3_1 1.3007541
#> 8               Duren Sawit       39  IDN.7.4.4_1 2.8430139
#> 9                    Gambir        1  IDN.7.2.2_1 2.7873936
#> 10         Grogolpetamburan       16  IDN.7.1.2_1 0.9117070
#> 11                Jagakarsa       31  IDN.7.3.2_1 1.5570230
#> 12               Jatinegara       35  IDN.7.4.5_1 1.1110417
#> 13               Johar Baru        8  IDN.7.2.3_1 1.8391780
#> 14                Kalideres       20  IDN.7.1.3_1 2.0554582
#> 15           Kebayoran Baru       29  IDN.7.3.4_1 1.3894347
#> 16           Kebayoran Lama       27  IDN.7.3.5_1 0.4958512
#> 17               Kebonjeruk       19  IDN.7.1.4_1 0.9636847
#> 18            Kelapa Gading       14  IDN.7.5.3_1 1.5283271
#> 19                Kemayoran        3  IDN.7.2.4_1 0.6448329
#> 20                Kembangan       22  IDN.7.1.5_1 0.7380678
#> 21 Kepulauan Seribu Selatan       43  IDN.7.6.1_1 3.8883591
#> 22   Kepulauan Seribu Utara       44  IDN.7.6.2_1 3.0302688
#> 23                     Koja       11  IDN.7.5.4_1 0.6520168
#> 24               Kramatjati       36  IDN.7.4.6_1 1.3174316
#> 25                  Makasar       40  IDN.7.4.7_1 0.8496886
#> 26         Mampang Prapatan       25  IDN.7.3.6_1 1.3940488
#> 27                 Matraman       33  IDN.7.4.8_1 0.6664353
#> 28                  Menteng        6  IDN.7.2.5_1 2.8176699
#> 29               Pademangan       13  IDN.7.5.5_1 1.1949798
#> 30                 Palmerah       21  IDN.7.1.6_1 0.2817538
#> 31                 Pancoran       30  IDN.7.3.7_1 0.2927823
#> 32             Pasar Minggu       26  IDN.7.3.8_1 1.4961706
#> 33                Pasarrebo       37  IDN.7.4.9_1 1.0703942
#> 34              Penjaringan        9  IDN.7.5.6_1 0.9902360
#> 35             Pesanggrahan       32  IDN.7.3.9_1 0.5700867
#> 36               Pulogadung       34 IDN.7.4.10_1 1.7423402
#> 37              Sawah Besar        2  IDN.7.2.6_1 2.6433461
#> 38                    Senen        4  IDN.7.2.7_1 2.2508194
#> 39               Setia Budi       24 IDN.7.3.11_1 2.0840550
#> 40                Tamansari       17  IDN.7.1.7_1 3.0029925
#> 41                  Tambora       18  IDN.7.1.8_1 0.9265753
#> 42               Tanahabang        7  IDN.7.2.8_1 1.1472453
#> 43            Tanjung Priok       10  IDN.7.5.7_1 3.6602417
#> 44                    Tebet       23 IDN.7.3.12_1 0.8911465

Dari indeks kerentanan sosial yang terbentuk, selanjutnya perlu untuk melakukan kategorisasi pengkategorian indeks menjadi 3, yaitu rendah, sedang, tinggi. Hal ini bertujuan untuk memudahkan insight yang dapat diambil. Pengkategorian indeks dilakukan dengan menggunakan Natural Breaks pada R menggunakan function getJenksBreaks() dari library BAMMtools.

function(df,var) -> nilai input dr function yang kita buat

Rendah a-b Sedang b-c Tinggi c-d

# Membagi SoVI ke dalam 3 kategori
natural_breaks <- function(df, var) {
  breaks <- BAMMtools::getJenksBreaks(df[[var]], k = 4)
  df %>% mutate("jenks" := 
              case_when(df[[var]] >= breaks[[1]] & df[[var]] < breaks[[2]] ~ 1,
                        df[[var]] >= breaks[[2]] & df[[var]] < breaks[[3]] ~ 2,
                        df[[var]] >= breaks[[3]] & df[[var]] <= breaks[[4]] ~ 3))
}
rentan_new <- natural_breaks(rentan_new, "SoVI")
rentan_new

#>                   Kecamatan Kode_shp         Kode      SoVI jenks
#> 1                    Cakung       38  IDN.7.4.1_1 2.8117822     3
#> 2             Cempaka Putih        5  IDN.7.2.1_1 1.9088870     2
#> 3                Cengkareng       15  IDN.7.1.1_1 3.2590887     3
#> 4                  Cilandak       28  IDN.7.3.1_1 0.4138422     1
#> 5                 Cilincing       12  IDN.7.5.1_1 1.4969243     2
#> 6                  Cipayung       42  IDN.7.4.2_1 1.9251900     2
#> 7                   Ciracas       41  IDN.7.4.3_1 1.3007541     2
#> 8               Duren Sawit       39  IDN.7.4.4_1 2.8430139     3
#> 9                    Gambir        1  IDN.7.2.2_1 2.7873936     3
#> 10         Grogolpetamburan       16  IDN.7.1.2_1 0.9117070     1
#> 11                Jagakarsa       31  IDN.7.3.2_1 1.5570230     2
#> 12               Jatinegara       35  IDN.7.4.5_1 1.1110417     1
#> 13               Johar Baru        8  IDN.7.2.3_1 1.8391780     2
#> 14                Kalideres       20  IDN.7.1.3_1 2.0554582     2
#> 15           Kebayoran Baru       29  IDN.7.3.4_1 1.3894347     2
#> 16           Kebayoran Lama       27  IDN.7.3.5_1 0.4958512     1
#> 17               Kebonjeruk       19  IDN.7.1.4_1 0.9636847     1
#> 18            Kelapa Gading       14  IDN.7.5.3_1 1.5283271     2
#> 19                Kemayoran        3  IDN.7.2.4_1 0.6448329     1
#> 20                Kembangan       22  IDN.7.1.5_1 0.7380678     1
#> 21 Kepulauan Seribu Selatan       43  IDN.7.6.1_1 3.8883591     3
#> 22   Kepulauan Seribu Utara       44  IDN.7.6.2_1 3.0302688     3
#> 23                     Koja       11  IDN.7.5.4_1 0.6520168     1
#> 24               Kramatjati       36  IDN.7.4.6_1 1.3174316     2
#> 25                  Makasar       40  IDN.7.4.7_1 0.8496886     1
#> 26         Mampang Prapatan       25  IDN.7.3.6_1 1.3940488     2
#> 27                 Matraman       33  IDN.7.4.8_1 0.6664353     1
#> 28                  Menteng        6  IDN.7.2.5_1 2.8176699     3
#> 29               Pademangan       13  IDN.7.5.5_1 1.1949798     2
#> 30                 Palmerah       21  IDN.7.1.6_1 0.2817538     1
#> 31                 Pancoran       30  IDN.7.3.7_1 0.2927823     1
#> 32             Pasar Minggu       26  IDN.7.3.8_1 1.4961706     2
#> 33                Pasarrebo       37  IDN.7.4.9_1 1.0703942     1
#> 34              Penjaringan        9  IDN.7.5.6_1 0.9902360     1
#> 35             Pesanggrahan       32  IDN.7.3.9_1 0.5700867     1
#> 36               Pulogadung       34 IDN.7.4.10_1 1.7423402     2
#> 37              Sawah Besar        2  IDN.7.2.6_1 2.6433461     3
#> 38                    Senen        4  IDN.7.2.7_1 2.2508194     3
#> 39               Setia Budi       24 IDN.7.3.11_1 2.0840550     2
#> 40                Tamansari       17  IDN.7.1.7_1 3.0029925     3
#> 41                  Tambora       18  IDN.7.1.8_1 0.9265753     1
#> 42               Tanahabang        7  IDN.7.2.8_1 1.1472453     1
#> 43            Tanjung Priok       10  IDN.7.5.7_1 3.6602417     3
#> 44                    Tebet       23 IDN.7.3.12_1 0.8911465     1

1
2

# Melihat threshold untuk masing-masing kategori 
getJenksBreaks(rentan_new$SoVI, k = 4)

#> [1] 0.2817538 1.1949798 2.2508194 3.8883591

Insight :

Nilai rendah(1) = nilai SoVI di antara nilai 0.2817538 - 1.1949798
Nilai sedang(2) = nilai SoVI di antara 1.1949798 - 2.2508194
Nilai tinggi(3) = nilai SoVI di antara 2.2508194 - 3.8883591

# Melihat jumlah kecamatan per masing-masing kategori
rentan_new %>% 
  group_by(jenks) %>% 
  summarise(Total_Kecamatan = n())

#> # A tibble: 3 × 2
#>   jenks Total_Kecamatan
#>   <dbl>           <int>
#> 1     1              18
#> 2     2              15
#> 3     3              11

Insight:

Di DKI Jakarta tingkat kerentanan sosial akibat COVID 19 lebih banyak berada di level rendah-sedang
Untuk kategori tinggi atau daerah kecamatan yang membutuhkan penanganan lebih berjumlah 11 kecamatan. Di antaranya kecamatan-kecamatan yang ada di bawah ini:

1
2
3

# insight bu Delvia
rentan_new %>% 
  filter(jenks == 3) # untuk nilai SoVI kategori tinggi

#>                   Kecamatan Kode_shp        Kode     SoVI jenks
#> 1                    Cakung       38 IDN.7.4.1_1 2.811782     3
#> 2                Cengkareng       15 IDN.7.1.1_1 3.259089     3
#> 3               Duren Sawit       39 IDN.7.4.4_1 2.843014     3
#> 4                    Gambir        1 IDN.7.2.2_1 2.787394     3
#> 5  Kepulauan Seribu Selatan       43 IDN.7.6.1_1 3.888359     3
#> 6    Kepulauan Seribu Utara       44 IDN.7.6.2_1 3.030269     3
#> 7                   Menteng        6 IDN.7.2.5_1 2.817670     3
#> 8               Sawah Besar        2 IDN.7.2.6_1 2.643346     3
#> 9                     Senen        4 IDN.7.2.7_1 2.250819     3
#> 10                Tamansari       17 IDN.7.1.7_1 3.002993     3
#> 11            Tanjung Priok       10 IDN.7.5.7_1 3.660242     3

Visualisasi Geospasial Kerentanan Sosial Kecamatan di Provinsi DKI Jakarta

Pemetaan dalam bentuk visual dilakukan untuk melihat persebaran indeks kerentanan sosial menurut kecamatan di DKI Jakarta. Visualisasi yang akan dilakukan dinamakan visualisasi geospasial, kita bisa memberikan informasi yang hanya dengan bantuan peta maka informasi tersebut dapat tersampaikan dengan baik.

Pada analisis ini akan dibuat visualisasi geocoding berupa choropleth.

Mengambil data spasial

Pertama-tama dilakukan penarikan data shapefile yang mengandung data geometri untuk kecamatan di DKI Jakarta melalui web https://hub.arcgis.com/datasets/18d008ec01aa41ab928c1959c451c6e7/explore?layer=0&location=-6.257734%2C106.896911%2C12.75

Membaca data geometri spasial

Untuk membaca data geometri spasial kita memerlukan sebuah package sf, dalam mempersiapkan package tersebut kita akan menjalankan code di bawah ini.

1
2
3

# Mendownload package sf dengan versi 1.0-8
# sf <- "http://cran.r-project.org/src/contrib/Archive/sf/sf_1.0-8.tar.gz"
# install.packages(sf, repos=NULL, type="source")

Langkah kedua, untuk membaca data .shp di R menggunakan function read_sf() pada packagesf. Class data sf itu merujuk kepada simple feature -> geometri file

Combining Data

Dalam melakukan penggabungan data, ada beberapa hal yang perlu dilakukan:

Mengubah tipe data OBJECTID menjadi character agar dapat dilakukan penggabungan data.

1
2
3

jakarta_rentan_new <- shapefile %>% 
  mutate(OBJECTID = as.character(OBJECTID))
glimpse(jakarta_rentan_new)

#> Rows: 44
#> Columns: 11
#> $ OBJECTID_1 <dbl> 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, …
#> $ OBJECTID   <chr> "43", "44", "1", "2", "3", "4", "5", "6", "7", "8", "9", "1…
#> $ NO_KAB     <chr> NA, NA, "71", "71", "71", "71", "71", "71", "71", "71", "72…
#> $ NO_KEC     <chr> NA, NA, "01", "02", "03", "04", "05", "06", "07", "08", "01…
#> $ KECAMATAN  <chr> "KEP. SERIBU SELATAN", "KEP. SERIBU UTARA", "GAMBIR", "SAWA…
#> $ WILAYAH    <chr> "KAB.ADM.KEP.SERIBU", "KAB.ADM.KEP.SERIBU", "JAKARTA PUSAT"…
#> $ Shape__Are <dbl> 1528586439, 3126051321, 7508915, 5358573, 7168174, 4350701,…
#> $ Shape__Len <dbl> 175577.253, 229963.261, 13230.033, 13639.388, 13630.925, 13…
#> $ Shape_Leng <dbl> 0.46211226, 0.85004576, 0.11956374, 0.12328821, 0.12317691,…
#> $ Shape_Area <dbl> 0.0003136000, 0.0005370378, 0.0006133681, 0.0004376953, 0.0…
#> $ geometry   <MULTIPOLYGON [°]> MULTIPOLYGON (((106.7359 -6..., MULTIPOLYGON (…

Menggabungkan data spasial yang dimiliki pada object jakarta_rentan_new dimana variabel OBJECTID akan dipasangkan dengan dataset rentan pada variabel Kode_shp
Melakukan reorder baris pada dataframe hasil gabungan menggunakan column name KECAMATAN

1
2
3

jakarta_rentan_clean <- jakarta_rentan_new %>% 
  left_join(rentan_new, by = c("OBJECTID" = "Kode_shp")) %>% 
  arrange(KECAMATAN) ##mengurutkan datanya berdasarkan variabel kecamatan

Mengubah object spasial geometri 3 dimensi yang kita miliki menjadi 2 dimensi

Dikarenakan output yang ingin ditampilkan adalah visualisasi geometri map 2 dimensi menggunakan leaflet, maka kita perlu mengubah dimensi object spasial geometri kita menjadi 2 dimensi. Di R, hal ini dapat dilakukan dengan bantuan function st_zm() dengan parameter drop = TRUE.

`1`	`jakarta_rentan_clean <- st_zm(jakarta_rentan_clean, drop = T, what = "ZM")`

Membuat Choropleth Map menggunakan `leaflet()`

Dalam penyajian data, pertama-tama dimasukkan data rentan dalam leaflet. Lalu, akan ditentukan pewarnaan yang sesuai untuk choropleth yang nanti akan dibentuk. Palette ditentukan dari kuning hingga merah (silahkan lihat cheatsheet yang sudah disediakan). Selain itu, dibuat juga popup pada tiap tempat. Referensi untuk melakukan pemilihan warna palette: https://www.datanovia.com/en/blog/the-a-z-of-rcolorbrewer-palette/

library(leaflet)
m <- leaflet(jakarta_rentan_clean)

# popup
popup.cont  <- paste("Region: ", jakarta_rentan_clean$KECAMATAN, "<br>",
                      "Value: ", jakarta_rentan_clean$SoVI, "<br>")

# mendefinisikan variabel yang menjadi dasar pewarnaan
col <- jakarta_rentan_clean$SoVI

# kategorisasi variabel
bins <- c(getJenksBreaks(jakarta_rentan_clean$SoVI, 4))

# melakukan pewarnaan
pal <- colorBin("YlOrRd", 
                domain = col, #variabel yang ingin dipetakan
                bins = bins )

Terakhir, kita akan membentuk peta dengan tiles yang disediakan oleh CartoDB menggunakan fungsi addProviderTiles. Lalu kita atur poligon dan legenda pada peta, serta mengubah nama tiap elemen pada legenda agar lebih masuk akal.

m %>% 
  addProviderTiles(providers$CartoDB.DarkMatter) %>% 
  addPolygons(fillColor = pal(col),
              weight = 1,
              opacity = 1,
              color = "black",
              dashArray = "3",
              fillOpacity = 0.8,
              popup = popup.cont) %>%
  addLegend("bottomright", 
            pal = pal,
            values = ~col,
            title = "Social Vulnerability Index",
            labFormat = labelFormat(digits = 2),
            opacity = 1)

Dive deeper Insight yang dapat ditarik dari pemetaan geospasial:

Pak Joe > Dilihat dari pemetaannya, hanya sedikit wilayah yang tingkat kerentanannya tinggi dan sebagian besar memiliki tingkat kerentanan rendah. Dilihat dari penyebarannya, daerah yang berdekatan dengan Bekasi memiliki tingkat kerentanan dari menengah hingga tinggi. Kecamatan Kep Seribu memiliki nilai tingkat kerentanan tertinggi sebesar 3.88. Luas wilayah tidak mempengaruhi tingkat kerentanan, perlu diketahui var. Apa yang memiliki pengaruh tertinggi yang digabungkan dengan masing-masing wilayah.
Pak Calvin

Semakin rendah jumlah titik banjirnya, SoVI-nya cenderung tinggi atau kerentanan sosial terhadap COVID nya semakin tinggi. Walau ini terlihat tersebar dengan cukup merata, namun titik banjir = 0 lebih banyak berada di Jenks = 3 (grup SoVI tinggi), dan sedikit sekali di Jenks = 1.
Angka informal dan positif_cov yang cukup paling tinggi lebih mungkin memiliki kerentanan / SoVI tinggi, kebanyakan dari kedua variabel ini yang tinggi berada di SoVI jenks 2 dan 3.
Selain itu, variabel-variabel lainnya tersebar dengan tidak terlalu menimbulkan pola yang terlihat jelas. Misalnya lansia, rawat, isoman, dan KK_perempuan yang sempat waktu FA mendapat angka korelasi tinggi, ternyata angkanya terlihat sangat beragam dan menyebar di sini, tidak terlalu mengumpul di grup Jenks tertentu.

Fuzzy C-Means Clustering terhadap Kerentanan Sosial

Setelah melakukan pembentukan indeks kerentanan sosial, selanjutnya clustering digunakan untuk mengetahui variabel yang berpengaruh pada masing-masing klaster. Supaya hasil klaster dapat diinterpretasikan, akan digunakan data rentan awal dengan variabel-variabel numerik yang mana sudah signifikan melalui pengujian KMO-MSA dan bartlett. Metode yang digunakan dalam melakukan klastering ini adalah Fuzzy C-Means Clustering.

Fuzzy C-means Methodology

Fuzzy C-means clustering merupakan suatu metode clustering yang hampir mirip seperti k-means clustering. Karena metode clustering ini mirip dengan k-means clustering, ada yang menyebut metode ini fuzzy k-means clustering.

Algoritma Fuzzy C-Means (FCM) pertama kali diperkenalkan oleh Dunn pada tahun 1973 dan dikembangkan kembali oleh Bezdek pada tahun 1981 sebagai metode clustering yang menggunakan konsep fuzzy.

Langkah pertama adalah menentukan banyaknya jumlah cluster pada data.
Dilanjutkan dengan mencari pusat cluster (centroid), yang akan menandai posisi rata-rata untuk masing-masing cluster.
Memperbarui setiap iterasi
Terakhir data dikelompokkan berdasarkan jarak antar data dengan setiap centroid.

Fuzzy c-means merupakan salah satu jenis soft clustering di mana dalam mengelompokan suatu data, setiap data bisa dimiliki lebih dari satu cluster. Contohnya sebuah tomat bisa dikelompokkan ke warna merah atau warna hijau jika dalam hard clustering, namun tomat bisa dikelompokkan ke warna merah dan hijau dalam fuzzy clustering. Tomat berwarna merah memiliki tingkatan yang sama dengan tomat berwarna hijau. Jika dimisalkan dengan angka dari 0 sampai 1, tomat warna merah 0.5 dan tomat warna hijau 0.5.

Kami melakukan pengolahan klaster menggunakan perangkat lunak R-Studio dengan package ppclust.

Data pre-processing

Langkah pertama mendefinisikan data yang akan digunakan dalam analisis FCM ini ke dalam object rentan_FCM dimana data ini berisi variabel-variabel bertipe data numerik.

1
2

rentan_FCM <- rentan[,1:13]
head(rentan_FCM)

#> # A tibble: 6 × 13
#>   Kecamatan  perem…¹ KK_pe…² lansia penga…³ wajib…⁴ banjir infor…⁵    RS puske…⁶
#>   <chr>        <dbl>   <dbl>  <dbl>   <dbl>   <dbl>  <dbl>   <dbl> <dbl>   <dbl>
#> 1 Cakung        49.2    18.1   6.22    2909    70.2     31    6.50     1       9
#> 2 Cempaka P…    50.1    23.0  10.9      544    75.1      0   13.9      2       4
#> 3 Cengkareng    49.3    18.5   6.59    2839    70.8     55    6.79     1      10
#> 4 Cilandak      50.4    20.8   9.31     584    73.7      2   10.9      3       6
#> 5 Cilincing     49.5    20.5   5.98    2746    74.4     38   11.3      1      10
#> 6 Cipayung      49.6    17.1   7.00    2388    71.2     15   14.3      0      11
#> # … with 3 more variables: positiv_cov <dbl>, rawat <dbl>, isoman <dbl>, and
#> #   abbreviated variable names ¹perempuan, ²KK_perempuan, ³pengangguran,
#> #   ⁴wajib_ktp, ⁵informal, ⁶puskesmas

Melakukan assign nilai dari kolom Kecamatan yang ada pada data rentan kita menjadi nama baris.

rentan_FCM <- rentan_FCM %>% 
  remove_rownames %>% 
  column_to_rownames(var="Kecamatan")
head(rentan_FCM)

#>               perempuan KK_perempuan    lansia pengangguran wajib_ktp banjir
#> Cakung          49.2096     18.05886  6.217495         2909  70.20066     31
#> Cempaka Putih   50.0816     23.00326 10.852327          544  75.12540      0
#> Cengkareng      49.3437     18.54968  6.585898         2839  70.84496     55
#> Cilandak        50.4252     20.79268  9.310706          584  73.73298      2
#> Cilincing       49.4534     20.51346  5.975218         2746  74.35980     38
#> Cipayung        49.6205     17.12363  7.000645         2388  71.17214     15
#>                informal RS puskesmas positiv_cov rawat isoman
#> Cakung         6.503218  1         9       11740    82    102
#> Cempaka Putih 13.910844  2         4        4998    50     17
#> Cengkareng     6.791649  1        10       13036   109    146
#> Cilandak      10.869546  3         6        8201    62    114
#> Cilincing     11.266610  1        10        7899    73     63
#> Cipayung      14.258273  0        11       10437    63    124

Normalisasi Min-Max

Normalisasi/standarisasi data diperlukan sebelum melakukan clustering untuk memudahkan kita melakukan identifikasi cluster. Apabila skala data untuk dilakukan clustering tidak sama, maka kita tidak akan dapat membandingkan setiap nilai dengan semua nilai lain dalam kumpulan data yang kita miliki.

Teknik penskalaan Min-Max merupakan salah satu metode normaliasi dengan melakukan perubahan linear terhadap data asli. Teknik normalisasi Min-Max dapat dilakukan dengan formula berikut:

$$ a’_i = \frac{a_i - min(a)}{max(a)-min(a)}$$ Selanjutnya, karena pada R tidak terdapat build-in function dari formula Min-Max ini, maka akan dibuat fungsi manual dengan minMax. Lalu melakukan pengaplikasian fungsi yang sudah dibuat pada data rentan_FCM kita.

# custom function to implement min max scaling
minMax <- function(x) {
  (x - min(x)) / (max(x) - min(x))
}
 
#Normalisasi data menggunakan fungsi custom 
rentan_scale <- as.data.frame(lapply(rentan_FCM, minMax))

#Menambahkan rownames kecamatan pada data yang sudah discale
rentan_FCM <- as.data.frame(rentan_scale, rentan$Kecamatan)

Pemilihan nilai k optimum

Pada analisis ini, akan digunakan metode silhouette untuk melakukan pemilihan nilai k (jumlah cluster) optimal.

Silhouette mengukur seberapa baik suatu observasi menjadi satu cluster atau memperkirakan nilai rata-rata jarak antar cluster. Nilai dari silhouette coefficient menginterpretasikan seberapa dekat suatu data pada satu cluster terhadap tetangga pada cluster yang lain. Semakin besar nilai silhouette (mendekati 1), maka semakin baik cluster tersebut.

Di mana akan dipilih nilai k pada saat nilai k ditambah, penurunan Total WSS (Whitin Sum of Squares) tidak terlalu drastis (dapat dikatakan sudah melandai). Secara statistik, dalam melakukan visualisasi elbow digunakan function fviz_nbclust() dari package factoextra.

library(cluster)
set.seed(123)
 
elbow <- 
fviz_nbclust(
  x = rentan_scale, #data
  FUNcluster = fanny, #fungsi clustering
  method = "silhouette") + labs(subtitle = "Metode Silhouette") #judul
elbow

Nilai k optimum dari data rentan_FCM didapatkan pada saat 2.

FCM menggunakan `fcm()`

Membuat object fcm.rentan dengan menggunakan fungsi fcm() dengan jumlah k sesuai banyaknya nilai k optimum yang didapat.

Parameter:

x: dataset
centers: banyaknya centroid (nilai k)

library(ppclust)
set.seed(123)
fcm.rentan <- fcm(x = rentan_FCM,
                  centers = 2)
summary(fcm.rentan)

#> Summary for 'fcm.rentan'
#> 
#> Number of data objects:  44 
#> 
#> Number of clusters:  2 
#> 
#> Crisp clustering vector:
#>  [1] 2 1 2 2 2 2 2 2 1 1 2 2 1 2 1 2 2 1 2 2 1 1 2 2 2 1 1 1 1 2 1 2 2 2 2 2 1 1
#> [39] 1 1 1 1 2 1
#> 
#> Initial cluster prototypes:
#>           perempuan KK_perempuan    lansia pengangguran  wajib_ktp    banjir
#> Cluster 1 0.4207626    0.2174588 0.0000000    0.5617909 0.04974973 0.0000000
#> Cluster 2 0.6848996    0.3799949 0.4628804    0.5716038 0.21393593 0.0338983
#>             informal        RS puskesmas positiv_cov     rawat    isoman
#> Cluster 1 0.01501601 0.1428571 0.9090909    0.553489 0.5140845 0.2857143
#> Cluster 2 0.01414969 0.2857143 0.8181818    1.000000 0.9788732 0.8285714
#> 
#> Final cluster prototypes:
#>           perempuan KK_perempuan    lansia pengangguran wajib_ktp    banjir
#> Cluster 1 0.5389905    0.6162944 0.5619571    0.2209669 0.4455888 0.1271250
#> Cluster 2 0.5210274    0.3989750 0.3105454    0.4079368 0.2562638 0.2963236
#>            informal        RS puskesmas positiv_cov     rawat    isoman
#> Cluster 1 0.4673842 0.2265458 0.3294159   0.3923736 0.3665275 0.2634173
#> Cluster 2 0.1986334 0.2812407 0.5169692   0.6422221 0.5513181 0.5059506
#> 
#> Distance between the final cluster prototypes
#>           Cluster 1
#> Cluster 2 0.4759771
#> 
#> Difference between the initial and final cluster prototypes
#>            perempuan KK_perempuan     lansia pengangguran  wajib_ktp    banjir
#> Cluster 1  0.1182279   0.39883559  0.5619571   -0.3408239 0.39583911 0.1271250
#> Cluster 2 -0.1638722   0.01898006 -0.1523350   -0.1636670 0.04232784 0.2624253
#>            informal           RS  puskesmas positiv_cov      rawat      isoman
#> Cluster 1 0.4523682  0.083688701 -0.5796750  -0.1611154 -0.1475570 -0.02229703
#> Cluster 2 0.1844837 -0.004473622 -0.3012126  -0.3577779 -0.4275551 -0.32262080
#> 
#> Root Mean Squared Deviations (RMSD): 1.015206 
#> Mean Absolute Deviation (MAD): 34.74745 
#> 
#> Membership degrees matrix (top and bottom 5 rows): 
#>               Cluster 1 Cluster 2
#> Cakung        0.2248502 0.7751498
#> Cempaka Putih 0.8585258 0.1414742
#> Cengkareng    0.2755678 0.7244322
#> Cilandak      0.4970376 0.5029624
#> Cilincing     0.2914759 0.7085241
#> ...
#>               Cluster 1 Cluster 2
#> Tamansari     0.7645847 0.2354153
#> Tambora       0.7258928 0.2741072
#> Tanahabang    0.8497685 0.1502315
#> Tanjung Priok 0.3348927 0.6651073
#> Tebet         0.7364816 0.2635184
#> 
#> Descriptive statistics for the membership degrees by clusters
#>           Size       Min        Q1      Mean    Median        Q3       Max
#> Cluster 1   20 0.5704471 0.6287925 0.7072370 0.7270529 0.7704006 0.8585258
#> Cluster 2   24 0.5029624 0.6715673 0.7042262 0.7164671 0.7553940 0.9221182
#> 
#> Dunn's Fuzziness Coefficients:
#> dunn_coeff normalized 
#>  0.6029278  0.2058555 
#> 
#> Within cluster sum of squares by cluster:
#>         1         2 
#>  9.563466 10.156608 
#> (between_SS / total_SS =  20.84%) 
#> 
#> Available components: 
#>  [1] "u"          "v"          "v0"         "d"          "x"         
#>  [6] "cluster"    "csize"      "sumsqrs"    "k"          "m"         
#> [11] "iter"       "best.start" "func.val"   "comp.time"  "inpargs"   
#> [16] "algorithm"  "call"

1
2

#Untuk menghitung jarak masing-masing individu terhadap pusat cluster
#fcm.rentan$d

Hasil dari fcm() :

Untuk mengetahui jumlah kecamatan untuk setiap cluster

`1`	`fcm.rentan$csize`

#>  1  2 
#> 20 24

Insight

Untuk cluster 1 ada sejumlah 20 kecamatan
Untuk cluster 2 ada sejumlah 24 kecamatan

Nilai untuk masing-masing kecamatan dan juga hasil cluster setiap kecamatannya.

`1`	`fcm.rentan$cluster`

#>  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 
#>  2  1  2  2  2  2  2  2  1  1  2  2  1  2  1  2  2  1  2  2  1  1  2  2  2  1 
#> 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 
#>  1  1  1  2  1  2  2  2  2  2  1  1  1  1  1  1  2  1

1
2
3

# membuat data frame gabungan antara rentan_FCM dengan hasil cluster
fcm_keep<-data.frame(rentan_FCM, fcm.rentan$cluster)
fcm_keep

#>                          perempuan KK_perempuan     lansia pengangguran
#> Cakung                   0.3464995  0.257657175 0.04693080   0.82858019
#> Cempaka Putih            0.6689966  0.646647181 0.67779574   0.10334253
#> Cengkareng               0.3960945  0.296270752 0.09707556   0.80711438
#> Cilandak                 0.7960723  0.472734676 0.46795975   0.11560871
#> Cilincing                0.4366656  0.450767265 0.01395345   0.77859552
#> Cipayung                 0.4984652  0.184079191 0.15352836   0.66881325
#> Ciracas                  0.5169200  0.307120134 0.25610534   0.49402024
#> Duren Sawit              0.6848996  0.379994918 0.46288044   0.57160380
#> Gambir                   0.4643663  0.810267949 1.00000000   0.15179393
#> Grogolpetamburan         0.7379711  0.640175890 0.74680051   0.30052131
#> Jagakarsa                0.5563815  0.286992779 0.19128910   0.30236124
#> Jatinegara               0.5151078  0.608250553 0.46869349   0.56117755
#> Johar Baru               0.3827434  0.602835796 0.48355587   0.19748543
#> Kalideres                0.4207626  0.217458793 0.00000000   0.56179086
#> Kebayoran Baru           0.7663745  0.606061960 0.62233733   0.08065011
#> Kebayoran Lama           0.5842672  0.444945652 0.40121807   0.17877952
#> Kebonjeruk               0.5802729  0.332212572 0.40789323   0.24471021
#> Kelapa Gading            1.0000000  0.536966368 0.97145883   0.07758356
#> Kemayoran                0.5076371  0.666277263 0.49483160   0.43207605
#> Kembangan                0.5665520  0.227503336 0.33098498   0.11039558
#> Kepulauan Seribu Selatan 0.0000000  0.000000000 0.07562168   0.00000000
#> Kepulauan Seribu Utara   0.3973520  0.009924764 0.06598185   0.02483901
#> Koja                     0.3856282  0.512254638 0.11573338   1.00000000
#> Kramatjati               0.5321203  0.479322129 0.27237959   0.38883778
#> Makasar                  0.5462850  0.278624962 0.27694544   0.30481447
#> Mampang Prapatan         0.4350383  0.491608209 0.28165521   0.04507820
#> Matraman                 0.6678501  0.699908247 0.52963631   0.50843300
#> Menteng                  0.5179186  0.727487357 0.83416971   0.15915363
#> Pademangan               0.3302637  0.561075943 0.37526201   0.05826434
#> Palmerah                 0.4393654  0.462745200 0.36869987   0.22109782
#> Pancoran                 0.6250971  0.503287615 0.33765364   0.09107636
#> Pasar Minggu             0.5450276  0.409313536 0.29890696   0.19871205
#> Pasarrebo                0.4938422  0.170310564 0.21261087   0.32505366
#> Penjaringan              0.4237213  0.504878996 0.51936571   0.19932536
#> Pesanggrahan             0.6430711  0.348132966 0.30645587   0.11683533
#> Pulogadung               0.6847147  0.600104789 0.48299608   0.38301135
#> Sawah Besar              0.4082252  0.816590504 0.69589161   0.18307268
#> Senen                    0.6108584  0.948075957 0.99188382   0.27292242
#> Setia Budi               0.6088613  0.751018436 0.48484996   0.12020852
#> Tamansari                0.6705499  1.000000000 0.97311418   0.15823367
#> Tambora                  0.3099597  0.639117672 0.52079335   0.38393131
#> Tanahabang               0.4301564  0.758500140 0.53646375   0.19595216
#> Tanjung Priok            0.4987241  0.492663228 0.31562520   0.87580497
#> Tebet                    0.7221051  0.689762495 0.54628073   0.22293775
#>                           wajib_ktp     banjir   informal        RS  puskesmas
#> Cakung                   0.00000000 0.52542373 0.01814115 0.1428571 0.54545455
#> Cempaka Putih            0.51684443 0.00000000 0.39627308 0.2857143 0.09090909
#> Cengkareng               0.06761805 0.93220339 0.03286447 0.1428571 0.63636364
#> Cilandak                 0.37071190 0.03389831 0.24102604 0.4285714 0.27272727
#> Cilincing                0.43649576 0.64406780 0.26129471 0.1428571 0.63636364
#> Cipayung                 0.10195493 0.25423729 0.41400807 0.0000000 0.72727273
#> Ciracas                  0.14525009 0.13559322 0.05311726 0.1428571 0.27272727
#> Duren Sawit              0.21393593 0.03389831 0.01414969 0.2857143 0.81818182
#> Gambir                   0.66548705 0.00000000 0.70812448 0.1428571 0.09090909
#> Grogolpetamburan         0.51810251 0.47457627 0.45608100 0.4285714 0.63636364
#> Jagakarsa                0.11974172 0.20338983 0.23868687 0.4285714 0.27272727
#> Jatinegara               0.26977340 0.54237288 0.25501448 0.2857143 0.81818182
#> Johar Baru               0.29897079 0.03389831 0.53446759 0.0000000 0.36363636
#> Kalideres                0.04974972 0.00000000 0.01501601 0.1428571 0.90909091
#> Kebayoran Baru           0.50283018 0.06779661 0.80261564 0.4285714 0.72727273
#> Kebayoran Lama           0.29749273 0.06779661 0.06371145 0.2857143 0.54545455
#> Kebonjeruk               0.23298989 0.00000000 0.17592762 0.5714286 0.45454545
#> Kelapa Gading            0.60690064 0.72881356 0.19133111 0.2857143 0.18181818
#> Kemayoran                0.31480828 0.44067797 0.27119979 0.1428571 0.45454545
#> Kembangan                0.16474356 0.00000000 0.17062707 0.1428571 0.54545455
#> Kepulauan Seribu Selatan 0.89030445 0.00000000 0.12321040 0.0000000 0.09090909
#> Kepulauan Seribu Utara   0.77275241 0.00000000 0.00000000 0.1428571 0.18181818
#> Koja                     1.00000000 0.28813559 0.19524772 0.4285714 0.45454545
#> Kramatjati               0.13909914 0.40677966 0.17963154 0.2857143 0.54545455
#> Makasar                  0.19855052 0.44067797 0.28931734 0.4285714 0.36363636
#> Mampang Prapatan         0.29789863 0.11864407 0.20563754 0.0000000 0.36363636
#> Matraman                 0.38994109 0.06779661 0.22649994 0.0000000 0.36363636
#> Menteng                  0.55560765 0.00000000 0.80798116 0.4285714 0.00000000
#> Pademangan               0.24567457 0.20338983 0.17419871 0.1428571 0.18181818
#> Palmerah                 0.24155332 0.00000000 0.21952901 0.2857143 0.63636364
#> Pancoran                 0.27205851 0.22033898 0.51407257 0.1428571 0.54545455
#> Pasar Minggu             0.22679113 0.27118644 0.09414903 0.2857143 0.63636364
#> Pasarrebo                0.13559586 0.27118644 0.17693053 0.1428571 0.27272727
#> Penjaringan              0.30135236 0.42372881 0.17655424 0.4285714 0.36363636
#> Pesanggrahan             0.16876389 0.18644068 0.12285780 0.1428571 0.27272727
#> Pulogadung               0.36422494 0.16949153 0.22195859 1.0000000 0.54545455
#> Sawah Besar              0.92895910 0.00000000 0.64538752 0.5714286 0.00000000
#> Senen                    0.20439747 0.00000000 0.61312886 0.1428571 0.36363636
#> Setia Budi               0.42413820 0.01694915 1.00000000 0.5714286 0.18181818
#> Tamansari                0.62403342 0.00000000 0.79566148 0.0000000 0.27272727
#> Tambora                  0.37136105 0.00000000 0.69774585 0.0000000 0.54545455
#> Tanahabang               0.33680315 0.16949153 0.53685872 0.0000000 0.27272727
#> Tanjung Priok            0.18412018 1.00000000 0.20422537 0.7142857 1.00000000
#> Tebet                    0.36855532 0.22033898 0.53528645 0.2857143 0.54545455
#>                          positiv_cov      rawat     isoman fcm.rentan.cluster
#> Cakung                   0.816213530 0.57746479 0.58285714                  2
#> Cempaka Putih            0.337616242 0.35211268 0.09714286                  1
#> Cengkareng               0.908213246 0.76760563 0.83428571                  2
#> Cilandak                 0.564988997 0.43661972 0.65142857                  2
#> Cilincing                0.543550792 0.51408451 0.36000000                  2
#> Cipayung                 0.723716902 0.44366197 0.70857143                  2
#> Ciracas                  0.659260311 0.51408451 0.61714286                  2
#> Duren Sawit              1.000000000 0.97887324 0.82857143                  2
#> Gambir                   0.252289345 0.17605634 0.24000000                  1
#> Grogolpetamburan         0.571874778 0.47183099 0.24571429                  1
#> Jagakarsa                0.819407965 0.56338028 1.00000000                  2
#> Jatinegara               0.618016611 0.38732394 0.49714286                  2
#> Johar Baru               0.279477532 0.17605634 0.10857143                  1
#> Kalideres                0.553489032 0.51408451 0.28571429                  2
#> Kebayoran Baru           0.434585078 0.38732394 0.56000000                  1
#> Kebayoran Lama           0.648328246 0.53521127 0.47428571                  2
#> Kebonjeruk               0.727408249 0.71830986 0.41714286                  2
#> Kelapa Gading            0.402995670 0.42957746 0.34285714                  1
#> Kemayoran                0.686945411 0.60563380 0.36000000                  2
#> Kembangan                0.627457940 0.47183099 0.43428571                  2
#> Kepulauan Seribu Selatan 0.000000000 0.00000000 0.00000000                  1
#> Kepulauan Seribu Utara   0.005324058 0.02112676 0.02285714                  1
#> Koja                     0.586569177 0.47887324 0.26857143                  2
#> Kramatjati               0.701000923 0.69014085 0.49714286                  2
#> Makasar                  0.604103074 0.41549296 0.59428571                  2
#> Mampang Prapatan         0.392489529 0.16901408 0.25142857                  1
#> Matraman                 0.437282601 0.33098592 0.21142857                  1
#> Menteng                  0.212962306 0.23239437 0.09714286                  1
#> Pademangan               0.340171790 0.34507042 0.29714286                  1
#> Palmerah                 0.530063179 0.38028169 0.24571429                  2
#> Pancoran                 0.438631362 0.50000000 0.35428571                  1
#> Pasar Minggu             0.797330872 0.64084507 0.75428571                  2
#> Pasarrebo                0.563214311 0.35915493 0.46857143                  2
#> Penjaringan              0.553134095 0.92253521 0.50285714                  2
#> Pesanggrahan             0.507844112 0.29577465 0.23428571                  2
#> Pulogadung               0.685667637 0.75352113 0.48571429                  2
#> Sawah Besar              0.271881877 0.33802817 0.14285714                  1
#> Senen                    0.300631788 0.25352113 0.21142857                  1
#> Setia Budi               0.371690211 0.43661972 0.33714286                  1
#> Tamansari                0.260168950 0.42957746 0.16000000                  1
#> Tambora                  0.397813587 0.48591549 0.24000000                  1
#> Tanahabang               0.396180876 0.45070423 0.16571429                  1
#> Tanjung Priok            0.869596081 1.00000000 0.51428571                  2
#> Tebet                    0.542485980 0.38028169 0.37714286                  1

# Inspeksi nilai-nilai kecamatan permasing-masing cluster
cluster_1 <- fcm_keep %>% 
  filter(fcm.rentan.cluster %in% 1)
cluster_1

#>                          perempuan KK_perempuan     lansia pengangguran
#> Cempaka Putih            0.6689966  0.646647181 0.67779574   0.10334253
#> Gambir                   0.4643663  0.810267949 1.00000000   0.15179393
#> Grogolpetamburan         0.7379711  0.640175890 0.74680051   0.30052131
#> Johar Baru               0.3827434  0.602835796 0.48355587   0.19748543
#> Kebayoran Baru           0.7663745  0.606061960 0.62233733   0.08065011
#> Kelapa Gading            1.0000000  0.536966368 0.97145883   0.07758356
#> Kepulauan Seribu Selatan 0.0000000  0.000000000 0.07562168   0.00000000
#> Kepulauan Seribu Utara   0.3973520  0.009924764 0.06598185   0.02483901
#> Mampang Prapatan         0.4350383  0.491608209 0.28165521   0.04507820
#> Matraman                 0.6678501  0.699908247 0.52963631   0.50843300
#> Menteng                  0.5179186  0.727487357 0.83416971   0.15915363
#> Pademangan               0.3302637  0.561075943 0.37526201   0.05826434
#> Pancoran                 0.6250971  0.503287615 0.33765364   0.09107636
#> Sawah Besar              0.4082252  0.816590504 0.69589161   0.18307268
#> Senen                    0.6108584  0.948075957 0.99188382   0.27292242
#> Setia Budi               0.6088613  0.751018436 0.48484996   0.12020852
#> Tamansari                0.6705499  1.000000000 0.97311418   0.15823367
#> Tambora                  0.3099597  0.639117672 0.52079335   0.38393131
#> Tanahabang               0.4301564  0.758500140 0.53646375   0.19595216
#> Tebet                    0.7221051  0.689762495 0.54628073   0.22293775
#>                          wajib_ktp     banjir  informal        RS  puskesmas
#> Cempaka Putih            0.5168444 0.00000000 0.3962731 0.2857143 0.09090909
#> Gambir                   0.6654871 0.00000000 0.7081245 0.1428571 0.09090909
#> Grogolpetamburan         0.5181025 0.47457627 0.4560810 0.4285714 0.63636364
#> Johar Baru               0.2989708 0.03389831 0.5344676 0.0000000 0.36363636
#> Kebayoran Baru           0.5028302 0.06779661 0.8026156 0.4285714 0.72727273
#> Kelapa Gading            0.6069006 0.72881356 0.1913311 0.2857143 0.18181818
#> Kepulauan Seribu Selatan 0.8903044 0.00000000 0.1232104 0.0000000 0.09090909
#> Kepulauan Seribu Utara   0.7727524 0.00000000 0.0000000 0.1428571 0.18181818
#> Mampang Prapatan         0.2978986 0.11864407 0.2056375 0.0000000 0.36363636
#> Matraman                 0.3899411 0.06779661 0.2264999 0.0000000 0.36363636
#> Menteng                  0.5556077 0.00000000 0.8079812 0.4285714 0.00000000
#> Pademangan               0.2456746 0.20338983 0.1741987 0.1428571 0.18181818
#> Pancoran                 0.2720585 0.22033898 0.5140726 0.1428571 0.54545455
#> Sawah Besar              0.9289591 0.00000000 0.6453875 0.5714286 0.00000000
#> Senen                    0.2043975 0.00000000 0.6131289 0.1428571 0.36363636
#> Setia Budi               0.4241382 0.01694915 1.0000000 0.5714286 0.18181818
#> Tamansari                0.6240334 0.00000000 0.7956615 0.0000000 0.27272727
#> Tambora                  0.3713610 0.00000000 0.6977458 0.0000000 0.54545455
#> Tanahabang               0.3368032 0.16949153 0.5368587 0.0000000 0.27272727
#> Tebet                    0.3685553 0.22033898 0.5352865 0.2857143 0.54545455
#>                          positiv_cov      rawat     isoman fcm.rentan.cluster
#> Cempaka Putih            0.337616242 0.35211268 0.09714286                  1
#> Gambir                   0.252289345 0.17605634 0.24000000                  1
#> Grogolpetamburan         0.571874778 0.47183099 0.24571429                  1
#> Johar Baru               0.279477532 0.17605634 0.10857143                  1
#> Kebayoran Baru           0.434585078 0.38732394 0.56000000                  1
#> Kelapa Gading            0.402995670 0.42957746 0.34285714                  1
#> Kepulauan Seribu Selatan 0.000000000 0.00000000 0.00000000                  1
#> Kepulauan Seribu Utara   0.005324058 0.02112676 0.02285714                  1
#> Mampang Prapatan         0.392489529 0.16901408 0.25142857                  1
#> Matraman                 0.437282601 0.33098592 0.21142857                  1
#> Menteng                  0.212962306 0.23239437 0.09714286                  1
#> Pademangan               0.340171790 0.34507042 0.29714286                  1
#> Pancoran                 0.438631362 0.50000000 0.35428571                  1
#> Sawah Besar              0.271881877 0.33802817 0.14285714                  1
#> Senen                    0.300631788 0.25352113 0.21142857                  1
#> Setia Budi               0.371690211 0.43661972 0.33714286                  1
#> Tamansari                0.260168950 0.42957746 0.16000000                  1
#> Tambora                  0.397813587 0.48591549 0.24000000                  1
#> Tanahabang               0.396180876 0.45070423 0.16571429                  1
#> Tebet                    0.542485980 0.38028169 0.37714286                  1

1
2
3

cluster_2 <- fcm_keep %>% 
  filter(fcm.rentan.cluster %in% 2)
cluster_2

#>                perempuan KK_perempuan     lansia pengangguran  wajib_ktp
#> Cakung         0.3464995    0.2576572 0.04693080    0.8285802 0.00000000
#> Cengkareng     0.3960945    0.2962708 0.09707556    0.8071144 0.06761805
#> Cilandak       0.7960723    0.4727347 0.46795975    0.1156087 0.37071190
#> Cilincing      0.4366656    0.4507673 0.01395345    0.7785955 0.43649576
#> Cipayung       0.4984652    0.1840792 0.15352836    0.6688132 0.10195493
#> Ciracas        0.5169200    0.3071201 0.25610534    0.4940202 0.14525009
#> Duren Sawit    0.6848996    0.3799949 0.46288044    0.5716038 0.21393593
#> Jagakarsa      0.5563815    0.2869928 0.19128910    0.3023612 0.11974172
#> Jatinegara     0.5151078    0.6082506 0.46869349    0.5611776 0.26977340
#> Kalideres      0.4207626    0.2174588 0.00000000    0.5617909 0.04974972
#> Kebayoran Lama 0.5842672    0.4449457 0.40121807    0.1787795 0.29749273
#> Kebonjeruk     0.5802729    0.3322126 0.40789323    0.2447102 0.23298989
#> Kemayoran      0.5076371    0.6662773 0.49483160    0.4320761 0.31480828
#> Kembangan      0.5665520    0.2275033 0.33098498    0.1103956 0.16474356
#> Koja           0.3856282    0.5122546 0.11573338    1.0000000 1.00000000
#> Kramatjati     0.5321203    0.4793221 0.27237959    0.3888378 0.13909914
#> Makasar        0.5462850    0.2786250 0.27694544    0.3048145 0.19855052
#> Palmerah       0.4393654    0.4627452 0.36869987    0.2210978 0.24155332
#> Pasar Minggu   0.5450276    0.4093135 0.29890696    0.1987121 0.22679113
#> Pasarrebo      0.4938422    0.1703106 0.21261087    0.3250537 0.13559586
#> Penjaringan    0.4237213    0.5048790 0.51936571    0.1993254 0.30135236
#> Pesanggrahan   0.6430711    0.3481330 0.30645587    0.1168353 0.16876389
#> Pulogadung     0.6847147    0.6001048 0.48299608    0.3830113 0.36422494
#> Tanjung Priok  0.4987241    0.4926632 0.31562520    0.8758050 0.18412018
#>                    banjir   informal        RS puskesmas positiv_cov     rawat
#> Cakung         0.52542373 0.01814115 0.1428571 0.5454545   0.8162135 0.5774648
#> Cengkareng     0.93220339 0.03286447 0.1428571 0.6363636   0.9082132 0.7676056
#> Cilandak       0.03389831 0.24102604 0.4285714 0.2727273   0.5649890 0.4366197
#> Cilincing      0.64406780 0.26129471 0.1428571 0.6363636   0.5435508 0.5140845
#> Cipayung       0.25423729 0.41400807 0.0000000 0.7272727   0.7237169 0.4436620
#> Ciracas        0.13559322 0.05311726 0.1428571 0.2727273   0.6592603 0.5140845
#> Duren Sawit    0.03389831 0.01414969 0.2857143 0.8181818   1.0000000 0.9788732
#> Jagakarsa      0.20338983 0.23868687 0.4285714 0.2727273   0.8194080 0.5633803
#> Jatinegara     0.54237288 0.25501448 0.2857143 0.8181818   0.6180166 0.3873239
#> Kalideres      0.00000000 0.01501601 0.1428571 0.9090909   0.5534890 0.5140845
#> Kebayoran Lama 0.06779661 0.06371145 0.2857143 0.5454545   0.6483282 0.5352113
#> Kebonjeruk     0.00000000 0.17592762 0.5714286 0.4545455   0.7274082 0.7183099
#> Kemayoran      0.44067797 0.27119979 0.1428571 0.4545455   0.6869454 0.6056338
#> Kembangan      0.00000000 0.17062707 0.1428571 0.5454545   0.6274579 0.4718310
#> Koja           0.28813559 0.19524772 0.4285714 0.4545455   0.5865692 0.4788732
#> Kramatjati     0.40677966 0.17963154 0.2857143 0.5454545   0.7010009 0.6901408
#> Makasar        0.44067797 0.28931734 0.4285714 0.3636364   0.6041031 0.4154930
#> Palmerah       0.00000000 0.21952901 0.2857143 0.6363636   0.5300632 0.3802817
#> Pasar Minggu   0.27118644 0.09414903 0.2857143 0.6363636   0.7973309 0.6408451
#> Pasarrebo      0.27118644 0.17693053 0.1428571 0.2727273   0.5632143 0.3591549
#> Penjaringan    0.42372881 0.17655424 0.4285714 0.3636364   0.5531341 0.9225352
#> Pesanggrahan   0.18644068 0.12285780 0.1428571 0.2727273   0.5078441 0.2957746
#> Pulogadung     0.16949153 0.22195859 1.0000000 0.5454545   0.6856676 0.7535211
#> Tanjung Priok  1.00000000 0.20422537 0.7142857 1.0000000   0.8695961 1.0000000
#>                   isoman fcm.rentan.cluster
#> Cakung         0.5828571                  2
#> Cengkareng     0.8342857                  2
#> Cilandak       0.6514286                  2
#> Cilincing      0.3600000                  2
#> Cipayung       0.7085714                  2
#> Ciracas        0.6171429                  2
#> Duren Sawit    0.8285714                  2
#> Jagakarsa      1.0000000                  2
#> Jatinegara     0.4971429                  2
#> Kalideres      0.2857143                  2
#> Kebayoran Lama 0.4742857                  2
#> Kebonjeruk     0.4171429                  2
#> Kemayoran      0.3600000                  2
#> Kembangan      0.4342857                  2
#> Koja           0.2685714                  2
#> Kramatjati     0.4971429                  2
#> Makasar        0.5942857                  2
#> Palmerah       0.2457143                  2
#> Pasar Minggu   0.7542857                  2
#> Pasarrebo      0.4685714                  2
#> Penjaringan    0.5028571                  2
#> Pesanggrahan   0.2342857                  2
#> Pulogadung     0.4857143                  2
#> Tanjung Priok  0.5142857                  2

Kluster 1 berisi 20, mencakup Kecamatan Cempaka Putih, Gambir, Grogolpetamburan, Johar Baru, Kebayoran Baru, Kelapa Gading, Kepulauan Seribu Selatan, Kepulauan Seribu Utara, Mampang Prapatan, Matraman, Menteng, Pademangan, Pancoran, Sawah Besar, Senen, Setia Budi, Tamansari, Tambora, Tanahabang, dan Tebet
Kluster 2 berisi 24 Kecamatan, mencakup Kecamatan Cakung, Cengkareng, Cilandak, Cilincing, Cipayung, Ciracas, Duren Sawit, Jagakarsa, Jatinegara, Kebayoran Lama, Kebon Jeruk, Kemayoran, Kembangan, Koja, Kramatjati, Makasar, Palmerah, Pasar Minggu, Pasarrebo, Penjaringan, Pesanggrahan, Pulogadung, dan Tanjung Priok.

Interpretasi: Cluster Profiling

Untuk dapat melakukan cluster profiling, dibutuhkan informasi pusat cluster (centroid) masing-masing variabel terhadap cluster yang dihasilkan. Dalam object fcm.rentan, kita dapat melakukan inspeksi menggunakan fcm.rentan$v yang mana merupakan matriks numerik yang berisi pusat masing-masing cluster (centroid).

Final cluster prototype

`1`	`fcm.rentan$v`

#>           perempuan KK_perempuan    lansia pengangguran wajib_ktp    banjir
#> Cluster 1 0.5389905    0.6162944 0.5619571    0.2209669 0.4455888 0.1271250
#> Cluster 2 0.5210274    0.3989750 0.3105454    0.4079368 0.2562638 0.2963236
#>            informal        RS puskesmas positiv_cov     rawat    isoman
#> Cluster 1 0.4673842 0.2265458 0.3294159   0.3923736 0.3665275 0.2634173
#> Cluster 2 0.1986334 0.2812407 0.5169692   0.6422221 0.5513181 0.5059506

Cluster Profiling:

Cluster 1 :
- Pada cluster ini dominan variabel yang terkait perempuan
- Variable yang lebih dominan pada cluster: perempuan, KK_perempuan, lansia, wajib_ktp, informal.
- Lansia pada kecamatan-kecamatan yang ada di cluster 1 memang memiliki persentase nilai yang tinggi pada data asli. Hal ini mengindikasikan, kecamatan-kecamatan yang ada di cluster 1 bisa dilakukan penanggulangan atau bantuan yang dikhususkan untuk lansia.
Cluster 2 :
- Pada cluster ini dominan variabel yang terkait var terkait lingkungan dan kesehatan
- Variable yang lebih dominan pada cluster: pengangguran, banjir, RS, puskesmas, positiv_cov, rawat, dan isoman
- Untuk variabel positiv_cov, rawat, dan juga isoman, variabel- variabel menjadi masalah di kecamatan-kecamatan yang ada pada cluster 2. hal ini dibuktikan dengan nilai jumlah positif covid, jumlah orang yang dirawat, serta jumlah orang melakukan isoman akibat COVID yang tinggi.
- keadaan positive_cov menyebabkan banyaknya warga yang harus diisoman dan di beberapa kecamatan terlihat lebih sedikit yang perlu dirawat

1
2
3

# pak calvin
rentan %>% 
  select(Kecamatan, positiv_cov, rawat, isoman)

#> # A tibble: 44 × 4
#>    Kecamatan        positiv_cov rawat isoman
#>    <chr>                  <dbl> <dbl>  <dbl>
#>  1 Cakung                 11740    82    102
#>  2 Cempaka Putih           4998    50     17
#>  3 Cengkareng             13036   109    146
#>  4 Cilandak                8201    62    114
#>  5 Cilincing               7899    73     63
#>  6 Cipayung               10437    63    124
#>  7 Ciracas                 9529    73    108
#>  8 Duren Sawit            14329   139    145
#>  9 Gambir                  3796    25     42
#> 10 Grogolpetamburan        8298    67     43
#> # … with 34 more rows

Berikut merupakan hasil visualisasi kluster dalam plot 2 dimensi.

# mengubah object ppclust menjadi tipe lain dari object cluster 
fcm.rentan2 <- ppclust2(fcm.rentan, "fanny")
# gunakan fungsi fviz_cluster() dari package factoextra
factoextra::fviz_cluster(fcm.rentan2, data = rentan_FCM, 
  ellipse.type = "convex",
  palette = "jco",
  repel = TRUE)

\> Untuk daerah yang terletak di antara 2 cluster atau memiliki membership degree yang fifty-fifty indikasinya daerah/kecamatan tersebut cenderung memiliki nilai yang cukup dominan/tinggi di variabel-varibel keseluruhan.

Implikasi Kebijakan

Peta kerentanan sosial juga dapat menjadi informasi yang baik dalam penanggulangan bencana. Integrasi dengan sistem peringatan yang dirancang untuk mencegah penularan penyakit dapat memberikan informasi yang berguna bagi wilayah yang memiliki risiko penyebaran COVID-19 yang tinggi. Permasalahan yang disebabkan oleh kependudukan, pekerjaan, kesehatan, dan pendidikan yang diakibatkan oleh COVID-19, dapat ditanggulangi dengan mengintegrasikan perencanaan penanggulangan resiko dengan program pengembangan lainnya.

Pada analisis ini ditemukan 2 faktor yang berpengaruh terhadap kerentanan sosial yaitu Faktor kesehatan, kependudukan/ Socio-health serta Faktor demografi, perempuan, dan pekerjaan Oleh karena itu, disarankan bahwa penanggulan dampak COVID-19 dapat diintegrasikan dengan program pengembangan lainnya yang sesuai dengan beberapa faktor yang telah disebutkan.

`1`	`knitr::include_graphics("assets/csr.PNG")`

Oleh karena itu, hasil analisis ini dapat menjadi masukan bagi pemerintah dalam hal ini Badan Penanggulangan Bencana Daerah (BPBD) DKI Jakarta, dimana hasilnya dapat digunakan untuk perencanaan program penanggulangan bencana kedepannya. Lebih jauh lagi berdasarkan hasil analisis ini dapat dilihat tingkat kerentanan yang berbeda pada masing-masing wilayah, Hal ini berarti kapasitas kesiapan dan respon terhadap COVID-19 berbeda beda. Bagi CSR perusahaan, hal ini dapat menjadi rujukan informasi yang bersesuaian agar lebih tepat sasaran dalam memberikan bantuan kepada daerah dengan karakteristik tertentu

Ketika informasi yang terbatas, penggunaan peta kerentanan sosial dapat dijadikan informasi untuk memprioritaskan terlebih dahulu daerah dengan tingkat kerentanan sosial yang tinggi.

Kesimpulan

Analisis ini merupakan upaya yang dilakukan untuk mendukung penanggulangan dampak yang diakibatkan oleh COVID-19 di DKI Jakarta. Pendekatan indeks dan klaster dilakukan untuk mengetahui bagaimana indeks kerentanan dan pengaruh variabel terhadap masing-masing wilayah kecamatan DKI Jakarta. Dari project ini, kami mendapatkan 2 faktor yang membentuk kerentanan sosial, yaitu Faktor kesehatan, kependudukan/ Socio-health serta Faktor demografi, perempuan, dan pekerjaan

Analisis ini menggunakan 12 variabel untuk mengukur kerentanan sosial. Hal ini dikarenakan keterbatasan data yang ada. Menggunakan data yang lebih lengkap dan terbaru maka pengukuran akan menjadi lebih baik. Akhirnya, poin penting bahwa project ini melakukan penilaian terhadap kerentanan sosial pada kecamatan di DKI Jakarta. Karakteristik wilayah yang berbeda beda menyebabkan penggunaan variabel yang berbeda juga, sehingga penilaian terhadap kerentanan sosial akan menjadi tantangan tersendiri.

Plot on Algoritma Technical Blog

Social Vulnerability Analysis to COVID-19 in DKI Jakarta Using Factor and Cluster Analysis

Preface

Pendahuluan

Library dan Setup

Tujuan Pembelajaran

Data Loading

Data Cleansing

Melakukan pengecekan missing value

Melakukan pengecekan duplikasi data

Exploratory Data Analysis

1. Gambaran umum variabel penyusun indeks kerentanan sosial

2. Melihat visualisasi antar beberapa variabel

Penyusun Kerentanan Sosial

Analisis Faktor / Factor Analysis dan Asumsi

Melakukan uji multikolinearitas pada seluruh variabel.

Kaiser–Meyer–Olkin (KMO) - Measure of Sampling Adequacy(MSA)

The Bartlett’s Test statistics

Membentuk Analisis Faktor

Pembentukan Faktor Umum dan Interpretasi

Indeks Kerentanan Sosial

Melakukan PCA pada data `rentan`

Visualisasi Geospasial Kerentanan Sosial Kecamatan di Provinsi DKI Jakarta

Mengambil data spasial

Membaca data geometri spasial

Combining Data

Mengubah object spasial geometri 3 dimensi yang kita miliki menjadi 2 dimensi

Membuat Choropleth Map menggunakan `leaflet()`

Fuzzy C-Means Clustering terhadap Kerentanan Sosial

Fuzzy C-means Methodology

Data pre-processing

Normalisasi Min-Max

Pemilihan nilai k optimum

FCM menggunakan `fcm()`

Interpretasi: Cluster Profiling

Implikasi Kebijakan

Kesimpulan

Reference

Plot on Algoritma Technical Blog

Social Vulnerability Analysis to COVID-19 in DKI Jakarta Using Factor and Cluster Analysis

Preface

Pendahuluan

Library dan Setup

Tujuan Pembelajaran

Topik: COVID-19 SOCIAL VULNERABILITY

Data Loading

Data Cleansing

Melakukan pengecekan missing value

Melakukan pengecekan duplikasi data

Exploratory Data Analysis

1. Gambaran umum variabel penyusun indeks kerentanan sosial

2. Melihat visualisasi antar beberapa variabel

Penyusun Kerentanan Sosial

Analisis Faktor / Factor Analysis dan Asumsi

Melakukan uji multikolinearitas pada seluruh variabel.

Kaiser–Meyer–Olkin (KMO) - Measure of Sampling Adequacy(MSA)

The Bartlett’s Test statistics

Membentuk Analisis Faktor

Pembentukan Faktor Umum dan Interpretasi

Indeks Kerentanan Sosial

Melakukan PCA pada data rentan

Visualisasi Geospasial Kerentanan Sosial Kecamatan di Provinsi DKI Jakarta

Mengambil data spasial

Membaca data geometri spasial

Combining Data

Mengubah object spasial geometri 3 dimensi yang kita miliki menjadi 2 dimensi

Membuat Choropleth Map menggunakan leaflet()

Fuzzy C-Means Clustering terhadap Kerentanan Sosial

Fuzzy C-means Methodology

Data pre-processing

Normalisasi Min-Max

Pemilihan nilai k optimum

FCM menggunakan fcm()

Interpretasi: Cluster Profiling

Implikasi Kebijakan

Kesimpulan

Reference

Melakukan PCA pada data `rentan`

Membuat Choropleth Map menggunakan `leaflet()`

FCM menggunakan `fcm()`