## REGRESI KUANTIL PENDEKATAN BOOTSTRAP UNTUK PEMODELAN KEMISKINAN

## DI PULAU JAWA
1 Gangga Anuraga, 2 Permadina K. Arieska
1 Program Studi Statistika, Fakultas Matematika dan Ilmu Pengetahuan Alam, Universitas PGRI Adi Buana Surabaya
2 Universitas Nahdhatul Ulama, Surabaya
Email : ganuraga@gmail.com

## ABSTRAK
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. Untuk itu, dalam rangka menunjang keberhasilan pelaksanaan program pembangunan terutama yang berkaitan dengan penanggulangan kemiskinan di Indonesia khususnya di Pulau Jawa, diperlukan suatu penelitian yang dapat mengetahui informasi mengenai faktor-faktor yang berpengaruh terhadap kemiskinan.Penelitian ini menggunakan regresi kuantil dengan pendekatan bootstrap , metode ini memiliki kelebihan dapat mengatasi masalah pencilan (outlier) dan heteroskedastisitas. Nilai kuantil yang digunakan dalam penelitian yaitu kuantil ke- τ = 0,25 ; τ = 0,5 ; τ = 0,75 dan τ = 0,99 dengan resampling 500 pada bootstrap standard error .Hasil regresi kuantil pada masing-masing kuantil dengan resampling 500 bootstrap standard error menunjukkan bahwa kuantil τ = 0,25 memiliki standar error yang relatif lebih kecil daripada kuantil τ = 0,5 , τ = 0,75 dan τ = 0,99 . Model regresi kuantil dalam penelitian ini mengindikasikan adanya multikolinieritas, sehingga dengan nilai pseudo R 2 sebesar 61,4% terdapat banyak variable yang tidak signifikan dalam model pada ke- τ = 0,25 .
Kata Kunci : Bootstrap standard error , Heteroskedastisitas, Kemiskinan, Regresi Kuantil

