## PEMODELAN KECEPATAN ANGIN RATA-RATA DI SUMENEP MENGGUNAKAN MIXTURE OF ANFIS
1 Syarifah Diana Permai, 2 Nur Iriawan, 3 Irhamah 1,2,3 Jurusan  Statistika, Fakultas Matematika dan Ilmu Pengetahuan Alam, Institut Teknologi Sepuluh November, Surabaya
Alamat e-mail : 1 syarifah.diana@ymail.com

## ABSTRAK
Angin memiliki dampak positif dan dampak negatif. Dampak positif dari angin yaitu memperlancar aktivitas pelayaran, membantu irigasi menggunakan kincir angin, pembangkit tenaga listrik dan lain sebagainya. Namun perubahan cuaca yang ekstrim akhir-akhir ini dapat menimbulkan angin kencang serta gelombang laut yang tinggi, sehingga menghambat pelayaran. Salah satu kabupaten yang mengalami dampak negatif ini adalah Sumenep, daerah pesisir di Pulau Madura. Oleh karena itu diperlukan pemodelan kecepatan angin rata-rata di Sumenep yang akurat. Dua metode yang diterapkan adalah Adaptive Neuro Fuzzy Inference System (ANFIS) dan Mixture of ANFIS. Mixture of ANFIS dilakukan melalui beberapa pembagian kelompok, yaitu dua, tiga, empat, lima dan enam kelompok. Evaluasi perbandingan kebaikan model dilakukan berdasarkan kriteria Root Mean Square Error (RMSE). Hasil evaluasi menunjukkan bahwa pemodelan kecepatan angin rata-rata menggunakan mixture of ANFIS dengan enam kelompok memiliki RMSE in sample dan out-sample lebih kecil daripada jumlah kelompok yang lain. Mixture of ANFIS memodelkan kecepatan angin rata-rata di Sumenep lebih baik dari ANFIS karena menghasilkan RMSE in dan out sample yang lebih kecil dari ANFIS.
Kata kunci : Kecepatan Angin, ANFIS, Mixture of ANFIS

