## PEMODELAN KEMATIAN BALITA MALNUTRISI DENGAN PENDEKATAN ZERO-INFLATED POISSON (ZIP) REGRESSION DI PROVINSI JAWA TENGAH
1 Prisca Shery Camelia, 2 Indah Manfaati Nur, 3 Moh. Yamin Darsyah
1,2,3 Program Studi S1 Statistika Fakultas Matematika dan Ilmu Pengetahuan Alam Universitas Muhammadiyah Semarang
Email : priscasherycamelia@gmail.com

## ABSTRAK
Kasus balita malnutrisi di Jawa Tengah merupakan salah satu kasus yang harus diperhatikan karena dapat menyebabkan kematian pada balita jika tidak ditangani dengan serius. Jumlah kematian balita malnutrisi untuk tiap kabupaten/kota di Provinsi Jawa Tengah sebagai variabel dependen mengikuti distribusi Poisson, analisis yang dapat digunakan adalah Metode Regresi Poisson. Metode regresi Poisson memiliki asumsi bahwa data harus equidispersi, namun setelah dilakukan pengujian overdispersi terbukti bahwa terdapat kondisi overdispersi pada model regresi Poisson, dimana nilai varians lebih besar dari pada mean. Metode yang dapat digunakan pada data yang mengalami overdispersi adalah Regresi Zero-Inflated Poisson (ZIP). Pemodelan dengan regresi ZIP menghasilkan model logit yang menunjukkan bahwa jumlah kematian balita karena malnutrisi di tiap kabupaten di Provinsi Jawa Tengah yang bernilai nol dipengaruhi oleh jumlah balita yang mengalami malnutrisi, jumlah pusat kesehatan masayarakat dan rata-rata kelengkapan imunisasi. Model log Regresi ZIP menjelaskan bahwa jumlah kematian balita akan berkurang jika jumlah pelayanan kesehatan masyarakat (puskesmas) tiap kab/kota di Jawa tengah bertambah jumlahnya.
Kata Kunci : Malnutrisi , Overdispersi, Regresi Poisson, Regresi Zero-Inflated Poisson