## PENDAHULUAN
Analisis regresi merupakan suatu teknik statistik yang digunakan untuk mendapatkan hubungan dan model matematis antara variabel respon (Y) dengan satu atau lebih variabel prediktor (X).Analisis regresi bertujuan untuk memperkirakan atau meramalkan nilai dari variabel respon apabila nilai dari variabel prdiktor sudah diketahui [7]. Dalam analisis regresi pendugaan parameter atau estimasi koefisien regresi dapat diduga dengan baik melalui metode kuadrat terkecil atau Ordinary Least Square (OLS).Metode OLS merupakan salah
satu metode yang populer dan telah sering digunakan untuk menduga parameter dalam analisis regresi klasik.Dalam regresi klasik,
[9] menyatakan bahwa untuk mendapatkan pendugaan yang tak bias atau Best Linear Unbiased Estimator (BLUE) harus memenuhi beberapa asumsi dalam
estimasinya. Asumsi yang harus dipenuhi dalam analisis regresi klasik diantaranya adalah tidak terjadi multikolinieritas atau korelasi antar variabel prediktor, residual atau error berdistribusi normal dengan rata- rata nol dan memiliki varian yang konstan, tidak terjadi otokorelasi dan heteroskedastisitas. [10] menyatakan
bahwa regresi klasik dengan metode OLS hanya difokuskan padasampel rata-rata atau didasarkan pada fungsi distribusi mean sebagai solusi dari meminimumkan jumlah kuadrat error . Dengan demikian, metode OLS menjadi sangat sensitif terhadap data yang memuat pencilan atau
outlier .Keberadaan pencilan atau outlier pada data juga dapat menyebabkan residual ( error ) semakin besar, sehingga varians pada residual tidak konstan, yang dalam regresi klasik dinamakan sebagai pelanggaran terhadap asumsi homogenitas. Regresi kuantil memberikan pandangan yang lebih komprehensif tentang hubungan antara
variabel respondan variabel prediktor .
Regresi kuantil diperkenalkan pertama kali oleh [11]. Metode ini merupakan pengembangan dari OLS pada regresi klasik, yaitu dengan menduga berbagai fungsi kuantil dari suatu distribusi Y sebagai fungsi dari X [13]. Teknik dan metode ini dapat mendeteksi hubungan yang lebih halus antara variable respon dengan prediktor dan sangat berguna jika distribusi data terdapat pencilan (outlier) , tidak homogen atau terjadi kasus heterokesdatisitas.
Regresi kuantil memberikan pendugaan parameter yang lebih stabil dan konsisten, yaitu dengan membatasi distribusi data yang memiliki pencilan (outlier) dan heterokesdatisitas .
[14] menyatakan bahwa pendugaan parameter dalam regresi kuantil diperoleh berdasarkan minimum jumlah mutlak dari residual (error) yang terboboti. [10] pendugaan parameter dalam regresi kuantil tidak dapat diperoleh secara analitik akan tetapi menggunakan algoritma berdasarkan pemrograman linear. Metode pendugaan parameter secara iterasi dengan pemrogramanlinear antara lain yaitu metode simplex , interior point dan smoothing. [12] metode simpleks (simplex) memberikan hasil estimasi yang konsisten pada data yang berukuran tidak terlalu besar (moderate sample) . [17] melakukan simulasi dengan
algoritma interior-point dan smoothing dengan sampel data berukuran besar yaitu n > 10 5 dan hasilnya menunjukkan estimasi yang konsisten. Selanjutnya untuk menghitung selang kepercayaan pada regresi kuantil dapat dilakukan dengan tiga pendekatan yaitu pendekatan sparsity , rank-score , dan resampling .
Metode sparsity digunakan bila residual (error) terdistribusi identik dan independen (i.i.d). Pendekatan rank- score , meng-hitung selang kepercayaan dengan menggunakan fungsi rank-score dan resampling dengan menggunakan
teknik bootstrap .
Berbagai penelitian terkait regresi kuantil, diantaranya adalah [15] melakukan penelitian tentang pemodelan tingkat pengangguran terbuka di indonesia, hasil kajian dengan estimasi sparsity menunjukkan bahwa hasil taksiran interval dari regresi kuantil lebih baik dan fleksibel. [21] membahas tentang model regresi kuantil pada kasus Indek Pembangunan Manusia (IPM). Selanjutnya [18] melakukan penelitian mengenai data temperatur suhu harian di Kota Sydney. Penelitian tersebut menyimpulkan bahwa pada kuantil 10%, 25%, 50%, 75%, dan 90% suhu kemarin ber-pengaruh terhadap suhu hari ini.
Kemudian pada bidang ekonomi regresi kuantil digunakan untuk pemodelan tingkat upah [6] yaitu menentukan pengaruh tingkat pendidikan, pengalaman danserikatkeanggotaanterhadap tingkat upah. Selanjutnya berkembang juga aplikasinya di ilmu lingkungan dan ekologi ([5]; [19]), serta di bidang kesehatan [1]. 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 [4]. Menelaah kemiskinan
secara multidimensional sangat diperlukan untuk perumuskan kebijakan pengentasan kemiskinan [20].
Berdasarkan uraian di atas, maka diperlukan suatu penelitian untuk mengetahui faktor-faktor yang mempengaruhi tingkat kemiskinan pada kabupaten / kota di Pulau Jawa, yang lebih lanjut dapat digunakan sebagai dasar kebijakan dalam usaha mengatasi kemiskinan. Untuk itu dalam penelitian ini, akan menggunakan pendekatan regresi kuantil pada kasus kemiskinan di Pulau Jawa. Metode ini memiliki kelebihan dapat mengatasi masalah pencilan (outlier) dan heteroskedastisitas. Sehingga
diharapkan memberikan taksiran parameter yang lebih stabil dan konsisten dengan menduga berbagai fungsi kuantil dari suatu distribusi Y sebagai fungsi dari X. Kajian akan difokuskan pada pendugaan selang kepercayaan menggunakan resampling
dengan teknik
bootstrap . Metode bootstrap telah dikembangkan oleh [8]
sebagai alat untuk membantu mengurangi ketidakandalan yang berhubungan
dengan kesalahan penggunaan distribusi normal dan penggunaannya. Bootstrap membuat data bayangan (pseudo data) dengan menggunakan informasi dari data asli dengan memperhatikan sifat-sifat dari data asli, sehingga data bayangan memiliki karakteristik yang sangat mirip dengan data asli.