## PENDAHULUAN
Cuaca yang terjadi di Indonesia dalam satu tahun terakhir ini sangatlah fluktuatif. Dampak negatif dari perubahan cuaca yang tak menentu ini adalah terjadinya musim hujan dengan curah hujan di atas normal, timbulnya angin kencang bahkan berpotensi terjadinya angin puting beliung, disertai hujan lebat dan petir serta gelombang laut yang tinggi [2]. Salah satu daerah yang terkena dampak negatif adalah Kabupaten Sumenep. Semenjak September 2010 terjadi musim hujan di Kabupaten Sumenep, arah angin dari arah timur ke barat dengan kecepatan
antara 5 hingga 35 kilometer per-jam, dan ini terjadi hingga bulan Maret bahkan April 2011. Akibat cuaca ekstrim yang terjadi di perairan Sumenep, membuat aktivitas pelayaran dari pelabuhan Kalianget ke Jangkar (Situbondo) menjadi lumpuh pada awal Januari 2011 [7]. Bahkan gelombang laut di perairan Sumenep bisa mencapai tinggi antara 1,75 sampai 6 meter, kecepatan angin 50 hingga 60 kilometer/jam, menyebabkan tingginya gelombang laut. Sehingga pelayaran dari pelabuhan Kalianget sempat ditunda beberapa waktu. Oleh karena itu, perlu adanya prediksi yang memiliki keakurasian tinggi terhadap kecepatan
angin agar aktivitas manusia seperti pelayaran, penerbangan, penggerak kincir angin dan aktivitas-aktivitas lain yang bergantung pada kecepatan angin menjadi lancar. Karena ketersediaan angin di alam tidak dapat disimpan dan tidak dapat dengan mudah diatur agar sesuai dengan kebutuhan manusia [6].
ANFIS merupakan singkatan dari Adaptive Neuro-Fuzzy Inference System atau Adaptive Network Based Fuzzy Inference System . ANFIS diartikan sebagai sistem inferensi fuzzy berbasis jaringan saraf adaptif. Parameter pada ANFIS, adalah fungsi keanggotaan premis dan konsekuensi. Proses dalam ANFIS adalah pengubahan parameter fungsi keanggotaan masukan dan keluaran menggunakan algoritma perambatan balik atau algoritma hybrid . Algoritma hybrid yaitu menggabungkan antara algoritma alur mundur dan algoritma alur maju [8]. Algoritma alur maju menggunakan metode kuadrat terkecil ( Least Squares Estimate ),
sedangkan algoritma alur mundur menggunakan
back-propogation .
Metode kuadrat terkecil digunakan untuk menentukan parameter konsekuensi, sedangkan langkah mundur digunakan untuk memperbaharui bobot premis. Arsitektur ANFIS model Sugeno order satu dapat dilihat pada Gambar 1 [5].
Metode Mixture of ANFIS diterapkan oleh [1] pada data CPU load . Hasil penelitian tersebut menunjukkan bahwa metode Mixture of ANFIS memprediksi CPU load lebih baik daripada Model AR dan Exponential Smoothing . Oleh karena itu tujuan penelitian ini adalah untuk menerapkan metode Mixture of ANFIS pada data kecepatan angin rata-rata di Sumenep. Selanjutnya melakukan perbandingan hasil peramalan metode ANFIS dengan Mixture of ANFIS.
Misalkan terdapat dua input, yaitu x 1
dan x 2 serta satu output , maka model fuzzy satu order sebagai berikut.
Rule  1:  If  x1  is  A1 and x2 is B1, then f1=p1 x1 + q1 x2 + r1 Rule  2:  If  x1  is  A2 and x2 is B2, then f2=p2 x1 + q2 x2 + r2 Rule  3:  If  x1  is  A1 and x2 is B2, then f3=p3 x1 + q3 x2 + r3 Rule  4:  If  x1  is  A2 and x2 is B1, then f4=p4 x1 + q4 x2 + r4 Secara umum apabila terdapat n input dan satu output , maka model ANFIS n+1 dimensi. Sebelum layer 1 yang terdapat pada Gambar 1 menunjukkan input yang digunakan dalam ANFIS yaitu dua input x 1 dan x 2 . Layer 1 menunjukkan proses fuzzifikasi dimana setiap simpul i dalam layer ini merupakan simpul aditif dengan fungsi berikut.
O 1,i = μ Ai (x 1 ) untuk i =1,2                (1) O 1,i = μ Bi-2 (x 2 ) untuk i =3,4              (2) O 1,i merupakan tingkat keanggotaan atau derajat keanggotaan dari himpunan fuzzy A ( A 1 , A 2 , B 1 , B 2 ) dengan x 1 dan x 2 adalah input untuk simpul ke- i . A i dan B i-2 adalah linguistic label . Sedangkan μ A merupakan fungsi keanggotaan untuk A yang sesuai dengan parameter fungsi keanggotaan untuk fungsi Gaussian, yaitu sebagai berikut.
µ 1
1 3 dimana a i , b i dan c i adalah parameter premis. Layer 2 atau rule layer, setiap simpul dalam layer ini merupakan fixed yang diberi label Π dan merepresentasikan setiap rule, sebagai berikut. , µ µ (4) , µ µ (5)
, µ µ (6)
, µ µ (7)
Pada layer 3 atau normalization layer dengan nilai simpul ke- i dalam layer ini merupakan ratio dari rule layer ke- i
dengan penjumlahan semua rule seperti rumus berikut.
, , 1,2,3,4 (8) Layer defuzzifikasi atau layer 4 mempunyai output dari setiap simpul yang disebut dengan nilai weighted consequent dengan rumus berikut. , (9)
dimana i = 1, 2, 3, 4 dan   merupakan output pada layer 3 sedangkan { p i , q i ,
r i } adalah parameter dari simpul ini yang disebut dengan consequent parameters . O verall output dari layer 5 yang disebut dengan summation layer yang merupakan fixed simpul dengan simbol Σ . Output pada layer ini merupakan penjumlahan seluruh output pada layer 4 yang dirumuskan berikut.
, ∑ ∑ ∑ , 1,2,3,4 (10)
Gambar 1 ANFIS Model Sugeno Order Satu Teknik optimasi menggunakan algoritma pembelajaran pada data training untuk mendapatkan parameter yang paling baik. Algoritma pembelajaran untuk ANFIS adalah algoritma hybrid yang meminimumkan error [10]. Satu input dan output pada data yang diolah dengan ANFIS, diiterasi menggunakan alur maju mulai layer satu sampai dengan layer terakhir dan melalui LSE didapatkan consequent parameter . Sedangkan perhitungan error didapatkan dari alur mundur sehingga diperoleh parameter premis atau parameter non linier [6].
Metode Mixture of ANFIS diterapkan oleh [1],  pada prediksi CPU load . Nilai prediksi dari metode Mixture of ANFIS mengkombinasikan antara mixture dari nilai prediksi ANFIS pada tiap cluster dengan probabilitas state
transition antar cluster dalam waktu ke- t . Probabilitas cluster state dari cluster j’ ke cluster j dimana j’ = 1, 2, ..., j didefinisikan sebagai berikut.
′ ′ ′ ′ 11
dimana probabilitas a j’j jika j ≠ j’ didefinisikan dengan teorema Bayes seperti berikut.
′ ′ ′ (12) menghitung output yang merupakan hasil prediksi Mixture of ANFIS melalui rumus berikut.
′
13 D imana
merupakan prediksi cluster ke- j .
N 1 N 2 N 3 N 4 A 2 B 1 B 2 Π 1 Π 2 Π 3 Π 4 ∑ A 1 2 3 4 1 x 1 x 2 Layer 1 Layer 3
Layer 2
Layer 4
Layer 5
Kajian mengenai kecepatan angin rata-rata pernah dilakukan oleh [5] terhadap kecepatan angin di Juanda dan
Sumenep yang menunjukkan bahwa metode AI ( Arificial Intelligence ) seperti Algoritma Genetika untuk identifikasi ARIMA menghasilkan MSE lebih kecil daripada correlogram . Penelitian kecepatan angin rata-rata di Sumenep menggunakan metode Adaptive Neuro Fuzzy Inference System (ANFIS) yang diperoleh kesimpulan bahwa peramalan menggunakan ANFIS menghasilkan akurasi yang lebih tinggi daripada ARIMA [3]. Hasil penelitian [1] menunjukkan bahwa metode Mixture of ANFIS memprediksi CPU load lebih baik daripada Model AR dan Exponential Smoothing . Oleh karena itu tujuan penelitian ini adalah untuk menerapkan metode Mixture of ANFIS pada data kecepatan angin rata-rata di Sumenep. Selanjutnya melakukan perbandingan hasil peramalan metode ANFIS dengan Mixture of ANFIS.

