HIERARCHICAL CLUSTERING MULTISCALE BOOTSTRAP UNTUK PENGELOMPOKAN KEMISKINAN DI JAWA TIMUR

## Gangga Anuraga
Program Studi Statistika Universitas PGRI Adi Buana (UNIPA) Surabaya Alamat e-mail : ganuraga @gmail.com

## ABSTRAK
Kemiskinan di Jawa Timur merupakan problematika hingga sekarang. Berbagai program telah direncanakan guna menanggulangi permasalahan tersebut. Agar program tersebut dapat terlaksana optimal, perlu strategi perencanaan yang mendalam terutama terkait pemetaan daerah dan karakteristik kemiskinannya. Metode hierarki complete linkage digunakan untuk mengelompokkan daerah-daerah berdasarkan kemiripan karakteristik. Penggunaan metode hierarki dengan teknik pengukuran jarak hanya memberikan satu solusi dalam penyelesaiannya, yaitu didasarkan pada ukuran kemiripan pada teknik jarak yang digunakan. Sehingga penelitian ini menggunakan pendekatan multiscale bootstrap untuk memberikan ukuran ketidaktentuan dalam metode pengelompokan klasik hierarchical clustering . Metode ini bekerja dengan pendekatan bootstrap resampling untuk setiap kelompok (cluster) .Pendekatan multiscale bootstrap dalam Hierarchical Cluster Analysis memiliki kestabilan pada saat iterasi B = 500. Hal itu dapat diketahui dari signifikansi nilai AU (P-values) yang mendekati 0.95 dengan nilai standar error bootstrap yang kecil.Terdapat lima kelompok (cluster) yang terbentuk dengan kemiripan satu dengan yang lain dan signifikan dengan nilai AU (Approximately Unbiased) p_value  0,95. Berdasarkan signifikansi Approximately unbiassed (AU) p-values  0,95, kelompok kelima yang terdiri dari Kabupaten Jember, Kabupaten Bondowoso, Kabupaten Situbondo, Kabupaten Probolinggo, Kabupaten Pasuruan dapat dikatakan sebagai daerah yang masih cukup tinggi persentase kemiskinannya dibandingkan dengan kelompok 1, 2, 3, dan 4.
Kata Kunci : Approximately Unbiassed (AU) p-values , clustering hierarki complete linkage, multiscale bootstrap , kemiskinan