## PENDAHULUAN
Angka Kematian Balita (AKABA) di Indonesia pada tahun 2011 lebih tinggi dari Angka Kematian Bayi (AKB) per 1000 kelahiran hidup. Kematian balita di Indonesia memberikan efek yang sangat signifikan terhadap perubahan akan kelangsungan hidup pada jaman selanjutnya, karena balita adalah generasi penerus bangsa. UNICEF menyebutkan, setiap tiga menit satu balita meninggal di Indonesia atau setara dengan 150.000 anak dalam setiap tahunnya. Angka Kematian Balita merupakan jumlah kematian anak usia 0-5 tahun selama satu tahun tertentu per 1000 anak umur yang sama pada pertengahan tahun itu
(termasuk kematian bayi). AKABA di Indonesia masih menjadi yang tertinggi di Asia yaitu 31 per 100.000 kelahiran. Salah satu penyebab adanya kejadian kematian pada balita adalah malnutrisi atau lebih dikenal dengan sebutan gizi buruk. Malnutrisi sering terjadi pada balita karena usia mereka yang masih rentan terhadap penyakit dan fungsi alat pencernaan yang belum sempurna, ini mengakibatkan banyak kasus malnutrisi yang terjadi pada balita serta akibat terburuknya adalah kematian pada balita. Menurut data terdapat 19,6% balita kekurangan gizi yang terdiri dari 5,7% balita dengan gizi buruk(malnutrisi) dan 13,9% berstatus gizi kurang serta 4,5% balita dengan gizi lebih. Berdasarkan
data dari [3] terdapat 2.978 kasus balita malnutrisi atau gizi buruk dengan 71 diantaranya meninggal dunia. Angka yang cukup tinggi mengingat pada tahun 2014 tercatat hanya 15 balita malnutrisi yang meninggal dengan 347 jumlah kejadian malnutrisi.
Jawa Tengah merupakan penyumbang kasus malnutrisi atau gizi buruk terbanyak ke-dua setelah Jawa Timur di Pulau Jawa. Kejadian balita
malnutrisi merupakan isu kesehatan yang harus selalu diwaspadai dan ditangani secara maksimal karena faktor terburuk dari kejadian balita malnutrisi dapat menyebabkan kematian pada balita jika tidak ditangani dengan benar.
Model Poisson banyak digunakan dalam berbagai bidang termasuk kesehatan masyarakat, epidemiologi,
sosiologi, psikologi, teknik, pertanian dan lainnya [1]. Metode regresi Poisson mensyaratkan adanya equidispersi yaitu suatu kondisi dimana nilai mean dan varian dari variabel respon bernilai sama, namun terkadang terjadi
kasus
overdispersi dan underdispersi .
Overdispersi yaitu jika terjadi kasus nilai varians lebih besar dari mean nya, sedangkan underdispers i terjadi jika nilai varians lebih kecil dari mean nya. Pada Regresi Poisson banyak data yang bernilai nol pada hasil observasi, untuk mengatasi hal tersebut maka banyak metode yang telah dikembangkan yaitu dengan metode Zero Inflated Poisson
Regression (ZIP) [4]. Regresi ZIP mampu mengendalikan overdispersi
dalam distribusi Poisson dan inflasi nilai 0 sehingga akurasi estimasi parameter dapat terjamin. Secara umum model regresi ZIP masih jarang digunakan untuk data count yang menunjukkan adanya inflasi akibat nilai 0 dan overdispersi. Model regresi ZIP lebih tepat diaplikasikan dari pada regresi Poisson untuk data yang mengandung lebih banyak kejadian 0. ZIP adalah model campuran untuk data diskrit
dengan banyak peristiwa yang bernilai 0.
Penulis belum menemukan penelitian maupun publikasi yang meneliti tentang kematian balita malnutrisi dengan menggunakan metode Regresi ZIP, kasus kematian balita malnutrisi merupakan masalah yang kompleks dan penyebab kejadian malnutrisi pada balita sangat bervariasi sehingga penulis ingin memodelkan Kematian Balita Malnutrisi di Provinsi Jawa Tengah dengan Pendekatan ZIP dan memperoleh faktor- faktor yang secara
signifikan mempengaruhi kasus kematian balita malnutrisi di Provinsi Jawa Tengah.

## METODE PENELITIAN

## Sumber Data dan Variabel Penelitian
Data yang digunakan dalam penelitian ini adalah data sekunder yang diperoleh dari Publikasi Provinsi Jawa Tengah Dalam Angka 2014 yang bersumber dari Badan Pusat Statistik Provinsi Jawa Tengah. Penelitian ini terdiri dari dari variabel dependen (Y) yaitu jumlah kematian balita malnutrisi di Jawa
Tengah pada tahun 2014, dan variabel independen (X) yaitu jumlah kejadian malnutrisi/gizi buruk balita tiap kab/kota(X 1 ), jumlah Pusat kesehatan (puskesmas) tiap kab/kota(X 2 ), Jumlah dokter tiap kab/kota(X 3 ), dan Rata-rata kelengkapan imunisasi tiap kab/kota(X 4 ).

## Metode Analisis
Langkah awal dalam penelitian ini adalah dengan melakukan uji kolinieritas pada variabel independen (X) dengan melihat nilai koefisien korelasi, nilai VIF dan nilai eigenvalue. Pengujian ini dilakukan untuk mengetahui bahwa variabel independen tidak memiliki hubungan dengan variabel independen lainnya. Kemudian, dilakukan pengujian distribusi data
menggunakan uji
Kolmogorov-Smirnov , uji ini bertujuan
untuk membuktikan bahwa bentuk
distribusi pada variabel angka kematian balita malnutrisi (Y) mengikuti distribusi Poisson.
Jika data mengikuti distribusi poisson maka dillanjutkan menggunakan analisis regresi Poisson, dalam analisis regresi Poisson dilakukan penaksiran parameter model regresi Poisson dan menentukan model yang paling fit terhadap data. Langkah selanjutnya adalah menghitung nilai devians untuk mengidentifikasi overdispersi pada data, jika terjadi overdispersi maka model regresi Poisson tidak layak digunakan.