## METODOLOGI PENELITIAN
Sumber Data dan Variabel Penelitian Penelitian ini menggunakan data sekunder yaitu data kecepatan angin yang diukur oleh Badan Meteorologi dan Geofisika Kalianget, Sumenep dari bulan Januari 2009 sampai dengan Desember 2010, sehingga total terdapat 730 data. Data training atau in sample yang digunakan yaitu data kecepatan angin rata-rata dari bulan Januari 2009 sampai dengan November 2010 sehingga terdapat 699 data in sample . Sedangkan data testing atau out sample adalah data bulan Desember 2010, sehingga terdapat 31 data.

## Metode Analisis
Tahapan analisis pada penelitian ini sebelum dilakukan analisis menggunakan Metode ANFIS terlebih dahulu data dianalisis dengan
menggunakan ARIMA. Input pada metode ANFIS berasal dari AR dari model terbaik ARIMA. Pada analisis ANFIS di penelitian ini ditentukan sebanyak dua membership funtion (mf) dengan tipe Gaussian sehingga diperoleh nilai RMSE in sample dan out of sample [9]. Selanjutnya dilakukan analisis dengan metode Mixture of ANFIS, dimana pada penelitian ini dilakukan mixture untuk dua cluster sampai dengan enam cluster sehingga akan diperoleh model terbaik Mixture of ANFIS dari nilai RMSE in sample dan out of sample terkecil.

## HASIL DAN PEMBAHASAN
Sebelum dilakukan pemodelan data, maka terlebih dahulu dilakukan pemeriksaan data secara deskriptif untuk mengetahui karakteristik data. Kemudian dilakukan analisis dengan ARIMA untuk menentukan input pada ANFIS. Selanjutnya dilakukan pemodelan data kecepatan angin di Sumenep menggunakan ANFIS dan Mixture of ANFIS. Hasil pemodelan dan peramalan pada kedua metode akan dibandingkan untuk memperoleh pemodelan dan peramalan terbaik.