## PENDAHULUAN
Analisis kelompok ( cluster analysis ) adalah salah satu analisis statistik yang bertujuan mengelompokkan obyek-obyek berdasarkan kesamaan karakteristik diantara obyek-obyektersebut. Obyek tersebut akan diklasifikasikan ke dalam satu atau beberapa kelompok ( cluster ) sehingga obyek-obyek yang berada dalam satu kelompok akan mempunyai kemiripan satu dengan yang lain. Dengan demikian
terdapat homogenitas (kesamaan) yang tinggi antar anggota
dalam satu kelompok ( within-cluster ) dan heterogenitas (perbedaan) yang tinggi antar kelompok yang satu dengan kelompok lainnya ( between-cluster ) [7]. Terdapat dua metode umum dalam algoritma cluster , yaitu hierarki dan non hierarki. Pada metode non-hierarki ( non hierarchical clustering ), jumlah kelompok ditentukan di awal sedangkan pada metode hierarki ( hierarchical clustering ), jumlah kelompok ditentukan di akhir analisis. Metode hierarki dibagi berdasarkan teknik pengukuran jarak antar kelompok yang digunakan untuk
mengelompokkan objek. Pengukuran jarak dalam metode hierarki antara lain single linkage yang menggunakan jarak minimum, complete linkage mengguna- kan jarak maksimum, average linkage menggunakan jarak rata-rata, dan median linkage yang menggunakan jarak median antar kelompok.
Penggunaan metode hierarki dengan teknik pengukuran jarak hanya memberikan satu solusi dalam penyelesaiannya, yaitu didasarkan pada ukuran kemiripan pada teknik jarak yang digunakan. Sehingga Nilai signifikansi statistik dan ukuran akurasi dalam pengelompokan tidak dapat ditemukan dengan metode klasik ini. [10] mengenalkan
pendekatan
multiscale bootstrap untuk memberikan ukuran ketidaktentuan dalam metode pengelompokan klasik hierarchical clustering . Metode ini bekerja dengan pendekatan bootstrap resampling untuk setiap kelompok (cluster) . Sehingga didapatkan ukuran validasi yaitu nilai p- values dalam membentuk kelompok (cluster) yang mempunyai kemiripan satu dengan yang lain. Terdapat dua nilai p- values dalam pendekatan bootstrap resampling yaitu probabability bootstrap (BP) value dan approximately unbiassed (AU) p-values. Multiscale bootstrap resampling digunakan untuk menghitung nilai Approximately unbiassed (AU) p- values, yang mana memberikan taksiran yang lebih baik dalam mengatasi bias. Beberapa penelitian tentang klasifikasi dengan pendekatan bootstrap antara lain adalah [9] memperkenalkan Approximately unbiassed (AU) p-values pada klasifikasi gen. [8] menggunakan pendekatan bootstrap untuk klasifikasi atau pengelompokan pada data microarray , yang digunakan untuk mengukur reliabilitas pengelompokan. Kemiskinan merupakan masalah multidimensi dan lintas sektor yang dipengaruhi oleh berbagai faktor yang saling berkaitan, antara lain tingkat
pendapatan, kesehatan, pendidikan, akses terhadap barang dan jasa, lokasi, geografis, gender, dan kondisi lingkungan (BAPPENAS, 2010). Menurut [2], pada bulan Maret 2012 penduduk miskin Jawa Timur sebanyak 5,071 juta (13,40 persen) atau turun 5,53 persen dibandingkan dengan jumlah penduduk miskin pada bulan Maret 2011 yang mencapai 5,356
juta jiwa (14,23 persen). Dengan wilayah yang luas sekitar ± 47.922 km 2 , Provinsi Jawa Timur terbagi atas 29 kabupaten dan 9 kota atau secara administratif terdapat 38 Kabupaten/Kota, sehingga dapat dikatakan berkembang dalam bidang ekonomi dan dapat menjadi potensi pembangunan. Namun demikian, dalam realitanya masih banyak masyarakat yang hidup di bawah persentase garis kemiskinan. Pada tahun 2011 persentase kemiskinan di Provinsi Jawa Timur masih berada diatas persentase garis kemiskinan tingkat nasional sebesar 13,33%. Berbagai program telah dicanangkan pemerintah baik pusat maupun daerah untuk menang-gulangi masalah kemiskinan mulai dari bidang pendidikan, pangan, kesehatan, perluasan kesempatan kerja, bantuan sarana dan prasarana pertanian, bantuan kredit usaha untuk masyarakat miskin dan lain sebagainya. Banyaknya program yang dicanangkan tentu saja perlu didukung oleh strategi perencanaan yang matang. Strategi perencanaan program bisa saja berbeda untuk tiap-tiap daerah.
Oleh karena itu, untuk mencapai terselenggaranya program pengentasan kemiskinan yang optimal, perlu dilakukan pengelompokan karakteristik kemiskinan di Kabupaten/Kota di Jawa Timur. Sehingga dilakukan penelitian dengan
menggunakan
Hierarchical clustering dengan pendekatan multiscale bootstrap
untuk pengelompokan kemiskinan di jawa timur tahun 2011, dimana variabel yang digunakan antara
lain indikator-indikator kemiskinan di sektor ekonomi, kesehatan, dan SDM.
Analisis kelompok (cluster analysis) digunakan untuk mengelompokkan objek-objek pengamatan menjadi beberapa kelompok berdasarkan pengukuran variabel-variabel yang diamati, sehingga objek dalam kelompok yang sama memiliki kemiripan karakteristik dan antarkelompok memiliki ketakmiripan karakteristik. Metode pengelompokan hierarki (Hierarchical clustering) memulai memulai pembentukan kelompok dengan jumlah kelompok sama dengan jumlah obyek (n) kemudian dua kelompok yang mempunyai jarak terdekat digabung menjadi satu kelompok sehingga jumlah kelompok menjadi (n-1). Pada penelitian ini digunakan metode complete linkage
karena metode ini menghasilkan perbedaan faktor-faktor yang mempengaruhi karakteristik kemiskinan pada masing-masing kelompok.
Selanjutnya jarak-jarak antar kelompok yang baru dihitung kembali. Prosedur ini diulang sampai terbentuk satu kelompok dengan anggota sebanyak n. Kriteria yang digunakan untuk menentukan jarak antar dua kelompok dalam penelitian ini adalah metode pautan lengkap ( complete linkage ) yaitu dengan memaksimumkan jarak antara kelompok. Jarak maksimum antara kelompok satu yang merupakan gabungan dari kelompok u dan v dengan kelompok wadalah sebagai berikut [7] :
    m a x , U W V W U V W d d d  (1)