## HASIL PENELITIAN

## Uji Kolinieritas
Pengujian kolinieritas bertujuan untuk mengetahui apakah variabel-variabel independent dalam keadaan yang saling tidak berkolinieritas. Ada tida cara yang dapat dilakukan untuk menguji
kolinieritas yaitu dengan melihat nilai koefisien korelasi, nilai VIF dan nilai eigenvalue.
Tabel 1. Koefisien Korelasi Model ࢞ ૝ ࢞ ૜ ࢞ ૛ ࢞ ૚ Correlations ݔ 4 1.00 ݔ 3 -.429 1.00 ݔ 2 -.251 -.243 1.00 ݔ 1 -.428 .384 -.296 1.00
Berdasarkan Korelasi dikatakan tinggi jika nilai korelasi antar variabel bernilai 0.95 atau lebih, karena pada hasil diatas tidak terdapat korelasi yang tinggi maka dapat disimpulkan bahwa tidak terdapat kolinieritas pada variabel independent.
Tabel 2. Nilai Variance Inflation Factor (VIF)
Variabel ࢞ ૚ ࢞ ૛ ࢞ ૜ ࢞ ૝ VIF 1.562 1.525 1.469 1.854
Berdasarkan Tabel 2 dapat disimpulkan bahwa tidak ada kondisi kolinieritas antar variabel independent karena nilai VIF tidak ada yang lebih dari 10. Cara ketiga yang dapat digunakan untuk mendeteksi adanya kondisi kolinieritas adalah nilai eigenvalue pada tiap variabel. Kondisi kolinieritas antar variabel independent ditunjukkan dengan besarnya nilai eigen minimum yang lebih kecil dari 0.05.
Nilai eigenvalue yang didapatkan adalah ߣ = (0.457 0.097 0.069 0.05 ) , dari nilai eigenvalue tersebut dapat diketahui bahwa tidak terdapat nilai yang kurang dari 0.05, dengan demikian berarti tidak terjadi kolinieritas pada variabel independen. Ketiga cara yang digunakan untuk menguji kolinieritas pada variabel independen dapat disimpulakan bahwa tidak terjadi kondisi kolinieritas pada variabel independen jadi semua variabel independen dapat dimasukkan kedalam analisis regresi Poisson.