## Analisis Deskriptif
Hasil analisis statistika deskriptif dari kecepatan angin rata-rata di Sumenep pada bulan Januari 2009 sampai Desember 2010 dapat dilihat pada Tabel 1.
Kecepatan angin rata-rata di Kalianget Sumenep selama dua tahun yaitu 2009 dan 2010 memiliki karakteristik yang berbeda setiap bulannya. Rata-rata dari kecepatan angin rata-rata setiap hari di Kalianget Sumenep pada tahun 2009 dan 2010 yang paling tinggi yaitu Bulan Agustus sebesar 8,355 knot dan 6,645 knot.
Tabel 1 Statistika Deskriptif Kecepatan Angin Rata-Rata Harian di Sumenep
Tahun Bulan Mean Varian Min Max
2009
Jan 7,097 10,49 2  15 Feb 7,143 24,349 2  22 Mar 2,093 1,09 1  5 Apr 3,867 4,189 1  7 Mei 4,387 3,512 2  9 Jun 6,467 2,12 4  10 Jul 7,742 4,398 4 12 Ags 8,355 2,437 5  13 Sep 7,1 3,472 4 11 Okt 7,258 3,798 3  12 Nov 5,767 9,702 2  11 Des 3,613 4,045 1 13 2010 Jan 6,065 14,262 2  15 Feb 2,929 0,884 2  5 Mar 2,871 0,783 2  5 Apr 2,533 0,74 1  5 Mei 3,871 3,583 1  8 Jun 4,8 3,752 1 8 Jul 5,935 6,262 2 11 Ags 6,645 2,903 4  10 Sep 5,233 4,668 1  9 Okt 4,516 4,125 2  10 Nov 3,8 1,2 2 6 Des 5,452 6,656 1  14
Pada tahun 2009, bulan Juni sampai Oktober memiliki rata-rata kecepatan angin rata-rata yang cukup tinggi, terlihat dari nilai mean yang bernilai antara 6 knot sampai 8 knot. Bulan Februari 2009 memiliki varians paling besar yaitu 24,349 yang di-akibatkan karena nilai minimum kecepatan angin rata-rata sebesar 2 knot, sedangkan nilai maksimumnya sebesar 22 knot. Kecepatan angin rata-rata yang cukup rendah pada tahun 2009 terjadi di bulan Maret sampai Mei dengan rata-rata sebesar 2 knot sampai 5 knot.
Tahun 2010 memiliki kecepatan angin rata-rata berkisar antara 2 knot sampai 7 knot. Maksimum kecepatan angin rata-rata sebesar 15 knot terjadi pada bulan Januari sehingga keragaman data pada bulan itu cukup tinggi yaitu sebesar 14,262. Kecepatan angin rata- rata minimum sebesar 2 knot terjadi bulan April. Pada bulan Februari sampai April kecepatan angin yang berhembus cukup stabil karena nilai varian kurang
dari satu dan rata-ratanya kurang lebih 2 knot dengan maksimum kecepatan angin rata-rata sebesar 5 knot. Perbedaan karakteristik kecepatan angin rata-rata tahun 2009 dan 2010 menunjukkan cuaca yang anomali.

