## GEOGRAPHICALLY WEIGHTED POISSON REGRESSION (GWPR) UNTUK PEMODELAN JUMLAH PENDERITA KUSTA DI JAWA TENGAH
1 Devy Noviani , 2 Rochdi Wasono, 3 Indah Manfaati Nur
1,2,3 Program Studi Statistika FMIPA Universitas Muhammadiyah Semarang Alamat e-mail : devygoo@gmail.com

## ABSTRAK
Regresi poisson adalah metode statistika yang digunakan untuk menganalisa hubungan antara variabel prediktor dan variabel respon dimana variabel respon berbentuk data cacahan atau counted data dan berdistribusi poisson. Data jumlah penderita kusta adalah salah satu contoh data yang asumsinya memenuhi distribusi poisson karena merupakan peristiwa yang jarang terjadi. Dalam penelitian ini peneliti ingin memodelkan jumlah penderita kusta di Provinsi Jawa Tengah dengan menggunakan pendekatan G eographically Weighted Poisson Regression atau GWPR. Penambahan pengaruh aspek spasial diharapkan mampu menghasilkan model terbaik dengan bias kecil dan menghasilkan model yang berbeda-beda ditiap wilayah. Pada pemodelan Pembobotan fungsi kernel yang digunakan untuk penaksiran parameter yaitu pembobotan kernel klasik kernel dan bisquare kernel, dari kedua pembobotan itu akan dipilih yang terbaik yaitu yang memiliki nilai AIC paling kecil. Pemodelan jumlah penderita kusta menggunakan GWPR dengan pembobot kernel gaussian adalah model yang terbaik diantara model regresi poisson dan GWPR dengan pembobot kernel bisquare karena memiliki nilai AIC terkecil. Dari hasil pemodelan menggunakan pembobot kernel gaussian diperoleh 22 kelompok wilayah berdasarkan variabel yag signifikan. Dari 35 kabupaten dan kota ternyata persentase rumah tangga ber-PHBS dan rata-rata lama sekolah penduduk menjadi faktor yang banyak berpengaruh untuk jumlah penderita kusta di Jawa Tengah.
Kata Kunci : Regresi Poisson, GWPR, Penderita Kusta.