## Model Regresi Poisson
Pengujian distribusi Poisson pada variabel dependen membuktikan bahwa variabel dependen mengikuti distribusi Poisson. Nilai T unuk uji Kolmogorov- Smirnov adalah adalah 0.145, dengan nilai ݓ ଴ . ଽଽ adalah 0.269. taraf signifikansi yang digunakan dalam pengujian distribusi dan pengujian selanjutnya adalah α = 0.1, karena penelitian ini merupakan penelitian sosial. Setelah melakukan uji distribusi Poisson maka langkah selanjutnya adalah melakukan pembentukan model regresi Poissson.
Pembentukan model Regresi Poisson dapat dilihat dari nilai devians pada setiap kemungkinan model yang
terbentuk. Penelitian ini menggunakan
Tabel 3. Nilai Estimasi Parameter dan Devian pada Tiap Kemungkinan Model Regresi Poisson
No Model Nilai Dugaan Parameter ࡰࢋ࢜࢏ࢇ࢔࢙ ࢼ ૙ ࢼ ૚ ࢼ ૛ ࢼ ૜ ࢼ ૝ 1 ߤ = exp( ߚ ଴ + ݔ ଵ ) 0.6362 0.0008 66.8563 2 ߤ = exp( ߚ ଴ + ݔ ଶ ) 0.5386 0.0016 67.0714 3 ߤ = exp( ߚ ଴ + ݔ ଷ ) 0.4824 0.0050 66.5888 4 ߤ = exp( ߚ ଴ + ݔ ସ ) 0.1103 0.0000 62.0203* 5 ߤ = exp( ߚ ଴ + ݔ ଵ + ݔ ଶ ) 0.5566 0.0007 0.0009 66.7984 6 ߤ = exp( ߚ ଴ + ݔ ଵ + ݔ ଷ ) 0.3794 0.0009 0.0055 66.0094 7 ߤ = exp( ߚ ଴ + ݔ ଵ + ݔ ସ ) 0.1036 0.1036 -0.0005 61.8416* 8 ߤ = exp( ߚ ଴ + ݔ ଶ + ݔ ଷ ) 0.4125 0.0009 0.0045 66.5290 9 ߤ = exp( ߚ ଴ + ݔ ଶ + ݔ ସ ) 0.3261 -0.0030 0.0000 61.5071* 10 ߤ = exp( ߚ ଴ + ݔ ଷ + ݔ ସ ) 0.1279 -0.0006 0.0000 62.0097* 11 ߤ = exp( ߚ ଴ + ݔ ଵ + ݔ ଶ + ݔ ଷ ) 0.4084 0.0010 -0.0004 0.0057 0.3029 65.9969 12 ߤ = exp( ߚ ଴ + ݔ ଵ + ݔ ଶ + ݔ ସ ) 0.3029 -0.0003 -0.0027 0.0000 61.4489* 13 ߤ = exp( ߚ ଴ + ݔ ଵ + ݔ ଷ + ݔ ସ ) 0.1533 -0.0007 -0.0019 0.0000 61.7583 14 ߤ = exp( ߚ ଴ + ݔ ଶ + ݔ ଷ + ݔ ସ ) 0.3331 -0.0030 -0.0003 0.0000 61.5042* 15 ߤ = exp( ߚ ଴ + ݔ ଵ + ݔ ଶ + ݔ ଷ + ݔ ସ ) 0.3216 -0.0004 -0.0026 -0.0012 0.0000 61.4192*
( Keterangan : * Model dengan nilai Devian terkecil )
empat variabel independent sehingga kombinasi kemungkinan Model regresi Poisson yang terbentuk adalah 15
kombinasi. Tabel 3 berikut adalah 15 kombinasi kemungkinan model regresi Poisson yang dapat dibentuk dan nilai estimasi serta nilai devians pada tiap kemungkinan Model Regrsi Poisson.
Bedasarkan model yang dibentuk pada tabel 3 diatas ada 7 model yang layak digunakan (tanda *) karena nilai devians untuk 7 model tersebut merupakan nilai
devians terkecil untk tiap-tiap kelompok model. Terdapat empat kelompok model berdasarkan jumlah variabel independent yang dimasukkan ke dalam model, yaitu kelompok model dengan satu variabel, dua variabel, tiga variabel dan empat variabel. Model-model yang layak
digunakan untuk pemodelan regresi Poisson , nilai devians dan ࣲ ( ௡ି௞ିଵ , ఈ ) ଶ selanjutnya ditunjukkan pada Tabel 4.
Tabel 4. Analisis Kesesuaian Model Regresi Poisson No. Model Devians ࢊ࢈ ढ ( ࢊ࢈ , ࢻ ) ૛ Devians/db
1. ߤ = exp( ߚ ଴ + ݔ ସ ) 62.0203 33 43.74518 1.879403
2.
ߤ = exp( ߚ ଴ + ݔ ଵ + ݔ ସ ) 61.8416 32 42.58475 1.93255
3.
ߤ = exp( ߚ ଴ + ݔ ଶ + ݔ ସ ) 61.5071 32 42.58475 1.922097
4.
ߤ = exp( ߚ ଴ + ݔ ଷ + ݔ ସ ) 62.0097 32 42.58475 1.937803
5.
ߤ = exp( ߚ ଴ + ݔ ଵ + ݔ ଶ + ݔ ସ ) 61.4489 31 41.42174 1.982223
6.
ߤ = exp( ߚ ଴ + ݔ ଶ + ݔ ଷ + ݔ ସ ) 61.5042 31 41.42174 1.984006
7.
ߤ = exp( ߚ ଴ + ݔ ଵ + ݔ ଶ + ݔ ଷ + ݔ ସ ) 61.4192 30 40.25602 2.047307
Berdasarkan ke-7 model Regresi
Poisson tersebut memiliki nilai Devians yang seluruhnya lebih besar dari nilai Chi-square tabel, yang artinya seluruh model signifikan pada tingkat signifikansi 0.1, sehingga dapat disimpulkan bahwa ke-7 model layak digunakan. Model-model
tersebut mengalami kondisi overdispersi karena berdasarkan nilai Devians yang dibagi dengan derajat bebasnya lebih besar dari 1, sehingga penggunaan Regresi Poisson tidak layak digunakan karena tidak memenuhi asumsi equidispersi . Model yang dapat digunakan adalah Model Regresi ZIP .