## METODE PENELITIAN

## Sumber Data dan Variabel Penelitian
Data yang digunakan pada penelitian ini adalah data sekunder berupa publikasi “Data dan
Informasi Kemiskinan Kabupaten / Kota Tahun 2011” yang merupakan hasil olah data triwulanan Survei Sosial Ekonomi Nasional (SUSENAS) yang diselenggarakan oleh Badan Pusat Statistik (BPS). Unit observasi dalam penelitian ini adalah
faktor-faktor
yang mempengaruhi
kemiskinan masing-masing provinsi,
yaitu Kabupaten / Kota di setiap Provinsi yang terletak di Pulau Jawa. Provinsi DKI Jakarta (6 kota), Jawa Barat (26 Kabupaten / Kota), Jawa Tengah (35 Kabupaten / Kota), Daerah Istimewa Yogyakarta (5 Kabupaten / Kota), Jawa
Timur (38 Kabupaten / Kota) dan Banten (8 Kabupaten / Kota). Enam Provinsi di Pulau Jawa tersebut selanjutnya digunakan sebagai unit analisis dalam regresi kuantil dengan pendekatan
bootstrap untuk mengetahui faktor-faktor yang mempengaruhi kemiskinan di Pulau Jawa.
Mengacu pada dasar-dasar analisis dan identifikasi kemiskinan BPS maka variabel-variabel yang digunakan dalam penelitian ini adalah indeks kedalaman kemiskinan (Y) sebagai variabel respon. Variabel prediktor yang terdiri dari indikator kesehatan antara lain Persentase perempuan pengguna alat KB di rumah tangga miskin (X 1 ), Persentase balita di rumah tangga miskin yang proses kelahirannya ditolong oleh tenaga kesehatan (X 2 ), Persentase balita di rumah tangga miskin yang telah diimunisasi (X 3 ), Persentase rumah tangga miskin dengan luas lantai perkapita ≤ 8 m 2 (X 4 ), Persentase rumah tangga miskin yang menggunakan air bersih (X 5 ), Persentase rumah tangga miskin yang menggunakan jamban sendiri / bersama (X 6 ), Persentase rumah tangga miskin yang mendapatkan pelayanan jamkesmas (X 7 ). Indikator SDM antara lain Persentase penduduk miskin usia 15 tahun keatas yang tidak tamat SD (X 8 ), Angka Melek Huruf (AMH) penduduk miskin usia 15-55 tahun (X 9 ), Angka Partisipasi Sekolah (APS) penduduk miskin usia 13-15 tahun (X 10 ), Indikator ekonomi antara lain Persentase penduduk miskin usia 15 tahun keatas yang tidak bekerja (X 11 ), Persentase penduduk miskin usia 15 tahun keatas yang bekerja di sektor
pertanian (X 12 ), Persentase rumah tangga yang pernah membeli beras raskin (X Persentase pengeluaran per kapita untuk non makanan (X 14 ).

## Metode Analisis
Langkah-langkah analisis tersebut dapat pula digambarkan dalam d alir seperti terdapat pada Gambar 1 :
Gambar 1. Diagram alir penelitian

## HASIL PENELITIAN
Pada bab ini akan diuraikan pembahasan pemodelan kemiskinan di Pulau Jawa dengan menggunakan regresi kuantil pendekatan bootstrap . Beberapa tahapan yang disajikan dalam bab ini diantaranya adalah analisis deskriptif dan variabel-variabel kemiskinan, pendugaan parameter dan interval prediksi dengan
regresi klasik OLS, melakukan pemeriksaan terhadap pencilan atau outlier da
), Persentase rumah tangga yang pernah membeli beras raskin (X 13 ), Persentase pengeluaran per kapita untuk
langkah analisis tersebut dapat pula digambarkan dalam diagram alir seperti terdapat pada Gambar 1 : Diagram alir penelitian akan diuraikan hasil dan kemiskinan di Pulau Jawa dengan menggunakan regresi kuantil
. Beberapa tahapan yang disajikan dalam bab ini diantaranya adalah analisis deskriptif dan scaterplot
variabel kemiskinan, pendugaan parameter dan interval prediksi dengan
regresi klasik OLS, melakukan pemeriksaan terhadap pencilan atau outlier dan
selanjutnya mendapatkan pemodelan dengan regresi kuantil pendekatan bootstrap Statistik Deskriptif Berdasarkan boxplot pada Gambar 2.diatas dapat diketahui bahwa terdapat keragaman yang cukup tinggi di beberapa indikator kemiskinan di Pulau Jawa.
terdapat kesenjangan antara di Pulau Jawa, sesuai dengan karakteristik wilayah itu sendiri.
Gambar 2. Boxplot Indikator Kemiskinan
Pulau Jawa Regresi Klasik Hasil regresi OLS pada diketahui bahwa terdapat beberapa indikator kemiskinan yang berpengaruh signifikan pada tingkat signifikansi 0,05. Idenstifikasi Pencilan / Outlier Berdasarkan pada Gambar 3.dan 4 dapat diketahui bahwa scatterplot variabel respon indeks kedalaman kemiskinan (Y) dengan 14 (empat belas) indikator kemiskinan (X menunjukkan bahwa ada beberapa data Kabupaten / Kota yang menyimpang dari sekumpulan data yang lainnya atau terdapat data yang jauh berbeda dibandingkan dengan keseluruhan data
Sehingga dapat disimpulkan secara visual berdasarkan pada scatterplot variabel respon indeks kedalaman kemiskinan (Y) dengan 14 (empat belas) indikator kemiskinan (X beberapa pengamatan yang
selanjutnya mendapatkan pemodelan dengan bootstrap .
Berdasarkan boxplot pada Gambar 2.diatas dapat diketahui bahwa terdapat keragaman yang cukup tinggi di beberapa indikator kemiskinan di Pulau Jawa. Artinya kesenjangan antara Kabupaten/ Kota Pulau Jawa, sesuai dengan karakteristik Boxplot Indikator Kemiskinan di Pulau Jawa Hasil regresi OLS pada Tabel 1.dapat diketahui bahwa terdapat beberapa indikator kemiskinan yang berpengaruh signifikan pada tingkat signifikansi ߙ =

## Outlier
Berdasarkan pada Gambar 3.dan 4 scatterplot antar variabel respon indeks kedalaman kemiskinan (Y) dengan 14 (empat belas) kemiskinan (X 1 -X 14 ) menunjukkan bahwa ada beberapa data Kabupaten / Kota yang menyimpang dari sekumpulan data yang lainnya atau data yang jauh berbeda ngkan dengan keseluruhan data.
Sehingga dapat disimpulkan secara visual scatterplot bahwa variabel respon indeks kedalaman kemiskinan (Y) dengan 14 (empat belas) kemiskinan (X 1 -X 14 ) terdapat beberapa pengamatan yang outlier.

## Ket : *Signifikan
Idenstifikasi Pencilan / Outlier
Gambar 3. Scatterplot antara Y dengan x 4 , x 5 , x 6
Gambar 4. Scatterplot antara Y dengan x 10 , x 11 , x 12 , x 13 , x . Deteksi Heterokesdatisitas Berdasarkan interval prediksi 95% pada Gambar 5.dan6., dapat diketahui bahwa terjadi kasus heteroskedastisitas. Hal itu dapat diketahui dari varians error konstan.Sebagai contoh, yaitu indikator persentase rumah tangga miskin yang
Outlier antara Y dengan x 1 , x 2 , x 3 ,
antara Y dengan x 7 , x 8 , x 9 , , x 14
Berdasarkan interval prediksi 95% pada Gambar 5.dan6., dapat diketahui bahwa terjadi kasus heteroskedastisitas. Hal itu error yang tidak konstan.Sebagai contoh, yaitu indikator ersentase rumah tangga miskin yang
menggunakan air bersih (X error yang membesar.
Gambar 5. Regresi Antara Y dengan X X 4 , X 5 dan X 6 berdasarkan Interval Prediksi 95%
Gambar 6. Regresi Antara Y dengan X
X 10 , X 11 , X 12 , X 13 Interval Prediksi
Regresi Kuantil Regresi kuantil dalam penelitian ini dilakukan pada kuantil ke- τ τ = 0,75 dan τ = 0,99 dengan 500 untuk masing-masing koefisien regresi pada data kemiskinan di Pulau Jawa.Berdasarkan pada tabel 2.dapat diketahui bahwa pada kuantil 0,25 dengan menggunakan bootstrap 500, standar error yang didapatkan relatif lebih kecil daripada kuantil
0,50 , τ = 0,75 dan τ = 0,99 kuantil dalam penelitian ini mengindika adanya multikolinieritas, hal tersebut dapat diketahuidari nilai pseudo R akan tetapi masih terdapat banyak variab yang tidak signifikan dalam model pada τ = 0,25
r bersih (X 5 ) dengan varians
Regresi Antara Y dengan X 1 , X 2 , X 3 , berdasarkan Interval Prediksi 95%
Regresi Antara Y dengan X 7 , X 8 , X 9 , dan X 14 berdasarkan Interval Prediksi 95%
Regresi kuantil dalam penelitian ini τ = 0,25 ; τ = 0,5 ;
dengan bootstrap masing koefisien regresi pada data kemiskinan di Pulau Berdasarkan pada tabel 2.dapat diketahui bahwa pada kuantil τ = dengan menggunakan resampling 500, standar error yang didapatkan relatif lebih kecil daripada kuantil τ = 99 . Model regresi kuantil dalam penelitian ini mengindikasikan hal tersebut dapat R 2 sebesar 61,4%, terdapat banyak variabel yang tidak signifikan dalam model pada ke-
Hasil Estimasi Regresi Kuantil pada Kuantil ke-0,25, 0,50, 0,75, 0,99 dengan
engan resampling 500

## KESIMPULAN
Hasil estimasi regresi kuantil dengan kuantil ke- ߬ = 0,25 pada resampling bootstrap 500 memiliki standard error yang relatif lebih kecil, dapat dijelaskan bahwa terdapat beberapa indikator yang bepengaruh signifikan terhadap indeks kedalaman kemiskinan, yaitu indikator persentase balita di rumah tangga miskin yang proses kelahirannya ditolong oleh tenaga kesehatan (X 2 ) berpengaruh positif dan signifikan terhadap model, indikator angka Melek Huruf (AMH) penduduk miskin usia 15-55 tahun (X 9 )
berpengaruh negatif dan signifikan
terhadap model, indikator Persentase penduduk miskin usia 15 tahun keatas yang tidak bekerja (X 11 ) berpengaruh positif dan signifikan terhadap model, indikator Persentase penduduk miskin usia 15 tahun keatas yang bekerja di sektor pertanian (X 12 ) berpengaruh positif dan signifikan terhadap model.
ܳ ଴ , ଶହ ( ݕ | ݔ ) = 8,4277 − 0,0097 ݔ ଵ + 0,0257 ݔ ଶ − 0,0109 ݔ ଷ + 0,0075 ݔ ସ − 0,0145 ݔ ହ + 0,0026 ݔ ଺ + 0,0084 ݔ ଻ − 0,0189 ݔ ଼ − 0,1164 ݔ ଽ − 0,0077 ݔ ଵ଴ + 0,0483 ݔ ଵଵ + 0,0375 ݔ ଵଶ + 0,0063 ݔ ଵଷ + 0,0353 ݔ ଵସ

## DAFTAR PUSTAKA
[1] Austin,P.,Tu,J.,Daly,P.,&Alter,D.,
2005, Theuse of quantile regression
inhealth care research: A case study examining gender differences in the timeliness of thrombolytictherapy, Statistics in
Medicine, 24, 791–816.
[2] Badan Pusat Statistik, 2012, Berita Resmi Statistik: Profil Kemiskinan di Indonesia Maret 2012 , BPS, Jakarta. [3] Badan Pusat Statistik, 2012, Penghitungan Dan Analisis
Kemiskinan Makro Indonesia 2012 , BPS, Jakarta. [4] BAPPENAS., 2010, Arah Kebijakan Penanggulangan Kemiskinan
[5] Cade,B.S.,Terrell,
J.W.,& Schroeder,R.L., 1999, Estimating effects of limiting factors with regression quantiles. Ecology, 80,
311–323. [6] Chamberlain,G., 1994, Quantile regression,censoring and the structure of wages.InC. Skins(Ed.), Advancesin
Econometrics, 171–209,
Cambridge,UK:Cambridge University Press.
[7] Drapper, N.R.,dan Smith, H., 1998, Applied Regression Analysis , John Wiley & Sons, Canada.
[8] Efron, B.dan Tibshirani, R.J., 1993, An Introduction to the Bootstrap , Chapman and Hall, Inc., New York.
[9] Gujarati,
D.N., 2004, Basic Econometric , Mc Graw Hill, New
York.
[10] Hao, L.,& Naiman, D. Q., 2007,
Quantile Regression . Sage Publications, Inc.
[11] Koenker,R. dan Basset, G., 1978,
Regression Quantiles, Econometrica .
January, 46:1, hal. 33-50.
[12] Koenker, R. dan D’Orey,V., 1993, Computing Regression Quantiles, . J.Roy . Statist.Soc.Ser. C
(Appl. Statist) . Vol. 43, pp. 410-414.
[13] Koenker, R. dan Hallock, K.F.,
2001, Quantile Regression. Journal of Economic Perspective. Vol.15,
No.4, hal. 143-156. Published by: AmericanEconomic Association.
[14] Koenker, R., 2005, Quantile Regression , Cambridge University Press, United Statesof America.
[15] Navianti, D.R., 2014, Regresi
Kuantil dengan Estimasi Metode Sparsity untuk Pemodelan Tingkat Pengangguran Terbuka, Tesis, Program Pasca sarjana Institut Teknologi Sepuluh Nopember, Surabaya. [16] Ngafiyah, A,N., 2014, Meta- Analisis dengan Pendekatan Two
Stage Structural Equation Modeling
(Studi Kasus Beberapa Penelitian
Faktor-faktor yang Mempengaruhi Kemiskinan di Pulau Jawa), Tesis, Program Pascasarjana Institut Teknologi Sepuluh Nopember, Surabaya. [17] Portnoy, S., 2003, Censored
Regression Quantiles, Journal of the American Statistical Association. Vol. 98, No. 464, hal.1001-1012.
[18] Rahmawati, R., 2011, Regresi
Kuantil (Studi Kasus pada Data Suhu Harian) , Skripsi, Program Studi Statistika FMIPA Undip [19] Scharf, F.S., Juanes, F., & Sutherland ,M., 1989, Inferringe cological relationships from the edges of scatter diagrams: Comparison of regression techniques, Ecology , 79, 448-460.
[20] Suryawati, C., 2005, Memahami
Kemiskinan Secara Multidimensional, Jurnal Manajemen Pelayanan Kesehatan (JMPK) , 8(3).
p.121-129.
[21] Wahyudi., V., E., 2014, Analisis
IPM Di Pulau Jawa Menggunakan Analisis Regresi Kuantil , Tesis, Program Pasca sarjana Institut Teknologi Sepuluh Nopember, Surab aya.