## PENDAHULUAN
Kusta disebut juga Morbus Hansen adalah penyakit menular yang disebabkan bakteri
Mycobacterium leprae . Penyakit ini dapat menyebabkan masalah yang kompleks, bukan hanya dari segi medis seperti cacat fisik tetapi juga sampai masalah sosial, ekonomi, budaya, keamanan dan ketahanan nasional. Bila tidak ditangani dengan cermat, kusta dapat menyebabkan cacat dan keadaan ini menjadi penghalang bagi pasien kusta dalam menjalani kehidupan bermasyarakat untuk memenuhi
kebutuhan sosial ekonominya [4].
Menurut WHO dari 2004 sampai 2012
Indonesia menduduki peringkat ke-tiga sebagai negara penyumbang penderita baru kusta terbesar di dunia setelah India dan Brasil.
Penelitian mengenai faktor-faktor dan penyebab penyakit kusta telah banyak dilakukan di Indonesia tapi
sangat terbatas penelitian yang mempertimbangkan aspek geografis antar wilayah. Analisis data secara spasial diperlukan agar dapat mengetahui faktor- faktor apa saja yang berpengaruh terhadap jumlah penderita kusta di
wilayah yang memiliki karakteristik berbeda satu sama lain. Dalam penelitian ini akan dimodelkan jumlah penderita kusta di Provinsi Jawa Tengah menggunakan regresi poisson dan GWPR dengan pembobot fungsi kernel gaussian dan kernel bisquare dengan begitu akan didapatkan variabel yang signifikan untuk pemodelan di setiap wilayah
dengan nilai penaksiran parameter yang berbeda-beda untuk tiap wilayah di Provinsi Jawa Tengah kemudian akan dipilih model terbaik untuk pemodelan jumlah penderita kusta dengan kriteria nilai AIC.
Regresi Poisson merupakan model regresi nonlinier dimana variabel respon mengikuti distribusi poisson. Peluang distribusi dinyatakan sebagai berikut :
,... 2 , 1 , 0 ; ! ) , (    y y e y f y    (1)
Model regresi poisson dapat dinyatakan dalam bentuk berikut ini.
         p j ij j ip p i i i x x x x 1 0 2 2 1 1 0 ... ) ln(        (2)
Penaksir parameter pada model regresi poisson dilakukan dengan metode MLE ( Maximum Likelihood Estimation ). Mencari persamaan likelihood dari fungsi peluang distribusi poisson, Subtitusi μ i terhadap y i dan mencari persamaan likelihood dalam bentuk ln.
                          n i i n i T i i n i T i y x y x L 1 1 1 ! exp ) exp( exp ln ) ( ln    (3)
              n i i n i T i i T i n i y x y x 1 1 1 ln exp   (4)
Persamaan di atas diturunkan terhadap T  yang merupakan bentuk vektor, karena dalam hal ini memiliki beberapa parameter.
           n i i i n i T i i T x y β x x L 1 1 exp ) ( ln   (5)
Persamaan di atas disamadengankan nol lalu diselesaikan denga iterasi
Newton-Raphson. Iterasi berhenti jika taksiran parameter konvergen.
GWPR merupakan model regresi linier lokal yang menghasilkan penaksir parameter model yang bersifat lokal untuk setiap titik atau lokasi dimana data tersebut dikumpulkan [3]. Model GWPR adalah :
              ij i i p j j i i i x v u v u , , exp 1 0    (6)
Keterangan:
μ i = nilai observasi variabel respon ke x ij = nilai observasi variabel prediktor j m pada pengamatan ke- i
u i = koordinat latitude pengamatan ke- i
v i = koordinat longitude prngamatan ke- i β = koefisien regresi Penaksir parameter pada model GWPR menggunakan MLE dengan
menambahkan faktor pembobot letak geografis pada fungsi ln-likelihoodnya sehingga didapatkan persamaan sebagai berikut:
               i i ij i i i n i T i i i i T i i i v u w y v u x y v u x v u L , ! ln , , exp , ln 1 *          (7)
Kemudian diferensiasikan persamaan terhadap   i i v u ,  dan hasilnya harus sama dengan nol.
        0 , , exp 1     i i ij n i i i T i i i i v u w v u x x x y  (8)
Kemudian diselesaikan dengan iterasi Newton-Raphson. Iterasi berhenti jika keadaan konvergen.Menurut [1] salah satu metode yang bisa digunakan untuk menentukan bandwidth optimum adalah metode Cross Validation (CV) yang didefinisikan sebagai berikut:
         n i i i h y y CV 1 2 ˆ (9)
Dimana   h y i  ˆ merupakan nilai penaksir untuk y i dengan radius h , tetapi pengamatan di titik i dikeluarkan dari proses penaksiran. Selanjutnya adalah
mencari nilai pembobot fungsi kernel dengan rumus sebagai berikut.
1. Kernel gaussian
                2 2 1 exp h d w ij ij (10)
2. Kernel bisquare
                      , 0 1 2 2 i ij ij h d w ,untuk i ij h d  (11)
Dimana w ij adalah nilai pembobot fungsi kernel dari lokasi pengamatan j untuk mengestimasi koefisien pada lokasi i dan d ij merupakan jarak eucliden antara lokasi i ditulis ( u i ,v i ) ke lokasi j dutulis ( u j ,v j ), h adalah nilai bandwidth optimum. Jarak eucliden dicari dengan rumus berikut.
    2 2 j i j i ij v v u u d     (12)
Penularan Kusta dapat disebabkan faktor orang, tempat, dan waktu.
Diagnosis kusta dapat ditegakkan dengan adanya kondisi sebagai berikut:
 Kelainan pada kulit (bercak) putih atau kemerahan disertai mati rasa.
 Penebalan saraf tepi yang disertai gangguan fungsi saraf berupa mati rasa dan kelemahan/kelumpuhan otot.
 Adanya kuman tahan asam di dalam kerokan jaringan kulit (BTA Positif).
Klasifikasi kusta menurut [5] disederhanakan menjadi dua yaitu Pausi
Basiler (PB) atau kusta kering dan Multi Basiler (MB) atau kusta basah. PB adalah tipe kusta yang tidak menular dengan tanda-tandanya antara lain bercak seperti panu dan bila disentuh masih terasa. Sedangkan MB adalah tipe kusta yang mudah menular ditandai dengan adanya bercak putih atau kemerahan yang menyebar merata di seluruh badan yang tidak terasa apabila disentuh disertai rusaknya organ tubuh seperti putusnya jari-jari tangan. Orang yang dalam kesehariannya melakukan kontak dan
tinggal serumah dengan penderita kusta jenis MB mempunyai resiko tertular lebih tinggi dari pada dengan penderita kusta jenis PB.

## METODE PENELITIAN

## Sumber Data dan Variabel Penelitian
Data yang digunakan adalah data dari buku Profil Kesehatan Provinsi Jawa Tengah tahun 2012, data Survey Sosial Ekonomi Nasional (SUSENAS) tahun 2012, dan data Badan Pusat Statistik (BPS) Provinsi Jawa Tengah (jateng.bps.go.id) dengan satuan observasi per kabupaten/kota di Jawa Tengah. Variabel yang digunakan adalah sebagai berikut
Tabel 1. Variabel Penelitian
Nama variabel
Y : Jumlah Penderita Kusta
X 1 : Persentase rumah tangga ber-perilaku hidup bersih sehat (PHBS)
X 2 : Persentase rumah sehat
X 3 : Banyaknya dokter
X 4 : Banyaknya puskesmas
X 5 : Persentase penduduk laki-laki
X 6 : Kepadatan penduduk
X 7 : Jumlah penyuluhan kesehatan kelompok
X 8 : Pengeluaran riil per kapita
X 9 : Persentase penduduk miskin
X 10 : Rata-rata lama sekolah penduduk

## Metode Analisis
Tahap-tahap penelitian adalah sebagai berikut: 1. Menentukan variabel-variabel prediktor untuk i ij h d 
2. Mendeskripsikan kabupaten/kota di
Jawa Tengah berdasarkan variabel penelitian. 3. Mengatasi kasus multikolinearitas, digunakan nilai VIF sebagai pendeteksi multikolinearitas.
4. Menganalisis model regresi poisson:
a. Mencari nilai estimasi parameter model regresi poisson.
b. Menguji signifikasi parameter model regresi poisson.
c. Mencari nilai AIC model regresi poisson.
5. Menguji efek heterogenitas spasial dengan uji Breusch-Pagan .
6. Menganalisis model GWPR: a. Menentukan latitude dan longitude b. Menghitung jarak euclidean .
c. Menentukan nilai bandwidth optimum dengan metode Cross Validation (CV).
d. Menghitung matriks pembobot fungsi kernel gaussian dan fungsi kernel bisquare .
e. Mencari nilai estimasi parameter model GWPR.
f. Meguji kesamaan antara model regresi poisson dengan model.
g. Menguji signifikasi parameter model GWPR secara parsial.
h. Mencari nilai AIC model GWPR
7. Pemilihan model terbaik.

## HASIL PENELITIAN

## Deskripsi Jumlah Penderita Kusta
Jumlah penderita kusta di Jawa Tengah tahun 2012 adalah 1513 jiwa yang tersebar diseluruh kabupaten/kota. Angka tertinggi terdapat di Kabupaten Tegal (170 jiwa) sedangkan angka terendah (0 jiwa) terdapat di Kabupaten Semarang, Kota Magelang, dan Salatiga. Dilihat dari gambar 1 wilayah yang masuk ke dalam jumlah penderita kusta kategori tinggi banyak terdapat di bagian utara yaitu Brebes, Pemalang, Pekalongan, dan Tegal artinya di
kabupaten tersebut terdapat penderita kusta antara 83-170 jiwa.
Gambar 1. Persebaran Jumlah Penderita Kusta
Persentase rumah tangga ber-PHBS (X 1 ) tertinggi adalah Kota Pekalongan sebesar 93,9% dan terendah adalah Kab. Banjarnegara sebesar 46,2%. Persentase rumah sehat (X 2 ) tertinggi terdapat di
Kab. Semarang sebesar 98,05% sedangkan persentase rumah sehat
terendah berada di Banjarnegara sebesar 36,12%. Banyaknya dokter (X 3 ) yang termasuk kategori tinggi adalah Kendal (88 jiwa) dan Kota Semarang (108 jiwa), sedangkan terendah adalah Kota Tegal (20 jiwa). Banyaknya puskesmas (X 4 ) terendah di Kota Magelang (5 unit) yang tertinggi di Banyumas (39 unit).
Persentase penduduk laki-laki (X 5 ) terendah adalah Kabupaten Pati 48,54% dan persentase penduduk laki-laki
tertinggi adalah Kabupaten Kendal 50,79%. Untuk variabel kepadatan penduduk (X 6 ) Surakarta adalah wilayah terpadat penduduknya karena Surakarta merupakan merupakan daerah tujuan urbanisasi dan kabupaten/kota lain berkategori sedang. Penyuluhan
kelompok (X 7 ) terbanyak dilakukan di Kab. Kendal (112.764 kali) dan paling sedikit dilakukan di Kab. Blora (66 kali) . Pengeluaran riil per kapita (X 8 ) terendah adalah Kab. Wonosobo sebesar 632,71 (ribu Rp) dan tertinggi adalah Kota Surakarta sebesar 658,92 (ribu Rp). Mayoritas kabupaten/kota di Jawa Tengah masuk di kategori rendah untuk rata-rata lama sekolah (X 10 ) yaitu sekitar
PATI BLORA CILACAP BREBES WONOGIRI GROBOGAN TEGAL DEMAK KEBUMEN JEPARA KENDAL BANYUMAS SRAGEN BOYOLALI BATANG REMBANG PEMALANG MAGELANG SEMARANG KLATEN PURWOREJO WONOSOBO BANJARNEGARA PEKALONGAN KUDUS TEMANGGUNG PURBALINGGA KARANGANYAR KOTA TEGAL SUKOHARJO KOTA SEMARANG KOTA SALATIGA KOTA SURAKARTA KOTA PEKALONGAN KOTA MAGELANG (jiwa) kategori : 0 - 38 (rendah) 39 - 82 (sedang) 83 - 170 (tinggi) N
6,07-7,23 tahun. Rata-rata lama sekolah tertinggi adalah di Kota Surakarta 10,5 tahun dan terendah adalah di Kabupaten Brebes 6,07 tahun.
Gambar 2. Scatterplot Jumlah Penderita Kusta dan Faktor-faktor yang mempengaruhinya
Pengaruh yang positif ditunjukan oleh persentase rumah tangga ber-PHBS (X 1 ), banyaknya dokter (X 3 ), banyaknya puskesmas (X 4 ), dan
persentase penduduk miskin (X 9 ). Variabel yang berpengaruh negatif atau variabel yang memiliki garis turun dari kiri atas ke kanan bawah untuk jumlah penderita kusta adalah persentase rumah sehat (X 2 ), Kepadatan penduduk (X 6 ). Pengeluaran riil per kapita (X 8 ), dan rata-rata lama sekolah (X 10 ).
Tabel 2. Statistik Deskriptif Variabel Penelitian
variabel Mean Var min max Y 43,40 2236,42 0 170 X 1 74,85 168,02 46,2 93,9 X 2 69,85 179,43 36,12 98,05 X 3 53,60 334,84 20 108 X 4 24,94 82,41 5 39 X 5 49,64 0,297 48,54 50,79 X 6 1986,06 5797913,94 472 11573 X 7 10565,3 369364342,37 66 112764 X 8 644,26 49,93 632,71 658,92 X 9 14,42 20,66 5,13 22,5 X 10 7,65 1,39 6,07 10,49 Pemodelan Jumlah Penderita Kusta

## Menggunakan Regresi Poisson
Hasil pengujian Multikolinearitas dapat disimpulkan tidak ada korelasi antar variabel prediktor karena semua variabel memiliki nilai VIF kurang dari
10. Tabel 3. VIF dari 10 Variabel prediktor Variabel VIF Variabel VIF X 1 1,983 X 6 3,529 X 2 2,189 X 7 1,626 X 3 3,210 X 8 2,667 X 4 3,053 X 9 2,405 X 5 1,894 X 10 4,381 Tabel 4. Estimasi Parameter Model Regresi Poisson Β Est Std.E Z hit P β 0 6,115 6,5359 0,936 0,374 β 1 0,0276 0,0033 8,431 0,000 β 2 -0,0091 0,0032 -2,890 0,004 β 3 0,0241 0,0028 8,468 0,000 β 4 -0,0349 0,0054 -6,479 0,000 β 5 0,1186 0,0751 1,578 0,109 β 6 0,0003 0,00003 10,90 0,000 β 7 -0,00001 0,000002 -2,702 0,008 β 8 -0,0022 0,0058 -0,379 0,731 β 9 --0,0118 0,0086 -1,377 0,179 β 10 -1,2151 0,06114 -19,875 0,000 Devians : 858,134; DB : 24; AIC : 880,134; α =5%

## 1. Uji Serentak Parameter Regresi Poisson
Uji secara serentak hasilnya tolak H 0 karena nilai devians = 858,134 lebih besar dari 36,4, atau D >   05 , 0 ; 24 2  ,
sehingga dapat disimpulkan paling sedikit ada satu 0  j  atau dengan kata lain minimal ada salah satu parameter yang berpengaruh secara signifikan terhadap model regresi poisson.
2. Uji Parsial Parameter Regresi Poisson Parameter yang signifikan terhadap model regresi poisson jumlah penderita kusta adalah , , , , , , 7 6 4 3 2 1       10 dan  karena Z hitung > Z α /2 . Jadi model regresi poisson untuk jumlah penderita kusta adalah sebagai berikut.
) 2151 , 1 00001 , 0 0003 , 0 0349 , 0 0241 , 0 0091 , 0 0276 , 0 115 , 6 exp( ˆ 10 7 6 4 3 2 1                
Untuk variabel rata-rata lama sekolah (X 10 ) jika ada pertambahan satu satuan maka akan menurunkan ln jumlah penderita kusta sebesar 1,2151. Hal ini juga berlaku untuk variabel persentase
Y 100 75 50 160 80 0 100 75 50 120 80 40 40 20 0 51 50 49 10000 5000 0 100000 50000 0 660 650 640 160 80 0 24 16 8 160 80 0 10 8 6 X1 X2 X3 X4 X5 X6 X7 X8 X9 X10 Scatterplot of Y vs X1; X2; X3; X4; X5; X6; X7; X8; X9; X10
rumah sehat, apabila bertambah satu persen maka akan menurunkan jumlah penderita kusta sebanyak exp(0,0091). Untuk variabel kepadatan penduduk (X 6 ) jika ada pertambahan satu penduduk/km 2 maka akan menambah jumlah penderita kusta sebesar exp(0,0003). Mengetahui ada tidaknya keragaman spasil pada model dilakukan pengujian Breusch- Pagan , hasilnya diperoleh nilai BP = 18,3979, p -value = 0,0486 pada taraf signifikasi 5%. Kesimpulannya adalah terdapat keragaman spasial antar wilayah.

## Pemodelan Jumlah Penderita Kusta Menggunakan GWPR
Bandwidth optimum untuk pembobot kernel gaussian adalah 63,4352 dengan
AICc minimum 396,188 sedangakan untuk pembobot bisquare didapatkan bandwidth yang berbeda-beda ditiap kabupaten /kota. Matriks diagonal untuk Kabupaten Cilacap yang terbentuk dari bandwidth 63,4352 ( gaussian ) dan
223,31 ( bisquare ) adalah sebagai berkut.
 
       


           0,313 0,172 0,020 0,026 0,006 0,093 0,334 0,028 0,264 0,172 0,142 0,219 0,096 0,00 0,006 0,002 0,002 0,000 0,000 000 , 0 0,001 0,002 0,003 0,003 0,007 0,020 0,020 0,097 0,236 0,220 0,495 0,469 0,689 0,888 00 , 1 , W 1 1 diag v u                        0,660 0,513 0,133 0,170 0,034 0,381 0,678 0,178 0,616 0,513 0,469 0,570 0,387 0,000 0,027 0,000 0,000 0,000 0,000 000 , 0 0,000 0,000 0,002 0,006 0,040 0,135 0,135 0,389 0,588 0,571 0,786 0,771 0,883 0,962 00 , 1 , W 1 1 diag v u

## Uji Kesamaan Model Regresi Poisson dan GWPR
Tabel 5. Uji Kesamaan Model
Model Devians Db Devians/ Db F hitung R.Poisson 858,13 24,00 35,756 GWPR ( Gaussian ) 262,68 10,39 25,272 1,415 GWPR ( Bisquare ) 631,408 19,797 31,894 1.121
Diketahui F hitung untuk model regresi poisson dan GWPR dengan pembobot kernel gaussian yaitu 1,41<F (0,05;24;10) yaitu 2,7 maka F hitung < F tabel disimpulkan gagal tolak H 0 artinya tidak terdapat perbedaan antara model regresi poisson dan model
GWPR dengan pembobot kernel gaussian. Sedangkan untuk model regresi poisson dan model GWPR dengan pembobot kernel bisquare didapatkan nilai F hitung 1,121 < F (0,05;24;19) yaitu 2,08 sehingga gagal tolak H 0 artinya tidak terdapat perbedaan antara model regresi poisson dan GWPR dengan pembobot kernel bisquare.

## Uji Parsial Parameter model GWPR
Untuk mengetahui faktor-faktor yang berpengaruh terhadap jumlah penderita kusta disetiap lokasi dilakukan uji parsial dan membandingkan nilai t hitung dengan nilai t tabel =2,06. Tolak H 0 jika t hitung >2,06 maka kesimpulannya parameter ke- j pada lokasi- i berpengaruh signifikan terhadap model. Kabupaten/kota di Jawa Tengah dikelompokan menjadi 21 kelompok menurut kesamaan variabel yang
signifikan jika dilakukan pemodelan
GWPR dengan pembobot kernel
gaussian . Pemodelan untuk Kabupaten Cilacap adalah :

 10 8 6 3 1 1 6774 , 1 0687 , 0 0007 , 0 0637 , 0 0576 , 0 811 , 70 exp X X X X X       
Model di atas menjelaskan bahwa jumlah penderita kusta akan bertambah sebesar exp(0,0576) jika X 1 bertambah 1% dengan syarat variabel lain konstan, hal yang sama juga berlaku untuk variabel X 3 dan X 6 . Sebaliknya jumlah penderita kusta akan berkurang exp(0,0687) jika variabel X 8 bertambah satu satuan dengan syarat variabel lain konstan, hal sama berlaku untuk variabel X 10 . Kelompok pertama yang memiliki 10 variabel yang signifikan adalah Kabupaten Semarang, Kabupaten Pekalongan, dan Kota Pekalongan.
Gambar 3. Pengelompokan kabupaten/kota berdasarkan variabel yang signifikan dengan pembobot kernel gaussian
Dari Gambar 2 kesamaan variabel yang signifikan tersebut terjadi di lokasi yang saling berdekatan yang ditandai dengan pola yang bebeda-beda untuk masing-masing kelompok. Persentase rumah tangga ber-PHBS (X 1 ) dan rata- rata lama sekolah penduduk (X 10 ) adalah variabel yang berpengaruh signifikan hampir di semua kabupaten/kota di Jawa Tengah.
Jika menggunakan pembobot kernel bisquare terdapat 9 kelompok kabupaten/ kota berdasarkan variabel yang signifikan. Pemodelan untuk Kabupaten Cilacap adalah sebagai berikut.
  10 8 6 4 3 1 1 619 , 1 044 , 0 0006 , 0 023 , 0 049 , 0 049 , 0 579 , 37 exp X X X X X X        
Model di atas menjelaskan bahwa jika ada pertambahan satu satuan dari variabel X 1 maka jumlah penderita kusta akan bertambah sebanyak exp(0,049) dengan syarat variabel lain konstan, hal yang sama berlaku untuk variabel X 3 dan X 6 . Sedangkan jika ada penambahan satu satuan dari variabel X 4 maka akan mengurangi jumlah penderita kusta sebanyak exp(0,023) dengan syarat variabel lain konstan, hal sama berlaku untuk variabel X 8 dan X 10 .
Gambar 4. Pengelompokan kabupaten/kota berdasarkan variabel yang signifikan dengan pembobot kernel bisquare Kelompok terbanyak yaitu kelompok kabupaten yang signifikan terhadap variabel persentase rumah tangga ber- PHBS (X1), persentase rumah sehat (X2), banyaknya dokter (X3), banyaknya puskesmas (X4), kepadatan penduduk (X6), persentase penduduk miskin (X8), dan rata-rata lama sekolah (X10) yang terdapat di Kabupaten Tegal dan beberapa kabupaten/ kota yang berada bagian Timur yang wilayahnya saling berdekatan. Variabel yang signifikan disemua kabupaten/kota adalah pesentase rumah tangga ber-PHBS (X 1 ), banyaknya dokter (X 3 ), kepadatan penduduk (X 6 ), dan rata-rata lama sekolah penduduk (X 10 ).

## Pemilihan Model Terbaik
Tabel 6. Perbandingan nilai AIC
Model Nilai AIC Regresi Poisson 880,134 GWPR ( gaussian ) 307,294 GWPR ( bisquare ) 659,165 Berdasarkan tabel model terbaik untuk pemodelan jumlah penderita kusta di Jawa Tengah tahun 2012 adalah model GWPR dengan pembobot kernel gaussian karena memiliki nilai AIC terkecil yaitu 307,294.
PATI BLORA CILACAP BREBES WONOGIRI GROBOGAN TEGAL DEMAK KEBUMEN JEPARA KENDAL BANYUMAS SRAGEN BOYOLALI BATANG REMBANG PEMALANG MAGELANG SEMARANG KLATEN PURWOREJO WONOSOBO BANJARNEGARA PEKALONGAN KUDUS TEMANGGUNG PURBALINGGA KARANGANYAR KOTA TEGAL SUKOHARJO KOTA SEMARANG KOTA SALATIGA KOTA SURAKARTA KOTA PEKALONGAN KOTA MAGELANG Variabel signifikan : X1, X2, X3, X4, X5, X6, X7, X8, X9, X10 X1, X2, X3, X4, X5, X6, X8, X10 X1, X2, X3, X5, X6, X10 X1, X2, X3, X5, X6, X7, X8, X9, X10 X1, X2, X3, X6, X7, X9, X10 X1, X2, X8, X10 X1, X3, X4, X5, X6, X8, X10 X1, X3, X4, X6, X7, X8, X10 X1, X3, X4, X6, X7, X8, X9, X10 X1, X3, X4, X8 X1, X3, X5, X6, X7, X8, X10
X1, X3, X5, X6, X8, X10 X1, X3, X6, X7, X8, X9, X10 X1, X3, X6, X8, X10 X1, X3, X6, X8, X9, X10 X1, X6, X7, X8, X9, X10 X1, X7, X10 X2, X3, X6, X7, X8, X10 X2, X3, X6, X9, X10 X2, X8, X9, X10 X7, X9, X10 PATI BLORA CILACAP BREBES WONOGIRI GROBOGAN TEGAL DEMAK KEBUMEN JEPARA KENDAL BANYUMAS SRAGEN BOYOLALI BATANG REMBANG PEMALANG MAGELANG SEMARANG KLATEN PURWOREJO WONOSOBO BANJARNEGARA PEKALONGAN KUDUS TEMANGGUNG PURBALINGGA KARANGANYAR KOTA TEGAL X1, X2, X3, X4, X6, X7, X10 X1, X2, X3, X4, X6, X7, X8, X10 X1, X2, X3, X4, X6, X7, X8, X9, X10 X1, X2, X3, X4, X6, X7,X9 ,X10 X1, X2, X3, X4, X6, X8, X10 X1, X2, X3, X4, X6, X8, X9, X10 X1, X2, X3, X4, X6, X9, X10 X1, X2, X3, X5, X6, X8, X10 X1, X3, X4, X6, X8, X10

## KESIMPULAN
Hasil pemodelan dengan regresi poisson variabel yang berpengaruh signifikan adalah persentase rumah tangga ber-PHBS, persentase rumah sehat, banyaknya dokter, banyaknya
puskesmas, kepadatan penduduk, jumlah penyuluhan kesehatan kelompok, dan rata-rata lama sekolah penduduk. Model regresi poisson dapat ditulis sebagai berikut.
) 2151 , 1 00001 , 0 0003 , 0 0349 , 0 0241 , 0 0091 , 0 0276 , 0 115 , 6 exp( ˆ 10 7 6 4 3 2 1       
        
Sedangkan untuk pemodelan GWPR dihasilkan 21 macam pemodelan dengan menggunakan bandwidth optimum 63,4352 untuk pembobot kernel gaussian
dan 9 macam pemodelan dengan
menggunakan bandwidth 223,31 untuk
pembobot kernel bisquare . Model terbaik untuk memodelkan jumlah penderita kusta di Jawa Tengah adalah
menggunakan pendekatan GWPR dengan
pembobot kernel gaussian.

## DAFTAR PUSTAKA
[1] Fortheringham, A.S., Brunsdon, C.,
dan Charlton, M. 2002. Geographically Weighted Regression: The Analysis of Spatially Varying Relationship. Chichester: John Wiley and Sons, ltd. [2] Hocking, R. 1996. Methods and
Application of Linear Models .
NewYork: John Wiley and Sons.
[3] Nakaya, T., Fortheringham, A.S.,
Brunsdon, C. dan Charlton, M. 2004. Geographically Weighted Poisson Regression for Disease Association Mapping. Statistics in Medichine.
24(17):2695-2717.
[4] Widoyono.(2008). Penyakit Tropis: Epidemiologi,
Penularan,
Pencegahan dan Pemberantasannya .
Penerbit Erlangga: Jakarta.
[5] World Health Organisation. 2013. Weekly epidemiological record Relevé épidémiologique hebdomadaire. Di Unduh di
http://www.who.int/entity/ wer/2013/wer8835.pdf?ua=1 diakses 14 Mei 2014.