Adapun rumus untuk menghitung jarak korelasi obyek i dan j pada k variabel adalah sebagai berikut
         1 1/2 2 1 1 1 1 dimana x adalah rata -rata dari variabel , , p i j ik jk k ij p p i j ik jk k k i i ip x x x x
d x x x x
x x             
 
    (2)
Metode bootstrap pertama kali
dikembangkan oleh Efron pada tahun
1979. Bootstrap adalah teknik yang dikembangkan untuk membuat inferensia statistik tertentu yang lebih sederhana dengan memerlukan komputer modern. Menurut Efron, bootstrap berbeda dengan inferensia statistika, metodenya sangat sederhana dan berdasar pada prosedur
resampling . Secara singkat algoritma bootstrap dapat dinyatakan sebagai berikut [4] :
1. Sampel data D didefinisikan sebagai data sampel berukuran n yang terdiri dari D i = D 1 , D 2 , ..., D n dimana D i merupakan vektor data pengamatan.
2. Sampel data diambil secara acak dengan pengembalian sebanyak n kali.
Diperoleh data sampel baru yang didefinisikan sebagai D *.
3. Langkah (2) dilakukan secara berulang sebanyak B kali sehingga didapatkan himpunan data bootstrap ( D * 1 , D * 2 , ..., D * B ) dimana setiap sampel bootstrap merupakan sampel acak yang saling independen. Multiscale bootstrap resampling digunakan untuk menghitung nilai probabability bootstrap (BP) value dan Approximately unbiassed (AU) p-values. Berikut ilustrasi penerapan Multiscale bootstrapresampling pada hierarchical clustering [10]:
1. Sampel data D didefinisikan sebagai data sampel berukuran n yang terdiri dari D i = D 1 , D 2 , ..., D n dimana D i merupakan vektor data pengamatan.
2. Sampel data diambil secara acak dengan pengembalian sebanyak n kali.
Diperoleh data sampel baru yang didefinisikan sebagai D *. 3. Langkah (2) dilakukan secara berulang sebanyak B kali sehingga didapatkan himpunan data Bootstrap ( D * 1 , D * 2 , ..., D * B ) dimana setiap sampel bootstrap merupa-kan sampel acak yang saling independen.

## METODE PENELITIAN

## Sumber Data dan Variabel Penelitian
Data yang digunakan pada penelitian ini merupakan publikasi Data dan Informasi Kemiskinan Kab/Kota Tahun 2011 hasil olah data triwulananSurvei Sosial Ekonomi Nasional (SUSENAS). Informasi yang didapatkan antara lain adalah 16 indikator kemiskinan di sektor ekonomi, kesehatan, dan SDM.
Variabel penelitian yang digunakan dalam penelitian ini adalah sebagai berikut :
Tabel 1. Variabel Penelitian

## Ekonomi
Persentase penduduk miskin usia 15 tahun keatas yang tidak bekerja (X 1 ) Persentase penduduk miskin usia 15 tahun keatas yang bekerja di sektor pertanian (X 2 ) Persentase rumah tangga yang pernah membeli beras raskin (X 3 ) Persentase pengeluaran per kapita untuk non makanan (X 4 )

## SDM
Persentasependudukmiskinusia 15 tahunkeatas yang tidaktamat SD (X 5 )
Angka Melek Huruf penduduk miskin usia 15- 55 tahun (X 6 )
Angka Partisipasi Sekolah penduduk miskin usia 13-15 tahun (X 7 )
Rata-rata lama sekolah (X 8 )

## Kesehatan
Persentase Perempuan pengguna alat KB di rumah tangga miskin (X 9 ) Persentase Balita di rumah tangga miskin yang proses kelahirannya ditolong oleh tenaga kesehatan(X 10 )
Persentase Balita di rumah tangga miskin yang telah diimunisasi (X 11 )
Persentase rumah tangga miskin dengan luas lantai perkapita ≤ 8 m 2 (X 12 )
Persentase rumah tangga miskin yang menggunakan air bersih (X 13 )
Persentase rumah tangga miskin yang menggunakan jamban sendiri/bersama (X 14 ) Persentase rumah tangga miskin yang mendapatkan pelayanan jamkesmas (X 15 ) Angka harapan hidup (X 16 )

