TUGAS AKHIR - SS141501
PEMODELAN FAKTOR-FAKTOR YANG MEMENGARUHI INDEKS PEMBANGUNAN KESEHATAN MASYARAKAT PROVINSI JAWA TIMUR MENGGUNAKAN PENDEKATAN REGRESI SEMIPARAMETRIK SPLINE MADE AYU DWI OCTAVANNY NRP 1313 100 012
Dosen Pembimbing Prof. Dr. Drs. I Nyoman Budiantara, M.Si. Co-Pembimbing Dr. Vita Ratnasari, S.Si, M.Si. PROGRAM STUDI S1 JURUSAN STATISTIKA FAKULTAS MATEMATIKA DAN ILMU PENGETAHUAN ALAM INSTITUT TEKNOLOGI SEPULUH NOPEMBER SURABAYA 2017
TUGAS AKHIR AKHIR -- SS141501 SS141501 TUGAS
PEMODELAN FAKTOR-FAKTOR YANG MEMENGARUHI PEMODELAN FAKTOR-FAKTOR YANG MEMENGARUHI INDEKS PEMBANGUNAN KESEHATAN MASYARAKAT INDEKS PEMBANGUNAN MANUSIA MENGGUNAKAN PROVINSI JAWA TIMUR MENGGUNAKAN PENDEKATAN REGRESI NONPARAMETRIK SPLINE DI JAWA TENGAH REGRESI SEMIPARAMETRIK SPLINE NI PUTU DERA YANTHI MADE AYU100 DWI040 OCTAVANNY NRP 1312 NRP 1313 100 012 Dosen Pembimbing Pembimbing Dosen Prof. Dr. I Nyoman Budiantara, Prof. Dr. IDrs. Nyoman Budiantara, M.Si. M.Si. Dr. Vita Ratnasari, S.Si, M.Si.
PROGRAM STUDI STUDI S1 S1 PROGRAM JURUSAN STATISTIKA JURUSAN STATISTIKA FAKULTAS MATEMATIKA MATEMATIKA DAN DAN ILMU ILMU PENGETAHUAN PENGETAHUAN ALAM ALAM FAKULTAS INSTITUT TEKNOLOGI SEPULUH NOPEMBER INSTITUT TEKNOLOGI SEPULUH NOPEMBER SURABAYA 2016 2017 SURABAYA
FINAL PROJECT - SS141501
MODELING FACTORS THAT INFLUENCE PUBLIC HEALTH DEVELOPMENT INDEX AT EAST JAVA USING SEMIPARAMETRIC SPLINE REGRESSION APPROACH
MADE AYU DWI OCTAVANNY NRP 1313 100 012
Supervisor Prof. Dr. Drs. I Nyoman Budiantara, M.Si. Dr. Vita Ratnasari, S.Si, M.Si.
UNDERGRADUATE PROGRAMME DEPARTMENT OF STATISTICS FACULTY OF MATHEMATICS AND NATURAL SCIENCE INSTITUT TEKNOLOGI SEPULUH NOPEMBER SURABAYA 2017
PEMODELAN FAKTOR-FAKTOR YANG MEMENGARUHI INDEKS PEMBANGUNAN KESEHATAN MASYARAKAT PROVINSI JAWA TIMUR MENGGUNAKAN PENDEKATAN REGRESI SEMIPARAMETRIK SPLINE Nama NRP Jurusan Dosen Pembimbing
: Made Ayu Dwi Octavanny : 1313 100 012 : Statistika : 1. Prof. Dr. Drs. I Nyoman Budiantara, M. Si. 2. Dr. Vita Ratnasari, S. Si, M. Si.
ABSTRAK Aspek kesehatan yang dilihat dari umur harapan hidup saat kelahiran belum dapat menjabarkan indeks kesehatan dengan baik sehingga diperlukan indikator yang dapat mengukur aspek kesehatan, yaitu Indeks Pembangunan Kesehatan Masyarakat (IPKM). IPKM dapat digunakan untuk mengetahui sejauh mana keberhasilan pembangunan kesehatan masyarakat suatu wilayah. Hampir semua kabupaten/kota di Provinsi Jawa Timur mengalami penurunan peringkat berdasarkan pengembangan IPKM 2013. Hal ini perlu menjadi perhatian bagi pemerintah untuk mengetahui faktor-faktor yang memengaruhi IPKM di Provinsi Jawa Timur sehingga pemerintah dapat melakukan usaha perbaikan kualitas kesehatan masyarakat. Pada penelitian ini, IPKM dan faktor-faktor yang memengaruhinya sebagian membentuk pola tertentu dan sebagian lagi tidak membentuk pola tertentu, sehingga digunakan metode regresi semiparametrik spline. Model terbaik didapatkan dari titik knot optimum berdasarkan nilai Generalized Cross Validation (GCV) minimum. Model regresi semiparametrik spline terbaik adalah dengan menggunakan tiga titik knot. Terdapat empat variabel yang signifikan, yaitu Angka Kematian Bayi, kepadatan penduduk, persentase rumah tangga berperilaku hidup bersih dan sehat, dan persentase rumah sehat. Nilai koefisien determinasi yang dihasilkan dari model ini adalah sebesar 92,99%. Kata kunci: Indeks Pembangunan Kesehatan Masyarakat, Jawa Timur, Regresi Semiparametrik, Spline Truncated
vii
(Halaman ini sengaja dikosongkan)
viii
MODELING FACTORS THAT INFLUENCE PUBLIC HEALTH DEVELOPMENT INDEX AT EAST JAVA USING SEMIPARAMETRIC SPLINE REGRESSION APPROACH Name NRP Department Supervisor
: Made Ayu Dwi Octavanny : 1313 100 012 : Statistics : 1. Prof. Dr. Drs. I Nyoman Budiantara, M. Si. 2. Dr. Vita Ratnasari, S.Si, M.Si. ABSTRACT
Health aspect which is seen from life expectancy has not been able to describe health index so well, so it needs an indicator that can measure the health aspect, namely Public Health Development Index (PHDI). PHDI can be used to determine the success of public health development of a region. Almost all districts/cities at East Java had decreased rank in PHDI in 2013. This development should be a concern for the government to know the factors that influence health index at East Java so that the government can do quality improvements in public health. In this research, PHDI and the factors that influence it partially forms a particular pattern and some does not form a specific pattern, so that the semiparametric spline regression is used. The best model is obtained from the optimum knots based on the minimum Generalized Cross Validation (GCV). The best semiparametric spline regression model is using three knots. There are four significant variables, namely Infant Mortality Rate, population density, the percentage of households with clean and healthy living behavior, and the percentage of healthy house. The coefficient of determination of this model is 92,99%. Keywords: East Java, Public Health Development Index, Semiparametric Regression, Spline Truncated
ix
(This page intentionally left blank)
x
KATA PENGANTAR Puji syukur penulis panjatkan kepada Ida Sang Hyang Widhi Wasa, Tuhan Yang Maha Esa atas limpahan rahmat-Nya sehingga penulis dapat menyelesaikan Tugas Akhir dengan judul “PEMODELAN FAKTOR-FAKTOR YANG MEMENGARUHI INDEKS PEMBANGUNAN KESEHATAN MASYARAKAT PROVINSI JAWA TIMUR MENGGUNAKAN PENDEKATAN REGRESI SEMIPARAMETRIK SPLINE”. Penulis menyadari bahwa dalam penyusunan Tugas Akhir ini tidak terlepas dari bantuan dan dukungan dari berbagai pihak. Oleh sebab itu, penulis mengucapkan terima kasih kepada: 1. Dr. Suhartono, selaku Ketua Jurusan Statistika FMIPA ITS. 2. Dr. Sutikno, S.Si, M.Si. selaku Ketua Prodi S1 Jurusan Statistika FMIPA ITS. 3. Prof. Dr. Drs. I Nyoman Budiantara, M.Si. dan Dr. Vita Ratnasari, S.Si, M.Si. selaku dosen pembimbing yang dengan sabar telah memberikan ilmu, waktu, dan pengarahan kepada penulis. 4. Dra. Madu Ratna, M.Si. dan Dr. Purhadi, M.Sc. selaku dosen penguji yang telah memberikan kritik dan saran demi perbaikan Tugas Akhir ini. 5. Dr. Ir. Setiawan, M.S. selaku dosen wali yang telah meluangkan waktu untuk membimbing penulis selama masa perkuliahan. 6. Seluruh bapak dan ibu dosen Statistika ITS yang telah memberikan ilmu, pengetahuan, dan bimbingan yang tak ternilai harganya, serta segenap seluruh karyawan Jurusan Statistika ITS atas bantuan terkait akademik. 7. Kedua orang tua penulis dan keluarga besar yang selalu memberikan doa, kasih sayang, bimbingan, semangat, dan motivasi serta kesabaran selama ini. 8. Teman seperjuangan spline: Krisna dan Dhira atas diskusi tugas akhir selama satu semester.
xi
9.
Teman perantauan terbaik dan sahabat-sahabat penulis: Alit, Fadhila, Halimah, Yulinda, Siska, Ajeng, Farah, Dina, atas semua keceriaan, persahabatan, dan motivasi selama ini. 10. Teman-teman TPKH-ITS, khususnya angkatan 2013 yang tidak bisa disebutkan satu persatu atas kebersamaannya selama kuliah di rantauan. 11. Seluruh keluarga besar Jurusan Statistika FMIPA ITSSurabaya, khususnya 24 atas kebersamaananya selama ini. 12. Serta pihak-pihak lain yang tidak dapat penulis sebutkan satu persatu. Penulis menyadari masih banyak kekurangan dalam laporan ini, maka saran dan kritik dari pembaca demi perbaikan sangat dibutuhkan oleh penulis. Semoga Tugas Akhir ini dapat bermanfaat bagi pembaca dan beberapa pihak.
Surabaya, Januari 2017
Penulis
xii
DAFTAR ISI Halaman HALAMAN JUDUL ................................................................ i COVER PAGE ........................................................................ iii LEMBAR PENGESAHAN ................................................... v ABSTRAK ............................................................................ vii ABSTRACT............................................................................. ix KATA PENGANTAR ........................................................... xi DAFTAR ISI ........................................................................ xiii DAFTAR TABEL ............................................................... xvii DAFTAR GAMBAR .......................................................... xix DAFTAR LAMPIRAN ....................................................... xxi BAB I PENDAHULUAN 1.1 Latar Belakang ............................................................... 1 1.2 Rumusan Masalah .......................................................... 4 1.3 Tujuan Penelitian ........................................................... 4 1.4 Manfaat Penelitian ......................................................... 4 1.5 Batasan Masalah ............................................................ 5 BAB II TINJAUAN PUSTAKA 2.1 Statistika Deskriptif ........................................................ 7 2.2 Pengujian Nonlinieritas .................................................. 7 2.3 Analisis Regresi ............................................................. 8 2.3.1 Regresi Parametrik ............................................... 9 2.3.2 Regresi Nonparametrik ...................................... 10 2.3.3 Regresi Semiparametrik ..................................... 10 2.4 Pemodelan Regresi Spline Truncated .......................... 11 2.5 Pemilihan Titik Knot Optimal ...................................... 13 2.6 Pengujian Parameter Model ......................................... 13 2.6.1 Pengujian Secara Serentak ................................. 14 2.6.2 Pengujian Secara Parsial .................................... 15 2.7 Kriteria Kebaikan Model Regresi ................................ 16 2.8 Pengujian Asumsi Residual .......................................... 17
xiii
2.8.1 Uji Asumsi Identik ............................................. 17 2.8.2 Asumsi Independen............................................ 17 2.8.3 Uji Asumsi Berdistribusi Normal ...................... 18 2.9 Indeks Pembangunan Kesehatan Masyarakat .............. 18 2.10 Kerangka Konsep Penelitian ........................................ 19 11 Penitian Sebelumnya ......................................................... 18 BAB III METODOLOGI PENELITIAN 3.1 Sumber Data ................................................................. 21 3.2 Variabel Penelitian ....................................................... 21 3.3 Langkah Penelitian ....................................................... 24 3.4 Diagram Alir ................................................................ 24 BAB IV ANALISIS DAN PEMBAHASAN 4.1 Karakteristik Indeks Pembangunan Kesehatan Masyarakat Provinsi Jawa Timur ................................. 27 4.2 Pemodelan Indeks Pembangunan Kesehatan Masyarakat Provinsi Jawa Timur Menggunakan Pendekatan Regresi Semiparametrik Spline .................................... 32 4.2.1 Pola Hubungan IPKM Provinsi Jawa Timur dengan Faktor-Faktor yang Diduga Memengaruhi 33 4.2.2 Pemilihan Titik Knot Optimum ......................... 37 1. Pola Berdasarkan Scatterplot ........................ 37 2. Pola Berdasarkan Uji RESET ....................... 43 4.2.3 Pemilihan Model Terbaik .................................. 47 4.2.4 Pemodelan IPKM Provinsi Jawa Timur dengan Menggunakan Titik Knot Optimum ................... 47 4.2.5 Pengujian Signifikansi Parameter Model Regresi Semiparametrik Spline ................................... 48 1. Pengujian Serentak ....................................... 48 2. Pengujian Individu ........................................ 49 4.2.6 Pemodelan IPKM Provinsi Jawa Timur dengan Empat Variabel .................................................. 51 1. Pemilihan Titik Knot dengan Satu Titik Knot .............................................................. 51
xiv
2. Pemilihan Titik Knot dengan Dua Titik Knot .............................................................. 52 3. Pemilihan Titik Knot dengan Tiga Titik Knot .............................................................. 53 4. Pemilihan Titik Knot dengan Kombinasi Titik Knot ......................................................... 54 4.2.7 Pemilihan Model IPKM Provinsi Jawa Timur Terbaik dengan Empat Variabel ........................ 55 4.2.8 Pemodelan IPKM Provinsi Jawa Timur Menggunakan Titik Knot Optimum dengan Empat Variabel .................................................. 55 4.2.9 Pengujian Signifikansi Parameter Model Regresi Semiparametrik Spline dengan Empat Variabel ...................................................................... 56 1. Pengujian Serentak ....................................... 56 2. Pengujian Individu ........................................ 57 4.2.10 Pengujian Asumsi Residual ............................... 58 1. Pengujian Asumsi Identik ............................. 58 2. Asumsi Independen ...................................... 59 3. Pengujian Asumsi Berdistribusi Normal ...... 59 4.2.11 Interpretasi Model Regresi Semiparametrik Spline ................................................................. 60 BAB V KESIMPULAN DAN SARAN 5.1 Kesimpulan ................................................................ 67 5.2 Saran .......................................................................... 68 DAFTAR PUSTAKA ........................................................... 69 LAMPIRAN .......................................................................... 73 BIODATA PENULIS ......................................................... 119
xv
(Halaman ini sengaja dikosongkan)
xvi
DAFTAR TABEL Tabel 2. 1 Tabel 3. 1 Tabel 4. 1 Tabel 4. 2 Tabel 4. 3 Tabel 4. 4 Tabel 4. 5 Tabel 4. 6 Tabel 4. 7 Tabel 4. 8 Tabel 4. 9 Tabel 4. 10 Tabel 4. 11 Tabel 4. 12 Tabel 4. 13 Tabel 4. 14 Tabel 4. 15 Tabel 4. 16 Tabel 4. 17
Halaman Analisis Ragam (ANOVA) ................................... 14 Variabel Penelitian ................................................ 21 Statistika Deskriptif Variabel Penelitian ............... 29 Uji RESET Antara IPKM dengan Faktor yang Memengaruhi ........................................................ 37 Nilai GCV dengan Satu Titik Knot Pola Berdasarkan Scatterplot ...................................................... 38 Nilai GCV dengan Dua Titik Knot Pola Berdasarkan Scatterplot ...................................................... 39 Nilai GCV dengan Tiga Titik Knot Pola Berdasarkan Scatterplot ...................................................... 41 Nilai GCV dengan Kombinasi Titik Knot Pola Berdasarkan Scatterplot ........................................ 42 Nilai GCV dengan Satu Titik Knot Pola Berdasarkan Uji RESET ..................................................... 43 Nilai GCV dengan Dua Titik Knot Pola Berdasarkan Uji RESET ..................................................... 44 Nilai GCV dengan Tiga Titik Knot Pola Berdasarkan Uji RESET ..................................................... 45 Nilai GCV dengan Kombinasi Titik Knot Pola Berdasarkan Uji RESET ....................................... 46 Perbandingan Nilai GCV Minimum ..................... 47 Hasil Pengujian Serentak ...................................... 48 Hasil Pengujian Individu....................................... 50 Nilai GCV dengan Satu Titik Knot Pada Model Empat Variabel ..................................................... 51 Nilai GCV dengan Dua Titik Knot Pada Model Empat Variabel ..................................................... 52 Nilai GCV dengan Tiga Titik Knot Pada Model Empat Variabel ..................................................... 53 Nilai GCV dengan Kombinasi Titik Knot Pada Model Empat Variabel .......................................... 54
xvii
Tabel 4. 18 Perbandingan Nilai GCV Minimum Pada Model Empat Variabel ..................................................... 55 Tabel 4. 19 Hasil Pengujian Serentak Pada Model Empat Variabel................................................................. 56 Tabel 4. 20 Hasil Pengujian Individu Pada Model Empat Variabel................................................................. 57 Tabel 4. 21 Hasil Pengujian Glejser ........................................ 58
xviii
DAFTAR GAMBAR Halaman Diagram Alir .................................................... 25 IPKM di Kabupaten/Kota Provinsi Jawa Timur ......................................................................... 28 Gambar 4. 2 Pola Hubungan IPKM dengan Persentase Penduduk Miskin ............................................. 32 Gambar 4. 3 Pola Hubungan IPKM dengan Angka Kematian Bayi ........................................................... 33 Gambar 4. 4 Pola Hubungan IPKM dengan Kepadatan Penduduk ......................................................... 34 Gambar 4. 5 Pola Hubungan IPKM dengan Angka KemaTian Ibu ............................................................ 34 Gambar 4. 6 Pola Hubungan IPKM dengan Persentase Rumah Tangga Berperilaku Hidup Bersih dan Sehat ................................................................ 35 Gambar 4. 7 Pola Hubungan IPKM dengan Persentase Rumah Sehat......................................................... 36 Gambar 4. 8 Plot ACF Residual ........................................... 59 Gambar 4. 9 Hasil Uji Kolmogorov-Smirnov ....................... 60 Gambar 4. 10 Peta Kepadatan Penduduk di Provinsi Jawa Timur ............................................................... 61 Gambar 4. 11 Peta Persentase Rumah Tangga Ber-PHBS di Provinsi Jawa Timur ........................................ 64 Gambar 4. 12 Peta Persentase Rumah Sehat di Provinsi Jawa Timur ............................................................... 65 Gambar 3. 1 Gambar 4. 1
xix
(Halaman ini sengaja dikosongkan)
xx
DAFTAR LAMPIRAN Lampiran 1. Lampiran 2. Lampiran 3. Lampiran 4. Lampiran 5.
Lampiran 6.
Lampiran 7. Lampiran 8. Lampiran 9. Lampiran 10. Lampiran 11. Lampiran 12. Lampiran 13. Lampiran 14. Lampiran 15.
Halaman Data IPKM Provinsi Jawa Timur dengan Faktor-Faktor yang Memengaruhi Tahun 2013 ..... 73 Program Pemilihan Titik Knot Optimal dengan Satu Titik Knot Menggunakan Software R ...... 75 Program Pemilihan Titik Knot Optimal dengan Dua Titik Knot Menggunakan Software R ...... 78 Program Pemilihan Titik Knot Optimal dengan Tiga Titik Knot Menggunakan Software R...... 81 Program Pemilihan Titik Knot Optimal dengan Kombinasi Titik Knot Pada Pola Berdasarkan Scatterplot Menggunakan Software R ............. 84 Program Pemilihan Titik Knot Optimal dengan Kombinasi Titik Knot Pada Pola Berdasarkan Uji RESET Menggunakan Software R ............ 90 Program Pengujian Parameter Menggunakan Software R ....................................................... 94 Program Uji Glejser dengan Tiga Titik Knot Menggunakan Software R................................ 98 Output Nilai GCV dengan Satu Titik Knot Pada Pola Berdasarkan Scatterplot ................ 100 Output Nilai GCV dengan Dua Titik Knot Pada Pola Berdasarkan Scatterplot ................ 101 Output Nilai GCV dengan Tiga Titik Knot Pada Pola Berdasarkan Scatterplot ................ 102 Output Nilai GCV dengan Kombinasi Titik Knot Pola Berdasarkan Scatterplot ................ 103 Output Nilai GCV dengan Satu Titik Knot Pada Pola Berdasarkan Uji RESET ............... 104 Output Nilai GCV dengan Dua Titik Knot Pada Pola Berdasarkan Uji RESET ............... 105 Output Nilai GCV dengan Tiga Titik Knot Pada Pola Berdasarkan Uji RESET ............... 106
xxi
Lampiran 16. Output Nilai GCV dengan Kombinasi Titik Knot Pada Pola Berdasarkan Uji RESET ...... 107 Lampiran 17. Output Nilai GCV dengan Satu Titik Knot Pada Model Empat Variabel .......................... 108 Lampiran 18. Output Nilai GCV dengan Dua Titik Knot Pada Model Empat Variabel .......................... 109 Lampiran 19. Output Nilai GCV dengan Tiga Titik Knot Pada Model Empat Variabel .......................... 110 Lampiran 20. Output Nilai GCV dengan Kombinasi Titik Knot Model Empat Variabel .......................... 111 Lampiran 21. Output Estimasi Parameter dan Uji Signifikansi .............................................................. 112 Lampiran 22. Output Residual ............................................. 114 Lampiran 23. Output Uji Gljeser.......................................... 115 Lampiran 24. Surat Pernyataan Data Sekunder.................... 117
xxii
BAB I PENDAHULUAN 1.1 Latar Belakang Tercapainya pembangunan sumber daya manusia ditandai dengan tingginya kualitas hidup yang dicapai masyarakat. Salah satu standar yang ditetapkan untuk mengukur sejauh mana keberhasilan pembangunan manusia adalah Indeks Pembangunan Manusia (IPM) atau Human Development Index (HDI). Dalam mengukur pencapaian pembangunan manusia, IPM dibentuk melalui pendekatan tiga dimensi dasar, yaitu lamanya hidup, pengetahuan, dan penghidupan yang layak (Badan Pusat Statistik, 2013). Dalam aspek pengetahuan, indeks pengetahuan dapat diukur dengan adanya program wajib belajar 9 tahun, sedangkan dalam aspek ekonomi, peningkatan ekonomi dapat digunakan untuk mengukur indeks ekonomi. Aspek kesehatan yang dilihat dari umur harapan hidup saat kelahiran belum dapat menjabarkan indeks kesehatan dengan baik. Oleh karena itu, diperlukan indikator yang dapat mengukur pembangunan kesehatan, yaitu Indeks Pembangunan Kesehatan Masyarakat (IPKM). IPKM merupakan indikator komposit yang menggambarkan kemajuan pembangunan kesehatan yang dirumuskan dari 24 indikator kesehatan. Indeks tersebut disusun oleh Badan Penelitian dan Pengembangan Kesehatan (Balitbangkes) Kementerian Kesehatan RI. IPKM didasarkan pada data Riset Kesehatan Dasar (Riskesdas), Survei Sosial Ekonomi Nasional (Susenas), dan Survei Potensi Desa (Podes). Kemudian dilakukan pengembangan model IPKM 2013 dengan menggunakan 30 indikator yang dapat digunakan sebagai rumus standar untuk IPKM yang akan datang (Tim Penyusun IPKM, 2014). Selain untuk menentukan peringkat provinsi dan kabupaten/kota dalam keberhasilan pembangunan kesehatan masyarakat, IPKM juga dapat dimanfaatkan sebagai bahan advokasi ke pemerintah sehingga sumber daya dan program kesehatan diprioritaskan serta sebagai salah satu kriteria penentuan
1
2 alokasi dana bantuan kesehatan pusat ke daerah (Kementerian Kesehatan RI, 2010). Pada tahun 2000, Pemerintah Indonesia turut serta dalam deklarasi Millenium Development Goals (MDGs). Kesepakatan tersebut memiliki tujuan untuk meningkatkan kesejahteraan masyarakat, baik untuk negaranya sendiri maupun masyarakat dunia. MDGs memiliki delapan tujuan yang beberapa diantaranya terkait dengan aspek kesehatan, yaitu menurunkan kematian anak, meningkatkan kesehatan ibu, dan pengendalian penyakit menular. Tujuan-tujuan tersebut terkandung ke dalam indikator penyusun IPKM, sehingga dengan melakukan peningkatan indeks kesehatan masyarakat, akan memengaruhi ketercapaian MDGs. Indeks kesehatan Indonesia menduduki peringkat ke 108 dari seluruh negara di dunia menurut Human Development Reports pada tahun 2013. Hampir semua kabupaten/kota di Provinsi Jawa Timur mengalami penurunan peringkat berdasarkan pengembangan IPKM 2013 (Tim Penyusun IPKM, 2014). Oleh karena itu, perlu diketahui faktor-faktor yang memengaruhi IPKM di Provinsi Jawa Timur sehingga pemerintah dapat melakukan upaya-upaya untuk meningkatkan IPKM sebagai usaha perbaikan kualitas kesehatan masyarakat. Metode yang dapat digunakan dalam menyelesaikan masalah pemodelan adalah analisis regresi. Analisis regresi adalah suatu metode statistika yang digunakan untuk mengetahui pola hubungan antara variabel respon dengan prediktor (Drapper & Smith, 1992). Terdapat tiga pendekatan dalam analisis regresi, yaitu regresi parametrik, regresi nonparametrik, dan regresi semiparametrik. Regresi parametrik digunakan jika bentuk kurva regresi diketahui, sedangkan regresi nonparametrik digunakan jika bentuk kurva regresi diasumsikan tidak diketahui (Budiantara, 2005). Regresi semiparametrik adalah gabungan dari regresi parametrik dan regresi nonparametrik, dimana sebagian bentuk kurva regresi diketahui dan sebagian lagi tidak diketahui. Dengan melihat scatterplot, pola data IPKM dan faktor-faktor yang memengaruhinya di Provinsi Jawa Timur sebagian mengikuti pola
3 tertentu dan sebagian lagi tidak mengikuti pola tertentu, sehingga digunakan regresi semiparametrik spline. Hubungan antara IPKM dengan variabel jumlah keluarga miskin dan Angka Kematian Bayi cenderung mengikuti pendekatan parametrik. Semakin meningkatnya jumlah keluarga miskin dan angka kematian bayi, maka IPKM akan semakin rendah, begitupun sebaliknya. Apabila persentase penduduk miskin dan Angka Kematian Bayi berkurang, maka IPKM akan meningkat. Sedangkan variabel nonparametrik adalah kepadatan penduduk, Angka Kematian Ibu, persentase rumah tangga berperilaku hidup bersih dan sehat, serta persentase rumah sehat. Semakin meningkatnya variabel tersebut, belum tentu akan meningkatkan maupun mengurangi IPKM. Sebaliknya, semakin berkurangnya variabel tersebut belum tentu pula akan meningkatkan atau mengurangi IPKM. Penelitian mengenai IPKM sebelumnya telah dilakukan oleh Riskiyanti (2010) di Provinsi Jawa Timur menggunakan regresi multivariat yang menghasilkan bahwa faktor yang berpengaruh adalah persentase persalinan yang dilakukan oleh tenaga medis dan persentase imunisasi lengkap. Berikutnya, penelitian oleh Prasetyo (2012) menghasilkan bahwa faktor-faktor yang memengaruhi data kesehatan Kabupaten Banyuwangi dengan regresi terboboti geografis adalah fasilitas kesehatan, kepadatan penduduk, dan jumlah keluarga miskin. Selanjutnya Maully (2014) melakukan penelitian mengenai faktor-faktor yang memengaruhi indeks kesehatan kabupaten dan kota di Provinsi Jawa Timur dengan regresi logistik biner. Faktor-faktor yang memengaruhi adalah persentase pertolongan pertama kelahiran pada ibu dan persentase bayi diberi imunisasi. Penelitian yang dilakukan oleh Diansuantari (2015) mengenai derajat kesehatan masyarakat Provinsi Bali dengan metode MARS menggunakan variabel Angka Kematian Bayi, jumlah kasus TB Paru, persentase balita gizi buruk, dan Angka Kematian Ibu. Akan tetapi, belum ada penelitian mengenai pemodelan Indeks Pembangunan Kesehatan Masyarakat (IPKM) di Provinsi Jawa Timur dengan metode analisis regresi semiparametrik spline.
4 1.2 Rumusan Masalah Berdasarkan latar belakang yang telah dijelaskan mengenai IPKM Provinsi Jawa Timur, maka permasalahan yang akan dibahas pada penelitian ini adalah karakteristik Indeks Pembangunan Kesehatan Masyarakat (IPKM) di Provinsi Jawa Timur dan pemodelan Indeks Pembangunan Kesehatan Masyarakat (IPKM) di Provinsi Jawa Timur menggunakan pendekatan regresi semiparametrik spline. 1.3 Tujuan Penelitian Tujuan yang ingin dicapai pada penelitian ini adalah sebagai berikut. 1. Mendeskripsikan karakteristik Indeks Pembangunan Kesehatan Masyarakat (IPKM) di Provinsi Jawa Timur. 2. Memodelkan faktor-faktor yang memengaruhi Indeks Pembangunan Kesehatan Masyarakat (IPKM) di Provinsi Jawa Timur menggunakan pendekatan regresi semiparametrik spline. 1.4 Manfaat Penelitian Manfaat yang dapat diperoleh dari hasil penelitian ini adalah sebagai berikut. 1. Memberikan informasi kepada pemerintah Provinsi Jawa Timur mengenai faktor-faktor yang memengaruhi Indeks Pembangunan Kesehatan Masyarakat (IPKM) sehingga dapat dijadikan pertimbangan dalam menentukan kebijakan maupun program-program yang terkait dengan pembangunan kesehatan masyarakat. 2. Mampu menerapkan ilmu regresi semiparametrik spline dalam bentuk nyata pada bidang sosial dan kependudukan. 3. Menjadi acuan dalam penelitian selanjutnya.
5 1.5 Batasan Masalah Dalam penelitian ini perlu diberikan batasan permasalahan agar penelitian yang akan dikerjakan lebih fokus dan sesuai dengan rentang waktu yang direncanakan. Batasan masalah dalam penelitian ini adalah sebagai berikut. 1. Data yang digunakan adalah data sekunder dari Badan Pusat Statistik (BPS) Provinsi Jawa Timur dan Profil Kesehatan Jawa Timur yang diterbitkan oleh Dinas Kesehatan Provinsi Jawa Timur tahun 2013. 2. Model spline yang digunakan adalah spline truncated linier. 3. Banyak knot yang digunakan adalah satu knot, dua knot, tiga knot, dan kombinasi knot. 4. Pemilihan titik knot menggunakan metode GCV (Generalized Cross Validation).
6
(Halaman ini sengaja dikosongkan)
BAB II TINJAUAN PUSTAKA 2.1 Statistika Deskriptif Statistika Deskriptif merupakan metode yang berkaitan dengan pengumpulan, penyusunan dan penyajian suatu gugus data (Walpole, 1995). Statistika Deskriptif hanya memberikan informasi mengenai data itu sendiri dan sama sekali tidak menarik inferensia atau kesimpulan apapun dari gugus data induknya yang lebih besar. Dengan Statistika Deskriptif, kumpulan data yang diperoleh akan tersaji dengan ringkas dan rapi serta dapat memberikan informasi inti dari kumpulan data yang ada. Pada penelitian ini, Statistika Deskriptif yang digunakan adalah scatterplot atau diagram pencar. Diagram pencar digunakan untuk menjelaskan hubungan antara dua variabel (Nasution, 2005). Hubungan antara dua variabel yang berbentuk garis lurus disebut dengan hubungan linier. Apabila plot cenderung ke atas maka data memiliki hubungan yang positif, sebaliknya plot yang cenderung ke bawah menggambarkan hubungan yang negatif. Jika plot menyebar secara acak tanpa condong ke arah tertentu, maka data tidak memiliki hubungan yang linier. 2.2 Pengujian Nonlinieritas Untuk mengetahui apakah dua variabel memiliki hubungan linier atau nonlinier, maka dapat dilakukan dengan pengujian nonlinieritas. Salah satu pengujian nonlinieritas adalah uji Ramsey’s RESET. RESET test pertama kali diperkenalkan oleh Ramsey pada tahun 1969 (Kim et al., 2004). Pengujian ini dilakukan untuk membantu pembuatan keputusan apakah kedua variabel memiliki hubungan linier atau nonlinier yang seringkali belum dapat diputuskan jika hanya berdasarkan scatterplot. Hipotesis yang digunakan dalam pengujian ini adalah sebagai berikut.
7
8 H0 : 1 0 atau model sesuai spesifikasi (model linier) H1 : 1 0 atau terdapat mispesifikasi pada model (model nonlinier) Berikut merupakan langkah-langkah pengujian RESET (Gujarati & Porter, 2009). 1. Melakukan regresi sederhana antara variabel respon dengan variabel prediktor sehingga diperoleh nilai koefisien determinasi pertama yang dinotasikan dengan R2old dengan model regresi sebagai berikut. yi 0 1 x1i 2 x2i ... p x pi i , i 1, 2,3,..., n (2.1) 2.
Melakukan regresi seperti langkah sebelumnya dengan tambahan prediktor, yaitu estimasi nilai respon yˆ dari langkah pertama yang dikuadratkan yˆ 2 sehingga diperoleh
3.
4.
koefesien determinasi kedua yang dinotasikan dengan R2new. Model yang didapatkan adalah sebagai berikut. yi 0 1 x1i 2 x2i ... p x pi 1 yˆi2 i , i 1, 2,3,..., n (2.2) Melakukan perhitungan nilai statistik uji F dengan menggunakan rumus sebagai berikut. 2 2 Rnew Rold m (2.3) F 2 1 Rnew n k dengan m merupakan banyaknya prediktor tambahan (dalam penelitian ini hanya satu, yaitu suku kuadratik) dan k merupakan banyaknya parameter dalam model yang baru. Apabila nilai F F ,m , n k atau P-value < α maka tolak H0 yang dapat disimpulkan bahwa model nonlinier.
2.3 Analisis Regresi Menurut Drapper & Smith (1992), analisis regresi adalah suatu metode Statistika yang digunakan untuk mengetahui pola hubungan antara variabel respon dengan prediktor. Terdapat tiga pendekatan dalam analisis regresi, yaitu pendekatan regresi
9 parametrik, regresi nonparametrik, dan regresi semiparametrik. Hubungan yang terjadi antara variabel tersebut tidak selalu memiliki pola parametrik seperti linear, kuadratik, kubik, dan lainnya, namun banyak pula terdapat pola nonparametrik (Budiantara, 2005). Apabila bentuk kurva regresi diketahui maka digunakan pendekatan regresi parametrik, namun jika kurva regresi memiliki bentuk yang tidak diketahui maka digunakan pendekatan regresi nonparametrik. Pendekatan regresi semiparametrik digunakan apabila bentuk kurva regresi sebagian diketahui dan sebagian tidak diketahui. 2.3.1 Regresi Parametrik Regresi parametrik didefinisikan sebagai regresi yang bentuk kurva regresi diasumsikan diketahui, misalnya linier, kuadratik, kubik, polinomial derajat p, eksponen, dan lain-lain (Budiantara, 2009). Misalkan terdapat data berpasangan xi , yi , hubungan antara variabel tersebut diasumsikan mengikuti model regresi sehingga bentuk umum modelnya adalah sebagai berikut. (2.4) yi f xi i , i 1, 2,3,..., n dimana, yi : variabel respon ke-i f : kurva regresi
i : error random ke-i yang diasumsikan memenuhi IIDN (0,σ2) Apabila terdapat p variabel prediktor dan satu variabel respon, maka regresi yang digunakan adalah parametrik linear berganda yang dapat ditulis sebagai berikut. yi 0 1 x1i 2 x2i ... p x pi i , i 1, 2,3,..., n (2.5) Persamaan (2.5) dapat ditulis dalam bentuk matriks sebagai berikut. (2.6) Y = Xβ + ε Variabel respon Y berupa vektor kolom berukuran n x 1, X merupakan matriks berukuran n x (p + 1) yang berisikan satu
10 kolom angka 1 dan p kolom data prediktor, serta error random ε yang merupakan vektor kolom berukuran n x 1. Untuk memperoleh estimasi parameter pada regresi parametrik, dapat dilakukan dengan metode Ordinary Least Square (OLS) dengan meminimumkan jumlah kuadrat error.
Estimasi parameter βˆ diberikan oleh: -1 βˆ X'X X'Y
(2.7)
dengan X diberikan oleh:
1 x11 1 x 12 X 1 x1n
x p1 x p 2 x pn
(2.8)
2.3.2 Regresi Nonparametrik Metode pendekatan regresi yang sesuai untuk pola data yang bentuk kurva regresinya tidak diketahui atau tidak terdapat informasi masa lalu yang lengkap tentang bentuk pola data adalah regresi nonparametrik. Menurut Eubank (1988), model regresi nonparametrik secara umum adalah sebagai berikut. (2.9) yi f xi i , i 1, 2,3,..., n dengan: yi : variabel respon ke-i f : kurva regresi yang diasumsikan bentuknya tidak diketahui
i : error random ke-i yang diasumsikan memenuhi IIDN (0,σ2) 2.3.3 Regresi Semiparametrik Regresi semiparametrik adalah gabungan dari regresi parametrik dan regresi nonparametrik (Ruppert et al., 2003). Selain itu, Budiantara (2009) mendefinisikan regresi semiparametrik sebagai model regresi yang bentuk kurva regresi sebagian diasumsikan diketahui polanya dan sebagian tidak diketahui
11 bentuk polanya. Misal terdapat data berpasangan
yi , xi , ti
dan
hubungan antara yi , xi , ti mengikuti regresi semiparametrik. Berikut merupakan model regresi semiparametrik (Budiantara, 2005). (2.10) yi x'iβ f ti i , i 1, 2,3,..., n Variabel respon yi berhubungan parametrik dengan variabel prediktor xi dan berhubungan nonparametrik dengan variabel prediktor ti . Persamaan (2.10) dapat ditulis dalam bentuk matriks sebagai berikut.
Y = Xβ + ε dengan X diberikan oleh: 1 1 X 1
x1
t1
x2
t2
t1 1 t2 1
xn
tn
tn 1
(2.11)
t1 K t2 K
(2.12)
tn K
2.4 Pemodelan Regresi Spline Truncated Spline merupakan salah satu model regresi nonparametrik dan semiparametrik. Model spline memiliki fleksibilitas yang tinggi dan memiliki kemampuan sangat baik untuk menangani data yang perilakunya berubah pada sub interval tertentu. Salah satu kelebihan pendekatan spline adalah model ini cenderung mencari sendiri estimasi data ke mana pola data tersebut bergerak. Kelebihan ini terjadi karena dalam spline terdapat titik-titik knot. Budiantara (2001) memberikan suatu basis dari ruang spline berorde p sebagai berikut.
1, x, x ,..., x , x 2
p
p 1
, x 2 ,..., x K p
p
(2.13)
12 Fungsi potongan (truncated) diberikan oleh:
x k
p
x k p , x k 0, x k
(2.14)
Pada regresi spline, fungsi regresi f dapat ditulis sebagai berikut. p
K
j 0
k 1
f x j x j p k x k
p
(2.15)
Sehingga model regresi spline dapat dituliskan sebagai berikut. (2.16) y f x atau dapat juga ditulis sebagai berikut. p
K
j 0
k 1
y j x j p k x k p
(2.17)
dimana: j : parameter model polinomial
x
: variabel prediktor
p k : parameter pada komponen truncated, k 1, 2,, K K
: banyaknya knot k : titik-titik knot Apabila terdapat lebih dari satu variabel prediktor dan satu variabel respon, maka regresi spline yang digunakan adalah regresi spline multivariabel. Misalkan yi merupakan variabel respon,
i 1, 2,..., n , sedangkan x1i , x2i ,..., x pi adalah variabel prediktor yang pola hubungannya dengan variabel respon tidak diketahui polanya. Hubungan x1i , x2i ,..., x pi , yi dapat dimodelkan dengan regresi nonparametrik multivariabel sebagai berikut.
yi f1 x1i f 2 x2i ... f p x pi i , i 1, 2,..., n (2.18) p
fl xli i l 1
(2.19)
13 Pada regresi spline, masing-masing fungsi f dapat ditulis sebagai berikut. p
K
j 0
k 1
f ( xli ) lj xlij l ( p k ) xli lk
p
(2.20)
Fungsi potongan (truncated) sebagai berikut.
xli lk
p
xli lk p , xli lk 0, xli lk
(2.21)
Error i diasumsikan identik, independen, dan berdistribusi normal. 2.5 Pemilihan Titik Knot Optimal Titik knot adalah titik perpaduan bersama dimana terdapat perubahan perilaku pola data pada interval yang berlainan (Budiantara, 2006). Model regresi spline terbaik tergantung pada titik knot yang optimal. Salah satu metode yang dapat digunakan untuk menentukan titik knot optimal adalah Generalized Cross Validation (GCV). Titik knot optimal diperoleh dari nilai GCV minimum. Secara umum, fungsi GCV adalah sebagai berikut. MSE 1 , 2 ,..., K (2.22) GCV 1 , 2 ,..., K 2 n 1trace I A 1 , 2 ,..., K
Dimana I adalah matriks indentitas, n adalah jumlah pengamatan, matriks X(XTX)-1XT, serta A 1 , 2 ,..., K merupakan
MSE 1 , 2 ,..., K diberikan oleh: n
MSE 1 , 2 ,..., K n 1 yi yˆi
2
(2.23)
i 1
2.6 Pengujian Parameter Model Pengujian parameter model digunakan untuk mengetahui apakah variabel prediktor berpengaruh atau tidak terhadap variabel respon. Tahap dalam pengujian parameter adalah pengujian secara
14 serentak dan pengujian secara parsial atau individu. Diberikan model regresi semiparametrik sebagai berikut. q
p
r
j 0
l 1
k 1
yi j x ji l til p k ti k i p
(2.24)
2.6.1 Pengujian Secara Serentak Pengujian secara serentak dilakukan untuk mengetahui apakah variabel prediktor secara bersama-sama berpengaruh terhadap variabel respon. Hipotesis yang digunakan dalam pengujian serentak adalah sebagai berikut. H0 : 1 2 ... q 1 2 ... p r 0 H1 : minimal terdapat satu j 0 atau l 0 , j 1, 2,..., q ,
l 1, 2,..., p r Dimana q adalah jumlah parameter komponen parametrik dan p+r merupakan jumlah parameter komponen nonparametrik. Statistik uji yang digunakan adalah sebagai berikut.
F
MSregresi
(2.25)
MSerror
Perhitungan nilai statistik uji F didapatkan dari Analisis Ragam (ANOVA) seperti yang ditunjukkan Tabel 2.1. Tabel 2. 1 Analisis Ragam (ANOVA)
Sumber Variasi
Derajat Bebas (df)
Regresi
q pr
Jumlah Kuadrat (SS)
Rataan Kuadrat (MS) n
n
yˆi y
2
yˆ y i 1
n q p r 1
n
n
y yˆ i 1
i
n
Total
n 1
y y i 1
i
2
i
2
y yˆ i 1
i
2
i
n q p r 1
-
MSregresi MSerror
q pr
i 1
Error
2
i
F
-
15 Daerah penolakan H0 adalah F F , q p r , n q p r 1 atau
P-value < α. Jika H0 ditolak maka dapat disimpulkan bahwa minimal terdapat satu parameter pada model regresi yang signifikan terhadap model. 2.6.2 Pengujian Secara Parsial Pengujian parsial atau individu dilakukan apabila pada pengujian parameter secara serentak didapatkan kesimpulan bahwa minimal terdapat satu parameter pada model regresi yang signifikan. Tujuan dari pengujian parsial adalah untuk mengetahui parameter mana yang berpengaruh dan tidak berpengaruh signifikan terhadap model regresi. Hipotesis yang digunakan pada pengujian parsial untuk komponen parametrik adalah sebagai berikut. H0 : j 0 H1 : j 0, j 1, 2,..., q Statistik uji yang digunakan untuk komponen parametrik adalah sebagai berikut. ˆ j (2.26) t , j 1, 2,..., q SE ˆ
j
merupakan akar dari elemen diagonal ke-j dari
Dengan SE ˆ j
2 matriks X'X ˆ , dimana ˆ 2 adalah Mean Square Error -1
(MSE) dan X diberikan oleh: x11 x 12 X x1n
x21 x22 x2 n
xq1 xq 2 xqn
(2.27)
Daerah penolakan untuk komponen parametrik adalah t t /2,n ( q p r ) 1 atau P-value < α.
16 Selanjutnya, hipotesis yang digunakan dalam pengujian parsial untuk komponen nonparametrik adalah sebagai berikut. H0 : l 0 H1 : l 0, l 1, 2,..., p r Statistik uji yang digunakan untuk komponen nonparametrik adalah sebagai berikut. ˆl (2.28) t , l 1, 2,..., p r SE ˆl Dengan SE ˆl merupakan akar dari elemen diagonal ke-l dari 2 matriks X'X ˆ , dimana ˆ 2 adalah Mean Square Error -1
(MSE) dan X diberikan oleh: t11 t X 12 t 1n
t11 k1 t12 k1
t11 kr t12 kr
t1n k1
t1n kr
t p2
t t
p 2 k1
t pn
t
pn
t p1
p1
k1
k1
t t
kr p 2 kr t k pn r p1
(2.29)
Daerah penolakan untuk komponen nonparametrik adalah t t /2,n ( q p r ) 1 atau P-value < α. 2.7 Kriteria Kebaikan Model Regresi Salah satu kriteria yang dapat digunakan untuk menentukan kebaikan model regresi adalah menggunakan koefisien determinasi (R2). Koefisien determinasi merupakan nilai proporsi keragaman total di sekitar nilai y yang dapat dijelaskan oleh model regresi (Drapper & Smith, 1992). Semakin tinggi nilai R2 yang dihasilkan suatu model, maka semakin baik pula variabel-variabel prediktor dalam model dalam menjelaskan variabilitas variabel respon. Berikut merupakan rumus dari koefisien determinasi (R2). n
R2
yˆ y
2
y y
2
i 1 n
i 1
i
i
(2.30)
17 2.8 Pengujian Asumsi Residual Setelah mendapatkan model terbaik dari regresi spline, perlu dilakukan pengujian asumsi residual untuk mengetahui apakah residual yang dihasilkan telah memenuhi asumsi identik, independen, dan berdistribusi normal (IIDN). 2.8.1 Uji Asumsi Identik Pengujian asumsi identik bertujuan untuk mengetahui homogenitas variansi residual. Pengujian asumsi identik dilakukan dengan uji Glejser. Uji Glejser dilakukan dengan cara meregresikan absolut residual dengan variabel prediktornya (Gujarati & Porter, 2009). Hipotesis yang digunakan dalam pengujian ini adalah sebagai berikut. H0 : 12 22 ... n2 2 H1 : minimal terdapat satu i2 2 , i 1, 2,..., n Statistik uji Glejser dirumuskan sebagai berikut.
eˆ n
F
i
i 1
e n
i 1
i
eˆi
ei 2
2
q p r
(2.31)
n q p r 1
Daerah penolakan adalah F F ,( q p r ) 1, n ( q p r )1 . Jika tolak H0 maka dapat disimpulkan bahwa terdapat indikasi adanya kasus heteroskedastisitas. 2.8.2 Asumsi Independen Asumsi independen dilakukan untuk mengetahui apakah terdapat autokorelasi atau tidak pada residual. Untuk mengetahui apakah terjadi autokorelasi secara visual dilakukan dengan melihat plot dari Autocorrelation Function (ACF) dari residual. Untuk mendapatkan nilai ACF digunakan rumus sebagai berikut. n
e
ˆ k t k 1
t k
e et e (2.32)
n
e e t 1
t
2
18 Asumsi independen dapat dideteksi dengan menggunakan interval konfidensi dengan rumus sebagai berikut. (2.33) t ,n 1 SE ˆ k k t ,n 1 SE ˆ k 2
2
Apabila terdapat ACF yang keluar dari interval konfidensi, maka diindikasikan adanya autokorelasi antar residual. Sebaliknya, jika semua nilai ACF berada di dalam batas signifikansi maka tidak terdapat kasus autokorelasi atau asumsi independen terpenuhi. 2.8.3 Uji Asumsi Berdistribusi Normal Untuk mengetahui apakah residual berdistribusi normal, secara visual dapat dilakukan dengan menggunakan normality plot. Apabila plot berada di sekitar garis regresi, maka dapat dikatakan bahwa residual berdistribusi normal. Sebaliknya, jika plot tidak berada di sekitar garis regresi, maka dapat dikatakan residual tidak berdistribusi normal. Selain dengan cara visual, untuk mengetahui apakah residual berdistribusi normal dapat dilakukan dengan pengujian Kolmogorov-Smirnov (Daniel, 1989). Hipotesis yang digunakan adalah sebagai berikut. H0 : F ( x) F0 ( x) atau residual berdistribusi normal H1 : F ( x) F0 ( x) atau residual tidak berdistribusi normal Nilai statistik uji Kolmogorov-Smirnov adalah sebagai berikut. D sup S ( x) F0 ( x) (2.34) x
F0(x) merupakan fungsi distribusi yang dihipotesiskan, sedangkan S(x) adalah fungsi peluang kumulatif yang dihitung dari data sampel. Daerah penolakan adalah D q1 pada tabel Kolmogorov-Smirnov atau P-value < α. Jika tolak H0, maka dapat disimpulkan bahwa residual tidak berdistribusi normal. 2.9 Indeks Pembangunan Kesehatan Masyarakat IPKM (Indeks Pembangunan Kesehatan Masyarakat) adalah indikator komposit yang menggambarkan kemajuan pembangunan kesehatan yang dirumuskan dari data kesehatan, yaitu Riset
19 Kesehatan Dasar (Riskesdas), Survei Sosial Ekonomi Nasional (Susenas), dan Survei Potensi Desa (Podes). IPKM pertama dirumuskan dari 24 indikator kesehatan yang kemudian dilakukan pengembangan model IPKM 2013 dengan menggunakan 30 indikator yang dapat digunakan sebagai rumus standar untuk IPKM yang akan datang. IPKM dapat dimanfaatkan untuk menentukan peringkat provinsi dan kabupaten/kota dalam keberhasilan pembangunan kesehatan masyarakat, sebagai bahan advokasi ke pemerintah sehingga sumber daya dan program kesehatan diprioritaskan, serta sebagai salah satu kriteria penentuan alokasi dana bantuan kesehatan pusat ke daerah (Kementerian Kesehatan RI, 2010). Indikator dalam IPKM 2013 ditentukan berdasarkan konsep determinan sosial kesehatan yang meliputi kesehatan perorangan, keluarga, masyarakat, dan sistem pelayanan kesehatan. Dimana aspek yang menjadi pertimbangan adalah prioritas program kesehatan nasional yang tertuang dalam rencana pembangunan jangka menengah dan panjang, komitmen untuk pembangunan kesehatan secara global, besaran masalah kesehatan yang menjadi masalah kesehatan utama secara nasional, pertimbangan secara referensi dan rekomendasi pelaksana program kesehatan, serta pertimbangan secara statistik mencakup aspek variasi data dan jumlah sampel untuk keterwakilan kabupaten/kota. Penentuan indikator dilakukan dengan melakukan pertemuan, konsultasi, dan diskusi dengan para pakar baik secara nasional maupun internasional serta para pengambil keputusan pada program kesehatan (Tim Penyusun IPKM, 2014). 2.10 Kerangka Konsep Penelitian Kerangka konsep merupakan kerangka hubungan antara konsep-konsep yang ingin diamati atau diukur melalui penelitian yang akan dilakukan (Notoatmodjo, 2002). Menurut Hendrik L. Blum, faktor yang memengaruhi kesehatan masyarakat atau determinan kesehatan adalah keturunan, lingkungan, fasilitas kesehatan, dan perilaku. Berdasarkan faktor determinan tersebut,
20 diperoleh aspek-aspek yang diduga memengaruhi indeks kesehatan, yaitu sebagai berikut. 1. Ekonomi Berdasarkan teori konsumsi dan ekonomi kesejahteraan menurut Pindyck & Rubinfeld (1998), untuk mencapai kesejahteraan tertentu individu akan mengonsumsi sejumlah barang dan jasa, dalam hal ini ditekankan dalam bentuk pelayanan kesehatan. Jika jumlah pendapatan meningkat, maka tingkat konsumsi akan meningkat pula. Sebaliknya jika pendapatan rendah maka individu cenderung kesulitan untuk memenuhi kebutuhannya termasuk kebutuhan dalam kesehatan sehingga akan berpengaruh pada ketercapaian indeks kesehatan. 2. Lingkungan dan Perilaku Menurut Blum (1974), faktor lingkungan dan perilaku berhubungan dengan tinggi rendahnya derajat kesehatan suatu masyarakat. Apabila lingkungan tempat tinggal suatu masyarakat sehat dan juga ditunjang dengan perilaku hidup bersih dan sehat, maka akan berpengaruh terhadap tingginya indeks kesehatan masyarakat. 3. Kependudukan Pertumbuhan penduduk yang tinggi berkaitan dengan sulitnya pelayanan kesehatan yang merata pada masyarakat karena terdapat banyak penduduk dalam suatu wilayah tertentu. Apabila pelayanan kesehatan yang didapatkan masyarakat tidak merata, maka dapat menyebabkan tidak tertolongnya masyarakat yang sedang mengalami sakit dan dapat mengakibatkan terjadinya kematian. Angka kematian merupakan indikator yang sensitif untuk menilai keberhasilan pembangunan pada sektor kesehatan. Menurut Afifah dkk (2008), Angka Kematian Bayi sering kali dianggap sebagai indikator umum dari pembangunan sosial dan secara khusus sebagai indikator kesehatan.
BAB III METODOLOGI PENELITIAN 3.1 Sumber Data Data yang digunakan pada penelitian ini adalah data sekunder yang diambil dari Badan Pusat Statistika Jawa Timur dan publikasi Profil Kesehatan Jawa Timur yang diterbitkan oleh Dinas Kesehatan Provinsi Jawa Timur. Data yang digunakan merupakan data pada tahun 2013. Banyak observasi adalah 38 kabupaten/kota yang ada di Provinsi Jawa Timur. 3.2 Variabel Penelitian Variabel respon yang digunakan pada penelitian ini adalah Indeks Pembangunan Kesehatan Masyarakat (IPKM) di 38 kabupaten/kota di Provinsi Jawa Timur. Variabel yang digunakan pada penelitian didapatkan dari variabel yang diduga berpengaruh dari penelitian sebelumnya. Berikut merupakan variabel penelitian yang digunakan pada penelitian ini. Tabel 3.1 Variabel Penelitian
Variabel y x1 x2 t1 t2 t3 t4
Keterangan IPKM Persentase Penduduk Miskin Angka Kematian Bayi Kepadatan Penduduk Angka Kematian Ibu Persentase Rumah Tangga Berperilaku Hidup Bersih dan Sehat Persentase Rumah Sehat
Berikut ini merupakan keterangan mengenai variabel penelitian yang digunakan. a. Variabel y merupakan variabel respon yang menyatakan Indeks Pembangunan Kesehatan Masyarakat. Indikator utama pembangunan kesehatan yang digunakan mencakup kesehatan balita, kesehatan reproduksi, pelayanan kesehatan, perilaku, penyakit tidak menular, penyakit menular, dan 21
22
b.
c.
kesehatan lingkungan. Data yang digunakan untuk menyusun IPKM 2013 adalah Riset Dasar Kesehatan (Riskesdas) 2013 dan Potensi Desa (Podes) 2011. Variabel x1 merupakan variabel yang menyatakan Persentase Penduduk Miskin. Persentase penduduk miskin atau yang disebut juga dengan Headcount Index merupakan persentase penduduk miskin yang berada di bawah Garis Kemiskinan (GK). Headcount Index secara sederhana mengukur proporsi yang dikategorikan miskin. Variabel x2 merupakan variabel yang menyatakan Angka Kematian Bayi. Angka Kematian Bayi (AKB) merupakan angka yang menunjukkan banyaknya kematian bayi usia 0 tahun dari setiap 1000 kelahiran hidup pada tahun tertentu atau dapat juga dikatakan sebagai probabilitas bayi meninggal sebelum mencapai usia satu tahun (dinyatakan dengan per seribu kelahiran hidup). AKB diperoleh dari rumus sebagai berikut. AKB
d.
Kematian bayi usia di bawah 1 tahun x1000 Kelahiran hidup
Variabel t1 merupakan variabel yang menyatakan Kepadatan Penduduk. Kepadatan Penduduk (KP) menunjukkan banyaknya jumlah penduduk untuk setiap kilometer persegi luas wilayah. Kepadatan penduduk dapat diperoleh dengan rumus sebagai berikut. KP
e.
Jumlah penduduk Luas wilayah (km 2 )
Variabel t2 merupakan variabel yang menyatakan Angka Kematian Ibu. Angka Kematian Ibu (AKI) merupakan banyaknya kematian perempuan pada saat hamil atau selama 42 hari sejak terminasi kehamilan tanpa memandang lama dan tempat persalinan, yang disebabkan karena kehamilannya atau pengelolaannya, dan bukan karena sebab-sebab lain, per 100.000 kelahiran hidup. AKI diperoleh dengan rumus sebagai berikut.
23 AKI
f.
g.
Kematian ibu dalam tahap kehamilan atau kelahiran x1000 Kelahiran hidup
Variabel t3 merupakan variabel yang menyatakan Persentase Rumah Tangga Berperilaku Hidup Bersih dan Sehat. Menurut Panduan Pembinaan dan Penilaian PHBS di Rumah Tangga yang diterbitkan oleh Pusat Promosi Kesehatan Kementerian Kesehatan Republik Indonesia, Persentase Rumah Tangga Berperilaku Hidup Bersih dan Sehat (Rumah Tangga BerPBHS) merupakan persentase rumah tangga yang memenuhi 10 indikator PHBS, yaitu persalinan ditolong oleh tenaga kesehatan, memberi bayi ASI ekslusif, menimbang balita setiap bulan, menggunakan air bersih, mencuci tangan dengan air bersih yang mengalir dan menggunakan sabun, menggunakan jamban sehat, memberantas jentik di rumah sekali seminggu, makan sayur dan buah setiap hari, melakukan aktivitas fisik setiap hari, dan tidak merokok di dalam rumah. Variabel t4 merupakan variabel yang menyatakan Persentase Rumah Sehat. Pembagian bobot penilaian rumah sehat meliputi bobot komponen rumah, bobot sarana sanitasi, dan bobot pada perilaku penghuni. Menurut Pedoman Teknis Penilaian Rumah Sehat Departemen Kesehatan Republik Indonesia tahun 2007, rumah dikatakan sehat apabila memenuhi kriteria, yaitu memenuhi kebutuhan psikologis antara lain privasi yang cukup, komunikasi yang sehat antar anggota keluarga dan penghuni rumah, adanya ruangan khusus untuk istirahat (ruang tidur) bagi masing-masing penghuni; memenuhi persyaratan pencegahan penularan penyakit antar penghuni rumah dengan penyediaan air bersih, pengelolaan tinja dan limbah rumah tangga, bebas vektor penyakit dan tikus, kepadatan hunian yang tidak berlebihan, cukup sinar matahari pagi, terlindungnya makanan dan minuman dari pencemaran dan penghawaan yang cukup; serta memenuhi persyaratan pencegahan terjadinya kecelakaan baik yang timbul karena pengaruh luar dan dalam rumah,
24 antara lain persyaratan garis sempadan jalan, konstruksi bangunan rumah, bahaya kebakaran dan kecelakaan di dalam rumah. 3.3 Langkah Penelitian Langkah analisis dalam penelitian ini adalah sebagai berikut. 1. Mengumpulkan data IPKM di Provinsi Jawa Timur beserta variabel-variabel yang berpengaruh. 2. Mendeskripsikan karakteristik data dari IPKM di Provinsi Jawa Timur beserta variabel-variabel yang berpengaruh. 3. Membuat scatterplot dan melakukan uji RESET untuk mengetahui pola data antara variabel IPKM dengan variabelvariabel prediktornya. 4. Mencari variabel prediktor komponen parametrik. 5. Mencari variabel prediktor komponen nonparametrik. 6. Memodelkan data dengan pendekatan spline dengan satu knot, dua knot, tiga knot, dan kombinasi knot. 7. Memilih titik knot optimal berdasarkan nilai Generalized Cross Validation (GCV) yang minimum. 8. Mendapatkan model regresi spline dengan titik knot optimal. 9. Menguji signifikansi parameter regresi spline secara serentak dan parsial. 10. Menghitung nilai koefisien determinasi (R2). 11. Menguji asumsi residual IIDN. Apabila residual model spline tidak memenuhi asumsi, maka dilakukan transformasi data. Selanjutnya memulai kembali dari langkah analisis (6). 12. Menginterpretasikan model dan menarik kesimpulan. 3.4 Diagram Alir Dari langkah analisis disajikan secara ringkas dalam diagram alir yang dijelaskan pada Gambar 3.1.
25
Mengumpulkan data IPKM di Jawa Timur beserta faktor-faktor yang diduga berpengaruh
Melakukan analisis karakteristik data dengan menggunakan statistika deskriptif
Mengidentifikasi pola data yang terbentuk dengan menggunakan scatterplot dan uji RESET antara IPKM dengan variabel prediktor
Menentukan komponen parametrik dan nonparametrik
Memilih titik knot optimal dengan melihat nilai GCV paling minimum
Melakukan pengujian signifikansi parameter
Transformasi data
Melakukan pemodelan IPKM di Provinsi Jawa Timur dengan metode regresi semiparametrik spline
Menghitung R2
Uji asumsi residual IIDN
Ya Interpretasi model dan membuat kesimpulan Gambar 3.1 Diagram Alir
Tidak
26
(Halaman ini sengaja dikosongkan)
BAB IV ANALISIS DAN PEMBAHASAN Pada bab ini akan dibahas mengenai karakteristik dari data Indeks Pembangunan Kesehatan Masyarakat (IPKM) di Provinsi Jawa Timur beserta faktor-faktor yang diduga memengaruhinya. Karakteristik dari data tersebut meliputi statistika deskriptif dan penentuan variabel yang termasuk komponen parametrik dan komponen nonparametrik. Selanjutnya adalah melakukan pemodelan menggunakan pendekatan regresi semiparametrik spline dengan menggunakan satu knot, dua knot, tiga knot, dan kombinasi knot. 4.1 Karakteristik Indeks Pembangunan Kesehatan Masyarakat Provinsi Jawa Timur IPKM merupakan indikator komposit yang menggambarkan kemajuan pembangunan kesehatan yang dirumuskan dari 30 indikator kesehatan dan disusun oleh Badan Penelitian dan Pengembangan Kesehatan (Balitbangkes) Kementerian Kesehatan RI. Gambar 4.1 menyajikan diagram batang IPKM di setiap kabupaten/kota yang berada di Provinsi Jawa Timur sesuai dengan data di Lampiran 1. Berdasarkan diagram batang yang telah diurutkan dari IPKM terendah hingga tertinggi pada Gambar 4.1, diketahui bahwa IPKM terendah di Provinsi Jawa Timur terdapat di Kabupaten Pamekasan dengan IPKM sebesar 0,5874. Indeks kesehatan di Kabupaten Pamekasan rendah disebabkan oleh rendahnya pengetahuan masyarakat tentang gizi bayi dan balita sehingga kasus gizi buruk dan gizi kurang masih banyak terjadi. Selain itu, masyarakat di Kabupaten Pamekasan masih kurang sadar akan pentingnya perilaku hidup bersih dan sehat. Sedangkan, Kota Madiun meraih IPKM tertinggi di Jawa Timur dengan nilai sebesar 0,79. Hal ini disebabkan oleh fasilitas kesehatan yang baik dari segi kualitas maupun kuantitas masyarakat di Kota Madiun telah memadai 27
28
Kabupaten/Kota
bahkan menjadi rujukan bagi rumah sakit di kabupaten sekitarnya. Selain itu angka harapan hidup Kota Madiun yang terus meningkat sejak tahun 2009 hingga tahun 2013 yang menandakan kesehatan penduduk yang semakin baik. Kota Madiun Kota Kediri Kota Blitar Kota Malang Kota Batu Kab. Nganjuk Kota Mojokerto Kab. Ponorogo Kota Surabaya Kab. Sidoarjo Kota Pasuruan Kab. Tulungagung Kab. Magetan Kab. Gresik Kab. Jombang Kab. Madiun Kota Probolinggo Kab. Mojokerto Kab. Kediri Kab. Ngawi Kab. Trenggalek Kab. Lamongan Kab. Pacitan Kab. Blitar Kab. Malang Kab. Banyuwangi Kab. Tuban Kab. Pasuruan Kab. Bojonegoro Kab. Sampang Kab. Lumajang Kab. Situbondo Kab. Probolinggo Kab. Jember Kab. Bangkalan Kab. Bondowoso Kab. Sumenep Kab. Pamekasan
0.79 0.783 0.7718 0.7588 0.7584 0.7556 0.749 0.744 0.7406 0.7395 0.7388 0.7343 0.7339 0.7298 0.727 0.7269 0.725 0.7246 0.7162 0.7132 0.7127 0.7116 0.7006 0.6948 0.6897 0.6878 0.6849 0.6778 0.6772 0.6643 0.6581 0.6517 0.6405 0.6391 0.6381 0.6083 0.6002 0.5874 0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
IPKM 2013
Gambar 4. 1 IPKM di Kabupaten/Kota Provinsi Jawa Timur
0.9
29 Terdapat 13 kabupaten/kota yang masih berada di bawah IPKM Indonesia (0,6879), yaitu Kabupaten Banyuwangi, Kabupaten Tuban, Kabupaten Pasuruan, Kabupaten Bojonegoro, Kabupaten Sampang, Kabupaten Lumajang, Kabupaten Situbondo, Kabupaten Probolinggo, Kabupaten Jember, Kabupaten Bangkalan, Kabupaten Bondowoso, Kabupaten Sumenep, dan Kabupaten Pamekasan. Walaupun sebagian besar kabupaten/kota Provinsi Jawa Timur sudah memiliki nilai yang lebih besar dari IPKM Indonesia, namun IPKM di Jawa Timur masih perlu ditingkatkan dengan memerhatikan faktor-faktor yang memengaruhi IPKM. Berikut ini adalah karakteristik dari keenam variabel yang diduga memengaruhi IPKM di Provinsi Jawa Timur. Tabel 4. 1 Statistika Deskriptif Variabel Penelitian
Variabel
Rata-rata
Varians
Minimum
Maksimum
y
0,70487
0,00254
0,5874
0,79
x1
12,54
27,136
4,77
27,08
x2
32,79
155,25
18,71
62,45
t1
1802
4664624
382
8551
t2
105,23
1856,33
30,82
212,71
t3
45,34
210,75
17,14
67,32
t4
38,85
548,91
1,02
81,03
Berdasarkan Tabel 4.1, variabel respon (y) merupakan IPKM di Provinsi Jawa Timur dengan rata-rata IPKM pada tahun 2013 sebesar 0,70487 dan varians sebesar 0,00254. IPKM tertinggi adalah sebesar 0,79 yang dicapai oleh Kota Madiun, sedangkan IPKM terendah adalah 0,5874 yang dicapai oleh Kabupaten Pamekasan. Variabel x1 merupakan variabel persentase penduduk miskin yang diduga memengaruhi IPKM. Pada Tabel 4.2 diketahui bahwa variabel persentase penduduk miskin memiliki rata-rata sebesar 12,54% dengan varians sebesar 27,136. Persentase tertinggi
30 terdapat di Kabupaten Sampang, yaitu sebesar 27,08%. Hal tersebut disebabkan oleh sumber daya alam dan infrastruktur yang kurang memadai. Selain itu, banyak penduduk di Kabupaten Sampang yang putus sekolah. Persentase penduduk miskin terendah, yaitu di Kota Batu dengan nilai sebesar 4,77%. Pemerintah Kota Batu telah menerapkan program pengentasan kemiskinan untuk mengurangi jumlah penduduk miskin sehingga persentase penduduk miskin semakin turun dari tahun ke tahun. Variabel x2 merupakan variabel Angka Kematian Bayi yang diduga memengaruhi IPKM. Berdasarkan Tabel 4.2, diketahui bahwa variabel Angka Kematian Bayi (AKB) memiliki nilai ratarata sebesar 32,79 dengan varians sebesar 155,25. AKB tertinggi terdapat pada Kabupaten Probolinggo dengan nilai sebesar 62,45. Hal tersebut disebabkan oleh masih banyak masyarakat yang bergantung pada tenaga non medis untuk membantu persalinan, seperti dukun. AKB terendah, yaitu 18,71 terdapat di Kota Blitar. Berbagai upaya telah dilakuan oleh pemerintah untuk menurunkan AKB di Kota Blitar dengan cara melakukan pelatihan tenaga kesehatan dan peningkatan mutu layanan kesehatan. Variabel t1 merupakan variabel kepadatan penduduk yang diduga memengaruhi IPKM. Dalam Tabel 4.2 diketahui bahwa variabel kepadatan penduduk memiliki rata-rata sebesar 1802 per km2 dengan varians sebesar 4664624. Nilai tertinggi terdapat di Kota Surabaya, yaitu sebesar 8551 per km2. Kepadatan penduduk di Surabaya disebabkan oleh tingginya tingkat urbanisasi untuk mencari pekerjaan atau menimba ilmu sehingga penduduk terus meningkat setiap tahun. Sedangkan kepadatan penduduk terendah terdapat di Kabupaten Pacitan dengan nilai sebesar 382 per km2. Kabupaten Pacitan memiliki wilayah yang sebagian besar merupakan hutan dan tanah ladang. Hanya 0,2% wilayahnya digunakan sebagai pemukiman penduduk. Variabel t2 merupakan variabel Angka Kematian Ibu (AKI) yang diduga memengaruhi IPKM. Berdasarkan Tabel 4.2, diketahui bahwa variabel AKI memiliki nilai rata-rata sebesar
31 105,23 dengan varians sebesar 1856,33. AKI tertinggi bernilai 212,71 yang terdapat di Kota Probolinggo. Masih banyak ibu di Kota Probolinggo yang menderita penyakit saat kehamilan sehingga menyebabkan kematian. Selain itu, usia ibu hamil yang masih muda dan melakukan persalinan dengan jasa dukun. Sedangkan nilai terendah yaitu 30,82 terdapat di Kota Batu. 95% persalinan dilakukan dengan pertolongan oleh tenaga kesehatan di Kota Batu. Variabel t3 merupakan variabel persentase rumah tangga berperilaku hidup bersih dan sehat yang diduga memengaruhi IPKM. Pada Tabel 4.2 variabel persentase rumah tangga ber-PBHS memiliki rata-rata sebesar 45,34% dengan varians sebesar 210,75. Persentase tertinggi terdapat di Kota Surabaya, yaitu sebesar 67,32%. Pemerintah Surabaya sejak tahun 2010 rutin melakukan kajian mengenai PHBS di rumah tangga untuk upaya percepatan capaian rumah tangga sehat. Sedangkan persentase terendah terdapat di Kabupaten Situbondo dengan nilai sebesar 17,14%. Hal tersebut disebabkan karena tingginya masyarakat yang merokok di dalam rumah, kurangnya kesadaran masyarakat akan jamban sehat, dan kurangnya pemberian ASI ekslusif. Variabel t4 merupakan variabel persentase rumah sehat yang diduga memengaruhi IPKM. Dalam Tabel 4.2 diketahui bahwa variabel persentase rumah sehat memiliki rata-rata sebesar 38,85% dengan varians sebesar 548,91. Persentase tertinggi terdapat di Kota Surabaya, yaitu sebesar 81,03%. Pemerintah Surabaya terus menggalakkan rumah sehat dengan melakukan sosialisasi program rumah sehat kepada masyarakat. Sedangkan persentase rumah sehat terendah terdapat di Kabupaten Blitar dengan nilai sebesar 1,02%. Masih terdapat banyak rumah yang tidak dibina menjadi rumah sehat di Kabupaten Blitar.
32 4.2 Pemodelan Indeks Pembangunan Kesehatan Masyarakat Provinsi Jawa Timur Menggunakan Pendekatan Regresi Semiparametrik Spline Pemodelan IPKM di Provinsi Jawa Timur menggunakan IPKM sebagai variabel respon dengan faktor-faktor yang diduga memengaruhinya menggunakan pendekatan regresi semiparametrik spline. Sebelumnya dilakukan analisis deksriptif dari variabel-variabel yang diduga berpengaruh terhadap IPKM di Provinsi Jawa Timur untuk menentukan komponen parametrik dan komponen nonparametrik. 4.2.1 Pola Hubungan IPKM Provinsi Jawa Timur dengan Faktor-Faktor yang Diduga Memengaruhi Salah satu analisis deskriptif dari variabel-variabel yang diduga berpengaruh terhadap IPKM di Provinsi Jawa Timur adalah membuat scatterplot antara variabel respon dengan variabel prediktor untuk mengetahui pola hubungan IPKM dengan masingmasing variabel prediktor. 0.80
IPKM 2013
0.75
0.70
0.65
0.60 5
10
15 20 % Penduduk miskin
25
30
Gambar 4. 2 Pola Hubungan IPKM dengan Persentase Penduduk Miskin
Pada Gambar 4.2 diketahui variabel IPKM dengan persentase penduduk miskin membentuk suatu pola tertentu
33 sehingga dilakukan variabel ini termasuk dalam komponen parametrik. Semakin tinggi persentase penduduk miskin, maka semakin rendah IPKM yang dicapai karena keterbatasan dalam mengakses fasilitas maupun layanan kesehatan. Sebaliknya, jika persentase penduduk miskin rendah, maka IPKM akan semakin tinggi. 0.80
IPKM 2013
0.75
0.70
0.65
0.60 20
30
40 Angka kematian bayi
50
60
Gambar 4. 3 Pola Hubungan IPKM dengan Angka Kematian Bayi
Dapat dilihat bahwa Gambar 4.3 merupakan pola hubungan IPKM dengan Angka Kematian Bayi (AKB) membentuk suatu pola tertentu, sehingga termasuk dalam komponen parametrik. Apabila AKB semakin tinggi, maka IPKM yang dicapai akan semakin rendah. Sebaliknya, jika AKB rendah maka IPKM yang dicapai akan semakin tinggi. Gambar 4.4 menunjukkan pola hubungan antara IPKM dengan kepadatan penduduk tidak membentuk suatu pola tertentu. Selain itu, pola tersebut cenderung berubah pada sub-sub interval tertentu sehingga termasuk dalam komponen nonparametrik. Semakin padat jumlah penduduk dalam suatu wilayah, tentu berkaitan dengan sulitnya pelayanan kesehatan yang merata.
34
0.80
IPKM 2013
0.75
0.70
0.65
0.60 0
1000
2000
3000 4000 5000 6000 Kepadatan penduduk
7000
8000
9000
Gambar 4. 4 Pola Hubungan IPKM dengan Kepadatan Penduduk
Berdasarkan Gambar 4.5, diketahui bahwa hubungan antara IPKM dengan Angka Kematian Ibu (AKI) pola tertentu, sehingga termasuk dalam komponen nonparametrik. Kematian dapat disebabkan karena tidak tertolongnya masyarakat akibat pelayanan kesehatan yang didapatkan masyarakat tidak merata. 0.80
IPKM 2013
0.75
0.70
0.65
0.60 50
100 150 Angka kematian ibu
200
Gambar 4. 5 Pola Hubungan IPKM dengan Angka Kematian Ibu
35 Dapat dilihat bahwa Gambar 4.6 merupakan pola hubungan IPKM dengan Persentase Rumah Tangga Berperilaku Hidup Bersih dan Sehat yang tidak membentuk suatu pola tertentu, sehingga termasuk ke dalam komponen nonparametrik. Apabila sebuah rumah tangga menerapkan perilaku hidup bersih dan sehat, maka rumah tangga tersebut cenderung akan terhindar dari sakit sehingga berpengaruh terhadap ketercapaian indeks kesehatan. 0.80
IPKM 2013
0.75
0.70
0.65
0.60 10
20
30 40 50 % Rumah tangga berPHBS
60
70
Gambar 4. 6 Pola Hubungan IPKM dengan Persentase Rumah Tangga Berperilaku Hidup Bersih dan Sehat
Gambar 4.7 menunjukkan pola hubungan antara IPKM dengan persentase rumah sehat tidak membentuk suatu pola tertentu. Kondisi lingkungan rumah yang tidak sehat dapat menimbulkan terjadinya berbagai jenis gangguan kesehatan. Apabila rumah sudah dijaga kesehatannya, tentunya masyarakat yang tinggal di rumah tersebut akan terhindar dari penyebaran penyakit.
36
0.80
IPKM 2013
0.75
0.70
0.65
0.60 0
10
20
30
40 50 % Rumah sehat
60
70
80
90
Gambar 4. 7 Pola Hubungan IPKM dengan Persentase Rumah Sehat
Selain menggunakan scatterplot, untuk mengetahui hubungan antara IPKM dengan faktor-faktor yang memengaruhinya dapat menggunakan uji RESET. Uji RESET dapat membantu dalam pembuatan keputusan apakah kedua variabel memiliki hubungan linier atau nonlinier yang seringkali belum dapat diputuskan jika hanya berdasarkan scatterplot. Berikut merupakan hasil pengujian RESET antara IPKM dengan faktor yang diduga memengaruhi yang disajikan pada Tabel 4.2. Berdasarkan Tabel 4.2, didapatkan kesimpulan yang berbeda dengan scatterplot yang telah dilakukan sebelumnya. Dengan menggunakan scatterplot, variabel yang termasuk dalam komponen parametrik adalah persentase penduduk miskin dan Angka Kematian Bayi, sedangkan variabel yang termasuk dalam komponen nonparametrik adalah kepadatan penduduk, Angka Kematian Ibu, persentanse rumah tangga ber-PHBS, dan persentase rumah sehat. Namun dengan menggunakan uji RESET, variabel yang termasuk dalam komponen parametrik adalah persentase penduduk miskin, Angka Kematian Bayi, Angka Kematian Ibu, dan persentase rumah sehat, sedangkan variabel
37 yang termasuk dalam komponen nonparametrik adalah kepadatan penduduk dan persentase rumah tangga ber-PHBS. Tabel 4. 2 Uji RESET Antara IPKM dengan Faktor yang Memengaruhi
Variabel x1 x2 t1 y t2 t3 t4
F 1,8377 0,77701 4,7216 0,055805 3,004 0,38616
P-value 0,1839 0,3841 0,03664 0,8146 0,09186 0,5384
Kesimpulan Gagal tolak, linier Gagal tolak, linier Tolak, nonlinier Gagal tolak, linier Tolak, nonlinier Gagal tolak, linier
Pada penelitian ini metode yang digunakan menggunakan pendekatan regresi semiparametrik spline karena terdapat komponen parametrik dan komponen nonparametrik. Berdasarkan scatterplot, terdapat dua komponen parametrik dan empat komponen nonparametrik. Sedangkan dengan menggunakan uji RESET didapatkan hasil bahwa terdapat empat komponen parametrik dan dua komponen nonparametrik. 4.2.2 Pemilihan Titik Knot Optimum Pada penelitian ini, regesi semiparametrik spline digunakan untuk memodelkan IPKM di Provinsi Jawa Timur dengan faktorfaktor yang memengaruhi. Model regresi semiparametrik spline terbaik didapatkan dari titik knot yang optimum. Untuk menentukan titik knot yang optimum, digunakan metode GCV (Generalized Cross Validation). Titik knot yang digunakan pada penelitian ini adalah satu knot, dua knot, tiga knot, dan kombinasi knot. Dari keempat titik knot tersebut, untuk membentuk model digunakan titik knot optimum yang diperoleh dari nilai GCV yang minimum. 1.
Pola Berdasarkan Scatterplot Berikut akan ditampilkan nilai GCV dengan satu titik knot, dua titik knot, tiga titik knot, dan kombinasi titik knot untuk pola berdasarkan scatterplot.
38 a.
Pemilihan Titik Knot dengan Satu Titik Knot Estimasi model regresi semiparametrik spline dengan satu titik knot pada IPKM Provinsi Jawa Timur pada pola berdasarkan scatterplot adalah sebagai berikut. yi 0 1 x1i 2 x2i 1t1i 2 t1i 1 3t2i 4 t2i 1 5t3i 1
1
6 t3i 1 7t4i 8 t4i 1 1
1
Tabel 4.3 menunjukkan bahwa nilai GCV pada satu titik knot dengan empat komponen nonparametrik yang paling minimum adalah 0,0008940 dengan titik knot optimum yaitu variabel kepadatan penduduk (t1) pada titik 1549; variabel Angka Kematian Ibu (t2) pada titik 56,8; variabel persentase rumah tangga ber-PHBS (t3) pada titik 24,31; serta variabel persentase rumah sehat (t4) pada titik 12,45. Tabel 4. 3 Nilai GCV dengan Satu Titik Knot Pola Berdasarkan Scatterplot
b.
GCV
t1
t2
t3
t4
0,0010628
548,71
34,53
18,16
2,65
0,0011558
715,43
38,24
19,19
4,29
0,0013117
882,14
41,96
20,21
5,92
0,0009344
1048,86
45,67
21,24
7,55
0,0009242
1215,57
49,38
22,26
9,18
0,0009001
1382,29
53,09
23,28
10,82
0,0008940 0,0009008
1549,00 1715,71
56,80 60,52
24,31 25,33
12,45 14,08
0,0009164
1882,43
64,23
26,36
15,72
0,0009310
2049,14
67,94
27,38
17,35
Pemilihan Titik Knot dengan Dua Titik Knot Setelah dilakukan pemilihan titik knot dengan satu titik knot, selanjutnya dilakukan pemilihan titik knot optimum menggunakan dua titik knot. Estimasi model regresi semiparametrik spline dengan dua titik knot pada IPKM Provinsi Jawa Timur pada pola berdasarkan scatterplot adalah sebagai berikut.
39 yi 0 1 x1i 2 x2i 1t1i 2 t1i 1 3 t1i 2 4t2i 1
1
5 t2i 1 6 t2i 2 7t3i 8 t3i 1 9 t3i 2 1
1
1
1
10t4i 11 t4i 1 12 t2i 2 1
1
Sepuluh nilai GCV minimum untuk model regresi semiparametrik spline dengan empat komponen nonparametrik menggunakan dua titik knot disajikan pada Tabel 4.4. Berdasarkan Tabel 4.4 diketahui bahwa nilai GCV minimum yang diperoleh dari pemodelan regresi semiparametrik spline dengan dua titik knot pada empat komponen nonparametrik adalah 0,0007568. Dengan titik knot optimum yaitu variabel kepadatan penduduk (t1) pada titik 382 dan 1048,86; variabel Angka Kematian Ibu (t2) pada titik 30,82 dan 45,67; variabel persentase rumah tangga ber-PHBS (t3) pada titik 17,14 dan 21,24; serta variabel persentase rumah sehat (t4) pada titik 1,02 dan 7,55. Tabel 4. 4 Nilai GCV dengan Dua Titik Knot Pola Berdasarkan Scatterplot
GCV 0,0008127 0,0008372 0,0007820 0,0007568 0,0007973 0,0008772
t1 382,00
t2 30,82
t3 17,14
t4 1,02
548,71
34,53
18,16
2,65
382,00
30,82
17,14
1,02
715,43
38,24
19,19
4,29
382,00
30,82
17,14
1,02
882,14
41,96
20,21
5,92
382,00
30,82
17,14
1,02
1048,86
45,67
21,24
7,55
382,00
30,82
17,14
1,02
1215,57
49,38
22,26
9,18
382,00
30,82
17,14
1,02
1382,29
53,09
23,28
10,82
40 Tabel 4. 4 Nilai GCV dengan Dua Titik Knot Pola Berdasarkan Scatterplot (lanjutan)
GCV 0,0008909 0,0008909 0,0008714 0,0008714
t1 382,00
t2 30,82
t3 17,14
t4 1,02
1549,00
56,80
24,31
12,45
382,00
30,82
17,14
1,02
1715,71
60,52
25,33
14,08
382,00
30,82
17,14
1,02
1882,43
64,23
26,36
15,72
382,00
30,82
17,14
1,02
2049,14
67,94
27,38
17,35
Nilai GCV dengan 2 titik knot untuk pola berdasarkan scatterplot bernilai lebih kecil dari satu titik knot sehingga model dengan dua titik knot lebih baik daripada dengan satu titik knot. Akan tetapi masih harus dicobakan dengan menggunakan tiga titik knot dan kombinasi titik knot untuk mendapatkan kemungkinan model yang lebih baik. c.
Pemilihan Titik Knot dengan Tiga Titik Knot Estimasi model regresi semiparametrik spline dengan tiga titik knot pada IPKM Provinsi Jawa Timur pada pola berdasarkan scatterplot adalah sebagai berikut. yi 0 1 x1i 2 x2i 1t1i 2 t1i 1 3 t1i 2 4 t1i 3 1
1
1
5t2i 6 t2i 1 7 t2i 2 8 t2i 3 9t3i 10 t3i 1 1
1
1
1
11 t3i 2 12 t3i 3 13t4i 14 t4i 1 15 t2i 2 1
1
1
1
16 t4i 3 1
Berikut ini merupakan lima nilai GCV minimum untuk model regresi semiparametrik spline dengan empat komponen nonparametrik menggunakan tiga titik knot yang disajikan dalam Tabel 4.5.
41 Tabel 4. 5 Nilai GCV dengan Tiga Titik Knot Pola Berdasarkan Scatterplot
GCV 0,0016185
0,0005271
0,0005665
0,0006434
0,0007953
t2
t3
t4
3049,43
t1
90,21
33,53
27,15
3216,14
93,92
34,55
28,78
8384,29
209,00
66,30
79,40
3049,43
90,21
33,53
27,15
3382,86
97,64
35,57
30,41
3549,57
101,35
36,60
32,04
3049,43
90,21
33,53
27,15
3382,86
97,64
35,57
30,41
3716,29
105,06
37,62
33,68
3049,43
90,21
33,53
27,15
3382,86
97,64
35,57
30,41
3883,00
108,77
38,65
35,31
3049,43
90,21
33,53
27,15
3382,86
97,64
35,57
30,41
4049,71
112,48
39,67
36,94
Berdasarkan Tabel 4.5 diketahui bahwa nilai GCV minimum yang diperoleh dari pemodelan regresi semiparametrik spline dengan empat komponen nonparametrik menggunakan tiga titik knot adalah 0,0005271. Dengan titik knot optimum yaitu variabel kepadatan penduduk (t1) pada titik 3049,43, 3382,86, dan 3549,57; variabel Angka Kematian Ibu (t2) pada titik 90,21, 97,64, dan 101,35; variabel persentase rumah tangga ber-PHBS (t3) pada titik 33,53, 35,57, dan 36,6; serta variabel persentase rumah sehat (t4) pada titik 27,15, 30,41, dan 32,04. Nilai GCV dengan tiga titik knot untuk pola berdasarkan scatterplot bernilai lebih kecil dari satu titik knot dan dua titik knot sebelumnya sehingga model dengan tiga titik knot lebih baik. Akan tetapi masih harus dicobakan lagi menggunakan kombinasi titik knot untuk mendapatkan kemungkinan model yang lebih baik.
42 d.
Pemilihan Titik Knot dengan Kombinasi Titik Knot Pemilihan titik knot optimum dengan satu titik knot, dua titik knot serta tiga titik knot telah dilakukan. Selanjutnya, dilakukan pemilihan titik knot optimum dengan kombinasi titik knot karena terdapat kemungkinan bahwa setiap pola data memiliki jumlah titik knot optimum yang berbeda-beda. Lima nilai GCV minimum untuk pola berdasarkan scatterplot yang disajikan dalam Tabel 4.6. Tabel 4. 6 Nilai GCV dengan Kombinasi Titik Knot Pola Berdasarkan Scatterplot
GCV
t1 1549,00
0,0007397 1549,00 0,0007409 1549,00 0,0005568 382,00 0,0009158
t3
t4
90,21
33,53
12,45
97,64
35,57
101,35 90,21
36,60 33,53
1,02
97,64
35,57
7,55
101,35 90,21
36,60 33,53
27,15
97,64
35,57
30,41
101,35 56,80
36,60 24,31
32,04 12,45
56,80
24,31
1,02
1048,86 382,00
0,0008876
t2
1048,86
7,55
Berdasarkan Tabel 4.6 diketahui bahwa nilai GCV minimum dari kombinasi knot untuk empat komponen nonparametrik adalah 0,0005568. Nilai GCV tersebut dihasilkan apabila digunakan kombinasi knot (1,3,3,3). Titik knot yang digunakan adalah 1549 untuk variabel kepadatan penduduk (t1); 90,21, 97,64, dan 101,35 untuk variabel Angka Kematian Ibu (t2); titik 33,53, 35,57, dan
43 36,6 untuk variabel persentase rumah tangga ber-PHBS (t3); serta 27,15, 30,41, dan 32,04 untuk variabel persentase rumah sehat (t4). 2.
Pola Berdasarkan Uji RESET Berikut akan ditampilkan nilai GCV dengan satu titik knot, dua titik knot, tiga titik knot, dan kombinasi titik knot untuk pola berdasarkan uji RESET. a. Pemilihan Titik Knot dengan Satu Titik Knot Estimasi model regresi semiparametrik spline dengan satu titik knot pada IPKM Provinsi Jawa Timur untuk pola berdasarkan uji RESET adalah sebagai berikut. 1 1 yi 0 1 x1i 2 x2i 3 x3i 4 x4i 1t1i 2 t1i 1 3t2i 4 t2i 1 Berdasarkan Tabel 4.7, diketahui bahwa nilai GCV pada satu titik knot dengan dua komponen nonparametrik yang paling minimum adalah 0,000837 dengan titik knot optimum yaitu variabel kepadatan penduduk (t1) pada titik 3382,86; serta variabel persentase rumah tangga ber-PHBS (t2) pada titik 35,57. Tabel 4. 7 Nilai GCV dengan Satu Titik Knot Pola Berdasarkan Uji RESET
b.
GCV
t1
t2
GCV
t1
t2
0,000844
2215,86
28,40
0,000842
3049,43
33,53
0,000844
2382,57
29,43
0,000839
3216,14
34,55
0,000845
2549,29
30,45
0,000846
2716,00
31,48
0,000837 0,000840
3382,86 3549,57
35,57 36,60
0,000846
2882,71
32,50
0,000843
3716,29
37,62
Pemilihan Titik Knot dengan Dua Titik Knot Setelah dilakukan pemilihan titik knot dengan satu titik knot, selanjutnya akan dilakukan pemilihan titik knot optimum menggunakan dua titik knot. Estimasi model regresi semiparametrik spline dengan dua titik knot pada IPKM Provinsi Jawa Timur untuk pola berdasarkan uji RESET adalah sebagai berikut.
44 yi 0 1 x1i 2 x2i 3 x3i 4 x4i 1t1i 2 t1i 1 3 t1i 2 1
1
4t2i 5 t2i 1 6 t2i 2 1
1
Tabel 4.8 menunjukkan bahwa nilai GCV minimum yang diperoleh dari pemodelan regresi semiparametrik spline dengan dua titik knot pada dua komponen nonparametrik adalah 0,0007515. Dengan titik knot optimum, yaitu variabel kepadatan penduduk (t1) pada titik 2549,29 dan 5126,71; serta variabel persentase rumah tangga ber-PHBS (t2) pada titik 30,45 dan 46,84. Tabel 4. 8 Nilai GCV dengan Dua Titik Knot Pola Berdasarkan Uji RESET
GCV
t1
t2
0,0007691
2549,29
30,45
4549,86
42,74
2549,29
30,45
4716,57
43,77
2549,29
30,45
4883,29
44,79
2549,29
30,45
5050,00
45,81
2549,29
30,45
0,0007605 0,0007548 0,0007520 0,0007515
5216,71
46,84
GCV 0,0007577 0,0007739 0,0008123 0,0008373 0,0008592
t1
t2
2549,29
30,45
5383,43
47,86
2549,29
30,45
5550,14
48,89
2549,29
30,45
5716,86
49,91
2549,29
30,45
5883,57
50,93
2549,29
30,45
6050,29
51,96
Nilai GCV dengan 2 titik knot baik untuk pola berdasarkan uji RESET bernilai lebih kecil dari satu titik knot sehingga model dengan dua titik knot lebih baik daripada dengan satu titik knot. Akan tetapi masih harus dicobakan dengan menggunakan tiga titik knot dan kombinasi titik knot untuk mendapatkan kemungkinan model yang lebih baik. c.
Pemilihan Titik Knot dengan Tiga Titik Knot Estimasi model regresi semiparametrik spline dengan tiga titik knot pada IPKM Provinsi Jawa Timur untuk pola berdasarkan uji RESET adalah sebagai berikut.
45 yi 0 1 x1i 2 x2i 3 x3i 4 x4i 1t1i 2 t1i 1 3 t1i 2 1
1
4 t1i 3 5t2i 6 t2i 1 7 t2i 2 8 t2i 3 1
1
1
1
Tabel 4. 9 Nilai GCV dengan Tiga Titik Knot Pola Berdasarkan Uji RESET
GCV
t1 548,714
0,0009725
0,0009741
0,0009743
0,0009749
0,0009747
t2
GCV
18,164
715,429
19,188
7550,714
61,176
548,714
18,164
715,429
19,188
7717,429
62,200
548,714
18,164
715,429
19,188
7884,143
63,224
548,714
18,164
715,429
19,188
8050,857
64,248
548,714
18,164
715,429
19,188
8217,571
65,272
0,0009727
0,0008793
0,0005639
0,0006155
0,0007454
t1
t2
548,714
18,164
715,429
19,188
8384,286
66,296
548,714
18,164
882,143
20,212
1048,857
21,236
548,714
18,164
882,143
20,212
1215,571
22,260
548,714
18,164
882,143
20,212
1382,286
23,284
548,714
18,164
882,143
20,212
1549,000
24,309
Tabel 4.9 menunjukkan bahwa nilai GCV minimum yang diperoleh dari pemodelan regresi semiparametrik spline dengan tiga titik knot pada dua komponen nonparametrik adalah 0,0005639. Dengan titik knot optimum, yaitu variabel kepadatan penduduk (t1) pada titik 548,714, 882,143, dan 1215,571; serta variabel persentase rumah tangga ber-PHBS (t2) pada titik 18,164, 20,212, dan 22,26. Nilai GCV dengan tiga titik knot untuk pola berdasarkan uji RESET bernilai lebih kecil dari satu titik knot dan dua titik knot sebelumnya sehingga model dengan tiga titik knot lebih baik. Akan
46 tetapi masih harus dicobakan lagi menggunakan kombinasi titik knot untuk mendapatkan kemungkinan model yang lebih baik. d.
Pemilihan Titik Knot dengan Kombinasi Titik Knot Setelah melakukan pemilihan titik knot optimum dengan satu titik knot, dua titik knot serta tiga titik knot, selanjutnya adalah pemilihan titik knot optimum dengan kombinasi titik knot karena terdapat kemungkinan bahwa setiap pola data memiliki jumlah titik knot optimum yang berbeda-beda. Berikut ini merupakan nilai GCV minimum untuk pola berdasarkan uji RESET yang disajikan dalam Tabel 4.10. Tabel 4. 10 Nilai GCV dengan Kombinasi Titik Knot Pola Berdasarkan Uji RESET
GCV
t1 3382,86
t2 35,57
0,0008372
3382,86
30,45 46,84
3382,86
0,0009511
18,16 20,21 22,26
t1 2549,29
0,0008954
0,0008894
0,0005278
GCV
0,0005562
t2 35,57
5216,71 2549,29
30,45
5216,71
46,84
2549,29
18,16
5216,71
20,21 22,26
Berdasarkan Tabel 4.10, diketahui bahwa bahwa nilai GCV minimum dari kombinasi knot untuk dua komponen nonparametrik adalah 0,0005278. Nilai GCV tersebut dihasilkan apabila digunakan kombinasi knot (1,3). Titik knot yang digunakan adalah 3382,86 untuk variabel kepadatan penduduk (t1); serta 18,164, 20,212, dan 22,26 untuk variabel persentase rumah tangga berPHBS (t2). Selanjutnya akan dilakukan pemilihan model terbaik berdasarkan nilai GCV yang paling minimum.
47 4.2.3 Pemilihan Model Terbaik Setelah mendapatkan nilai GCV minimum pada pemilihan titik knot optimum dengan satu titik knot, dua titik knot, tiga titik knot, dan kombinasi knot, selanjutnya dilakukan pemilihan model terbaik dengan membandingkan nilai GCV. Tabel 4. 11 Perbandingan Nilai GCV Minimum
Keterangan Pola Berdasarkan Scatterplot
Pola Berdasarkan Uji RESET
Jumlah Knot
GCV Minimum
Satu Titik Knot
0,0008940
Dua Titik Knot
0,0007568
Tiga Titik Knot Kombinasi Knot (1,3,3,3)
0,0005271 0,0005568
Satu Titik Knot
0,000837
Dua Titik Knot
0,0007515
Tiga Titik Knot
0,0005639
Kombinasi Knot (1,3)
0,0005278
Berdasarkan Tabel 4.11 diketahui bahwa pemodelan yang menghasilkan nilai GCV paling minimum adalah pemodelan regresi semiparametrik spline pada pola berdasarkan scatterplot dengan menggunakan tiga titik knot pada empat variabel nonparametrik. Oleh karena itu, diputuskan bahwa model terbaik yang akan dipilih adalah model regresi semiparametrik spline dengan menggunakan tiga titik knot untuk masing-masing variabel komponen nonparametrik, yaitu variabel kepadatan penduduk (t1), Angka Kematian Ibu (t2), persentase rumah tangga ber-PHBS (t3), dan persentase rumah sehat (t4). 4.2.4 Pemodelan IPKM Provinsi Jawa Timur dengan Menggunakan Titik Knot Optimum Untuk melakukan pemodelan IPKM Provinsi Jawa Timur, maka digunakan titik knot optimum yang berdasarkan GCV minimum. Estimasi model untuk regresi semiparametrik spline adalah sebagai berikut.
48 yˆ 0,94873 0, 0083 x1 0, 00262 x2 0, 00001t1 0, 00008 t1 3049, 43 1
0, 00003 t1 3382,86 0, 00008 t1 3549,57 0, 00046t2 1
1
0, 00342 t2 90, 21 0, 01149 t2 97, 64 0, 00929 t2 101,35 1
1
1
0, 00801t3 0, 07391 t3 33,53 0,142 t3 35,57 0, 07571 t3 36, 6 1
1
1
0, 00023t4 0, 0596 t4 27,15 0,1876 t4 30, 41 0,12836 t4 32, 04 1
1
1
4.2.5 Pengujian Signifikansi Parameter Model Regresi Semiparametrik Spline Pengujian parameter dilakukan untuk mengetahui apakah parameter yang telah didapatkan dari hasil pemodelan dengan regresi semiparametrik spline berpengaruh secara signifikan terhadap variabel IPKM atau tidak. Pengujian dilakukan dengan dua tahap, yaitu pengujian parameter secara serentak dan individu. Apabila pada pengujian secara serentak menghasilkan bahwa parameter berpengaruh signifikan terhadap IPKM maka dilanjutkan pada pengujian secara individu untuk mengetahui pengaruh signifikansi dari tiap-tiap parameter terhadap IPKM. 1. Pengujian Serentak Pengujian secara serentak dilakukan pada parameter model regresi terhadap variabel IPKM secara bersama-sama atau serentak. Hipotesis yang digunakan untuk pengujian secara serentak adalah sebagai berikut. H0 : 1 2 1 2 ... 16 0 H1 : minimal terdapat satu j 0 atau l 0 , j 1,2, l 1,2,...,16 Hasil pengujian secara serentak ditampilkan pada Tabel 4.12. Tabel 4. 12 Hasil Pengujian Serentak
Sumber Regresi
df 18
SS 0,088347
MS 0,004908
Error
19
0,005548
0,000292
Total
37
0,093895
F 16,80864
P-value 0,00
49 Berdasarkan Tabel 4.12, diketahui bahwa diperoleh nilai F sebesar 16,80864 dan P-value sebesar 0. Dengan taraf signifikansi (α) sebesar 5% maka didapatkan keputusan tolak H0 karena nilai F lebih besar dari F(0,05;18,19), yaitu 16,80864 > 2,18 dan P-value < α, sehingga dapat disimpulkan bahwa minimal terdapat satu parameter pada model yang signifikan. Nilai R2 yang diperoleh adalah 94,09%. Hal ini menunjukkan bahwa model mampu menjelaskan keragaman IPKM di Provinsi Jawa Timur sebesar 94,09%, sedangkan sisanya dijelaskan oleh variabel lain. 2. Pengujian Individu Hasil dari pengujian serentak adalah minimal terdapat satu parameter pada model yang signifikan, sehingga dapat dilanjutkan pada pengujian parameter secara individu. Pengujian individu dilakukan untuk mengetahui parameter mana saja yang berpengaruh signifikan terhadap model regresi. Berikut merupakan perumusan hipotesis untuk komponen parametrik. H0 : j 0 H1 : j 0, j 1, 2 Sedangkan hipotesis yang digunakan untuk komponen nonparametrik adalah sebagai berikut. H0 : l 0 H1 : l 0, l 1, 2,...,16 Hasil dari pengujian parameter model regresi secara individu disajikan pada Tabel 4.13. Berdasarkan Tabel 4.13, dengan taraf signifikansi (α) sebesar 5%, terdapat empat variabel yang signifikan dan dua variabel yang tidak signifikan. Dari 19 parameter pada model regresi semiparametrik spline, terdapat 10 parameter yang signifikan dan 9 parameter yang tidak signifikan terhadap model. Empat variabel yang signifikan adalah variabel Angka Kematian Bayi (x2), kepadatan penduduk (t1), persentase rumah tangga ber-PHBS (t3), dan persentase rumah sehat (t4). Sedangkan variabel yang tidak signifikan adalah persentase
50 penduduk miskin (x1) dan Angka Kematian Ibu (t2). Berdasarkan hasil yang telah diperoleh, terdapat dua variabel yang tidak signifikan, maka dilakukan pemodelan kembali dengan menghilangkan variabel yang tidak signifikan tersebut, sehingga terdapat empat variabel prediktor, yang terdiri dari satu variabel komponen parametrik dan empat variabel komponen nonparametrik. Tabel 4. 13 Hasil Pengujian Individu
Variabel Konstan
Parameter 0
Estimasi 0,94873
t 12,699
Keputusan Tolak
x1
1
-0,00083
-0,863
Gagal tolak
x2
2
-0,00262
-6,883
Tolak
1
0,00001
1,695
Gagal tolak
2
0,00008 -0,00003
1,418 -2,580
Gagal tolak Tolak
-0,00008
-1,899
Gagal tolak
0,00046
1,393
Gagal tolak
-0,00342
-0,935
Gagal tolak
0,01149 -0,00929
1,147 -1,375
Gagal tolak Gagal tolak
-0,00801
-2,690
Tolak
0,07391
3,010
Tolak
-0,14200
-3,415
Tolak
0,07571
3,349
Tolak
0,00023
0,314
Gagal tolak
t1
t2
t3
t4
3 4
5 6 7
8 9 10 11 12 13 14
-0,05960
-3,15
Tolak
15
0,18760
3,360
Tolak
16
-0,12836
-3,414
Tolak
Kesimpulan Signifikan Tidak Signifikan Signifikan
Signifikan
Tidak Signifikan
Signifikan
Signifikan
51 4.2.6 Pemodelan IPKM Provinsi Jawa Timur dengan Empat Variabel Pengujian individu menghasilkan dua variabel tidak signifikan sehingga dilakukan pemodelan kembali dengan menghilangkan variabel yang tidak signifikan tersebut dari model. Pemodelan dilakukan menggunakan satu knot, dua knot, tiga knot, dan kombinasi knot. Dari keempat titik knot tersebut, untuk membentuk model digunakan titik knot optimum yang diperoleh dari nilai GCV yang minimum. 1.
Pemilihan Titik Knot dengan Satu Titik Knot Berikut merupakan sepuluh nilai GCV minimum dengan tiga komponen nonparametrik menggunakan satu titik knot yang disajikan pada Tabel 4.14. Tabel 4.14 menunjukkan bahwa nilai GCV pada satu titik knot dengan tiga komponen nonparametrik yang paling minimum adalah 0,0008246 dengan titik knot optimum yaitu variabel kepadatan penduduk (t1) pada titik 1549; variabel persentase rumah tangga ber-PHBS (t3) pada titik 24,31; dan variabel persentase rumah sehat (t4) pada titik 12,45. Tabel 4. 14 Nilai GCV dengan Satu Titik Knot Pada Model Empat Variabel
GCV
t1
t3
t4
0,0010184
548,71
18,16
2,65
0,0011310
715,43
19,19
4,29
0,0012847
882,14
20,21
5,92
0,0008746
1048,86
21,24
7,55
0,0008488
1215,57
22,26
9,18
0,0008272
1382,29
23,28
10,82
0,0008246 0,0008248
1549,00 1715,71
24,31 25,33
12,45 14,08
0,0008269
1882,43
26,36
15,72
0,0008329
2049,14
27,38
17,35
52 2.
Pemilihan Titik Knot dengan Dua Titik Knot Setelah dilakukan pemilihan titik knot dengan satu titik knot, selanjutnya dilakukan pemilihan titik knot optimum menggunakan dua titik knot. Lima nilai GCV minimum untuk model regresi semiparametrik spline dengan tiga komponen nonparametrik menggunakan dua titik knot disajikan pada Tabel 4.15. Tabel 4. 15 Nilai GCV dengan Dua Titik Knot Pada Model Empat Variabel
GCV 0,0008971 0,0008971 0,0006411 0,0007975 0,0008803
t1 2215,86
t3 28,40
t4 18,98
8384,29
66,30
79,40
2215,86
28,40
18,98
8551,00
67,32
81,03
2382,57
29,43
20,61
2549,29
30,45
22,25
2382,57
29,43
20,61
2716,00
31,48
23,88
2382,57
29,43
20,61
2882,71
32,50
25,51
Berdasarkan Tabel 4.15 diketahui bahwa nilai GCV minimum yang diperoleh dari pemodelan regresi semiparametrik spline dengan dua titik knot pada tiga komponen nonparametrik adalah 0,0006411. Dengan titik knot optimum yaitu variabel kepadatan penduduk (t1) pada titik 2382,57 dan 2549,29; variabel persentase rumah tangga ber-PHBS (t3) pada titik 29,43 dan 31,48; serta variabel persentase rumah sehat (t4) pada titik 20,61 dan 22,25. Nilai GCV dengan 2 titik knot bernilai lebih kecil dari satu titik knot sehingga model dengan dua titik knot lebih baik daripada dengan satu titik knot. Akan tetapi masih harus dicobakan dengan menggunakan tiga titik knot dan kombinasi titik knot untuk mendapatkan kemungkinan model yang lebih baik.
53 3.
Pemilihan Titik Knot dengan Tiga Titik Knot Nilai GCV minimum untuk model regresi semiparametrik spline dengan tiga komponen nonparametrik menggunakan tiga titik knot yang disajikan dalam Tabel 4.16. Tabel 4. 16 Nilai GCV dengan Tiga Titik Knot Pada Model Empat Variabel
GCV 0,0012818
0,0011922
0,0004716
0,0004339
0,0004478
t1
t3
t4
2716,00
31,48
23,88
3216,14
34,55
28,78
8217,57
65,27
77,76
2716,00
31,48
23,88
3216,14
34,55
28,78
8384,29
66,30
79,40
2716,00
31,48
23,88
3382,86
35,57
30,41
3549,57
36,60
32,04
2716,00
31,48
23,88
3382,86
35,57
30,41
3716,29
37,62
33,68
2716,00
31,48
23,88
3382,86
35,57
30,41
3883,00
38,65
35,31
Tabel 4.16 menunjukkan bahwa nilai GCV minimum yang diperoleh dari pemodelan regresi semiparametrik spline dengan tiga komponen nonparametrik menggunakan tiga titik knot adalah 0,0004339. Dengan titik knot optimum yaitu variabel kepadatan penduduk (t1) pada titik 2716, 3382,86, dan 3716,29; variabel persentase rumah tangga ber-PHBS (t3) pada titik 31,48, 35,57, dan 37,62; serta variabel persentase rumah sehat (t4) pada titik 23,88, 30,41, dan 33,68. Nilai GCV dengan tiga titik knot bernilai lebih kecil dari satu titik knot dan dua titik knot sebelumnya sehingga model dengan
54 tiga titik knot lebih baik. Akan tetapi masih harus dicobakan lagi menggunakan kombinasi titik knot untuk mendapatkan kemungkinan model yang lebih baik. 4.
Pemilihan Titik Knot dengan Kombinasi Titik Knot Selanjutnya dilakukan pemilihan titik knot optimum dengan kombinasi titik knot karena terdapat kemungkinan bahwa setiap pola data memiliki jumlah titik knot optimum yang berbeda-beda. Nilai GCV minimum untuk pola berdasarkan scatterplot yang disajikan dalam Tabel 4.17. Tabel 4. 17 Nilai GCV dengan Kombinasi Titik Knot Pada Model Empat Variabel
GCV 0,0009235
0,0005420
0,0009060
t1
t3
t4
2382,57
31,48
20,61
2549,29
35,57
22,25
2382,57
37,62 31,48
23,88
2549,29
35,57
30,41
2716,00
37,62 24,31
33,68 12,45
24,31
20,61
3382,86 3716,29 2716,00
0,0009493
3382,86 3716,29 2716,00
0,0004438
22,25 24,31
23,88
3382,86
30,41
3716,29
33,68
Berdasarkan Tabel 4.17 diketahui bahwa nilai GCV minimum dari kombinasi knot untuk empat komponen nonparametrik adalah 0,0004438. Nilai GCV tersebut dihasilkan apabila digunakan kombinasi knot (3,1,3). Titik knot yang digunakan adalah 2716,
55 3382,86, dan 3716,29 untuk variabel kepadatan penduduk (t1); titik 24,31 untuk variabel persentase rumah tangga ber-PHBS (t3); dan 23,88, 30,41, dan 33,68 untuk variabel persentase rumah sehat (t4). 4.2.7 Pemilihan Model IPKM Provinsi Jawa Timur Terbaik dengan Empat Variabel Setelah mendapatkan nilai GCV minimum pada pemilihan titik knot optimum dengan satu titik knot, dua titik knot, tiga titik knot, dan kombinasi knot, selanjutnya dilakukan pemilihan model terbaik dengan membandingkan nilai GCV. Tabel 4. 18 Perbandingan Nilai GCV Minimum Pada Model Empat Variabel
Jumlah Knot
GCV Minimum
Satu Titik Knot
0,0008246
Dua Titik Knot
0,0006411
Tiga Titik Knot Kombinasi Knot (3,1,3)
0,0004339 0,0004438
Berdasarkan Tabel 4.18 diketahui bahwa pemodelan yang menghasilkan nilai GCV paling minimum adalah pemodelan regresi semiparametrik spline menggunakan tiga titik knot pada tiga variabel nonparametrik, yaitu variabel kepadatan penduduk (t1), persentase rumah tangga ber-PHBS (t3), dan persentase rumah sehat (t4). 4.2.8 Pemodelan IPKM Provinsi Jawa Timur Menggunakan Titik Knot Optimum dengan Empat Variabel Untuk melakukan pemodelan IPKM Provinsi Jawa Timur, maka digunakan titik knot optimum yang berdasarkan GCV minimum. Estimasi model untuk regresi semiparametrik spline adalah sebagai berikut. yˆ 0, 73408 0, 00246 x2 0, 00006t1 0, 00091 t1 2716 0, 0024 t1 3382,86 1
1
0, 00156 t1 3716, 29 0, 00057t3 0, 00832 t3 31, 48 0, 02661 t3 35,57 1
1
1
0, 01873 t3 37, 62 0, 00108t4 0, 02051 t4 23,88 0, 061 t4 30, 41 1
0, 04178 t4 33, 68 1
1
1
56 4.2.9 Pengujian Signifikansi Parameter Model Regresi Semiparametrik Spline dengan Empat Variabel Pengujian parameter dilakukan untuk mengetahui apakah parameter yang telah didapatkan dari hasil pemodelan dengan regresi semiparametrik spline dengan empat variabel berpengaruh secara signifikan terhadap variabel IPKM atau tidak. Pengujian dilakukan dengan dua tahap, yaitu pengujian parameter secara serentak dan individu. 1. Pengujian Serentak Pengujian secara serentak dilakukan pada parameter model regresi terhadap variabel IPKM secara bersama-sama atau serentak. Hipotesis yang digunakan untuk pengujian secara serentak adalah sebagai berikut. H0 : 2 1 2 ... 12 0 H1 : minimal terdapat satu 2 0 atau l 0 , l 1, 2,...,12 Hasil pengujian secara serentak ditampilkan pada Tabel 4.19. Tabel 4. 19 Hasil Pengujian Serentak Pada Model Empat Variabel
Sumber Regresi
df 13
SS 0,087318
MS 0,006717
Error
24
0,006577
0,000274
Total
37
0,093895
F 24,50918
P-value 0,00
Berdasarkan Tabel 4.19, diketahui bahwa diperoleh nilai F sebesar 24,50918 dan P-value sebesar 0. Dengan taraf signifikansi (α) sebesar 5% maka didapatkan keputusan tolak H0 karena nilai F lebih besar dari F(0,05;13,24), yaitu 16,80864 > 2,155 dan P-value < α, sehingga dapat disimpulkan bahwa minimal terdapat satu parameter pada model yang signifikan. Nilai R2 yang diperoleh adalah 92,99%. Hal ini menunjukkan bahwa model mampu menjelaskan keragaman IPKM di Provinsi Jawa Timur sebesar 92,99%, sedangkan sisanya dijelaskan oleh variabel lain.
57 2. Pengujian Individu Pengujian serentak menghasilkan bahwa minimal terdapat satu parameter pada model yang signifikan, sehingga dapat dilanjutkan pada pengujian parameter secara individu. Berikut merupakan perumusan hipotesis untuk komponen parametrik. H0 : 2 0 H1 : 2 0 Sedangkan hipotesis yang digunakan untuk komponen nonparametrik adalah sebagai berikut. H0 : l 0 H1 : l 0, l 1, 2,...,12 Hasil dari pengujian parameter model regresi secara individu disajikan pada Tabel 4.20. Tabel 4. 20 Hasil Pengujian Individu Pada Model Empat Variabel
Variabel Konstan
Parameter 0
Estimasi 0,73408
t 13,8592
Keputusan Tolak
Kesimpulan Signifikan
x2
2
-0,00246
-8,4150
Tolak
Signifikan
1
0,00006
3,9135
Tolak
2
-0,00091 0,00240
-3,0295 3,0541
Tolak Tolak
-0,00156
-3,1215
Tolak
-0,00057
-0,2756
Gagal tolak
0,00832
1,0270
Gagal tolak
-0,02661
-1,9891
Gagal tolak
0,01873
2,2049
Tolak
0,00108
1,6401
Gagal tolak
-0,02051
-5,7326
Tolak
0,06100
6,3390
Tolak
-0,04178
-6,2259
Tolak
t1
t3
t4
3 4 5 6 7 8 9 10 11 12
Signifikan
Signifikan
Signifikan
58 Berdasarkan Tabel 4.20, dengan taraf signifikansi (α) sebesar 5%, semua variabel prediktor berpengaruh signifikan terhadap model. Dari 14 parameter pada model regresi semiparametrik spline, terdapat 10 parameter yang signifikan dan 4 parameter yang tidak signifikan. Variabel yang signifikan adalah variabel Angka Kematian Bayi (x2), kepadatan penduduk (t1), persentase rumah tangga ber-PHBS (t3), dan persentase rumah sehat (t4). 4.2.10 Pengujian Asumsi Residual Pengujian asumsi residual dilakukan untuk mengetahui apakah residual yang dihasilkan dari model regresi telah memenuhi asumsi identik, independen, dan berdistribusi normal (IIDN) atau tidak. Apabila suatu model regresi dengan kriteria model terbaik dan parameter signifikan namun tidak memenuhi asumsi IIDN, maka model regresi tidak layak untuk menggambarkan hubungan antara variabel prediktor dengan variabel respon. 1. Pengujian Asumsi Identik Pengujian asumsi identik pada residual digunakan untuk mengetahui apakah terjadi kasus heteroskedastisitas atau variansi residual dari model harus homogen. Pengujian asumsi residual identik dilakukan dengan menggunakan uji Glejser. Hasil uji Glejser akan ditampilkan pada Tabel 4.21 sebagai berikut. Tabel 4. 21 Hasil Pengujian Glejser
Sumber
df
SS
MS
F
P-value
Regresi
13
0,000723
0,000056
0,934
0,536
Error
24
0,001429
0,000060
Total
37
0,002151
Berdasarkan Tabel 4.21 diketahui bahwa nilai F yang dihasilkan adalah 0,934 dan P-value sebesar 0,536. Dengan taraf signifikansi (α) sebesar 5% maka didapatkan keputusan gagal tolak H0 karena nilai F < F(0,05;13,24), yaitu 0,934 < 2,155 dan P-value > α, sehingga dapat disimpulkan bahwa tidak terjadi heteroskedastisitas atau
59 dengan kata lain variansi antar residual sama. Hal ini berarti bahwa asumsi residual identik telah terpenuhi. 2. Asumsi Independen Residual yang independen berarti bahwa tidak terjadi autokorelasi antar residual. Salah satu cara untuk mengetahui ada atau tidaknya autokorelasi antar residual adalah dengan menggunakan Plot Autocorrelation Function (ACF). Apabila ada autokorelasi yang keluar dari batas atas maupun batas bawah interval konfidensi maka dapat disimpulkan bahwa terdapat autokorelasi antar residual. 1.0 0.8
Autocorrelation
0.6 0.4 0.2 0.0 -0.2 -0.4 -0.6 -0.8 -1.0 1
5
10
15
20 Lag
25
30
35
Gambar 4. 8 Plot ACF Residual
Gambar 4.8 menunjukkan bahwa tidak terlihat adanya nilai autokorelasi yang keluar batas interval konfidensi, sehingga dapat disimpulkan bahwa tidak terdapat autokorelasi antar residual. Hal ini mengindikasikan bahwa asumsi independen pada residual model telah terpenuhi. 3. Pengujian Asumsi Distribusi Normal Pengujian asumsi distribusi normal dilakukan untuk mengetahui apakah residual telah berdistribusi normal atau tidak. Pengujian ini dilakukan dengan uji Kolmogorov-Smirnov. Hasil uji Kolmogorov-Smirnov ditampilkan pada Gambar 4.9.
60
99
Mean StDev N KS P-Value
95 90
-1.48832E-13 0.01333 38 0.070 >0.150
80
Percent
70 60 50 40 30 20 10 5
1
-0.04
-0.03
-0.02
-0.01 0.00 residual
0.01
0.02
0.03
Gambar 4. 9 Hasil Uji Kolmogorov-Smirnov
Berdasarkan Gambar 4.9 diketahui bahwa nilai KolmogorovSmirnov sebesar 0,070 dimana nilai ini lebih kecil dibandingkan q(1-α) pada taraf signifikansi (α) sebesar 5%, yaitu 0,215 dan Pvalue yang dihasilkan lebih besar dari 0,150 sehingga gagal tolak H0. Hal ini menunjukkan bahwa residual model regresi semiparametrik spline telah memenuhi asumsi distribusi normal. 4.2.11 Interpretasi Model Regresi Semiparametrik Spline Setelah melakukan pengujian parameter model regresi dan semua asumsi residual telah terpenuhi, maka selanjutnya adalah menginterpretasi model regresi yang telah didapatkan. Model regesi yang terbentuk menggunakan titik knot optimum, yakni tiga titik knot yang ditunjukkan pada persamaan sebagai berikut. 1 1 yˆ 0, 73408 0, 00246 x2 0, 00006t1 0, 00091 t1 2716 0, 0024 t1 3382,86 0, 00156 t1 3716, 29 0, 00057t3 0, 00832 t3 31, 48 0, 02661 t3 35,57 1
1
1
0, 01873 t3 37, 62 0, 00108t4 0, 02051 t4 23,88 0, 061 t4 30, 41 1
1
1
0, 04178 t4 33, 68 1
Berikut merupakan interpretasi dari model regresi semiparametrik spline yang dilakukan terhadap variabel-variabel yang berpengaruh signifikan.
61 1.
Dengan mengasumsikan variabel prediktor selain x2 konstan, maka pengaruh variabel Angka Kematian Bayi terhadap IPKM di Provinsi Jawa Timur dapat ditulis sebagai berikut. yˆ 0, 73408 0, 00246 x2
2.
Apabila terjadi kenaikan Angka Kematian Bayi sebanyak satu persen, maka IPKM di Provinsi Jawa Timur akan turun sebesar 0,00246. Dengan mengasumsikan variabel prediktor selain t1 konstan, maka pengaruh variabel kepadatan penduduk terhadap IPKM di Provinsi Jawa Timur dapat ditulis sebagai berikut. ; t1 2716 0, 0000595t1 2, 46857 0, 00085t ; 2716 t1 3382,86 1 yˆ 5, 64928 0, 00155t1 ;3382,86 t1 3716, 29 0,15408 0, 00001t1 ; t1 3716, 29 Pengelompokan kabupaten/kota di Jawa Timur berdasarkan kepadatan penduduk disajikan pada Gambar 4.10 sebagai berikut.
Gambar 4. 10 Peta Kepadatan Penduduk di Provinsi Jawa Timur
Untuk kepadatan penduduk yang kurang dari 2716, apabila terjadi kenaikan pada kepadatan penduduk sebesar satu
62 satuan, maka akan terjadi penurunan pada IPKM sebesar 0,0000595. Kabupaten/kota di Provinsi Jawa Timur yang termasuk interval ini adalah Kabupaten Pacitan, Kabupaten Ponorogo, Kabupaten Trenggalek, Kabupaten Tulungagung, Kabupaten Blitar, Kabupaten Kediri, Kabupaten Malang, Kabupaten Lumajang, Kabupaten Jember, Kabupaten Banyuwangi, Kabupaten Bondowoso, Kabupaten Situbondo, Kabupaten Probolinggo, Kabupaten Pasuruan, Kabupaten Mojokerto, Kabupaten Jombang, Kabupaten Nganjuk, Kabupaten Madiun, Kabupaten Magetan, Kabupaten Ngawi, Kabupaten Bojonegoro, Kabupaten Tuban, Kabupaten Lamongan, Kabupaten Gresik, Kabupaten Bangkalan, Kabupaten Sampang, Kabupaten Pamekasan, Kabupaten Sumenep, dan Kota Batu. Untuk kepadatan penduduk yang berada di antara 2716 hingga 3382,86, apabila terjadi kenaikan pada kepadatan penduduk sebesar satu satuan, maka akan terjadi kenaikan pada IPKM sebesar 0,00085. Hanya Kabupaten Sidoarjo yang berada pada interval ini. Selanjutnya, untuk kepadatan penduduk yang berada di antara 3382,86 hingga 3716,29, apabila terjadi kenaikan pada kepadatan penduduk sebesar satu satuan, maka akan terjadi kenaikan pada IPKM sebesar 0,00155. Tidak ada kabupaten/kota yang berada pada interval ini. Untuk kepadatan penduduk yang lebih dari 3716,29, apabila terjadi kenaikan pada kepadatan penduduk sebesar satu satuan, maka akan terjadi penurunan pada IPKM sebesar 0,00001. Kabupaten/kota di Jawa Timur yang berada pada interval ini adalah Kota Kediri, Kota Blitar, Kota Malang, Kota Probolinggo, Kota Pasuruan, Kota Mojokerto, Kota Madiun, dan Kota Surabaya. Pemodelan ini kurang sesuai dengan teori, dimana apabila kepadatan meningkat seharusnya IPKM juga meningkat. Ketidaksesuaian ini disebabkan kemungkinan kurangnya jumlah fasilitas kesehatan dan tenaga kesehatan yang
63
3.
memadai baik secara kualitas maupun kuantitas untuk daerah yang memiliki penduduk yang padat, seperti Kabupaten Sidoarjo, Kota Kediri, Kota Blitar, Kota Malang, Kota Probolinggo, Kota Pasuruan, Kota Mojokerto, Kota Madiun, dan Kota Surabaya. Dengan mengasumsikan variabel prediktor selain t3 konstan, maka pengaruh variabel persentase rumah tangga berperilaku hidup bersih dan sehat terhadap IPKM di Provinsi Jawa Timur dapat ditulis sebagai berikut. ; t3 31, 48 0, 00057t3 0, 262 0, 00776t ;31, 48 t3 35,57 3 yˆ ;35,57 t3 37, 62 0, 68436 0, 01885t3 0, 02027 0, 00012t3 ; t3 37, 62 Untuk persentase rumah tangga ber-PHBS yang kurang dari 31,48, apabila terjadi kenaikan pada persentase rumah tangga ber-PHBS sebesar satu persen, maka akan terjadi penurunan pada IPKM sebesar 0,00057. Kabupaten/kota di Provinsi Jawa Timur yang termasuk interval ini adalah Kabupaten Trenggalek, Kabupaten Bondowoso, Kabupaten Situbondo, Kabupaten Probolinggo, Kabupaten Sampang, Kabupaten Pamekasan, dan Kota Batu. Untuk persentase rumah tangga ber-PHBS yang berada di antara 31,48 hingga 35,57, apabila terjadi kenaikan pada persentase rumah tangga ber-PHBS sebesar satu persen, maka akan terjadi kenaikan pada IPKM sebesar 0,00776. Hanya Kabupaten Ponorogo yang berada pada interval ini. Selanjutnya, untuk persentase rumah tangga ber-PHBS yang berada di antara 33,57 hingga 37,62, apabila terjadi kenaikan pada persentase rumah tangga ber-PHBS sebesar satu persen, maka akan terjadi penurunan pada IPKM sebesar 0,01885. Kabupaten Tulungagung, Kabupaten Nganjuk, dan Kota Malang berada pada interval ini. Untuk persentase rumah tangga ber-PHBS yang lebih dari 37,62, apabila terjadi kenaikan pada
64 persentase rumah tangga ber-PHBS sebesar satu persen, maka akan terjadi penurunan pada IPKM sebesar 0,00012. Kabupaten/kota di Jawa Timur yang berada pada interval ini adalah Kabupaten Pacitan, Kabupaten Blitar, Kabupaten Kediri, Kabupaten Malang, Kabupaten Lumajang, Kabupaten Jember, Kabupaten Banyuwangi, Kabupaten Pasuruan, Kabupaten Sidoarjo, Kabupaten Mojokerto, Kabupaten Jombang, Kabupaten Madiun, Kabupaten Magetan, Kabupaten Ngawi, Kabupaten Bojonegoro, Kabupaten Tuban, Kabupaten Lamongan, Kabupaten Gresik, Kabupaten Bangkalan, Kabupaten Sumenep, Kota Kediri, Kota Blitar, Kota Probolinggo, Kota Pasuruan, Kota Mojokerto, Kota Madiun, dan Kota Surabaya. Secara visual, pengelompokan kabupaten/kota di Jawa Timur berdasarkan persentase rumah tangga ber-PHBS disajikan pada Gambar 4.11 sebagai berikut.
Gambar 4. 11 Peta Persentase Rumah Tangga Ber-PHBS di Provinsi Jawa Timur
4.
Dengan mengasumsikan variabel prediktor selain t4 konstan, maka pengaruh variabel persentase rumah sehat terhadap IPKM di Provinsi Jawa Timur dapat ditulis sebagai berikut.
65 ; t4 23,88 0, 00108t4 0, 48968 0, 01942t ; 23,88 t4 30, 41 4 yˆ 1,36523 0, 04157t4 ;30, 41 t4 33, 68 ; t4 33, 68 0, 04204 0, 00021t4 Pengelompokan kabupaten/kota di Jawa Timur berdasarkan persentase rumah sehat disajikan pada Gambar 4.12.
Gambar 4. 12 Peta Persentase Rumah Sehat di Provinsi Jawa Timur
Untuk persentase rumah sehat yang kurang dari 23,88, apabila terjadi kenaikan pada persentase rumah sehat sebesar satu persen, maka akan terjadi kenaikan pada IPKM sebesar 0,00108. Kabupaten/kota di Provinsi Jawa Timur yang termasuk interval ini adalah Kabupaten Trenggalek, Kabupaten Blitar, Kabupaten Kediri, Kabupaten Malang, Kabupaten Bondowoso, Kabupaten Situbondo, Kabupaten Probolinggo, Kabupaten Pasuruan, Kabupaten Jombang, Kabupaten Ngawi, Kabupaten Sampang, Kabupaten Sumenep, Kota Kediri, dan Kota Blitar. Untuk persentase rumah sehat yang berada di antara 23,88 hingga 30,41, apabila terjadi kenaikan pada persentase rumah sehat sebesar satu persen, maka akan terjadi penurunan pada IPKM
66 sebesar 0,01942. Kabupaten Lumajang dan Kota Probolinggo berada pada interval ini. Selanjutnya, untuk persentase rumah sehat yang berada di antara 30,41 hingga 33,68, apabila terjadi kenaikan pada persentase rumah sehat sebesar satu persen, maka akan terjadi kenaikan pada IPKM sebesar 0,04157. Hanya Kabupaten Pamekasan yang berada pada interval ini. Untuk persentase rumah sehat yang lebih dari 33,68, apabila terjadi kenaikan pada persentase rumah sehat sebesar satu persen, maka akan terjadi penurunan pada IPKM sebesar 0,00021. Kabupaten/kota di Jawa Timur yang berada pada interval ini adalah Kabupaten Pacitan, Kabupaten Ponorogo, Kabupaten Tulungagung, Kabupaten Jember, Kabupaten Banyuwangi, Kabupaten Sidoarjo, Kabupaten Mojokerto, Kabupaten Nganjuk, Kabupaten Madiun, Kabupaten Magetan, Kabupaten Bojonegoro, Kabupaten Tuban, Kabupaten Lamongan, Kabupaten Gresik, Kabupaten Bangkalan, Kota Malang, Kota Pasuruan, Kota Mojokerto, Kota Madiun, Kota Surabaya, dan Kota Batu.
BAB V KESIMPULAN DAN SARAN 5.1 Kesimpulan Berdasarkan hasil analisis dan pembahasan yang telah dilakukan, dapat disimpulkan beberapa hal sebagai berikut. 1. Pada tahun 2013, IPKM terendah di Provinsi Jawa Timur terdapat di Kabupaten Pamekasan dengan nilai sebesar 0,5874. Sedangkan, Kota Madiun meraih IPKM tertinggi di Jawa Timur dengan nilai sebesar 0,79. Berdasarkan IPKM Indonesia (0,6879), terdapat 13 kabupaten/kota di Provinsi Jawa Timur yang masih berada di bawah nilai tersebut. Persentase penduduk miskin tertinggi terdapat di Kabupaten Sampang sedangkan terendah berada di Kota Batu. Angka Kematian Bayi (AKB) tertinggi terdapat pada Kabupaten Probolinggo sedangkan AKB terendah berada di Kota Blitar. Kepadatan penduduk tertinggi terdapat di Kota Surabaya dan Kabupaten Pacitan merupakan kabupaten dengan kepadatan penduduk terendah. Angka Kematian Ibu (AKI) tertinggi terdapat di Kota Probolinggo, sedangkan AKI terendah terdapat di Kota Batu. Persentase rumah tangga berperilaku hidup bersih dan sehat tertinggi terdapat di Kota Surabaya, sedangkan persentase terendah terdapat di Kabupaten Situbondo. Nilai tertinggi pada persentase rumah sehat terdapat di Kota Surabaya sedangkan persentase rumah sehat terendah terdapat di Kabupaten Blitar. 2. Model regresi semiparametrik spline terbaik dalam pemodelan IPKM di Provinsi Jawa Timur adalah dengan menggunakan tiga titik knot. Model ini memiliki empat variabel yang signifikan, yaitu Angka Kematian Bayi, kepadatan penduduk, persentase rumah tangga berperilaku hidup bersih dan sehat, dan persentase rumah sehat. Model regesi semiparametrik spline yang diperoleh adalah sebagai berikut. 67
68 yˆ 0, 73408 0, 00246 x2 0, 00006t1 0, 00091 t1 2716 0, 0024 t1 3382,86 1
1
0, 00156 t1 3716, 29 0, 00057t3 0, 00832 t3 31, 48 0, 02661 t3 35,57 1
1
1
0, 01873 t3 37, 62 0, 00108t4 0, 02051 t4 23,88 0, 061 t4 30, 41 1
1
1
0, 04178 t4 33, 68 1
Nilai koefisien determinasi yang dihasilkan dari model ini adalah sebesar 92,99%. Hal ini menunjukkan bahwa model tersebut mampu menjelaskan keragaman nilai IPKM di Provinsi Jawa Timur sebesar 92,99%, sedangkan sisanya dijelaskan oleh variabel lain. 5.2 Saran Saran yang dapat diberikan oleh penulis untuk pemerintah adalah sebaiknya memerhatikan faktor-faktor yang berpengaruh terhadap IPKM baik dari segi ekonomi, lingkungan, perilaku, dan kependudukan. Selain itu, dapat membuat program-program maupun kebijakan yang terkait dengan kesehatan sehingga dapat meningkatkan IPKM di Provinsi Jawa Timur.
DAFTAR PUSTAKA Afifah Tin, Djaja, S., & Irianto, J. (2008). Trend and Disparity of Infant Mortality Rate (IMR), Child Mortality Rate (CMR), and Under Five Mortality Rate (U5MR) and Social Economics Status in Indonesia in 1998, 2001, and 2003. Jurnal Ekologi Kesehatan, 7 (3), 838-848. Badan Pusat Statistik. (2013). Indeks Pembangunan Manusia 2013. Jakarta: Badan Pusat Statistik. Blum, H. (1974). Planning for Health; Development and Application of Social Change Theory. New York: Human Sciences Press. Budiantara, I. (2001). Estimasi Parametrik dan Nonparametrik untuk Pendekatan Kurva Regresi. Seminar Nasional Statistika V, Jurusan Statistika, FMIPA, ITS Surabaya. Budiantara, I. (2005). Model Keluarga Spline Polinomial Truncated Dalam Regresi Semiparametrik. Berkala Ilmiah MIPA, 15 (3). Budiantara, I. (2006). Model Spline dengan Knots Optimal. Jurnal Ilmu Dasar, FMIPA Universitas Jember, 7, 77-85. Budiantara, I. (2009). Spline dalam Regresi Nonparametrik dan Semiparametrik: Sebuah Pemodelan Statistika Masa Kini dan Masa Mendatang, Pidato Pengukuhan Untuk Jabatan Guru Besar pada Jurusan Statistika, Institut Teknologi Sepuluh Nopember (ITS) Surabaya. Surabaya: ITS Press. Daniel, W. (1989). Statistika Nonparametrik Terapan. Diterjemahkan oleh: Alex Tri Kantjono W. Jakarta: PT Gramedia. Diansuantari, N. (2015). Analisis Derajat Kesehatan Masyarakat Provinsi Bali dengan Menggunakan Metode Metode Multivariate Adaptive Regression Splines (MARS). Tugas Akhir S1: Jurusan Matematika, Universitas Udayana. Draper, N. R., & Smith, H. (1992). Analisis Regresi Terapan. Diterjemahkan oleh: Bambang Sumantri. Jakarta: Gramedia Pustaka Utama. 69
70 Eubank, R. (1988). Spline Smoothing and Nonparametric Regression. New York: Marcel Dekker Inc. Gujarati, D., & Porter, D. (2009). Basic Econometrics, 5th Edition. New York: The McGraw-Hill Companies, Inc. Kementerian Kesehatan RI. (2010). Indeks Pembangunan Kesehatan Masyarakat. Jakarta: Kementerian Kesehatan Republik Indonesia. Kim, T., Lee, Y., & Newbold, P. (2004). Spurious Nonlinear Regressions in Econometrics. Nottingham NG7 2RD, UK: School of Economics, University of Nottingham. Maully, A. (2014). Faktor-Faktor yang Mempengaruhi Indeks Kesehatan Kabupaten dan Kota di Provinsi Jawa Timur. Tugas Akhir D3: Jurusan Statistika, Institut Teknologi Sepuluh Nopember. Nasution, M. (2005). Manajemen Mutu Terpadu (Total Quality Management) Edisi Kedua. Bogor: Ghalia Indonesia. Notoatmodjo, S. (2002). Metodologi Penelitian Kesehatan. Jakarta: PT. Rineka Cipta. Pindyck, R., & Rubinfeld, D. (1998). Econometric Models and Economic Forecasts, 4th Ed. New York: McGraw-Hill. Prasetyo, D. (2012). Pemodelan Data Kesehatan Kabupaten Banyuwangi dengan Regresi Terboboti Geografis. Tugas Akhir S1: Departemen Statistika, Institut Pertanian Bogor. Riskiyanti, R. (2010). Analisis Regresi Multivariat Berdasarkan Faktor-Faktor yang Memengaruhi Derajat Kesehatan di Provinsi Jawa Timur. Tugas Akhir S1: Jurusan Statistika, Institut Teknologi Sepuluh Nopember. Ruppert, D., Wand, M., & Carroll, R. (2003). Semiparametric Regression. United States of America: Cambridge University Press. Tim Penyusun IPKM. (2014). Indeks Pembangunan Kesehatan Masyarakat. Jakarta: Badan Penelitian dan Pengembangan Kesehatan.
71 Walpole, R. (1995). Pengantar Statistika. Diterjemahkan oleh: Bambang Sumantri. Jakarta: PT. Gramedia Pustaka Utama.
72
(Halaman ini sengaja dikosongkan)
LAMPIRAN Lampiran 1. Data IPKM Provinsi Jawa Timur dengan FaktorFaktor yang Memengaruhi Tahun 2013 Kab/Kota Kab. Pacitan Kab. Ponorogo Kab. Trenggalek Kab. Tulungagung Kab. Blitar Kab. Kediri Kab. Malang Kab. Lumajang Kab. Jember Kab. Banyuwangi Kab. Bondowoso Kab. Situbondo Kab. Probolinggo Kab. Pasuruan Kab. Sidoarjo Kab. Mojokerto Kab. Jombang Kab. Nganjuk Kab. Madiun Kab. Magetan Kab. Ngawi Kab. Bojonegoro Kab. Tuban Kab. Lamongan Kab. Gresik Kab. Bangkalan Kab. Sampang Kab. Pamekasan Kab. Sumenep Kota Kediri Kota Blitar Kota Malang Kota Probolinggo Kota Pasuruan Kota Mojokerto Kota Madiun Kota Surabaya Kota Batu
y 0.7006 0.744 0.7127 0.7343 0.6948 0.7162 0.6897 0.6581 0.6391 0.6878 0.6083 0.6517 0.6405 0.6778 0.7395 0.7246 0.727 0.7556 0.7269 0.7339 0.7132 0.6772 0.6849 0.7116 0.7298 0.6381 0.6643 0.5874 0.6002 0.783 0.7718 0.7588 0.725 0.7388 0.749 0.79 0.7406 0.7584
x1 16.73 11.92 13.56 9.07 10.57 13.23 11.48 12.14 11.68 9.61 15.29 13.65 21.21 11.26 6.72 10.99 11.17 13.60 12.45 12.19 15.45 16.02 17.23 16.18 13.94 23.23 27.08 18.53 21.22 8.23 7.42 4.87 8.55 7.60 6.65 5.02 6.00 4.77
x2 22.12 25.83 20.8 21.4 23.12 26.83 29.46 36.92 55.42 32.56 52.28 53.82 62.45 49.74 23.36 23.99 27.05 30.46 30.64 22.29 25.83 38.24 32.86 33.25 22.65 53.69 51.72 49 47.48 23.3 18.71 22.84 23.13 38.38 21.38 22.62 22.48 27.91
t1 382 606 544 871 644 998 721 564 718 439 484 397 650 1043 2838 1079 1091 795 596 879 587 527 573 674 981 716 743 1032 508 4129 4112 7644 3998 5060 6190 5121 8551 981
t2 144.22 102.03 103.49 114.44 96.65 135.97 89.31 143.59 101.3 142.11 206.44 192.35 65.93 112.36 72.82 133.95 89.72 156.66 110.83 89.2 100.19 108.57 71.31 90.9 112.16 60.2 110.63 98.14 58.87 74.88 49.48 149.78 212.71 54.91 48.1 44.57 119.15 30.82
t3 55.82 34.61 28.02 36.9 43.05 53.06 56.25 38.36 63.92 40.98 19.07 17.14 22.9 41.98 59.81 45.18 51.42 35.78 46.05 59.34 40.51 55.49 58.84 59.27 66.54 56.69 23.98 21.13 55 52.49 38.65 37.09 57.46 39.65 55.16 65.48 67.32 22.42
t4 40.57 66.18 21.24 62.9 1.02 4.48 20.88 25.61 44.22 62.41 17.07 17.1 21.06 15.73 46.18 72.42 21.46 40.61 42.21 47.04 16.94 36.84 61.61 71.01 69.09 45.27 17.29 31.19 3.56 10.39 20.49 35.53 27.15 64.64 78.19 78.39 81.03 37.31
73
74 Keterangan: y = Indeks Pembangunan Kesehatan Masyarakat x1 = Persentase Penduduk Miskin x2 = Angka Kematian Bayi t1 = Kepadatan Penduduk t2 = Angka Kematian Ibu t3 = Persentase Rumah Tangga Ber-PHBS t4 = Persentase Rumah Sehat
75 Lampiran 2. Program Pemilihan Titik Knot Optimal dengan Satu Titik Knot Menggunakan Software R gcv1=function(para) { data=read.table("d://data.txt",header=FALSE) data=as.matrix(data) p=length(data[,1]) q=length(data[1,]) m=ncol(data)-para-1 dataA=data[,(para+2):q] F=matrix(0,nrow=p,ncol=p) diag(F)=1 nk= length(seq(min(data[,2]),max(data[,2]),length.out=50)) knot1=matrix(ncol=m,nrow=nk) for (i in (1:m)) { for (j in (1:nk)) { a=seq(min(dataA[,i]),max(dataA[,i]),length.out=50) knot1[j,i]=a[j] } } a1=length(knot1[,1]) knot1=knot1[2:(a1-1),] aa=rep(1,p) data1=matrix(ncol=m,nrow=p) data2=data[,2:q] a2=nrow(knot1) GCV=rep(NA,a2) Rsq=rep(NA,a2) for (i in 1:a2) { for (j in 1:m) { for (k in 1:p) { if(data[k,(j+para+1)]
76 } } mx=cbind(aa,data2,data1) mx=as.matrix(mx) C=pinv(t(mx)%*%mx) B=C%*%(t(mx)%*%data[,1]) yhat=mx%*%B SSE=0 SSR=0 for (r in (1:p)) { sum=(data[r,1]-yhat[r,])^2 sum1=(yhat[r,]-mean(data[,1]))^2 SSE=SSE+sum SSR=SSR+sum1 } Rsq[i]=(SSR/(SSE+SSR))*100 MSE=SSE/p A=mx%*%C%*%t(mx) A1=(F-A) A2=(sum(diag(A1))/p)^2 GCV[i]=MSE/A2 } GCV=as.matrix(GCV) Rsq=as.matrix(Rsq) cat("======================================","\n") cat("Nilai Knot dengan Spline linear 1 knot","\n") cat("======================================","\n") print (knot1) cat("=======================================","\n") cat("Rsq dengan Spline linear 1 knot","\n") cat("=======================================","\n") print (Rsq) cat("=======================================","\n") cat("HASIL GCV dengan Spline linear 1 knot","\n") cat("=======================================","\n") print (GCV) s1=min(GCV) print(max(Rsq))
77 cat("======================================","\n") cat("HASIL GCV terkecil dengan Spline linear 1 knot","\n") cat("======================================","\n") cat(" GCV =",s1,"\n") write.csv(GCV,file="d:/coba/output GCV1.csv") write.csv(Rsq,file="d:/coba/output Rsq1.csv") write.csv(knot1,file="d:/coba/output knot1.csv") }
78 Lampiran 3. Program Pemilihan Titik Knot Optimal dengan Dua Titik Knot Menggunakan Software R gcv2=function(para) { data=read.table("d://data.txt",header=FALSE) data=as.matrix(data) p=length(data[,1]) q=length(data[1,]) m=ncol(data)-para-1 F=matrix(0,nrow=p,ncol=p) dataA=data[,(para+2):q] diag(F)=1 nk= length(seq(min(data[,2]),max(data[,2]),length.out=50)) knot=matrix(ncol=m,nrow=nk) for (i in (1:m)) { for (j in (1:nk)) { a=seq(min(dataA[,i]),max(dataA[,i]),length.out=50) knot[j,i]=a[j] } } z=(nk*(nk-1)/2) knot2=cbind(rep(NA,(z+1))) for (i in (1:m)) { knot1=rbind(rep(NA,2)) for ( j in 1:(nk-1)) { for (k in (j+1):nk) { xx=cbind(knot[j,i],knot[k,i]) knot1=rbind(knot1,xx) } } knot2=cbind(knot2,knot1) } knot2=knot2[2:(z+1),2:(2*m+1)]
79 aa=rep(1,p) data2=matrix(ncol=(2*m),nrow=p) data1=data[,2:q] a1=length(knot2[,1]) GCV=rep(NA,a1) Rsq=rep(NA,a1) for (i in 1:a1) { for (j in 1:(2*m)) { if(mod(j,2)==1) b=floor(j/2)+1 else b=j/2 for (k in 1:p) { if (data1[k,b]
80 Rsq=as.matrix(Rsq) cat("======================================","\n") cat("Nilai Knot dengan Spline linear 2 knot","\n") cat("======================================","\n") print (knot2) cat("======================================","\n") cat("Rsq dengan Spline linear 2 knot","\n") cat("======================================","\n") print (Rsq) cat("======================================","\n") cat("HASIL GCV dengan Spline linear 2 knot","\n") cat("======================================","\n") print (GCV) s1=min(GCV) cat("======================================","\n") cat("HASIL GCV terkecil dengan Spline linear 2 knot","\n") cat("======================================","\n") cat(" GCV =",s1,"\n") write.csv(GCV,file="d:/coba/output GCV2.csv") write.csv(Rsq,file="d:/coba/output Rsq2.csv") write.csv(knot2,file="d:/coba/output knot2.csv") }
81 Lampiran 4. Program Pemilihan Titik Knot Optimal dengan Tiga Titik Knot Menggunakan Software R gcv3=function(para) { data=read.table("d://data.txt",header=FALSE) data=as.matrix(data) p=length(data[,1]) q=length(data[1,]) m=ncol(data)-para-1 F=matrix(0,nrow=p,ncol=p) dataA=data[,(para+2):q] diag(F)=1 nk= length(seq(min(data[,2]),max(data[,2]),length.out=50)) knot=matrix(ncol=m,nrow=nk) for (i in (1:m)) { for (j in (1:nk)) { a=seq(min(dataA[,i]),max(dataA[,i]),length.out=50) knot[j,i]=a[j] } } knot=knot[2:(nk-1),] a2=nrow(knot) z=(a2*(a2-1)*(a2-2)/6) knot1=cbind(rep(NA,(z+1))) for (i in (1:m)) { knot2=rbind(rep(NA,3)) for ( j in 1:(a2-2)) { for (k in (j+1):(a2-1)) { for (g in (k+1):a2) { xx=cbind(knot[j,i],knot[k,i],knot[g,i]) knot2=rbind(knot2,xx) }
82 } } knot1=cbind(knot1,knot2) } knot1=knot1[2:(z+1),2:(3*m+1)] aa=rep(1,p) data1=matrix(ncol=(3*m),nrow=p) data2=data[,(para+2):q] a1=length(knot1[,1]) GCV=rep(NA,a1) Rsq=rep(NA,a1) for (i in 1:a1) { for (j in 1:ncol(knot1)) { b=ceiling(j/3) for (k in 1:p) { if (data2[k,b]
83 A1=(F-A) A2=(sum(diag(A1))/p)^2 GCV[i]=MSE/A2 } GCV=as.matrix(GCV) Rsq=as.matrix(Rsq) cat("======================================","\n") cat("Nilai Knot dengan Spline linear 3 knot","\n") cat("======================================","\n") print (knot1) cat("======================================","\n") cat("Rsq dengan Spline linear 3 knot","\n") cat("======================================","\n") print (Rsq) r=max(Rsq) print (r) cat("======================================","\n") cat("HASIL GCV dengan Spline linear 3 knot","\n") cat("======================================","\n") print (GCV) s1=min(GCV) cat("======================================","\n") cat("HASIL GCV terkecil dengan Spline linear 3 knot","\n") cat("======================================","\n") cat(" GCV =",s1,"\n") write.csv(GCV,file="d:/coba/output GCV3.csv") write.csv(Rsq,file="d:/coba/output Rsq3.csv") write.csv(knot1,file="d:/coba/output knot3.csv") }
84 Lampiran 5. Program Pemilihan Titik Knot Optimal dengan Kombinasi Titik Knot Pada Pola Berdasarkan Scatterplot Menggunakan Software R gcvkom=function(para) { data=read.table("d://data.txt",header=FALSE) data=as.matrix(data) p1=length(data[,1]) q1=length(data[1,]) v=para+2 F=matrix(0,nrow=p1,ncol=p1) diag(F)=1 x1=read.table("d:/x1.txt") x2=read.table("d:/x2.txt") x3=read.table("d:/x3.txt") x4=read.table("d:/x4.txt") n2=nrow(x1) a=matrix(nrow=4,ncol=3^4) m=0 for (i in 1:3) for (j in 1:3) for (k in 1:3) for (l in 1:3) { m=m+1 a[,m]=c(i,j,k,l) } a=t(a) GCV=matrix(nrow=nrow(x1),ncol=3^4) for (i in 1:3^4) { for (h in 1:nrow(x1)) { if (a[i,1]==1) { gab=as.matrix(x1[,1]) gen=as.matrix(data[,v]) aa=matrix(nrow=nrow(x1)*nrow(data),ncol=1)
85 for (j in 1:1) for (w in 1:nrow(data)) { if (gen[w,j]
86 } else if (a[i,2]==2) { gab=as.matrix(x2[,2:3] ) gen=as.matrix(cbind(data[,(v+1)],data[,(v+1)])) bb=matrix(nrow=nrow(x1)*nrow(data),ncol=2) for (j in 1:2) for (w in 1:nrow(data)) { if (gen[w,j]
87 gen=as.matrix(cbind(data[,(v+2)],data[,(v+2)])) cc=matrix(nrow=nrow(x1)*nrow(data),ncol=2) for (j in 1:2) for (w in 1:nrow(data)) { if (gen[w,j]
88 if (gen[w,j]
89 if (a[i,1]==2) sp=x1[,2:3] else sp=x1[,4:6] if (a[i,2]==1) spl=x2[,1] else if (a[i,2]==2) spl=x2[,2:3] else spl=x2[,4:6] if (a[i,3]==1) splin=x3[,1] else if (a[i,3]==2) splin=x3[,2:3] else splin=x3[,4:6] if (a[i,4]==1) spline=x4[,1] else if (a[i,4]==2) spline=x4[,2:3] else spline=x4[,4:6] kkk=cbind(sp,spl,splin,spline) cat("============================================== ==","\n") print(i) print(kkk) print(Rsq) } write.csv(GCV,file="d:/coba/output GCV kombinasi.csv") }
90 Lampiran 6. Program Pemilihan Titik Knot Optimal dengan Kombinasi Titik Knot Pada Pola Berdasarkan Uji RESET Menggunakan Software R gcvkom=function(para) { data=read.table("d://data.txt",header=FALSE) data=as.matrix(data) p1=length(data[,1]) q1=length(data[1,]) v=para+2 F=matrix(0,nrow=p1,ncol=p1) diag(F)=1 x1=read.table("d:/x1.txt") x2=read.table("d:/x2.txt") n2=nrow(x1) a=matrix(nrow=2,ncol=3^2) m=0 for (i in 1:3) for (j in 1:3) { m=m+1 a[,m]=c(i,j) } a=t(a) GCV=matrix(nrow=nrow(x1),ncol=3^2) for (i in 1:3^2) { for (h in 1:nrow(x1)) { if (a[i,1]==1) { gab=as.matrix(x1[,1]) gen=as.matrix(data[,v]) aa=matrix(nrow=nrow(x1)*nrow(data),ncol=1) for (j in 1:1)
91 for (w in 1:nrow(data)) { if (gen[w,j]
92 { if (gen[w,j]
93 SSR=0 for (r in 1:nrow(data)) { sum=(data[r,1]-yhat[r,])^2 sum1=(yhat[r,]-mean(data[,1]))^2 SSE=SSE+sum SSR=SSR+sum1 } Rsq=(SSR/(SSE+SSR))*100 MSE=SSE/p1 A=mx%*%C%*%t(mx) A1=(F-A) A2=(sum(diag(A1))/p1)^2 GCV[h,i]=MSE/A2 } if (a[i,1]==1) sp=x1[,1] else if (a[i,1]==2) sp=x1[,2:3] else sp=x1[,4:6] if (a[i,2]==1) spl=x2[,1] else if (a[i,2]==2) spl=x2[,2:3] else spl=x2[,4:6] kkk=cbind(sp,spl) cat("========================================= =======","\n") print(i) print(kkk) print(Rsq) } write.csv(GCV,file="d:/coba/output GCV kombinasi.csv") }
94 Lampiran 7. Program Pengujian Parameter Menggunakan Software R uji=function(alpha,para) { data=read.table("d://data.txt") knot=read.table("d://knot.txt") data=as.matrix(data) knot=as.matrix(knot) ybar=mean(data[,1]) m=para+2 p=nrow(data) q=ncol(data) dataA=cbind(data[,m],data[,m],data[,m],data[,m+1],data[,m+1] ,data[,m+1],data[,m+2],data[,m+2],data[,m+2],data[,m+3],data [,m+3],data[,m+3]) dataA=as.matrix(dataA) satu=rep(1,p) n1=ncol(knot) data.knot=matrix(ncol=n1,nrow=p) for (i in 1:n1) { for(j in 1:p) { if(dataA[j,i]
95 cat("=======================================","\ n") print (B) n1=nrow(B) yhat=mx%*%B res=data[,1]-yhat SSE=sum((data[,1]-yhat)^2) SSR=sum((yhat-ybar)^2) SST=SSR+SSE MSE=SSE/(p-n1) MSR=SSR/(n1-1) Rsq=(SSR/(SSR+SSE))*100 #uji F (uji serentak) Fhit=MSR/MSE pvalue=pf(Fhit,(n1-1),(p-n1),lower.tail=FALSE) if (pvalue<=alpha) { cat("------------------------------------","\n") cat("Kesimpulan hasil uji serentak","\n") cat("------------------------------------","\n") cat("Tolak Ho yakni minimal terdapat 1 prediktor yang signifikan","\n") cat("","\n") } else { cat("------------------------------------","\n") cat("Kesimpulan hasil uji serentak","\n") cat("------------------------------------","\n") cat("Gagal Tolak Ho yakni semua prediktor tidak berpengaruh signifikan","\n") cat("","\n") }
96 #uji t (uji individu) thit=rep(NA,n1) pval=rep(NA,n1) SE=sqrt(diag(MSE*(pinv(t(mx)%*%mx)))) cat("------------------------------------","\n") cat("Kesimpulan hasil uji individu","\n") cat("------------------------------------","\n") thit=rep(NA,n1) pval=rep(NA,n1) for (i in 1:n1) { thit[i]=B[i,1]/SE[i] pval[i]=2*(pt(abs(thit[i]),(p-n1),lower.tail=FALSE)) if (pval[i]<=alpha) cat("Tolak Ho yakni prediktor signifikan dengan pvalue",pval[i],"\n") else cat("Gagal tolak Ho yakni prediktor tidak signifikan dengan pvalue",pval[i],"\n") } thit=as.matrix(thit) cat("=======================================","\ n") cat("nilai t hitung","\n") cat("=======================================","\ n") print (thit) cat("------------------------------------","\n") cat("Analysis of Variance","\n") cat("======================================","\n ") cat("Sumber df SS MS Fhit","\n") cat("Regresi ",(n1-1)," ",SSR," ",MSR,"",Fhit,"\n") cat("Error ",p-n1," ",SSE,"",MSE,"\n") cat("Total ",p-1," ",SST,"\n") cat("======================================","\n ") cat("s=",sqrt(MSE)," Rsq=",Rsq,"\n")
97 cat("pvalue(F)=",pvalue,"\n") write.csv(res,file="d:/coba/output uji residual.csv") write.csv(pval,file="d:/coba/output uji pvalue.csv") write.csv(mx,file="d:/coba/output uji mx.csv") write.csv(yhat,file="d:/coba/output uji yhat.csv") }
98 Lampiran 8. Program Uji Glejser dengan Tiga Titik Knot Menggunakan Software R glejser=function(alpha,para) { data=read.table("d:/data.txt") knot=read.table("d:/knot.txt") res=read.table("d:/residual.txt") data=as.matrix(data) knot=as.matrix(knot) res=abs(res) res=as.matrix(res) rbar=mean(res) m=para+2 p=nrow(data) q=ncol(data) dataA=cbind(data[,m],data[,m],data[,m],data[,m+1],data[,m+1],data[, m+1],data[,m+2],data[,m+2],data[,m+2],data[,m+3],data[,m+3],data[, m+3]) dataA=as.matrix(dataA) satu=rep(1,p) n1=ncol(knot) data.knot=matrix(ncol=n1,nrow=p) for (i in 1:n1) { for(j in 1:p) { if (dataA[j,i]
99 SSE=sum((res-yhat)^2) SSR=sum((yhat-rbar)^2) SST=SSR+SSE MSE=SSE/(p-n1) MSR=SSR/(n1-1) Rsq=(SSR/SST)*100 #uji F (uji serentak) Fhit=MSR/MSE pvalue=pf(Fhit,(n1-1),(p-n1),lower.tail=FALSE) if (pvalue<=alpha) { cat("------------------------------------","\n") cat("Kesimpulan hasil uji serentak","\n") cat("------------------------------------","\n") cat("Tolak Ho yakni minimal terdapat 1 prediktor yang signifikan atau terjadi heteroskedastisitas","\n") cat("","\n") } else { cat("------------------------------------","\n") cat("Kesimpulan hasil uji serentak","\n") cat("------------------------------------","\n") cat("Gagal Tolak Ho yakni semua prediktor tidak berpengaruh signifikan atau tidak terjadi heteroskedastisitas","\n") cat("","\n") } cat("Analysis of Variance","\n") cat("======================================","\n") cat("Sumber df SS MS Fhit","\n") cat("Regresi ",(n1-1)," ",SSR," ",MSR,"",Fhit,"\n") cat("Error ",p-n1," ",SSE,"",MSE,"\n") cat("Total ",p-1," ",SST,"\n") cat("======================================","\n") cat("s=",sqrt(MSE)," Rsq=",Rsq,"\n") cat("pvalue(F)=",pvalue,"\n") }
100 Lampiran 9. Output Nilai GCV dengan Satu Titik Knot Pada Pola Berdasarkan Scatterplot No
GCV
t1
t2
t3
t4
1
0.0010628
548.7143
34.53204
18.16408
2.652857
2
0.0011558
715.4286
38.24408
19.18816
4.285714
3
0.0013117
882.1429
41.95612
20.21224
5.918571
4
0.0009344
1048.857
45.66816
21.23633
7.551429
5
0.0009242
1215.571
49.3802
22.26041
9.184286
6
0.0009001
1382.286
53.09224
23.28449
10.81714
7
0.0008940
1549
56.80429
24.30857
12.45
8
0.0009008
1715.714
60.51633
25.33265
14.08286
9
0.0009164
1882.429
64.22837
26.35673
15.71571
10
0.0009310
2049.143
67.94041
27.38082
17.34857
⋮ 40
⋮ 0.0009678
⋮ 7050.571
⋮ 179.3016
⋮ 58.10327
⋮ 66.33429
41
0.0009736
7217.286
183.0137
59.12735
67.96714
42
0.0009781
7384
186.7257
60.15143
69.6
43
0.0009850
7550.714
190.4378
61.17551
71.23286
44
0.0009957
7717.429
194.1498
62.19959
72.86571
45
0.0009994
7884.143
197.8618
63.22367
74.49857
46
0.0010091
8050.857
201.5739
64.24776
76.13143
47
0.0010351
8217.571
205.2859
65.27184
77.76429
48
0.0009683
8384.286
208.998
66.29592
79.39714
101 Lampiran 10. Output Nilai GCV dengan Dua Titik Knot Pada Pola Berdasarkan Scatterplot No
GCV
1
0.0008127
2
0.0008372
3
0.0007820
⋮
⋮
596
0.0008271
597
0.0008271
598
0.0008271
599
0.0008271
600
0.0008791
⋮
⋮
1224
0.0009317
1225
0. 0009317
t1
t2
t3
382
30.82
17.14
1.02
548.71
34.532
18.164
2.6529
382
30.82
17.14
1.02
715.43
38.244
19.188
4.2857
382
30.82
17.14
1.02
882.14
41.956
20.212
5.9186
⋮
⋮
⋮
t4
⋮
2716
82.789
31.477
23.88
2882.7
86.501
32.501
25.513
2716
82.789
31.477
23.88
3049.4
90.213
33.525
27.146
2716
82.789
31.477
23.88
3216.1
93.925
34.549
28.779
2716
82.789
31.477
23.88
3382.9
97.637
35.573
30.411
2716
82.789
31.477
23.88
3549.6
101.35
36.598
32.044
⋮
⋮
⋮
⋮
8217.6
205.29
65.272
77.764
8551
212.71
67.32
81.03
8384.3
209
66.296
79.397
8551
212.71
67.32
81.03
102 Lampiran 11. Output Nilai GCV dengan Tiga Titik Knot Pada Pola Berdasarkan Scatterplot No 1
2
3 ⋮ 8929
8930
GCV 0.0010694
0.0011480
0.0010397 ⋮ 0.0015183
0.0016064
⋮
⋮
17295
0.0010942
17296
0.0010771
t1
t2
t3
t4
548.71
34.532
18.164
2.6529
715.43
38.244
19.188
4.2857
882.14
41.956
20.212
5.9186
548.71
34.532
18.164
2.6529
715.43
38.244
19.188
4.2857
1048.9
45.668
21.236
7.5514
548.71
34.532
18.164
2.6529
715.43
38.244
19.188
4.2857
1215.6
49.38
22.26
9.1843
⋮
⋮
2215.9
⋮
71.6524
⋮
28.405
18.981
2549.3
79.077
30.453
22.247
8050.9
201.57
64.248
76.131
2215.9
71.6524
28.405
18.981
2549.3
79.077
30.453
22.247
8217.6
205.29
65.272
77.767
⋮
⋮
⋮
⋮
7884.1
197.8618
63.224
74.499
8217.6
205.29
65.272
77.764
8384.3
209
66.296
79.397
8050.9
201.5738
64.248
76.131
8217.6
205.29
65.272
77.764
9394.3
209
66.296
79.397
103 Lampiran 12. Output Nilai GCV dengan Kombinasi Titik Knot Pola Berdasarkan Scatterplot
1
Kombinasi Knot 1111
0.0008940
2
1112
0.0008698
No
GCV
t1
t2
t3
t4
1549
56.804
24.309
12.45
1549
56.804
24.309
1.02 7.5514 27.146
3
1113
0.0006033
1549
56.804
24.309
30.411 32.044
⋮
⋮
⋮
40
2221
0.0009594
⋮ 382
⋮
80
81
2223
⋮
3332
3333
⋮
30.82
17.14
45.668
21.236
30.82
17.14
1048.9
45.668
21.236
⋮
⋮
⋮
3049.4
90.213
33.525
3382.9
97.637
35.473
3549.6 3049.4 3382.9 3549.6
101.35 90.213 97.637 101.35
36.598 33.525 35.473 36.598
1048.9 382
41
⋮
0.0006062
⋮
0.0007063
0.0005271
⋮ 12.45 27.146 30.411 32.044 ⋮ 1.02 7.5514 27.146 30.411 32.044
104 Lampiran 13. Output Nilai GCV dengan Satu Titik Knot Pada Pola Berdasarkan Uji RESET No
GCV
t1
t2
1
0.001011508
548.7143
18.16408
2
0.001190013
715.4286
19.18816
3
0.001364926
882.1429
20.21224
4
0.000900162
1048.857
21.23633
5
0.000889297
1215.571
22.26041
6
0.000862348
1382.286
23.28449
7
0.000846766
1549
24.30857
8
0.000842988
1715.714
25.33265
9
0.000842146
1882.429
26.35673
10
0.000843269
2049.143
27.38082
⋮
⋮
⋮
⋮
40
0.000886388
7050.571
58.10327
41
0.000887215
7217.286
59.12735
42
0.000887743
7384
60.15143
43
0.00089082
7550.714
61.17551
44
0.000894301
7717.429
62.19959
45
0.000894312
46
0.000895608
7884.143 8050.857
63.22367 201.5739
47
0.000901015
8217.571
205.2859
48
0.000902826
8384.286
208.998
105 Lampiran 14. Output Nilai GCV dengan Dua Titik Knot Pada Pola Berdasarkan Uji RESET No
GCV
1
0.000827087
2
0.00088176
3
0.00088176
⋮
⋮
596
0.000807918
597
0.000823275
598
0.000823647
599
0.000830385
600
0.000836337
⋮
⋮
1224
0.000827087
1225
0.000827087
t1
t2 382
17.14
548.71
18.164
382
17.14
715.43
19.188
382
17.14
882.14
20.212
⋮
⋮
2716
31.477
2882.7
32.501
2716
31.477
3049.4
33.525
2716
31.477
3216.1
34.549
2716
31.477
3382.9
35.573
2716
31.477
3549.6
36.598
⋮
⋮
8217.6
65.272
8551
67.32
8384.3
66.296
8551
67.32
106 Lampiran 15. Output Nilai GCV dengan Tiga Titik Knot Pada Pola Berdasarkan Uji RESET No 1
2
3 ⋮ 8929
8930
GCV 0.000966366
0.000954141
0.000845779 ⋮ 0.001105164
0.001121472
⋮
⋮
17295
0.000953282
17296
0.000951113
t1
t2
548.71
18.164
715.43
19.188
882.14
20.212
548.71
18.164
715.43
19.188
1048.9
21.236
548.71
18.164
715.43
19.188
1215.6
22.26
⋮
⋮
2215.9
28.40489
2549.3
30.453
8050.9
64.248
2215.9
28.40489
2549.3
30.453
8217.6
65.272
⋮
⋮
7884.1
63.2236
8217.6
65.272
8384.3
66.296
8050.9
64.2477
8217.6
65.272
8384.3
66.296
107 Lampiran 16. Output Nilai GCV dengan Kombinasi Titik Knot Pada Pola Berdasarkan Uji RESET No
Kombinasi Knot
GCV
t1
t2
1
11
0.000837175
3382.9
35.573
2
12
0.000889416
3382.9
30.453 46.838 18.164
3
13
0.000527796
3382.9
20.212 22.26
4
21
0.00089545
5
22
0.000951058
2540.3 5216.7 2540.3
30.453
5216.7
46.838
2540.3 6
23
35.573
0.00055619
18.164 20.212
5216.7
22.26
548.71 7
31
0.000989711
882.14
35.573
1215.6 548.71 8
32
0.001060501
882.14 1215.6
9
33
0.000563922
30.453 46.838
548.71
18.164
882.14
20.212
1215.6
22.26
108 Lampiran 17. Output Nilai GCV dengan Satu Titik Knot Pada Model Empat Variabel No
GCV
t1
t3
t4
1
0.001018423
548.7143
18.16408
2.652857
2
0.001130984
715.4286
19.18816
4.285714
3
0.001284738
882.1429
20.21224
5.918571
4
0.000874602
1048.857
21.23633
7.551429
5
0.000848849
1215.571
22.26041
9.184286
6
0.000827162
1382.286
23.28449
10.81714
7
0.000824607
1549
24.30857
12.45
8
0.000824776
1715.714
25.33265
14.08286
9
0.000826942
1882.429
26.35673
15.71571
10
0.000832915
2049.143
27.38082
17.34857
⋮
⋮
⋮
⋮
⋮
40
0.000878616
7050.571
58.10327
66.33429
41
0.00088872
7217.286
59.12735
67.96714
42
0.00090022
7384
60.15143
69.6
43
0.000915251
7550.714
61.17551
71.23286
44
0.000931081
7717.429
62.19959
72.86571
45
0.000935496
46
0.000942301
7884.143 8050.857
63.22367 64.24776
74.49857 76.13143
47
0.000948987
8217.571
65.27184
77.76429
48
0.000888626
8384.286
66.29592
79.39714
109 Lampiran 18. Output Nilai GCV dengan Dua Titik Knot Pada Model Empat Variabel No
GCV
1
0.000850043
2
0.000850043
3
0.000850043
⋮
⋮
596
0.000922619
597
0.00092253
598
0.000923752
599
0.000923941
600
0.00092323
⋮
⋮
1224
0.000850043
1225
0.000850043
t1
t3
t4
382
17.14
1.02
548.71
18.164
2.6529
382
17.14
1.02
715.43
19.188
4.2857
382
17.14
1.02
882.14
20.212
5.9186
⋮
⋮
⋮
2716
31.477
23.88
2882.7
32.501
25.513
2716
31.477
23.88
3049.4
33.525
27.146
2716
31.477
23.88
3216.1
34.549
28.779
2716
31.477
23.88
3382.9
35.573
30.411
2716
31.477
23.88
3549.6
36.598
32.044
⋮
⋮
⋮
8217.6
65.272
77.764
8551
67.32
81.03
8384.3
66.296
79.397
8551
67.32
81.03
110 Lampiran 19. Output Nilai GCV dengan Tiga Titik Knot Pada Model Empat Variabel No 1
2
3 ⋮ 8929
8930
GCV 0.000990136
0.001034584
0.00083531 ⋮ 0.001173373
0.001172964
⋮
⋮
17295
0.000937542
17296
0.000941117
t1
t3
t4
548.71
18.164
2.6529
715.43
19.188
4.2857
882.14
20.212
5.9186
548.71
18.164
2.6529
715.43
19.188
4.2857
1048.9
21.236
7.5514
548.71
18.164
2.6529
715.43
19.188
4.2857
1215.6
22.26
9.1843
⋮
⋮
⋮
2215.9
28.405
18.981
2549.3
30.453
22.247
8050.9
64.248
76.131
2215.9
28.405
18.981
2549.3
30.453
22.247
8217.6
65.272
77.767
⋮
⋮
⋮
7884.1
63.224
74.499
8217.6
65.272
77.764
8384.3
66.296
79.397
8050.9
64.248
76.131
8217.6
65.272
77.764
8384.3
66.296
79.397
111 Lampiran 20. Output Nilai GCV dengan Kombinasi Titik Knot Model Empat Variabel
1
Kombinasi Knot 111
0.000824607
2
112
0.000853992
⋮
⋮
⋮
14
222
0.000963595
No
15
⋮
26
27
223
⋮
332
333
GCV
t1
t3
t4
1549
24.309
1549
24.309
⋮
⋮
12.45 20.614 22.247 ⋮
2382.571
29.429
20.614
2549.286
30.453
22.247
2382.571
29.429
0.000567568
⋮
0.0007132
0.0007063
23.88 30.411
2382.571
29.429
⋮
⋮
2716
31.477
3382.857
35.573
3716.286 2716 3382.857 2716
37.621 31.477 35.573 31.477
33.677 ⋮ 20.614 20.614 23.88 30.411 33.677
112 Lampiran 21. Output Estimasi Parameter dan Uji Signifikansi ======================================= Estimasi Parameter ======================================= [,1] [1,] 7.340836e-01 [2,] -2.456648e-03 [3,] 5.953191e-05 [4,] -9.089387e-04 [5,] 2.399663e-03 [6,] -1.561649e-03 [7,] -5.670136e-04 [8,] 8.322865e-03 [9,] -2.660555e-02 [10,] 1.873022e-02 [11,] 1.081512e-03 [12,] -2.050594e-02 [13,] 6.099663e-02 [14,] -4.178350e-02 -----------------------------------Kesimpulan hasil uji serentak -----------------------------------Tolak Ho yakni minimal terdapat 1 prediktor yang signifikan -----------------------------------Kesimpulan hasil uji individu -----------------------------------Tolak Ho yakni prediktor signifikan dengan pvalue 5.998151e-13 Tolak Ho yakni prediktor signifikan dengan pvalue 1.274644e-08 Tolak Ho yakni prediktor signifikan dengan pvalue 0.0006554105 Tolak Ho yakni prediktor signifikan dengan pvalue 0.005785281 Tolak Ho yakni prediktor signifikan dengan pvalue 0.005455059 Tolak Ho yakni prediktor signifikan dengan pvalue 0.004641598 Gagal tolak Ho yakni prediktor tidak signifikan dengan pvalue 0.7852209 Gagal tolak Ho yakni prediktor tidak signifikan dengan pvalue 0.3146465
113 Gagal tolak Ho yakni prediktor tidak signifikan dengan pvalue 0.058204 Tolak Ho yakni prediktor signifikan dengan pvalue 0.03729008 Gagal tolak Ho yakni prediktor tidak signifikan dengan pvalue 0.1140177 Tolak Ho yakni prediktor signifikan dengan pvalue 6.598624e-06 Tolak Ho yakni prediktor signifikan dengan pvalue 1.489498e-06 Tolak Ho yakni prediktor signifikan dengan pvalue 1.960266e-06 ======================================= nilai t hitung ======================================= [,1] [1,] 13.8592275 [2,] -8.4149600 [3,] 3.9135257 [4,] -3.0294833 [5,] 3.0541156 [6,] -3.1214947 [7,] -0.2755914 [8,] 1.0270307 [9,] -1.9891109 [10,] 2.2048902 [11,] 1.6401463 [12,] -5.7326276 [13,] 6.3389662 [14,] -6.2259009 -----------------------------------Analysis of Variance ====================================== Sumber df SS MS Fhit Regresi 13 0.08731821 0.006716785 24.50918 Error 24 0.006577245 0.0002740519 Total 37 0.09389545 ====================================== s= 0.01655451 Rsq= 92.99514 pvalue(F)= 1.050539e-10
114 Lampiran 22. Output Residual No
Residual
No
Residual
1
-0.00829
20
-0.00237
2
0.00299
21
0.01434
3
-0.00975
22
-0.00155
4
-0.01653
23
-0.00417
5
0.00341
24
0.01951
6
0.01030
25
-0.00614
7
-0.01060
26
-0.01202
8
0.01369
27
0.00794
9
-0.00625
28
-0.00709
10
0.00401
29
-0.02457
11
-0.03381
30
0.01393
12
0.01742
31
-0.02132
13
0.01135
32
-0.00625
14
0.01201
33
0.00354
15
0.00000
34
-0.00109
16
-0.01573
35
-0.01507
17
-0.00246
36
0.01808
18
0.00573
37
0.00819
19
0.02538
38
0.01324
115 Lampiran 23. Output Uji Gljeser -----------------------------------Kesimpulan hasil uji serentak -----------------------------------Gagal Tolak Ho yakni semua prediktor tidak berpengaruh signifikan atau tidak terjadi heteroskedastisitas Analysis of Variance ====================================== Sumber df SS MS Fhit Regresi 13 0.0007226654 5.558965e-05 0.9338176 Error 24 0.001428707 5.952945e-05 Total 37 0.002151372 ====================================== s= 0.007715533 Rsq= 33.59091 pvalue(F)= 0.5356631
116
(Halaman ini sengaja dikosongkan)
117 Lampiran 24. Surat Pernyataan Data Sekunder SURAT PERNYATAAN Saya yang bertanda tangan di bawah ini, mahasiswa Jurusan Statistika FMIPA ITS: Nama NRP
: Made Ayu Dwi Octavanny : 1313 100 012
Menyatakan bahwa data yang digunakan dalam Tugas Akhir/Thesis ini merupakan data sekunder yang diambil dari penelitian/buku/Tugas Akhir/Thesis/publikasi lainnya, yaitu: Sumber
: Website Badan Pusat Statistik Provinsi Jawa Timur Profil Kesehatan Provinsi Jawa Timur Indeks Pembangunan Kesehatan Masyarakat
Keterangan: 1. Tabel Persentase Penduduk Miskin Menurut Kabupaten/Kota Provinsi Jawa Timur Tahun 2013 2. Angka Kematian Bayi per 1000 Kelahiran Hidup Provinsi Jawa Timur Tahun 2013 3. Tabel Kepadatan Penduduk Pertengahan Tahun Menurut Kabupaten/Kota 2013 4. Angka Kematian Ibu per 100.000 Kelahiran Hidup Provinsi Jawa Timur Tahun 2013 5. Tabel Persentase Rumah Tangga Berperilaku Hidup Bersih dan Sehat Menurut Kabupaten/ Kota Provinsi Jawa Timur Tahun 2013 6. Tabel Persentase Rumah Sehat Menurut Kabupaten/ Kota Provinsi Jawa Timur Tahun 2013 7. Tabel IPKM Provinsi Jawa Timur 2013 Surat pernyataan ini dibuat dengan sebenarnya. Apabila terdapat pemalsuan data maka saya siap menerima sanksi sesuai aturan yang berlaku. Mengetahui Pembimbing Tugas Akhir
(Prof. Dr. I Nyoman Budiantara, M.Si) NIP. 19650603 198903 1 003
Surabaya, 23 Januari 2017
(Made Ayu Dwi Octavanny) NRP 1313 100 012
118
(Halaman ini sengaja dikosongkan)
BIODATA PENULIS Made Ayu Dwi Octavanny dilahirkan di Denpasar, 20 Oktober 1995. Anak bungsu dari I Wayan Suja dan Ni Made Adi Sri Masyuni ini memiliki hobi jalan-jalan dan makan. Penulis telah menempuh pendidikan formal di SD 5 Saraswati Denpasar, SMPN 1 Denpasar, dan SMA Negeri 1 Denpasar. Selanjutnya, penulis melanjutkan ke jenjang perguruan tinggi yaitu di Jurusan Statistika ITS pada tahun 2013 melalui jalur SNMPTN undangan. Semasa perkuliahan, penulis mengikuti organisasi di luar perkuliahan, yaitu Sekretaris Departemen Pengembangan Sumber Daya Mahasiswa TPKH-ITS periode kepengurusan 2014/2015, Bendahara Umum TPKH-ITS periode kepengurusan 2015/2016, dan anggota Divisi Pers HIMASTA-ITS selama dua periode kepengurusan, yaitu 2014-2016. Selain itu, penulis aktif dalam beberapa kepanitian di dalam kampus, baik menjadi panitia inti maupun anggota sie. Penulis pernah mendapatkan beasiswa dari PPA (Peningkatan Prestasi Akademik) pada tahun ke-3 di bangku perkuliahan. Apabila pembaca ingin berdiskusi maupun memberikan kritik dan saran mengenai tugas akhir ini, penulis dapat dihubungi melalui: Email :
[email protected] Telepon: +6281 999 661 663
119
120
(Halaman ini sengaja dikosongkan)