## Model Regresi ZIP
Pengujian parameter secara serentak untuk semua alternative model
membuktikan bahwa semua model alternatif layak digunakan.
Tabel 6. Pemilihan Model Terbaik Variabel Model regresi ZIP AIC ࢞ ૝ 138.2 ࢞ ૚ , ࢞ ૝ 136.9* ࢞ ૛ , ࢞ ૝ 141.3 ࢞ ૜ , ࢞ ૝ 140.8 ࢞ ૚ , ࢞ ૛ , ࢞ ૝ 138.0* ࢞ ૛ , ࢞ ૜ , ࢞ ૝ 140.6 ࢞ ૚ , ࢞ ૛ , ࢞ ૜ , ࢞ ૝ 146.7
( keterangan * Model dengan nilai AIC terkecil)
Model dengan nilai AIC terkecil adalah model regresi ZIP dengan dua variabel yaitu x 1 dan x 4 , namun jika dilihat jumlah parameter yang signifikan maka model ini tidak sesuai karena hanya satu parameter yang signifikan dari enam parameter. Model selanjutnya yang memiliki nilai AIC terkecil adalah model dengan tiga variabel independen yaitu x 1 , x 2 dan x 4 . Parameter yang signifikan .adalah ߛ ଴ , ߛ ଵ , ߛ ଶ , ߛ ସ , ߚ ଴ , ߚ ସ . Berikut adalah model regresi ZIP terbaik : log( ߤ ௜ ) = 1.0993 − 0.00006 ݔ ଵ ೔ − 0.00812 ݔ ଶ ೔ + 0.00003 ݔ ସ ೔
dan ݈݋݃݅ݐ ( ߱ ௜ ) = 97.822 − 2.1306 ݔ ଵ ೔ − 1.9009 ݔ ଶ ೔ + 0.0041 ݔ
ସ ೔ dimana ݔ ଵ menyatakan jumlah balita penderita malnutrisi, ݔ ଶ menyatakan jumlah pusat kesehatan masyarakat (puskesmas), dan ݔ ସ menyatakan rata- rata kelengkapan imunisasi. Model logit regresi ZIP menjelaskan bahwa peluang variabel dependen bernilai nol dipengaruhi oleh jumlah balita penderita malnutrisi dan rata-rata kelengkapan imunisasi tiap kabupaten. Sedangkan model log menjelaskan bahwa jumlah kematian balita malnutrisi akan berkurang sebesar 0.00006 sama dengan 1.00006 jika jumlah balita penderita malnutrisi bertambah sebesar satu satuan, dengan syarat bahwa variabel independen yang lain adalah konstan, jumlah kematian balita malnutrisi akan berkurang sebesar 0.00812 sama dengan 1.00812 jika jumlah pelayanan kesehatan masyarakat (puskesmas) bertambah sebesar satu satuan, dengan syarat bahwa variabel independen yang lain adalah konstan dan jumlah kematian balita malnutrisi akan bertambah sebesar
0.00003 sama dengan 1.00003 jika jumlah rata-rata kelengkapan imunisasi bertambah sebesar satu satuan, dengan syarat bahwa variabel independen yang lain adalah konstan.