## Langkah-langkah Penelitian
Metode dan tahapan analisis yang akan digunakan dalam mencapai tujuan penelitian adalah sebagai berikut :
a. Sampel data D didefinisikan sebagai data sampel berukuran n yang terdiri dari D i = D 1 , D 2 , ..., D n dimana D i merupakan vektor data pengamatan.
b. Sampel data diambil secara acak dengan pengembalian sebanyak n kali. Diperoleh data sampel baru yang didefinisikan sebagai D *. c. Langkah (2) dilakukan secara berulang sebanyak B kali
sehingga didapatkan himpunan data Bootstrap ( D * 1 , D * 2 , ..., D * B ) dimana setiap sampel bootstrap merupa-kan sampel acak yang saling independen.
d. Untuk setiap kelompok (cluster) dari sampel bootstrap D * B , nilai probababilitybootstrap (BP) value adalah
  1 1 ( ),dimana . adalahinversdari fungsi normalstandar. z BP    
e. Approximately unbiassed (AU) p- values adalah 1 ( ) AU v c     
dimana v dan c adalah estimasi parameter dari setiap kelompok (cluster) . Approximately unbiassed (AU) p-values memiliki bias yang lebih kecil dibandingkan dengan probababilitybootstrap (BP) value .

## HASIL PENELITIAN

## Kemiskinan di JawaTimur
Pada Gambar 1 diketahui bahwa daerah kantong kemiskinan masih seputar Sampang (30.21%), Bangkalan (26.22 persen), Sumenep (23.10 persen), Pamekasan (20.94 persen),Ngawi (16.74
persen), Bojonegero (17.47 persen), Tuban (18.78 persen),Lamongan (17.41), Pacitan (18.13 persen), Kota Probolinggo
(17.74 persen), Probolinggo (23.48
persen), dan Bondowoso (16.66 persen).
Gambar 1. Persebaran persentase penduduk miskin di Jatim

## Hierarchical Clustering Multiscale Bootstrap
Berikut disajikan analisis dan pembahasan metode hierarki complete linkage pendekatan multiscale bootstrap .
Gambar 2. Dendrogram dengan nilai AU/BP (B=500) pada masing-masing Kabupaten / Kotamadya di Provinsi Jawa Timur
Gambar 2 menunjukkan dendrogram beserta nilai Approximately unbiassed
(AU) p-values
dan probabability bootstrap (BP) value untuk 38
Kabupaten/Kota di Propinsi Jawa Timur. Kelompok dengan AU (approximately unbiased)  0,95 ditunjukkan dalam satu kelompok (cluster) dan dimana dapat diinterpretasikan sebagai kelompok yang signifikan dengan tingkat signifikansi 0.05. Berdasarkan signifikansi Approximately unbiassed (AU) p-values kelompok pertama terdiri dari Kabupaten Pacitan dan Trenggalek. Kelompok kedua adalah Kabupaten Sidoarjo, Kota. Madiun, Kota. Batu. Kelompok ketiga
K a b S a m p a n g K a b P a m e k a s a n K a b S u m e n e p K a b J e m b e r K a b P ro b o li n g g o K a b P a s u ru a n K a b B o n d o w o s o K a b S it u b o n d o K o ta M a la n g K o ta B li ta r K o ta K e d ir i K o ta M o jo k e rt o K o ta S u ra b a y a K o ta P ro b o li n g g o K o ta P a s u ru a n K a b S id o a rj o K o ta M a d iu n K o ta B a tu K a b B a n g k a la n K a b T u lu n g a g u n g K a b M a g e ta n K a b L a m o n g a n K a b P o n o ro g o K a b M a d iu n K a b N g a w i K a b K e d ir i K a b B a n y u w a n g i K a b J o m b a n g K a b T u b a n K a b B li ta r K a b M a la n g K a b M o jo k e rt o K a b G re s ik K a b P a c it a n K a b T re n g g a le k K a b B o jo n e g o ro K a b L u m a ja n g K a b N g a n ju k 0 .0 0 .2 0 .4 0 .6 0 .8 Cluster dendrogram with AU/BP values (%) Cluster method: complete Distance: correlation H e ig h t 90 99 90 100 95 91 95 98 88 89 90 92 70 97 98 95 93 96 81 95 0 90 58 95 86 71 96 93 96 85 90 84 87 93 91 92 au 35 35 29 65 32 17 31 24 14 30 21 9 50 43 8 2 5 10 23 17 1 15 0 1 20 0 17 6 9 0 1 1 1 8 0 1 bp 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 27 28 29 30 31 32 33 34 35 36 edge #
terdiri dari Kabupaten Ponorogo, Kabupaten Tulungagung, Kabupaten Blitar, Kabupaten Kediri, Kabupaten Malang, Kabupaten Banyuwangi, Kabupaten Mojokerto, Kabupaten Jombang, Kabupaten Madiun, Kabupaten
Magetan, Kabupaten Ngawi, Kabupaten Tuban, Kabupaten Lamongan, Kabupaten Gresik. Kelompok keempat terdiri dari Kota Kediri, Kota Blitar, Kota Malang, Kota Mojokerto. Dan kelompok kelima yaitu terdiri dari Kabupaten Jember,
Kabupaten Bondowoso, Kabupaten Situbondo, Kabupaten Probolinggo,
Kabupaten Pasuruan.
Gambar 3. AU P-values dengan standard error bootstrap 500
Kestabilan proses dalam membentuk kelompok (cluster) yang mempunyai
kemiripan satu dengan yang lain dan signifikan didapatkan pada replikasi B = 500, dimana nilai standar error yang kecil dengan nilai AU/BP value yang mendekati nilai lebih besar dari 0,95.

