PENDEKATAN COPULA UNTUK PENYUSUNAN PETA KERAWANAN PUSO TANAMAN PADI DI JAWA TIMUR DENGAN INDIKATOR EL-NINO SOUTHERN OSCILLATION (ENSO) Pratnya Paramitha O.1, Sutikno2 dan Heri Kuswanto3 1 Mahasiswa Jurusan Statistika, ITS, Surabaya 2,3 Dosen Pembimbing, Jurusan Statistika, ITS, Surabaya
[email protected];
[email protected];
[email protected] Abstrak
Makalah ini membahas identifikasi pola hubungan antara indikator ENSO dan luas puso padi di beberapa kabupaten di Jawa Timur dengan metode Copula. Hasil identifikasi digunakan untuk menyusun peta kerawanan puso. Metode Copula digunakan karena pola hubungan antara indikator ENSO dan luas puso padi tidak mengikuti distribusi normal. Penelitian dilakukan di 5 kabupaten, yaitu : Bojonegoro, Lamongan, Jember, Ngawi, dan Banyuwangi. Hasil penelitian menunjukkan analisis korelasi pearson tidak bisa digunakan untuk mengidentifikasi pola hubungan indikator ENSO dan luas puso padi agregat di Jawa Timur dan kabupaten lokasi penelitian. Kesimpulan lain menunjukkan struktur dependensi antara semua indikator ENSO dan luas puso padi agregat di Jawa Timur serta Banyuwangi mengikuti Copula Gumbel. Di Bojonegoro, struktur dependensi anomali SST Nino 1+2, SST Nino 3, SST Nino 4 dan luas puso padi mengikuti Copula Gumbel, SST Nino 3.4 mengikuti Copula Frank, sedangkan di Lamongan dan Jember tidak mengikuti Copula apapun. Di Ngawi, struktur dependensi anomali SST Nino 4, SST Nino 3.4 dan luas puso padi mengikuti Copula Gumbel. Berdasarkan analisis Copula disimpulkan bahwa Kabupaten Bojonegoro, Banyuwangi dan Ngawi rawan puso padi, sedangkan Kabupaten Lamongan dan Jember tidak rawan puso padi. Kata Kunci : Luas Puso Padi, ENSO, Copula
1. PENDAHULUAN Jawa Timur merupakan salah satu provinsi di Indonesia yang sebagian besar penduduknya bermatapencaharian sebagai petani. Produksi padi berkaitan erat dengan keragaman iklim terutama kejadian ekstrem El-Nino dan La-Nina yang seringkali disebut El-Nino Southern Oscillation (ENSO) (Naylor dkk., 2001). Berdasarkan data pada tahun 2004, jumlah lahan puso lahan padi di Jawa Timur seluas 4.570 Ha. Upaya penanggulangan yaitu pemberian subsidi bibit dan pupuk. Upaya lain untuk meminimalkan dampak kerugian produksi pertanian akibat iklim adalah identifikasi wilayah rawan puso akibat anomali iklim, terutama lumbung padi, selanjutnya disusun peta kerawanan puso. Penelitian tentang keterkaitan produksi padi dan indikator ENSO, diantaranya Naylor, Falcon, Wada and Rochberg (2002); Irawan (2006). Naylor dkk. (2002) menyimpulkan ada keterkaitan luas tanaman padi dengan indikator ENSO di Pulau Jawa. Irawan (2006) menyimpulkan pengaruh kejadian ElNino dan La-Nina terhadap produksi pangan. Penelitian tersebut masih bersifat agregat lokasi, keterkaitan yang lebih spesifik lokasi belum banyak dilakukan, termasuk di Jawa Timur. Metode identifikasi pola hubungan produksi pertanian dan indikator ENSO sering menggunakan korelasi Pearson dan pemodelannya dengan metode regresi. Penggunaan kedua metode tersebut harus memenuhi asumsi data mengikuti distribusi normal. Jika asumsi terlanggar maka kesahihan kesimpulan berkurang. Asumsi kenormalan akan sering terlanggar dengan adanya data ekstrem, pola hubungan sulit terdeteksi, meskipun terdapat keterkaitan cukup erat antara kejadian ekstrem dengan produksi pertanian. Pendekatan copula adalah metode yang tepat mengatasi masalah ini. Copula merupakan salah satu metode statistika untuk mengetahui hubungan dua variabel atau lebih, dimana kasusnya distribusi yang dimiliki bisa normal maupun tidak (Sklar, 1959). Penelitian tentang Copula antara lain Embrechts, Lindskog dan McNeil (2001) menerapkan Copula untuk manajemen resiko; Salvadori, Michele, Passoni dan Vezzoli (2007); Villarini, Serinaldi, dan Krajewski. (2008); Schӧlzel dan Friederichs (2008); Syahrir (2011) menerapkan Copula pada bidang klimatologi, hidrologi dan meteorologi. Hasil penelitian menunjukkan metode Copula merupakan solusi yang tepat dalam kondisi asumsi-asumsi kenormalan terlanggar. Dengan menggunakan Copula dapat diketahui hubungan luas puso padi di Jawa Timur dan indikator ENSO, selanjutnya dapat disusun peta kerawanan puso di Jawa Timur. Penelitian ini dilakukan di kabupaten lumbung padi Jawa Timur, yaitu Bojonegoro, Lamongan, Jember, Banyuwangi dan Ngawi.
1
2. TINJAUAN PUSTAKA 2.1 Korelasi Korelasi dibedakan menjadi korelasi untuk statistik parametrik dan nonparamerik. Korelasi untuk statistik parametrik adalah korelasi Pearson, sedangkan korelasi untuk statistik nonparametrik adalah korelasi Spearman dan Tau Kendall. Penjelasan dari masing-masing korelasi tersebut sebagai berikut. a. Korelasi Pearson ( r ) Karl Pearson mengembangkan koefisien korelasi pada tahun 1896, korelasi ini disebut Korelasi Pearson (Dorey dan Joubert, 2007). Menurut Samian (2008) korelasi ini masuk kategori statistik parametrik sehingga ada syarat yang harus dipenuhi. Syarat tersebut adalah data yang digunakan berskala interval atau rasio dan sebaran data mengikuti distribusi normal. Formula korelasi Pearson adalah. n
s r XY s X sY
( X i 1
n
( X i 1
i
i
X )(Yi Y )
X)
2
n
(Y Y ) i 1
i
; -1 r 1
(1)
2
dengan X dan Y merupakan dua variabel yang akan dihitung korelasinya, n menunjukkan banyaknya amatan. Hipotesis : H0 : = 0 (tidak ada korelasi antara X dan Y) dan H1 : ≠ 0 (ada korelasi antara X dan Y) Statistik Uji : t
r n2
~T 1 r2 Tolak H0 jika t t ,dengan v = n-2 atau tolak H0 jika p-value < α. 2
(2)
,v
b. Korelasi Spearman ( rs ) Pada perhitungan korelasi Spearman, data yang dianalisis tidak harus memenuhi asumsi berdistribusi normal. Data yang digunakan untuk korelasi Spearman berskala ordinal. Data merupakan sampel acak yang terdiri atas n pasangan hasil pengamatan (Daniel, 1989). Hipotesis : H0 : = 0 (tidak ada korelasi antara X dan Y) dan H1 : ≠ 0 (ada korelasi antara X dan Y) n2
Statistik Uji : t rs
6
1 rs
n
d
i
2
~ T , untuk n > 30
(3)
2 n
2
n
[ R( X i ) R(Yi )]2 n(n 1) i 1 i 1 R(Xi) adalah peringkat dari Xi dan R(Yi) adalah peringkat dari Yi Tolak H0 jika t t ,dengan v = n-2 atau tolak H0 jika p-value < α. (Naga, 2005).
dengan rs 1
i 1 2
2
;
di
(4)
,v
c. Korelasi Tau Kendall ( ˆ ) Koefisien korelasi Tau Kendall didapatkan berdasarkan peringkat pada hasil amatan, dan memiliki nilai dari -1 hingga +1, data yang digunakan berskala ordinal (Daniel, 1989). Hipotesis : H0 : = 0 (tidak ada korelasi antara X dan Y) dan H1 : ≠ 0 (ada korelasi antara X dan Y) 3ˆ n(n 1) Statistik Uji : Z (5) ~ N (0,1) , untuk n > 10 2(2n 5) S ; S P Q (6) n(n 1) / 2 dimana P adalah banyaknya pasangan berurutan wajar dan Q merupakan banyaknya pasangan berurutan terbalik. Tolak H0 jika Z Z / 2 atau p-value < α. (Naga, 2005).
dengan ˆ
2
2.2 Uji Kolmogorov-Smirnov Uji Kolmogorov-Smirnov digunakan untuk menguji normalitas suatu data. Hipotesis : H0 : F(x) = F0(x) untuk semua nilai x dan H1: F(x) ≠ F0(x) untuk sekurang-kurangnya sebuah nilai x Statistik Uji : D= supx |S(x)- F0(x)| (7) dimana S(x) = (banyaknya nilai pengamatan dalam sampel yang kurang dari/sama dengan x)/n. Tolak H0 jika D > Dα , dimana Dα didapatkan dari tabel statistik Uji Kolmogorov-Smirnov (Daniel, 1989). 2.3 Copula dan Variabel Random Multivariat Menurut Schölzel dan Friederichs (2008), apabila terdapat vektor random x dengan dimensi m yang memiliki fungsi distribusi kumulatif marginal FX1 , FX 2 ,..., FX m dengan domain R yang tidak turun, yaitu FX () 0 dan FX () 1, maka sesuai Teorema Sklar (1959), distribusi bersamanya adalah. 1
1
FX ( x) C X ( FX1 ( x1 ), FX 2 ( x2 ),..., FX m ( xm )) (8) dimana CX:[0,1]x… x[0,1] [0,1] merupakan distribusi bersama dari transformasi variabel random U j FX j ( X j ) untuk j = 1,2,…,m, dengan Uj berdistribusi marginal Uniform.
Apabila distribusi marginal FX ( X j ) kontinu, maka CX adalah unik (Nelsen, 2005) dan persamaan j
untuk CX dapat dituliskan sebagai berikut.
u1
um
0
0
C X (u1 ,..., u m ) ... c X (u1 ,..., u m )du1 ...du m
(9)
dimana CX adalah Copula dan cX adalah persamaan densitas Copula. 2.4 Keluarga dari Copula 1. Copula Elips Copula Elips merupakan Copula dari distribusi elips. Terdapat dua tipe Copula yang termasuk dalam keluarga Copula Elips, yaitu Copula Normal dan Copula-t. a. Copula Normal Vektor random x dengan dimensi m yang memiliki fungsi distribusi kumulatif marginal FX , FX ,..., FX , dengan U j FX ( X j ) ~ U (0,1) , maka untuk masing-masing komponen vektor random bisa ditransformasikan ke distribusi normal standar sebagai berikut. 1
2
m
j
Z j FN(10,1) ( FX j ( X j )) ~ N (0,1)
(10) dimana j=1,2,…,m dan diasumsikan bahwa Z = (Z1,Z2,…,Zm) mengikuti distribusi normal standar multivariat N(0,∑), dengan pdf FN(o,∑) dan matriks covarians ∑. Fungsi Copula Normal adalah. C X (u1 , u 2 ...,u m ) FN (0,) ( FN(10,1) (u1 ), FN(10,1) (u 2 ),...,FN(10,1) (u m )) (11) T
Densitas Copula untuk Copula Normal adalah. c X (u1 , u 2 ,...,u m )
f N ( 0, ) ((FN(10,1) (u1 ), FN(10,1) (u 2 ),..., FN(10,1) (u m ))
m j 1
( f N ( 0.1) ( FN1( 0,1) (u j )))
(12)
b. Copula-t Menurut Schӧlzel dan Friederichs (2008), fungsi Copula untuk Copula-t adalah sebagai berikut. C X (u1 , u 2 ...,u m ) Ft (v,) ( Ft (v1) (u1 ), Ft (v1) (u 2 ),..., Ft (v1) (u m ))
(13) dimana Ft(v) adalah distribusi t univariat dengan derajat bebas v (pada kasus bivariat, nilai v = n-2) dan matriks covarians ∑, mempunyai parameter sebanyak 1+m(m-1)/2. 2. Copula Archimedean Copula Archimedean memiliki bentuk distribusi yang mempunyai tail dependensi yang menunjukkan probabilitas kondisional daerah ekstrem.Persamaan untuk Copula Archimedean adalah. C (u1 , u 2 ,...,u m ) 1{ (u1 ) (u 2 ) ... (u m )} (14) disebut Generator Copula C dan memiliki satu parameter (Kpanzou, 2007). Berbagai macam keluarga dari Archimedean Copula dapat dilihat pada Tabel 1. 3
Tabel 1. Keluarga Copula Archimedean
Generator (u )
Keluarga
Copula Bivariat C (u1 , u 2 )
Clayton
u
Gumbel
( log(u)) , [1, )
1 exp [( log( u1 )) ( log( u 2 )) ]
Frank
eu 1 , R \ 0 log e 1
(eu1 1)(eu2 1) log1 e 1
1
, (0, )
(u1
u2
1)
1
1
(Sumber : Kpanzou, 2007)
Copula Clayton mempunyai tail dependensi di bagian bawah, Copula Frank tidak mempunyai tail dependensi, dan Copula Gumbel mempunyai tail dependensi di bagian atas. 2.5 Estimasi Parameter Copula Menggunakan Maximum Likelihood Estimation Menurut Mikosch (2006), estimasi parameter Copula didapatkan dengan Maximum Likelihood Estimation. Maximum Likelihood Estimation untuk Copula ini diperoleh dari memaksimumkan fungsi log likelihood. Dalam kasus bivariat, fungsi densitas f dapat dituliskan dalam bentuk Copula, yaitu. 2
f (x )
f ( x1 , x 2 ) c{F1 ( x1 ), F2 ( x 2 )}
i 1
i
; x1 , x 2 R
i
(15)
dengan metode Maximum Likelihood Estimation, fungsi likelihood L dan loglikelihood berikut. L
n
f (x j 1
ln L
( j) ( j) 1 , x2 )
n
n
c{F ( x j 1
1
ln c{F ( x j 1
1
2
f (x
( j) ( j) 1 ), F2 ( x2 )}
( j) ( j) 1 ),F2 ( x 2 )}
2
i 1
i
n
ln f i 1 j 1
(16)
( x i( j ) )
(17)
( j) i )
i
Penyelesaian persamaan (17) tidak dapat memperoleh bentuk closed form, sebagai alternatif estimasi parameter Copula Archimedean dilakukan dengan pendekatan Tau Kendall (Nelsen, 2005). Estimasi parameter untuk Copula Archimedean dengan pendekatan Tau Kendall sebagai berikut. Keluarga Clayton Gumbel Frank
Tabel 2. Estimasi Parameter Archimedean Copula
Estimasi
C maka 2 C 1 C 2
1
1 1 maka G 1 G
1 4(1 D1 ( F )) / F , Dk (x) fungsi Debye, Dk ( x) k k x
(Sumber : Embrecth dkk., 2001)
x
0
uk du eu 1
2.6 El-Nino Southern Oscillation (ENSO) ENSO (El-Nino Southern Oscillation) merupakan fenomena peningkatan SPL (Suhu Permukaan Laut)/ SST (Sea Surface Temperature) yang sebelumnya dingin, mengakibatkan perairan yang subur dan kaya ikan menjadi sebaliknya. ENSO dibagi menjadi ENSO hangat (El-Nino) dan ENSO dingin (LaNina). Fenomena El-Nino biasanya disertai perubahan sirkulasi di atmosfer dikenal sebagai osilasi selatan, sehingga El-Nino juga disebut ENSO (El-Nino Southern Oscillation). Berdasarkan NOAA National Weather Service, terdapat empat indikator ENSO, yaitu anomali SST Nino 1+2, SST Nino 3, SST Nino 4 dan SST Nino 3.4 yang masing-masing didapatkan dari Nino 1+2, Nino 3, Nino 4 dan Nino 3.4 sesuai dengan letak Nino tersebut. 3.
METODE PENELITIAN Data yang digunakan adalah data sekunder, yang diperoleh dari Direktorat Perlindungan Tanaman Indonesia dan NOAA National Weather Service yang diaskses melalui http://www.cpc.ncep.noaa.gov/. Data dari Direktorat Perlindungan Tanaman adalah data luas puso tanaman padi di Jawa Timur tahun 1995-2004, meliputi 5 kabupaten, yaitu Kabupaten Bojonegoro, Lamongan, Jember, Banyuwangi dan Ngawi. Data indikator ENSO meliputi anomali SST Nino 1+2, SST Nino 3, SST Nino 4 dan SST Nino 3.4 pada tahun yang sama, dengan periode bulanan. 4
Langkah-langkah analisis data yang digunakan dalam penelitian untuk mencapai tujuan yang diharapkan adalah sebagai berikut. 1. Membuat scatterplot untuk kedua variabel X (indikator ENSO) dan Y (luas puso padi), yaitu. a. X1(anomali SST Nino 1+2) dengan Y b. X2(anomali SST Nino 3) dengan Y c. X3(anomali SST Nino 4) dengan Y d. X4(anomali SST Nino 3.4) dengan Y 2. Menghitung korelasi Pearson, Tau Kendall, dan Spearman untuk variabel-variabel pada langkah (2) 3. Menguji kenormalan data menggunakan histogram dan uji Kolmogorov-Smirnov 4. Membuat scatterplot transformasi [0,1] untuk kedua variabel X (indikator ENSO) dan Y (luas puso padi), yaitu : a. X1(anomali SST Nino 1+2) dengan Y b. X2(anomali SST Nino 3) dengan Y c. X3(anomali SST Nino 4) dengan Y d. X4(anomali SST Nino 3.4) dengan Y 5. Mengidentifikasi struktur dependensi dari scatterplot yang diperoleh 6. Fitting Copula untuk struktur dependensi dari scatterplot yang diperoleh 7. Melakukan Goodness of fit test 8. Melakukan Pemilihan Copula Terbaik 9. Membuat Peta Kerawanan Puso Tanaman Padi 4. HASIL DAN PEMBAHASAN 4.1 Hubungan Luas Puso Padi Agregat di Jawa Timur dan Indikator ENSO dengan Korelasi Pearson, Spearman dan Tau Kendal
0
1 ANOM1.2
(a)
2
3
4
2500
2500
-1
0
1
2
3
-1.0
-0.5
0.0
ANOM4
(b)
1000
Puso.Padi.Jatim
500 0
-1.5
ANOM3
1500
2000
2000 1500 500 0
0
0
-1
1000
Puso.Padi.Jatim
1500 500
1000
Puso.Padi.Jatim
1500 1000 500
Puso.Padi.Jatim
2000
2000
2500
2500
4.1.1 Identifikasi Pola Hubungan Menggunakan Scatterplot Struktur dependensi antara masing-masing variabel anomali SST Nino 1+2, SST Nino 3, SST Nino 4, SST Nino 3.4 dan luas puso padi agregat di Jawa Timur digambarkan dalam bentuk scatterplot. Scatterplot hasil analisis disajikan pada Gambar 1.
(c)
0.5
1.0
-2
-1
0
1
2
ANOM3.4
(d)
Gambar 1. Scatterplot antara Luas Puso Padi di Jawa Timur dan Anomali SST : Nino 1+2 (a) Nino 3 (b) Nino 4 (c) Nino 3.4 (d)
Gambar 1 menunjukkan bahwa pola hubungan antara luas puso padi agregat Jawa Timur dan SST adalah acak, beberapa amatan berdekatan dan beberapa amatan yang lain nampak berjauhan. Untuk mengidentifikasi hubungan antara variabel-variabel tersebut tidak cukup menggunakan scatterplot saja dikarenakan plot-plotnya juga tidak membentuk pola yang spesifik, sehingga sulit untuk dijelaskan. 4.1.2 Analisis Korelasi Analisis hubungan luas puso padi agregat Jawa Timur dan indikator ENSO dengan korelasi Pearson, Spearman dan Tau Kendall disajikan pada Tabel 4.
5
Tabel 4. Nilai Korelasi, Statistik Uji, P-value antara Anomali SST dan Luas Puso Padi Agregat Jawa Timur No SST Metode Nilai Korelasi Statistik Uji P-Value 0,0466 t = 0,5063 0,6136 Nino 1+2 Pearson 1 0,0719 S = 267273,2 0,4351 Spearman 0,0541 z = 0,8622 0,3886 Tau Kendall 0,1546 t = 1,6998 0,0918 Nino 3 Pearson 2 0,1418 S = 247141,8 0,1223 Spearman 0,0921 z = 1,4696 0,1417 Tau Kendall 0,1181 t = 1,2921 0,1989 Nino 4 Pearson 3 0,0819 S = 264384,1 0,3736 Spearman 0,0519 z = 0,8281 0,4076 Tau Kendall 0,1541 t = 1,6944 0,0928 Nino 3.4 Pearson 4 0,1498 S = 244839,8 0,1025 Spearman 0,0931 z = 1,4855 0,1374 Tau Kendall
Tabel 4 menunjukkan bahwa semua indikator ENSO (anomali SST Nino 1+2, anomali SST Nino 3, anomali SST Nino 4 dan anomali SST Nino 3.4) tidak berkorelasi dengan luas puso padi agregat di Jawa Timur. Hal ini ditunjukkan oleh p-value semua metode baik Pearson, Spearman maupun Tau Kendall lebih dari α=0.05, sehingga gagal tolak H 0 (H0 : tidak ada korelasi antar variabel).
(c)
20
frequency
30
frequency
0
0
0
10
20
10
20
40
frequency
60
40
30
50
(b)
80
(a)
40
4.1.3 Pengujian Kenormalan Data Pengujian kenormalan data dilakukan menggunakan histogram dan uji Kolmogorov-Smirnov. Histogram masing-masing variabel SST dan luas puso padi agregat Jawa Timur disajikan pada Gambar 2.
1500
2000
2500
-2
-1
2
3
4
-2
5
-1
(e)
1
2
3
4
DataJatim$ANOM3
frequency
10
20
0
0
0
5
10
frequency
1
DataJatim$ANOM1.2
30
(d)
0
25
1000
DataJatim$Puso.Padi.Jatim
20
500
15
0
-2.0
-1.5
-1.0
-0.5
0.0
0.5
1.0
1.5
-2
DataJatim$ANOM4
-1
0
1
2
3
DataJatim$ANOM3.4
Gambar 2. Histogram : Luas Puso Padi Jawa Timur (a) Anomali SST Nino 1+2 (b) Anomali SST Nino 3 (c) Anomali SST Nino 4 (d) Anomali SST Nino 3.4 (e)
Gambar 2 (a) menunjukkan bahwa distribusi luas puso padi agregat di Jawa Timur tidak berdistribusi normal, karena bentuk histogram cenderung turun dan memiliki kemiringan (skewness) ke kanan. Demikian juga anomali SST Nino 1+2, SST Nino 3, SST Nino 4, dan SST Nino 3.4 menunjukkan bahwa tidak memiliki distribusi normal. Hal ini ditunjukkan oleh histogram tidak simetris, memiliki kemiringan dan memiliki tail panjang (Gambar 2 b-e). Selanjutnya dilakukan uji Kolmogorov-Smirnov, seperti yang disajikan pada Tabel 5. Tabel 5. Nilai KS (D) dan P-value Uji Kolmogorov-Smirnov LuasPuso Padi Jawa Timur dan Anomali SST No. Variabel Nilai KS (D) P-Value 1 2 3 4 5
Luas Puso Padi Agregat Jawa Timur Anomali SST Nino 1+2 Anomali SST Nino 3 Anomali SST Nino 4 Anomali SST Nino 3.4
0,262 0,202 0,168 0,126 0,083
< 0,01 < 0,01 < 0,01 < 0,01 0,046
Berdasarkan Tabel 5 dapat disimpulkan bahwa p-value dari masing-masing variabel kurang dari α = 0,05 dan nilai D>D0,05 = 0,124 kecuali untuk Anomali SST Nino 3.4, sehingga disimpulkan Tolak H 0 yaitu dengan H0 : Data mengikuti distribusi Normal, maka data Luas Puso Padi Jawa Timur, Anomali SST Nino 1+2, Anomali SST Nino 3, dan Anomali SST Nino 4 tidak mengikuti distribusi normal. Untuk Anomali 6
SST Nino 3.4, nilai D
Identifikasi hubungan Luas Puso Padi Agregat di Jawa Timur dan Indikator ENSO dengan Pendekatan Copula
0.0
0.2
0.4
0.6
0.8
0.0
1.0
1.0 0.8
1.0
0.6
0.8
v
0.6
0.4
v
0.2
0.4 0.2
0.2
0.2
0.4
0.4
v
v
0.6
0.6
0.8
0.8
1.0
1.0
4.2.1 Scatterplot Hasil Transformasi Variabel pada Domain Uniform [0,1] Sebelum melakukan analisis Copula, terlebih dahulu data yang digunakan harus ditransformasi ke domain Uniform [0,1]. Scatterplot untuk masing-masing anomali SST Nino 1+2, SST Nino 3, SST Nino 4, SST Nino 3.4 dan luas puso padi agregat di Jawa Timur setelah terlebih dahulu ditransformasi ke domain [0,1] disajikan pada Gambar 3.
0.2
0.4
0.6
0.8
1.0
0.0
0.2
0.4
u
u
(a)
(b)
(c)
0.6
0.8
1.0
0.0
0.2
0.4
0.6
0.8
1.0
u
u
(d)
Gambar 3. Scatterplot Transformasi [0,1] antara Luas Puso Padi di Jawa Timur dengan Anomali SST : Nino 1+2 (a) Nino 3 (b) Nino 4 (c) Nino 3.4 (d)
Gambar 3 menjelaskan bahwa beberapa amatan berdekatan, sehingga diduga antara variabel tersebut ada hubungan meskipun hubungannya kecil. Selanjutnya dilakukan analisis Copula untuk mengidentifikasi apakah ada struktur dependensi khusus. 4.2.2 Estimasi Parameter Copula Pendekatan Copula yang digunakan adalah Copula Clayton, Gumbel dan Frank karena variabel yang dianalisis tidak mengikuti distribusi normal. Estimasi parameter Copula Clayton, Gumbel dan Frank dengan pendekatan Tau Kendall untuk masing-masing variabel disajikan pada Tabel 6. Tabel 6. Estimasi Parameter Copula dengan Pendekatan Tau Kendall untuk SST Nino dan Luas Puso Padi Agregat di Jawa Timur No.
Variabel
1
SST Nino 1+2 dan Luas Puso Padi SST Nino 3 dan Luas Puso Padi SST Nino 4 dan Luas Puso Padi SST Nino 3.4 dan Luas Puso Padi
2 3 4
Jenis Copula Clayton Gumbel Frank Clayton Gumbel Frank Clayton Gumbel Frank Clayton Gumbel Frank Normal
Parameter
P-Value
0,1143 1,0571 0,4876 0,2029 1,1015 0,8349 0,1096 1,0548 0,4685 0,2053 1,1027 0,8439 0,1457
0,3951 0,0000 0,3709 0,1281 0,0000 0,0983 0,3877 0,0000 0,3643 0,1174 0,0000 0,0887 0,0821
Tabel 6 menunjukkan struktur dependensi masing-masing anomali SST Nino 1+2, SST Nino 3, SST Nino 4, SST Nino 3.4 dan luas puso padi agregat di Jawa Timur semuanya mengikuti Copula Gumbel, karena p-value untuk estimasi parameter Copula tersebut kurang dari α = 0,05. Langkah selanjutnya adalah melakukan model fitting Copula Gumbel dengan MLE (Maximum Likelihood Estimation) dengan hasil pada Tabel 7. 7
Tabel 7. Model Fitting Copula Gumbel dengan MLE (Jawa Timur)
No.
Variabel
Jenis Copula
Estimasi
Log Likelihood
P-Value
1
SST Nino 1+2 dan Luas Puso Padi
Gumbel
1,0459
0,3219
0,0000
2
SST Nino 3 dan Luas Puso Padi
Gumbel
1,1033
1,5055
0,0000
3
SST Nino 4 dan Luas Puso Padi
Gumbel
1,0459
0,2746
0,0000
4
SST Nino 3.4 dan Luas Puso Padi
Gumbel
1,1050
1,4996
0,0000
0.0
0.2
0.4
(a)
0.6
0.8
1.0
1.0
0.8
0.8
v
0.6
0.6
0.0
1.0
0.4
v
0.2 0.0
0.2 0.0
0.0
0.0
0.2
0.2
0.4
0.4
0.4
v
v
0.6
0.6
0.8
0.8
1.0
1.0
Tabel 7 memperlihatkan bahwa p-value model fitting Copula Gumbel kurang dari α = 0,05. Hal ini menunjukkan estimasi parameter Copula Gumbel menggunakan pendekatan Tau Kendall signifikan, sehingga Copula Gumbel yang terpilih. Scatterplot rank Copula Gumbel disajikan pada Gambar 4.
u
0.2
0.4
0.6
0.8
0.0
1.0
u
(b)
0.2
0.4
0.6
0.8
0.0
1.0
u
(c)
0.2
0.4
0.6
0.8
1.0
u
(d)
Gambar 4. Scatterplot rank Copula Gumbel antara Luas Puso Padi di Jawa Timur dengan Anomali SST : Nino 1+2 (a) Nino 3 (b) Nino 4 (c) Nino 3.4 (d).
0.0
0.2
0.4
0.6 u
(a)
0.8
1.0
1.0 0.8
1.0
0.6
0.8
v
0.6
0.4
v
0.2
0.4
0.0
0.2 0.0
0.0
0.0
0.2
0.2
0.4
0.4
v
v
0.6
0.6
0.8
0.8
1.0
1.0
Gambar 4 memperlihatkan Scatterplot rank dari Copula Gumbel dengan parameter θ sesuai dengan estimasi dengan pendekatan tau kendall. Bentuk rank Copula Gumbel yang tail dependensinya pada bagian atas kurang begitu jelas terlihat, karena terkendala jumlah sampel yang kecil, yaitu 120. Dengan menggunakan parameter θ dari Copula Gumbel, dapat dibangkitkan data sebanyak 2500 yang bentuk Copula Gumbelnya akan terlihat lebih jelas pada Gambar 5.
0.0
0.2
0.4
0.6
0.8
1.0
0.0
u
(b)
0.2
0.4
0.6
0.8
1.0
0.0
0.2
0.4
0.6
u
u
(c)
(d)
0.8
1.0
Gambar 5. Scatterplot rank Copula Gumbel bangkitan data (n = 2500) sesuai dengan parameter : θ = 1,0571 (a) θ = 1,1015 (b) θ = 1,0548 (c) θ = 1,1027 (d).
4.3
Hubungan Luas Puso Padi di Kabupaten Lumbung Padi Jawa Timur dan Indikator ENSO dengan Korelasi Pearson, Spearman dan Tau Kendal
4.3.1 Identifikasi Pola Hubungan Menggunakan Scatterplot Scatterplot indikator ENSO dan luas puso padi yang disajikan adalah di salah satu Kabupaten lumbung padi, yaitu Bojonegoro, karena kabupaten yang lain memiliki hasil yang hampir sama. Scatterplot yang dihasilkan dapat disajikan pada Gambar 6.
8
1
2
3
-1
4
0
1
2
3
(a)
1500 500 0
-1.5
-1.0
ANOM3
ANOM1.2
1000
Puso.Padi.Bojonegoro
1500 1000
Puso.Padi.Bojonegoro
500 0
0
0
2000
2000
2000 1500
Puso.Padi.Bojonegoro
500
1000
2000 1500 1000
Puso.Padi.Bojonegoro
500 0
-1
-0.5
0.0
0.5
1.0
-2
(b)
-1
0
1
2
ANOM3.4
ANOM4
(c)
(d)
Gambar 6. Scatterplot antara Luas Puso Padi di Bojonegoro dengan Anomali SST : Nino 1+2 (a) Nino 3 (b) Nino 4 (c) Nino 3.4 (d).
Gambar 6 memperlihatkan bahwa plot-plot memiliki kerapatan, namun ada beberapa yang terpencar. Hal tersebut menimbulkan dugaan adanya hubungan antar variabel. Kesimpulan tentang hubungan antar variabel sulit diperoleh jika hanya melihat bentuk scatterplot saja, karena tidak menunjukkan pola spesifik. Kejadian sama juga terjadi di Lamongan, Jember, Banyuwangi dan Ngawi. 4.3.2 Analisis Korelasi Analisis korelasi untuk anomali SST Nino 1+2, SST Nino 3, SST Nino 4, SST Nino 3.4 dan luas puso padi di kabupaten lumbung padi Jawa Timur disajikan pada Tabel 8-12. Tabel 8. Nilai Korelasi, Statistik Uji, P-value antara Anomali SST dan Luas Puso Padi di Bojonegoro No.
Variabel
1
SST Nino 1+2 dan Luas Puso Padi
2
SST Nino 3 dan Luas Puso Padi
3
SST Nino 4 dan Luas Puso Padi
4
SST Nino 3.4 dan Luas Puso Padi
Jenis Korelasi Pearson Spearman Tau Kendall Pearson Spearman Tau Kendall Pearson Spearman Tau Kendall Pearson Spearman Tau Kendall
Nilai Korelasi
P-Value
-0,0743 0,0631 0,0487
0,4202 0,4935 0,4920
0,0008 0,2098 0,1609
0,9934 0,0215 0,0232
0,0600 0,1742 0,1328
0,5150 0,0571 0,0610
0,0385 0,2474 0,1897
0,6760 0,0064 0,0074
Tabel 9. Nilai Korelasi, Statistik Uji, P-value antara Anomali SST dan Luas Puso Padi di Lamongan
Tabel 10. Nilai Korelasi, Statistik Uji, P-value antara Anomali SST dan Luas Puso Padi di Jember No.
Variabel
1
SST Nino 1+2 dan Luas Puso Padi
2
SST Nino 3 dan Luas Puso Padi
3
SST Nino 4 dan Luas Puso Padi
4
SST Nino 3.4 dan Luas Puso Padi
Jenis Korelasi Pearson Spearman Tau Kendall Pearson Spearman Tau Kendall Pearson
Nilai Korelasi -0,0729 -0,0329 -0,0248
Spearman Tau Kendall Pearson Spearman Tau Kendall
0,4288 0,7215 0,7351 0,4015 0,3878 0,3720
-0,0960
0,2967
-0,0236 -0,0192
0,7979 0,7942
-0,1068 -0,0877 -0,0715
0,2458 0,3406 0,3299
Variabel
1
SST Nino 1+2 dan Luas Puso Padi
2
SST Nino 3 dan Luas Puso Padi
3
SST Nino 4 dan Luas Puso Padi
4
SST Nino 3.4 dan Luas Puso Padi
Jenis Korelasi Pearson Spearman Tau Kendall Pearson Spearman Tau Kendall Pearson Spearman Tau Kendall Pearson Spearman Tau Kendall
Nilai Korelasi
P-Value
-0,0801 -0,0842 -0,0709
0,3847 0,3606 0,3131
-0,0981 -0,0682 -0,0473
0,2866 0,4591 0,5019
-0,1386 -0,1152 -0,0859
0,1310 0,2103 0,2224
-0,1271 -0,0938 -0,0677
0,1666 0,3083 0,3356
Tabel 11. Nilai Korelasi, Statistik Uji, P-value antara Anomali SST dan Luas Puso Padi di Banyuwangi
P-Value
-0,0773 -0,0795 -0,0656
No.
9
No.
Variabel
1
SST Nino 1+2 dan Luas Puso Padi
2
SST Nino 3 dan Luas Puso Padi
3
SST Nino 4 dan Luas Puso Padi
4
SST Nino 3.4 dan Luas Puso Padi
Jenis Korelasi Pearson Spearman Tau Kendall Pearson Spearman Tau Kendall Pearson Spearman Tau Kendall Pearson Spearman Tau Kendall
Nilai Korelasi
P-Value
0,3359 0,2809 0,2267
0,0002 0,0019 0,0021
0,3699 0,2882 0,2328
0,0321 0,0014 0,0016
0,1679
0,0667
0,1921 0,1535
0,0355 0,0375
0,3298 0,2697 0,2149
0,0002 0,0028 0,0036
Tabel 12. Nilai Korelasi, Statistik Uji, P-value antara Anomali SST dan Luas Puso Padi di Ngawi No.
Variabel
1
SST Nino 1+2 dan Luas Puso Padi
2
SST Nino 3 dan Luas Puso Padi
3
SST Nino 4 dan Luas Puso Padi
4
SST Nino 3.4 dan Luas Puso Padi
Jenis Korelasi Pearson Spearman Tau Kendall Pearson Spearman Tau Kendall Pearson Spearman Tau Kendall Pearson Spearman Tau Kendall
Nilai Korelasi
P-Value
-0,1121 -0,1796 -0,1436 -0,0004 0,0446 0,0355 0,1152 0,1156 0,0931 0,0494 0,1145
0,2229 0,0497 0,0515 0,9967 0,6284 0,6304 0,2100 0,2085 0,2069 0,5921 0,2132
0,0920
0,2119
Tabel 8 memperlihatkan bahwa anomali SST Nino 3 dan SST Nino 3.4 mempunyai korelasi dengan luas puso padi Bojonegoro, karena nilai korelasi Spearman dan Tau Kendall signifikan. Tabel 9 dan 10 menunjukkan bahwa nilai korelasi antar variabel tidak ada yang signifikan untuk Kabupaten Lamongan dan Jember. Tabel 11 memperlihatkan semua nilai korelasi signifikan kecuali korelasi Pearson antara anomali SST Nino 4 dan luas puso padi di Banyuwangi karena p-value > α = 0,05. Tabel 12 menunjukkan bahwa hanya korelasi Spearman untuk SST Nino 1+2 dan luas puso padi di Ngawi yang nilainya signifikan. Pola hubungan SST Nino dan luas puso padi Kabupaten lumbung padi Jawa Timur juga tidak bisa langsung disimpulkan melalui uji korelasi ini, dikarenakan sebelumnya saat pembuatan scatterplot terdapat plot yang sangat berdekatan satu sama lain.
1000
1500
100 80 60
frequency
20
2000
0
200
400
600
800
1000
1200
0
20
LumbungPadi$Puso.Padi..Lamongan
40
60
80
120
LumbungPadi$Puso.Padi.Jember
(e) frequency
40
60 0
0
20
20
40
frequency
80
80
100
(d)
100
500
LumbungPadi$Puso.Padi.Bojonegoro
60
0
0
0
0
20
20
40
40
60
frequency
80
80 60 40
frequency
(c)
100
(b)
100
(a)
120
4.3.3 Pengujian Kenormalan Data Pengujian kenormalan untuk indikator ENSO, sudah dilakukan sebelumnya yaitu pada Gambar 2 dan Tabel 5 yang menghasilkan indikator ENSO tidak mengikuti distribusi normal. Histogram luas puso padi di Kabupaten lumbung padi Jawa Timur disajikan pada Gambar 7.
0
50
100
150
200
250
0
300
100
200
300
400
500
LumbungPadi$Puso.Padi..Ngawi
LumbungPadi$Puso.Padi.Banyuwangi
Gambar 7. Histogram Luas Puso Padi : Bojonegoro (a) Lamongan (b) Jember (c) Banyuwangi (d) Ngawi (e)
Gambar 7 menunjukkan bahwa luas puso padi di lima Kabupaten lumbung padi Jawa Timur diduga sebarannya tidak mengikuti distribusi normal karena terlihat bahwa histogram tidak simetris, luas puso padi rata-rata kebanyakan berada pada rentang 0-50 ha. Untuk lebih jelasnya, dilakukan pengujian kenormalan data menggunakan uji Kolmogorov-Smirnov yang disajikan pada Tabel 13. No. 1 2 3 4 5
Tabel 13. Nilai KS (D) dan P-value Uji Kolmogorov-Smirnov Luas Puso Padi di Kabupaten Lumbung Padi Jawa Timur dan Anomali SST Variabel Nilai KS (D) P-Value Luas Puso Padi Bojonegoro 0,404 < 0,01 Luas Puso Padi Lamongan 0,410 < 0,01 Luas Puso Padi Jember 0,474 < 0,01 Luas Puso Padi Banyuwangi Luas Puso Padi Ngawi
0,482 0,491
10
< 0,01 < 0,01
Tabel 13 menunjukkan bahwa luas puso padi di Bojonegoro, Lamongan, Jember, Banyuwangi dan Ngawi tidak mengikuti distribusi normal, karena didapatkan p-value < α = 0,05 dan nilai D >D0,05 = 0,124 sehingga diambil keputusan tolak H 0 (H0 : Data mengikuti distribusi normal). Sama seperti analisis untuk luas puso padi agregat di Jawa Timur sebelumnya, hasil analisis korelasi juga tidak cukup bisa mengidentifikasi pola hubungan variabel-variabel di lima Kabupaten lumbung padi. Oleh karena itu, dilakukan analisis Copula untuk melihat ada tidaknya struktur dependensi yang khusus antara luas puso padi di lima Kabupaten lumbung padi tersebut dengan masing-masing indikator ENSO. 4.4
Identifikasi hubungan Luas Puso Padi di Kabupaten Lumbung Padi di Jawa Timur dan Indikator ENSO dengan Pendekatan Copula
0.2
0.4
0.6
0.8
1.0
1.0
1.0
0.9
0.9
v
0.6 0.5 0.4
0.4 0.0
0.2
u
0.4
0.6
0.8
1.0
0.0
0.2
0.4
u
(a)
0.7
0.8
0.8 0.6 0.5
0.5 0.4
0.4 0.0
0.7
v
0.7 0.6
v
0.7 0.5
0.6
v
0.8
0.8
0.9
0.9
1.0
1.0
4.4.1 Scatterplot Hasil Transformasi Variabel pada Domain Uniform [0,1] Transformasi semua variabel ke Uniform [0,1] dilakukan sebelum analisis Copula. Scatterplot pada Gambar 8 adalah scatterplot transformasi [0,1] untuk indikator ENSO dan luas puso padi di Bojonegoro.
0.6
0.8
1.0
0.0
0.2
(b)
0.4
0.6
0.8
1.0
u
u
(c)
(d)
Gambar 8. Scatterplot Transformasi [0,1] antara Luas Puso Padi di Bojonegoro dengan Anomali SST : Nino 1+2 (a) Nino 3 (b) Nino 4 (c) Nino 3.4 (d).
Gambar 8 memperlihatkan bahwa setelah ditransformasi ke domain Uniform [0,1], masih terlihat plot-plot yang rapat satu sama lain, maka diduga antar variabel terdapat hubungan. Hal yang sama juga terjadi pada empat kabupaten yang lain, yaitu Lamongan, Jember, Banyuwangi dan Ngawi. 4.4.2 Estimasi Parameter Copula Indikator ENSO dan luas puso padi di Kabupaten lumbung padi Jawa Timur tidak mengikuti distribusi normal, sehingga Copula yang digunakan adalah Copula Clayton, Gumbel dan Frank. Estimasi parameter Copula dengan pendekatan Tau Kendall pada disajikan pada Tabel 14-18. Tabel 14. Estimasi Parameter Copula dengan Pendekatan Tau Kendall untuk SST Nino dan Luas Puso Padi di Bojonegoro No. 1
2
3
4
Jenis Copula
Parameter
P-Value
SST Nino 1+2 dan Luas Puso Padi SST Nino 3 dan Luas Puso Padi SST Nino 4 dan Luas Puso Padi
Clayton
0,1024
0,5965
Gumbel
1,0512
0,0000
Frank
0,4392
0,5793
Clayton Gumbel Frank
0,3837
0,0949
1,1918
0,0000
1,4799
0,0565
Clayton
0,3064
0,1462
Gumbel
1,1532
0,0000
Frank
1,2130
0,1035
SST Nino 3.4 dan Luas Puso Padi
Clayton Gumbel Frank Normal
0,4682
0,0489
1,2341 1,7589 0,2936
0,0000 0,0220 0,0123
Variabel
Tabel 15. Estimasi Parameter Copula dengan Pendekatan Tau Kendall untuk SST Nino dan Luas Puso Padi di Lamongan No. 1
2
3
4
11
Jenis Copula
Parameter
P-Value
SST Nino 1+2 dan Luas Puso Padi SST Nino 3 dan Luas Puso Padi SST Nino 4 dan Luas Puso Padi
Clayton
-0,1325
0,3374
Gumbel
-
-
Frank
-0,6414
0,3743
Clayton Gumbel Frank
-0,0902 -0,4260
0,5688 0,5877
Clayton
-0,1582
0,3224
-
-
Frank
-0,7777
0,3679
SST Nino 3.4 dan Luas Puso Padi
Clayton Gumbel Frank Normal
-0,1269 -0,6119 -0,1062
0,4232 0,4566 0,4515
Variabel
Gumbel
Tabel 16. Estimasi Parameter Copula dengan Pendekatan Tau Kendall untuk SST Nino dan Luas Puso Padi di Jember No. 1
2 3
4
Variabel SST Nino 1+2 dan Luas Puso Padi SST Nino 3 dan Luas Puso Padi SST Nino 4 dan Luas Puso Padi SST Nino 3.4 dan Luas Puso Padi
Jenis Copula
Parameter
P-Value
Clayton
-0,0485
0,7821
Gumbel Frank Clayton Gumbel Frank Clayton Gumbel Frank Clayton Gumbel Frank Normal
-0,2238 -0,1231 -0,5922 -0,0376 -0,1726 -0,1335 -0,6465 -0,1121
Tabel 17. Estimasi Parameter Copula dengan Pendekatan Tau Kendall untuk SST Nino dan Luas Puso Padi di Banyuwangi No. 1
0,7874 0,4889 0,5189 0,8469 0,8498 0,4456 0,4802 0,4747
2 3
4
Variabel SST Nino 1+2 dan Luas Puso Padi SST Nino 3 dan Luas Puso Padi SST Nino 4 dan Luas Puso Padi SST Nino 3.4 dan Luas Puso Padi
Jenis Copula
Parameter
P-Value
Clayton
0,5864
0,0502
Gumbel Frank Clayton Gumbel Frank Clayton Gumbel Frank Clayton Gumbel Frank Normal
1,2932 2,1304 0,6068 1,3034 2,1926 0,3626 1,1813 1,4082 0,5476 1,2738 2,0107 0,3313
0,0000 0,0201 0,0478 0,0000 0,0184 0,1612 0,0000 0,1113 0,0703 0,0000 0,0327 0,0165
Tabel 18. Estimasi Parameter Copula dengan Pendekatan Tau Kendall untuk SST Nino dan Luas Puso Padi di Ngawi No.
Variabel
1
SST Nino 1+2 dan Luas Puso Padi
2
SST Nino 3 dan Luas Puso Padi
3
SST Nino 4 dan Luas Puso Padi
4
SST Nino 3.4 dan Luas Puso Padi
Jenis Copula Clayton Gumbel Frank Clayton Gumbel Frank Clayton Gumbel Frank Clayton Gumbel Frank Normal
Parameter -0,2512 -1,3147 0,0736 1,0368 0,3197 0,2054 1,1027 0,8439 0,2027 1,1014 0,8341 0,1441
P-Value 0,1299 0,2005 0,7192 0,0000 0,7099 0,3659 0,0000 0,3256 0,3567 0,0000 0,3167 0,3066
Tabel 14 menunjukkan bahwa untuk variabel anomali SST Nino 1+2, SST Nino 3, SST Nino 4 dan Luas Puso Padi di Bojonegoro, parameter Copula yang signifikan adalah Copula Gumbel karena mempunyai p-value < α = 0,05, sedangkan untuk anomali SST Nino 3.4 dan luas puso padi di Bojonegoro, parameter semua Copula (Clayton, Gumbel, Frank dan Normal) signifikan, maka struktur dependensi masing-masing variabel mengikuti Copula yang terpilih. Pada tabel 15 dan 16, yaitu untuk kabupaten Lamongan dan Jember tidak terdapat parameter Copula apapun yang signifikan, kemudian Copula Gumbel pada tabel tersebut tidak dapat diestimasi karena hasil perhitungan menunjukkan bahwa nilai θ dari Copula Gumbel kurang dari 1, padahal seharusnya berapa pada [1, ), hal ini disebabkan nilai korelasi Tau Kendall untuk SST Nino dan luas puso padi di Lamongan dan Jember bernilai negatif. Struktur dependensi variabel di Lamongan dan Jember tidak mengikuti Copula. Tabel 17 yang memperlihatkan estimasi parameter untuk Banyuwangi, menghasilkan bahwa untuk anomali SST Nino 1+2 dan luas puso padi, Copula yang parameternya signifikan adalah Gumbel dan Frank. Untuk anomali SST Nino 3 dan luas puso padi, struktur dependensinya mengikuti Copula Clayton, Gumbel dan Frank, karena semuanya signifikan, SST Nino 4 dengan luas puso padi, yang signifikan adalah parameter Copula Gumbel, sedangkan SST Nino 3.4 dengan luas puso padi, parameter yang signifikan adalah parameter Copula Gumbel, Frank dan Normal. Tabel 18 memperlihatkan bahwa untuk anomali SST Nino 1+2 dengan luas puso padi di Ngawi, tidak ada parameter Copula yang signifikan, selain itu juga terjadi kasus yang sama dengan Lamongan dan Jember, parameter Copula Gumbel tidak dapat diestimasi. Untuk anomali SST Nino 3, SST Nino 4, SST Nino 3.4 dan luas puso padi di Ngawi, Copula Gumbel mempunyai parameter yang signifikan. Model fitting Copula dilakukan dengan memilih nilai loglikelihood terbesar untuk memilih Copula mana yang paling sesuai. Hasil model fitting Copula masing-masing Kabupaten disajikan pada Tabel 19.
12
Tabel 19. Model Fitting Copula dengan MLE (Bojonegoro) No. 1 2 3
Variabel SST Nino 1+2 dan Luas Puso Padi SST Nino 3 dan Luas Puso Padi SST Nino 4 dan Luas Puso Padi
4
SST Nino 3.4 dan Luas Puso Padi
Jenis Copula Gumbel Gumbel Gumbel Clayton Gumbel Frank Normal
Estimasi 1,0448 1,1675 1,0900 0,4682 1,1860 1,9779 0,3151
Log Likelihood 0,2302 2,5203 0,6346 2,2968 2,8609 3,7583 3,2318
P-Value 0,0000 0,0000 0,0000 0,0228 0,0000 0,0044 0,0022
0.0
0.2
0.4
0.6
0.8
1.0
1.0
1.0
0.8
0.8
0.6
0.6
v
v
0.4
0.4
0.2
0.2 0.0
0.2
u
(a)
0.4
(b)
0.6
0.8
1.0
0.0
0.0
0.0
0.0
0.2
0.2
0.4
0.4
v
v
0.6
0.6
0.8
0.8
1.0
1.0
Tabel 19 memperlihatkan bahwa struktur dependensi anomali SST Nino 1+2, SST Nino 3, SST Nino 4 dan luas puso padi di Bojonegoro mengikuti Copula Gumbel yang merupakan Copula Archimedean dengan tail dependensi pada bagian atas. Apabila anomali SST Nino 1+2, SST Nino 3 dan SST Nino 4 mengalami kenaikan, maka luas puso padi di Bojonegoro juga akan mengalami penambahan. Terdapat korelasi antara anomali SST Nino 1+2, SST Nino 3, SST Nino 4 dan luas puso padi pada titiktitik ekstrim atas. Struktur dependensi anomali SST Nino 3.4 dan luas puso padi di Bojonegoro mengikuti Copula Frank, yaitu Copula Archimedean yang tidak mempunyai tail dependensi. Masing-masing Copula terpilih di Bojonegoro yang digambarkan sesuai parameternya dalam bentuk scatterplot disajikan pada Gambar 9.
0.0
u
0.2
0.4
0.6
0.8
1.0
0.0
u
(c)
(d)
0.2
0.4
0.6
0.8
1.0
u
Gambar 9. Scatterplot rank Copula : Gumbel θ = 1,051201 (a) Gumbel θ = 1,191833 (b) Gumbel θ = 1,153196 (c) Frank θ = 1,758931 (d)
Model fitting Copula untuk Kabupaten Banyuwangi disajikan pada Tabel 20. Tabel 20. Model Fitting Copula dengan MLE (Banyuwangi) No.
Variabel
1
SST Nino 1+2 dan Luas Puso Padi
2
SST Nino 3 dan Luas Puso Padi
3
SST Nino 4 dan Luas Puso Padi
4
SST Nino 3.4 dan Luas Puso Padi
Jenis Copula Gumbel Frank Clayton Gumbel Frank Gumbel Gumbel Frank Normal
Estimasi
Log Likelihood
P-Value
1,2582 2,7042 0,4113 1,2856 2,8017 1,1089 1,3003 2,7975 0,4441
4,7976 4,3886 1,4919 5,5790 4,7160 0,7251 6,1100 4,4583 5,5282
0,0000 0,0005 0,0587 0,0000 0,0003 0,0000 0,0000 0,0004 6,02x10-7
Tabel 20 memperlihatkan bahwa Copula yang terpilih untuk mendeskripsikan struktur dependensi antara semua indikator ENSO (SST Nino 1+2, SST Nino 3, SST Nino 4, SST Nino 3.4) dan luas puso padi di Banyuwangi adalah Copula Gumbel, karena nilai loglikelihood pada Copula tersebut terbesar dan mempunyai parameter yang signifikan. Korelasi indikator ENSO dengan luas puso padi Jawa Timur lebih tinggi pada titik maksimum. Jika anomali SST Nino 1+2, SST Nino 3, SST Nino 4 dan SST Nino 3.4 tinggi, maka luas puso padi di Banyuwangi akan bertambah. Berdasarkan hasil dari estimasi Tau Kendall sebelumnya, model fitting Copula pada Kabupaten Ngawi dapat disajikan pada Tabel 21.
1
Tabel 21. Model Fitting Copula dengan MLE (Ngawi) Jenis Log Variabel Estimasi Copula Likelihood SST Nino 4 dan Luas Puso Padi Gumbel 1,0606 0,2933
2
SST Nino 3.4 dan Luas Puso Padi
No.
Gumbel
13
1,0044
0,0012
P-Value 0,0000 0,0000
Tabel 21 menunjukkan bahwa Copula Gumbel mewakili struktur dependensi antara anomali SST Nino, SST Nino 3.4 dan luas puso padi di Ngawi. Hal ini berarti apabila anomali SST Nino 4 dan SST Nino 3.4 naik, luas puso padi di Ngawi juga akan bertambah jumlahnya. 4.5
Peta Kerawanan Puso Padi Jawa Timur Hasil analisis untuk luas puso padi agregat di Jawa Timur, didapatkan bahwa Copula yang mewakili struktur dependensi antara indikator ENSO dan luas puso padi agregat di Jawa Timur adalah Copula Gumbel, yaitu Copula Archimedean yang memiliki tail dependensi pada bagian atas. Hal ini menunjukkan bahwa jika dianalisis secara agregat, semua daerah di Jawa Timur rawan terkena puso padi. Namun, dalam penelitian kali ini lebih difokuskan pada lima Kabupaten lumbung padi Jawa Timur. Berdasarkan hasil analisis di lima Kabupaten lumbung padi Jawa Timur, dapat disusun peta kerawanan puso yang disajikan pada Gambar 10. Peta Kerawanan Puso Padi
LA M O NG AN BOJ ON E GO R O NG AW I
Jatim_.shp Raw an Puso Padi Bukan Wilayah Penelitian Tidak Raw an Puso Padi JEM B ER BAN YU W AN GI
Gambar 10. Peta Kerawanan Puso Padi Kabupaten Lumbung Padi di Jawa Timur terhadap Indikator ENSO
Gambar 10 memperlihatkan bahwa dari lima Kabupaten lumbung padi di Jawa Timur, terdapat tiga Kabupaten yang rawan puso padi terhadap indikator ENSO, yaitu Bojonegoro, Banyuwangi dan Ngawi (warna merah pada peta). Hal tersebut ditunjukkan oleh analisis menggunakan Copula sebelumnya, yang menghasilkan bahwa struktur dependensi indikator ENSO dan luas puso padi di Bojonegoro, Banyuwangi dan Ngawi mengikuti Copula tertentu, maka mengindikasikan adanya hubungan antara indikator ENSO dan luas puso padi di tiga Kabupaten tersebut. Lamongan dan Jember tidak rawan puso padi terhadap indikator ENSO (warna hijau pada peta), sesuai dengan analisis Copula sebelumnya yang menunjukkan bahwa struktur dependensi indikator ENSO dan luas puso padi di Lamongan dan Jember tidak mengikuti Copula manapun, sehingga tidak ada hubungan antara indikator ENSO dan luas puso padi di dua Kabupaten lumbung padi tersebut. 5.
KESIMPULAN DAN SARAN Hasil analisis korelasi pearson tidak bisa mengidentifikasi pola hubungan antara indikator ENSO dan luas puso padi agregat di Jawa Timur maupun luas puso padi di lima kabupaten lumbung padi Jawa Timur. Struktur dependensi anomali SST Nino 1+2, SST Nino 3, SST Nino 4, SST Nino 3.4 dan luas puso padi agregat di Jawa Timur dan di Kabupaten Banyuwangi mengikuti Copula Gumbel. Pada Bojonegoro, anomali SST Nino 1+2, SST Nino 3, SST Nino 4 dan luas puso padi mengikuti Copula Gumbel, sedangkan SST Nino 3.4 mengikuti Copula Frank. Pada Kabupaten Lamongan dan Jember tidak mengikuti Copula manapun. Pada Kabupaten Ngawi, indikator ENSO yang berpengaruh pada luas puso padi adalah anomali SST Nino 4 dan SST Nino 3.4, dan struktur dependensinya mengikuti Copula Gumbel. Berdasarkan hasil penelitian, menunjukkan bahwa SST Nino yang selalu berpengaruh terhadap luas puso padi adalah SST Nino 3.4. Hal ini menunjukkan bahwa SST Nino 3.4 merupakan SST Nino yang paling berpengaruh pada wilayah Indonesia. Peta kerawanan puso padi disusun berdasarkan hasil analisis Copula. Kabupaten Bojonegoro, Banyuwangi dan Ngawi merupakan Kabupaten yang rawan puso padi. dan Kabupaten Lamongan dan Jember tidak rawan puso padi. Saran yang dapat diberikan adalah analisis menggunakan Copula lebih baik lagi apabila dilakukan dengan jumlah sampel yang besar agar bentuk scatterplot rank Copula lebih terlihat jelas. Oleh karena itu, untuk penelitian selanjutnya disarankan untuk menggunakan data dengan jumlah n lebih besar agar hasil yang didapatkan lebih baik lagi dan dapat dilakukan di lokasi lain yang berbeda. 14
6. DAFTAR PUSTAKA Daniel, W. W. (1989). Statistika Nonparametrik Terapan. George State University. Jakarta : Gramedia. Dorey, M. dan Joubert, P. (2007). Modelling Copulas : An Overview. University of Grenwich and University of Witwatersrand. Embrechts, P., Lindskog, F., dan McNeil, A. (2001). Modelling Dependence with Copulas and Application to Risk Management. Department of Mathematics, ETHZ CH-8092 Zürich, Switzerland. Embrechts, P. (2009). Copulas : A personal View. Department of Mathematics, ETH Zürich, Switzerland. Genest, C. dan Nešlehová, J. (2007). Modelling Count Data With Copula : Should We ?. Tartu. Genest, C. dan Nešlehová, J. (2010). Copulas : Introduction to the Theory and Implementation in R with Applications in Finance and Insurance. Universite Laval and McGill University. Irawan, B. (2006). Fenomena Anomali Iklim El-Nino dan La-Nina : Kecenderungan Jangka Panjang dan Pengaruhnya terhadap Produksi Pangan. Pusat Analisis dan Kebijakan Pertanian Bogor. Kpanzou, T.A. (2007). Copulas in Statistics. African Institute for Mathematical Sciences (AIMS). Mikosch, T. (2008). Copulas: Tales and Facts. Laboratory of Actuarial Mathematics, University of Copenhagen, Universitetsparken 5, DK-2100, Denmark. Naga, D.S. (2005). Statistika Terapan Bab 13A Nonparametrik : Data Peringkat I,
, diakses pada tanggal 23 Oktober 2011, pukul 18:35 WIB. Naylor, R.L., Falcon, W., Wada, N., dan Rochberg, D. (2001). Using El Niño -Southern Oscillation Climate Data To Predict Rice Production In Indonesia. Climatic Change 50, 255–265. Naylor, R.L., Falcon W., Wada N., dan Rochberg, D. (2002). Using El Niño-Southern Oscillation Climate Data To Improve Food Policy Planning In Indonesia. Bulletin of Indonesian Economic Studies 38: 75-91. Nelsen, R. B. (2005). An Introduction to Copulas : Second Edition. New York : Springer. Salvadori, G., Michele, C.D., Passoni, G., dan Vezzoli, R. (2007). A Multivariate Model of Sea Storms Using Copulas. Coastal Engineering, Volume 54, Issue 10, Pages 734-751. Samian. (2008). Korelasi Sederhana Pearson, , diakses pada 29 September 2011, pukul 19:53 WIB. Schӧlzel, C. dan Friederichs, P. (2008). Multivariate Non-Normally Distributed Random Variables In Climate Research – Introduction to The Copula Approach. Nonlin. Processes Geophys., 15, 761–772. Serinaldi, F. (2009). Comment on "Multivariate Non-Normally Distributed Random Variables In Climate Research – Introduction to The Copula Approach" by C. Schӧlzel and P. Friederichs. Nonlin. Processes Geophys., 16, 331–332. Setiawan, A. (2005). Oseanografi, Awal Kehidupan Bermula di Laut : ENSO, , diakses pada 27 September 2011, pukul 20:29 WIB. Sklar, A. (1959), Fonctions de répartition à n dimensions et leurs marges. Publ. Inst. Statist. Univ. Paris 8: 229–231 Syahrir, I. (2011). Estimasi Parameter Copula dan Aplikasinya pada Klimatologi. Thesis Jurusan Statistika FMIPA ITS. Surabaya : ITS. Villarini, G., Serinaldi, F., dan Krajewski, W. F. (2008). Modeling Radar-Rainfall Estimation Uncertainties Using Parametric and Non-Parametric Approaches. Advances in Water Resources, Volume 31, Issue 12, Pages 1674-1686. Yan, J. (2006). Enjoy The Joy Copula. Department of Statistics and Actuarial Science, University of Iowa, 241 Schaeffer Hall, Iowa City, IA 52242, USA. Yan, J. (2007). Enjoy The Joy of Copulas with a Package Copula. American Statistical Association.
15