## Pemodelan ANFIS
Sebelum analisis menggunakan ANFIS, maka terlebih dahulu data dianalisis dengan ARIMA untuk menentukan input pada ANFIS. Identifikasi stasioneritas data dalam varian menggunakan Box-Cox plot dimana dipeorleh nilai rounded value sebesar 0,5. Hal ini mengartikan bahwa data kecepatan angin rata-rata di Kalianget Sumenep tidak stasioner terhadap varian, sehingga data ditransformasi akar. Selanjutnya identifikasi stasioneritas terhadap mean dengan menggunakan ACF plot menunjukkan bahwa nilai ACF turun sedikit demi sedikit atau turun lambat sampai dengan lag ke- 30. Hal ini menunjukkan bahwa data transformasi tidak stasioner terhadap mean, maka data yang telah ditransformasi tesebut perlu di- differencing sebanyak satu kali. Setelah itu data hasil transformasi dan differencing diidentifikasi modelnya meng-gunakan ACF dan PACF plot . Hasil ACF dan PACF plot disajikan pada Gambar 1 dan Gambar 2.
Berdasarkan plot ACF dan PACF maka dipeorleh beberapa model dugaan untuk ARIMA. Selanjutnya dilakukan pemeriksaan signifikansi parameter terhadap semua kemungkinan model ARIMA dan pemeriksaan white noise serta normal.
Hasil pengujian signifikansi parameter terhadap semua kemungkinan model ARIMA dan pemeriksaan white noise serta normal diperoleh model ARIMA (2,1,1) dan ARIMA ([1,3],1,1), sedangkan yang memenuhi semua asumsi adalah ARIMA ([1,2,4],1,1).
AIC 658,
SBC pada
RMS 1,800 terse ([1,2 kemu O pemo Hasil RMS
meng dan ke-3. mem adala funct terda yaitu berik linier dan b
Gambar 2 A
Gambar 3 P
yang pal 1931 pada A minimum s a model AR SE in samp 0534 dim ebut terdapa 2,4],1,1), se ungkinan in Oleh kare odelan terh l analisis SE paling ggunakan  i differencing . (Z * t-1 d bership fun ah dua bua tion nya  y apat dua pa u mean (c) d kut nilai pa r dan linier backpropog
ACF Plot Data
PACF Plot Da ling minim
ARIMA([1 sebesar 676
RIMA(2,1,1 ple paling  k mana nilai at pada mo ehingga terd nput AR pad ena itu, hadap beb ANFIS  m kecil  p input data g pada lag dan Z * t-3 ). nction yan ah dan tipe yaitu Gaus arameter pa dan standar arameter-pa r setelah al gation .
a Differencing
ata Differencin
mum sebes ,2,4],1,1) d 6,224 terdap 1) sedangk kecil bernil i minimu odel ARIM dapat bany da ANFIS. dilakuk berapa inpu menunjukk pada ANF transforma ke-1 dan l Banyakny ng digunak e membersh ssian, mak ada Gaussia rt deviasi ( σ arameter no lgoritma LS
g ng
sar an pat an lai um MA ak an ut. an IS asi ag ya an hip ka an, σ ). on SE
A di an m tr Z ep pa se se
P da N m en pe (F m m A m M
pe m st
Statistika,
Tabel 2 P I MF Par
1 2 Tabel 3 I Aturan 1 -0 2 -0 3 - 4 2 Model inyatakan s + +
Berdasark ngin rata-ra menggunaka ransformasi Z * t-3 diperole poch seban ada data in edangkan o ebesar 3,713
Pemodelan M Data kec ari bulan Ja November menjadi dua nam kelom engelompok FCM). Sete masing-masi menggunaka
ANFIS tiap menjadi satu
Mixture of A
Proses embagian masing-masi tasioneritas,
Vol. 1, No. 2
arameter Non Input Z * t-1 dan
rameter
c - σ 0 c 2 σ 0
3 Parameter L Input Z * t-1 dan p 0,4953 -0 0,1175 -0, -3,243 1 2,648 3 ANFIS  y eperti berik kan hasil pr ata di Kali an ANFIS d dan differ eh nilai erro nyak 170. n sample s ut sample d 321.
Mixture of epatan ang anuari 2009 2010 a, tiga, em mpok mengg kkan Fu elah dikelo ing kelom an ANFIS. p kelompo u, inilah y ANFIS. dilakukan kelompok ing kelomp
, identif
2, November
n Linier ANFI n Z * t-3 Input Z * t-1 Z 1,779 -1 0,7677 1 2,543 2 0,7261 1 inier ANFIS n Z * t-3 q 0,1376 -0, ,03373 0, 1,743 5 3,862 -7 yang terb kut. rediksi kece ianget Sum dengan inpu rencing Z * t- or minimum RMSE res sebesar 3,8 diperoleh R
f ANFIS
gin di Sum 9 sampai de dikelompo mpat, lima gunakan  m uzzy C-M mpokkan, mpok dian Kemudian, ok digabun yang dinam dimulai dengan  F ok diidenti fikasi m
S
Z * t-3 1,404 ,489 ,336 ,811 r ,1791 1639 5,71 7,575 entuk
+ epatan menep t data
1 dan m pada sidual 81142 RMSE menep engan okkan a dan metode Means maka nalisis hasil ngkan makan dari FCM, fikasi model
ARIM terba predi perhi
dan p H angin deng Tabe D kemu Mixtu anali clust meto yaitu seban kecep data. ANF a. K Se
data meng mene ident 1 me
cox d telah
mean
ARIM MA, pemi aik, penen iksi ANFIS itungan pro prediksi mix Hasil pengelo n rata-rata gan dua ke el 4 dan Tab Tabel 4 Uru wa 4 . 6 6 6 Tabel 5 Uru wa . 6 6 6
Data pada udian dia ure of AN isis Mixture ter. Pembag ode FCM m u kecepatan nyak 411 d patan angi Selanjutn
FIS pada ma Kelompok 1 ebelum dia kelomp ggunakan entukan inp tifikasi stasi enggunakan dan ACF m h stasioner n. Selanjut MA mengg
ilihan mod ntuan inp
S. Selanjutn obabilitas sta xture of AN ompokan da a menggun elompok ad bel 5. Data Kelomp utan aktu Data 2 2 3 3 4 3 ... ... 97 3 98 4 99 4 Data Kelomp utan aktu Data 1 7 8 6 9 8 ... ... 79 6 80 6 81 6 Tabel 4 d analisis m NFIS. Berik e of ANFIS
gian dua cl menghasilkan n angin ren data. Kelom in tinggi t nya dilaku asing-masin
analisis den pok 1 ARIMA put pada A ioneritas pa n time seri menunjukkan terhadap tnya identi gunakan AC
del ARIM
put ANFI nya dilakuk ate transitio
NFIS. ata kecepat nakan FC dalah sepe pok 1 a pok 2 a dan Tabel menggunak kut ini ha S dengan du luster deng n kelompok ndah terdap mpok 2 yai terdapat 28 ukan analis ng kelompok ngan ANFI dianalis A untu ANFIS. Ha ada kelompo es plot , bo n bahwa da varians d fikasi mod CF dan PAC
MA IS, an on an M erti
5 an sil ua an k 1 pat itu 88 sis k.
IS, sis uk sil ok ox- ata an del CF pl 4 nu ad S si
B
di an di pa si no
V W *N va di T se ke 1)
lot yang dis . Gamba Gambar Hasil iden unjukkan pl dalah lag k edangkan  p ignifikan  a Berdasarkan ilakukan  p nalisis iperoleh arameter ignifikan da ormal seper Tabel
Model Konstanta φ 1 θ 1
Var estimate
Std error AIC SBC
White Noise
Normal Nilai yang be alue Berdasark iperoleh mo
Tabel 6, m ebagai inpu ecepatan an ). Jika men
Vol. 1, No. 2
sajikan pad r 4 ACF Plot r 5 PACF Plot
ntifikasi mo lot ACF lag ke-1, 2, 3, 4 pada plot P
adalah lag hasil  t pemodelan menggunak ARIMA( pada m an white no rti Tabel 6. l 6 ARIMA K AR 3,316 0,52839 0,77996 1,033348 1,016537 1182,839 1194,895 Ya Tidak erada di dalam kan anal odel ARIM maka AR y ut dalam AN ngin rata-ra nggunakan
da Gambar 3 Kelompok 1 t Kelompok 1 odel ARIMA g-lag yang k 4 dan lag k PACF lag g ke-1 da tersebut, ARIMA. kan AR (1,0,1), model ter oise, tetapi Kelompok 1 RIMA(1,0,1) (< 0,0 (< 0,0 (< 0,0 8 7 9 5 (0,03
m kurung ada lisis AR A (1,0,1) se yang digun NFIS adalah ata lag ke-1 MF type G r 2013
3 dan
A me- keluar ke-15. yang an 2. maka Hasil RIMA yang rsebut tidak 0001) 0001)
0001)
382) alah p- RIMA, eperti nakan h data 1 (Zt- Gauss
maka mean anali meng linier Tab M T At A meng param linier sebag
b. K B meng cox p
kecep varia
1, se mem - 1,6
harus mode Gam
PAC
a terdapat n (c) dan st isis ANFIS ghasilkan  p r yang disaj
bel 7 Paramet MF 1
2
Tabel 8 Param turan 1 2 Analisis AN g-hasilkan meter non r. Model A
gai berikut.
Kelompok 2 erdasarkan ggunakan ti pada kelomp patan angin ans, karena elain itu nila muat angka 1
62 dan -0,5 s ditransf el ARIMA mbar 5 dan
CF kelompok
Gambar 5 A
dua para tandart devi S dengan parameter  n jikan di Tab
ter Non Linier Parameter C Σ C Σ meter Linier K Param P 0,2191 -8,985 NFIS pada masing-ma linier da ANFIS kelom identifikasi ime series p pok 2 menu n tidak stasio rounded va ai lower dan
1 dimana ni
2. Oleh ka formasi .
A dapat 6 yaitu pl k 2.
ACF Plot Kel
ameter yai iasi ( σ ). Ha input Zt non linier d bel 7 dan 8.
r Kelompok 1 Input Z t-1 1,191 1,571 5,619 0,5147 Kelompok 1 meter q 2,457 48,78 kelompok asing emp an paramet mpok 1 yai
i stasionerit plot dan bo unjukkan da oner terhad alue sebesar n upper tid ilainya anta arena itu da Identifika dilihat pa lot ACF d lompok 2
itu sil t-1 an
1 pat ter itu
tas ox- ata ap r - ak ara ata asi da an ad S la ka ya in si ([ A se pe pa T M V St A SB W N N R m ya A ya A la da
da
Statistika, Gambar
Pada plot dalah lag edangkan ag-lag ke-1, arena itu ad ang nantin nput pada ignifikansi [6],0,4),  A ARIMA ( emua param emilihan  m ada kriteri Tabel 9. Tabel Model A ( Variance 0 td error 0 AIC - BC - White Noise Normal Ya MSE 0 Berdasark model terbai ang pali ARIMA([6], ang digun ANFIS adal ag ke-6 atau an 11 adala an linier yan Tabel 10 Par MF 1 2
Vol. 1, No. 2 r 6 PACF Plot
t ACF lag ke-1, 2, plot PACF 6 dan 18 y da beberapa nya dapat ANFIS.  B parameter ARIMA [1,6,8],0,0) meter signi model ARI ia model l 9 ARIMA K ARIMA ([6],0,4) 0,00054 ,023126 1345,96 1323,98 Ya a (>0,15) Y 0,38851 kan krite ik, didapatk ing min 0,4). Oleh akan sebag lah data ha u Xt-6 . Beri ah hasil para ng optimum
rameter Non L
Parameter C Σ C Σ
t Kelompok 2 yang signi 3 dan F menunju ang cut off . a model AR dipilih me Berdasarkan model AR ([1,6],0,0) menunju ifikan, seh IMA didas terbaik se Kelompok 2 ARIMA ([1,6],0,0) 0,000542 0,023288 -1344,99 -1334 Ya Ya (>0,15) 0,39325 eria pemi kan nilai R imum a karena itu, gai input asil transfor ikut ini Tab ameter non
m.
Linier Kelomp r Input 0,038 0,07 0,14 0,054
2 ifikan ke-4. ukkan Oleh RIMA enjadi n uji RIMA dan ukkan ingga arkan eperti ARIMA ([1,6,18], 0,0) 0,000535 0,023129 -1347,61 -1332,96 Tidak Ya (>0,15) 0,39011 ilihan RMSE adalah , data pada rmasi bel 10 linier pok 2 t X t-6 863 757 409 444
Tabel 11 Parameter Linier Kelompok 2 Aturan Parameter P Q 1 -2,597 0,2724 2 -0,1287 0,2414 Satu input dan dua mebership function menghasilkan dua aturan yang dengan empat parameter linier seperti pada Tabel 11. Model ANFIS untuk kelompok 2 dapat ditulis seperti persamaan berikut. 2,597 0,272 0,129 0,241
c. Mixture of ANFIS Dua Kelompok
Pada saat pengelompokan data kecepatan angin rata-rata di Sumenep menggunakan FCM membagi data menjadi dua kelompok. Namun ketika data prediksi masing-masing kelompok di- mix , maka data dikembalikan berdasarkan urutan waktu yang sesuai. Inilah proses mixture pada dua kelompok. Langkah selanjutnya, setelah data kedua kelompok di- mix , maka dilakukan perhitungan probabilitas state transition untuk mendapatkan probabilitas cluster state sehingga diperoleh prediksi dari mixture of ANFIS. Berikut ini Tabel 12 menyajikan hasil perhitungan probabilitas cluster state.
Tabel 12 Probabilitas Cluster State Dua
Kelompok P(CS t = 1) 0,70712 P(CS t = 2) 0,46854
Berdasarkan Tabel 12 dapat dilihat bobot pada kelompok 1 sebesar 0,70712 dan bobot untuk kelompok 2 adalah 0,46856. Selanjutnya diperoleh nilai prediksi mixture of ANFIS dengan cara mengalikan nilai prediksi masing- masing kelompok dengan nilai bobot masing-masing kelompok.
Tabel 13 Prediksi Kecepatan Angin Rata-Rata Menggunakan Mixture of ANFIS Dua Kelompok Urutan Waktu Kelompok Prediksi Mixture of ANFIS 1 2  3,626257697 2 1  2,341284377 3 1  2,024176253 4 1  2,231855649 5 1  2,231855649 ... ... ... 695 1 2,56135089 696 1  2,024176253 697 1  2,231855649 698 1  2,231855649 699 1 2,56135089
Setelah diperoleh nilai prediksi mixture of ANFIS sepeti pada Tabel 13, selanjutnya dihitung residual hasil mixture of ANFIS dan diperoleh Nilai RMSE in sample pada residual mixture dua kelompok sebesar 3,3271.
d. Perbandingan Mixture of ANFIS dengan 2, 3, 4, 5 dan 6 Cluster Analisis Mixture of ANFIS dengan 3, 4, 5 dan 6 cluster dilakukan dengan cara yang sama seperti mixter 2 cluster tersebut. Berikut ini hasil rangkuman deskriptif dari residual dan RMSE pada tiap mixture dengan cluster yang berbeda-beda.
Tabel 14 Statistika Deskriptif Residual Mixture Cluster Mean  Min  Maks RMSE
2 2,4041 - 1,7157 18,494 3,3271 3 1,9054 -1,358 16,246 3,1313 4 1,9283 - 0,5742 15,469 2,9924 5 1,5063 - 0,6035 20,86 2,5258 6 1,08 - 1,4059 20,272 2,3883
Berdasarkan perbandingan residual yang dihasilkan dari metode mixture of ANFIS dengan 2, 3, 4, 5 dan 6 cluster seperti Tabel 14 menunjukkan bahwa rata-rata residual paling mendekati nol adalah mixture enam cluster. Selain itu dari hasil kriteria pemilihan model yaitu RMSE menunjukkan bahwa nilai
RMSE paling kecil adalah mixture
enam cluster dengan nilai RMSE sebesar 2,3883.

## Perbandingan ANFIS dan Mixture of ANFIS
Berdasarkan pembahasan sebelumnya yaitu analisis kecepatan angin rata-rata di Sumenep menggunakan ANFIS dan Mixture of ANFIS diperoleh RMSE pada masing- masing metode. Tabel 15 berikut merupakan hasil perbandingan metode ANFIS dan Mixture of ANFIS berdasarkan kriteria RMSE.
Tabel 15 Perbandingan ANFIS dan Mixture Metode Kriteria ANFIS Mixture of ANFIS RMSE in sample 3,81142 2,3883 RMSE out of sample 3,71321 2,6777 Tabel 15 menunjukkan bahwa nilai RMSE baik in sample ataupun out sample metode Mixture of ANFIS dengan enam cluster lebih kecil daripada ANFIS dengan input Zt-1 dan Zt-3. Hal ini menunjukkan bahwa metode Mixture of ANFIS lebih akurat dalam memprediksi kecepatan angin rata-rata di Sumenep.

## KESIMPULAN
Berdasarkan analisis dan pembahasan, analisis kecepatan angin rata-rata di Sumenep menggunakan ANFIS dengan dua input Z t-1 dan Z t-3 , dua membership function dengan tipe Gauss menghasilkan empat aturan dan 20 parameter. RMSE residual pada data in sample sebesar 3,81142. Nilai error minimum di-peroleh dengan epoch sebanyak 170. Nilai RMSE minimum pada analisis kecepatan angin rata-rata di Sumenep menggunakan mixture of ANFIS sebesar 2,3883 diperoleh ketika banyaknya cluster enam. Nilai residual dengan frekuensi paling banyak bernilai
disekitar nol. Sehingga hasil analisis kecepatan angin rata-rata dengan mixture of ANFIS lebih akurat daripada menggunakan ANFIS, karena nilai RMSE in sample dan out sample nya lebih kecil.
Penelitian selanjutnya yang dapat dikembangkan adalah dengan menggunakan banyaknya membership function yang berbeda serta type membership function yang lain, sehingga diperoleh pemodelan menggunakan mixture of ANFIS yang bisa dibandingkan dengan hasil pada penelitian ini. Forecasting menggunakan Mixture of ANFIS dapat dilakukan dengan langkah,yaitu diawali dengan meramalkan cluster untuk kecepatan angin pada waktu ke t+k , misalkan cluster j. Apabila pada waktu sebelumnya berada dalam cluster yang sama atau dari tahun ke tahun berada dalam cluster yang tetap yaitu cluster j , maka yang digunakan untuk meramalkan kecepatan angin pada waktu ke t+k menggunakan model cluster j . Sedangkan bila terjadi perubahan cluster pada tahun-tahun sebelumnya, maka peramalan kecepatan angin menggunakan probabilitas state transition .

## DAFTAR PUSTAKA
[1] Benhammadi, F., Bey, K. B., Mokhtari, A., dan Gessoum, Z. 2010. Mixture of ANFIS System For CPU Load Prediction In Metacomputing Environment . Future Generation Computer System hal. 1003- 1011.
[2] Daryono. 2010. Cuaca Ekstrim Dampak La Nina . [http:// daryonobmkg.wordpress.com] diakses pada tanggal 1 Maret 2011 Pukul 10.05 [3] Faulina, R. 2011. Adaptive Neuro- Fuzzy Inference System untuk

## Peramalan Kecepatan Angin Rata-
Rata Harian di Sumenep . Laporan Tugas Akhir Program Sarjana Jurusan Statistika FMIPA. Surabaya : ITS.
[4] Irhamah, Prasetyo, D.D., dan Fithriasari, K. 2010. Pengembangan
Model Ramalan Kecepatan Angin Menggunakan Hybrid Time Series dan Algoritma Genetika . Penelitian Produktif ITS – 2010. Surabaya : ITS. [5] Jang, R., Jyh-Shing. 1993. ANFIS:
Adaptive Network Based Fuzzy Inference System . Transactions on Systems, Man and Cybernetics Vol. 23 No. 3.
[6]  James, N., dan Castellanos, F. 2009. Average Hourly Wind Speed Forecasting with ANFIS . 11-th Americas Conference on Wind Engineering. San Juan-Puerto Rico. [7] Kabar bisnis. 2011. Cuaca ekstrem,
Pelayaran Kalianget-Jangkar lumpuh . [http://www.kabarbisnis.com/umum /transportasi/2817377- Cuaca_ekstrem__Pelayaran_ Kalianget _Jangkar_lumpuh.html] diakses pada tanggal 9 Februari 2011 Pukul 12.40.
[8] Kusumadewi, S., dan Hartati, S. 2006. Neuro-Fuzzy Integrasi Sistem
Fuzzy dan Jaringan Syaraf . Yogyakarta : Graha Ilmu. [9] Naba, E. A. 2009. Belajar Cepat
Fuzzy Logic Menggunakan Matlab .
Yogyakarta : Andi.
[10] Widodo, T. S. 2005. Sistem Neuro Fuzzy untuk Pengolahan Informasi, Pemodelan, dan Kendali . Yogyakarta : Graha Ilmu.