## KESIMPULAN
Dari hasil penelitian dapat disimpulkan bahwa terdapat lima kelompok (cluster) yang terbentuk dengan kemiripan satu dengan yang lain dan signifikan dengan nilaiAU (Approximately Unbiased) p_value  0,95. Kestabilan proses dalam membentuk
kelompok (cluster) tercapai pada replikasi B = 500. Berdasarkan signifikansi Approximately unbiassed (AU) p-values  0,95, kelompok kelima yang terdiri dari Kabupaten Jember, Kabupaten Bondowoso, Kabupaten Situbondo, Kabupaten Probolinggo, Kabupaten Pasuruan dapat dikatakan sebagai daerah yang masih cukup tinggi persentase kemiskinannya dibandingkan dengan kelompok 1, 2, 3, dan 4.

## DAFTAR PUSTAKA
[1] Anuraga, G., & Otok, B. W., 2014, Spasial Structural Equation Modeling-Partial Least Square untuk Pemodelan Kemiskinan di Jawa Timur, Statistika [2] BPS, 2012, Data dan Informasi Kemiskinan Kabupaten/Kota 2011 , Badan Pusat Statistik.Jakarta: [3] Badan Pusat Statistik, 2012, Perhitungan dan Analisis Kemiskinan Makro Indonesia 2012
[4] Efron, B. dan Tibshirani, R.J., 1993, An Introduction to the Bootstrap . Chapman Hall, New York.
[5] Hair, Jr. F., Black, W.C., Babin, B.J., dan Anderson, R.E., 2010
, Multivariate Data Analysis , Pearson Prentice Hall.
[6] Hastie, T., Tibshirani, R., dan
Friedman, J., 2001, The Elements of
Statistical Learning: Data Mining,
Inference, and Prediction , Springer- Verlag, New York.
[7] Johnson, R.A. dan Wichern, D.W., 2007, Applied Multivariate Statistical Analysis , Prentice Hall, New Jersey
[8] Kerr,M.K. and Churchill,G.A., 2001,
Bootstrapping cluster analysis: assessing the reliability of
0.0 0.2 0.4 0.6 0.8 1.0 0 .0 0 .5 1 .0 1 .5 2 .0 p-value vs standard error plot AU p-value S ta n d a rd E rr o r
conclusions
from microarray experiments , USA.
[9] Shimodaira,H., 2004, Approximately unbiased tests of regions using multistepmultiscale bootstrap resamplin g. Ann. Stat., 32, 2616–
2641.
[10] Suzuki, R. dan Shimodaira, H.,
2004, An application of multiscale
bootstrap resampling to hierarchical clustering of microarray: How accurate are these clusters, In proceedings by the fifteenth International Conference on Genome
Informatics (GIW 2004), hal. 34