## KESIMPULAN
1. Model Regresi Poisson tidak memehuni asumsi equidispersi, maka digunakan model lain, model yang diusulkan adalah model Regresi Zero-Inflated Poisson (ZIP) , model Regresi ZIP salah satu kabupaten/kota di provinsi jawa Tengah yaitu Kabupaten Grobogan adalah sebagai berikut :
log( ߤ ହ ) = 1.0993 − 0.00006 ݔ ଵ ఱ − 0.00812 ݔ ଶ ఱ + 0.00003 ݔ ସ ఱ dan
݈݋݃݅ݐ ( ߱ ହ ) = 97.822 − 2.1306 ݔ ଵ ఱ − 1.9009 ݔ ଶ ఱ + 0.0041 ݔ ସ ఱ
Model logit regresi ZIP untuk Kabupaten Grobogan menjelaskan bahwa peluang jumlah kematian
balita malnutrisi di Kabupaten Grobogan bernilai nol dipengaruhi oleh jumlah balita penderita malnutrisi di kabupaten Grobogan dan rata-rata kelengkapan imunisasi di Kabupaten Grobogan. Sedangkan model log menjelaskan bahwa jumlah kematian balita malnutrisi di Kabupaten Grobogan akan berkurang jika jumlah pusat kesehatan masyarakat (Puskesmas) di Kabupaten Grobogan mengalami peningkatan, kondisi ini sesuai
dengan dugaan umum karena jika jumlah Puskesmas di kabupaten Grobogan jumlahnya meningkat maka balita yang menderita malnutrisi dapat dirawat di puskesmas dan jumlah kematian
balita malnutrisi di Kabupaten Grobogan akan berkurang jika jumlah balita penderita malnutrisi di Kabupaten Grobogan mengalami peningkatan, kondisi ini tidak sesuai dengan dugaan umum, hal ini disebabkan karena sesuai dengan analisis sebelumnya karena jika di Kabupaten Grobogan jumlah Puskesmas meningkat maka balita yang mengalami malnutrisi akan segera ditangani oleh tim medis di Puskesmas sehingga walaupun jumlah balita penderita malnutrisi meninggkan tidak akan menambah jumlah kematian balita malnutrisi. Sebaliknya jumlah rata-rata kelengkapan imunisasi di Kabupaten Grobogan mengalami kondisi ini tidak sesuai dengan dugaan umum,
hal ini disebabkan oleh jumlah rata- rata kelengkapan imunisasi di kabupaten Grobogan yang semakin meningkat, kondisi ini terjadi karena
kematian balita malnutrisi tidak sepenuhnya dapat dicegah dengan pemberian imunisasi namun juga kebersihan lingkungan serta asupan gizi diperlukan.
2. Faktor-faktor yang secara signifikan mempengaruhi kejadian kematian balita malnutrisi di Jawa Tengah adalah jumlah balita penderita malnutrisi di tiap kabupaten/kota di Provinsi Jawa tengah dan jumlah pusat kesehatan
masyarakat
(Puskesmas) di tiap kabupaten/kota di
Provinsi Jawa Tengah.

## DAFTAR PUSTAKA
[1] Bohning, D., Dietz, E., Schlattmann, P.2012. Zero Inflated Count Model and Their Application in Public Health and Social Science . Paper dalam http://www.ipn.unikiel.de
(sitasi tanggal 18 Desember 2015 pukul 20.20 WIB).
[2] Hocking, R., 1996, Methods and Applications of Linier Models , New York : John Wiley & Sons.
[3] http:// jateng.bps.go.id
[4] Khoshgoftaar, T.m., Gao,K, dan Szabo, R.M.,2004. Comparing software fault prediction of pure and zero-inflated Poisson regression models. Internasional Jurnal of System Science 36,11 : 705-715.