APLIKASI MODEL HIDDEN MARKOV DISKRET UNTUK MENDETEKSI PENYALAHGUNAAN KARTU KREDIT
PUTRI UTARI
DEPARTEMEN MATEMATIKA FAKULTAS MATEMATIKA DAN ILMU PENGETAHUAN ALAM INSTITUT PERTANIAN BOGOR BOGOR 2012
ABSTRAK PUTRI UTARI. Aplikasi Model Hidden Markov Diskret untuk Mendeteksi Penyalahgunaan Kartu Kredit. Dibimbing oleh BERLIAN SETIAWATY dan NGAKAN KOMANG KUTHA ARDANA. Saat ini penggunaan kartu kredit sudah semakin meningkat. Hal tersebut menyebabkan penyalahgunaan terhadap kartu kredit juga meningkat. Jika rincian transaksi kartu kredit merupakan proses observasi dan penyebab besarnya transaksi tidak diamati secara langsung dan diasumsikan membentuk suatu rantai Markov, maka pasangan dari proses observasi dan penyebabnya dapat dimodelkan dengan model Hidden Markov. Model Hidden Markov dicirikan oleh parameter-parameternya yaitu matriks peluang state transisi, matriks peluang dari proses observasi, dan vektor peluang state awal. Penduga parameter dalam karya ilmiah ini diestimasi dengan menggunakan algoritme Rabiner yang meliputi algoritme forward-backward, algoritme Viterbi, dan algoritme Baum-Welch. Model Hidden Markov diskret diaplikasikan pada data transaksi suatu kartu kredit untuk mendeteksi penyalahgunaan terhadap kartu kredit. Proses deteksi kartu kredit dilakukan dengan menghitung peluang observasi untuk setiap transaksi baru. Jika perbandingan selisih dari peluang observasi transaksi baru dengan transaksi sebelumnya lebih besar atau sama dengan nilai ambang batas maka transaksi baru tersebut terdeteksi sebagai penyalahgunaan. Untuk mempermudah mencari penduga parameter model Hidden Markov diskret, dibuat suatu program komputasi dengan menggunakan Mathematica 7.0. Dari hasil yang diperoleh, model Hidden Markov diskret dapat memodelkan transaksi kartu kredit dengan cukup baik. Akurasi yang diperoleh dari hasil deteksi mencapai 77%. Kata kunci:
transaksi kartu kredit, deteksi penyalahgunaan, model Hidden Markov diskret, algoritme Rabiner.
ABSTRACT PUTRI UTARI. Application of Discrete Hidden Markov Model to Credit Card Fraud Detection. Supervised by BERLIAN SETIAWATY and NGAKAN KOMANG KUTHA ARDANA. Today the use of credit cards has increased which causes the increase of credit card fraud. If credit card transactions are observation process and the cause of magnitude transaction is assumed to form a Markov chain, then credit card transactions can be modeled by the Hidden Markov model. The Hidden Markov model is characterized by its parameters which are state transition probability matrix, observation symbol probability matrix, and initial state probability vector. Parameters of model are estimated using Rabiner algorithm that consists of forward-backward algorithm, Viterbi algorithm, and Baum-Welch algorithm. The discrete Hidden Markov model is applied to data transactions of a credit card to detect the fraud of credit card. Detection process of credit card is done by counting probability of observation for each new transaction. If the ratio of difference of the new transaction observation probability with previous transaction is greater than or equal to the threshold value then the new transaction is detected as fraud. The estimation of parameters are implemented on computational program by using Mathematica version 7.0. From the results obtained, it can be shown that discrete Hidden Markov model can model credit card transactions quite well. The accuracy of detection is 77%. Keywords: credit card transactions, fraud detection, discrete Hidden Markov model, Rabiner algorithm.
APLIKASI MODEL HIDDEN MARKOV DISKRET UNTUK MENDETEKSI PENYALAHGUNAAN KARTU KREDIT
PUTRI UTARI
Skripsi Sebagai salah satu syarat untuk memperoleh gelar Sarjana Sains pada Departemen Matematika
DEPARTEMEN MATEMATIKA FAKULTAS MATEMATIKA DAN ILMU PENGETAHUAN ALAM INSTITUT PERTANIAN BOGOR BOGOR 2012
Judul Skripsi : Aplikasi Model Hidden Markov Penyalahgunaan Kartu Kredit Nama : Putri Utari NIM : G54080058
Diskret
untuk
Mendeteksi
Menyetujui Pembimbing I
Pembimbing II
Dr. Berlian Setiawaty, MS. NIP. 19650505 198903 2 004
Ir. Ngakan Komang Kutha Ardana, M.Sc. NIP. 19640823 198903 1 001
Mengetahui Ketua Departemen Matematika
Dr. Berlian Setiawaty, MS. NIP. 19650505 198903 2 004
Tanggal Lulus : ............................................................
PRAKATA Puji syukur penulis panjatkan kepada Allah SWT atas segala rahmat dan karunia-Nya sehingga penulis dapat menyelesaikan karya ilmiah ini. Penyusunan karya ilmiah ini tidak terlepas dari dukungan dan bantuan dari berbagai pihak. Pada kesempatan kali ini, penulis juga ingin mengucapkan terima kasih yang sebesar-besarnya kepada: 1. Bapak dan Ibu tersayang, terima kasih atas kasih sayang, didikan, nasihat, semangat serta doa yang tiada henti-hentinya buat penulis. 2. Dr. Berlian Setiawaty, MS. selaku dosen pembimbing I dan Ir. N. K. Kutha Ardhana, M.Sc. selaku pembimbing II. Terima kasih atas waktu, ilmu yang diberikan dan kesabarannya dalam membimbing penulis. 3. Dr. Ir. I Wayan Mangku, M.Sc. selaku dosen penguji. Terima kasih atas waktu dan ilmu yang sangat bermanfaat bagi penulis. 4. Kakakku tersayang (Wisnu Mahadi), nenekku tersayang, terima kasih atas doa, semangat dan dukungan yang diberikan kepada penulis. 5. Semua dosen Departemen Matematika, terima kasih atas ilmu yang diberikan. 6. Seluruh staf pegawai Departemen Matematika, terima kasih atas bantuannya dalam memperlancar administrasi akademik bagi penulis. 7. Tante Ria, Om Iman, Tante Nita, Om Son, Tante Wiwik, dan Tante Euis. Terima kasih atas doa, semangat dan dukungan yang diberikan kepada penulis. 8. Sepupu-sepupuku tersayang: Mas Heru, Mba Lita, Rahma, Hana, Raja, Dimas, Dinda, Dara, Fahrezi, Fahira. Terima kasih atas doa, semangat, dan dukungan yang diberikan kepada penulis. 9. Guru-guru SMAN 2 Tangerang, SMPN 5 Tangerang, SDN 3 Poris Gaga, dan TK Islam AlFathir. Berkat jasa-jasa kalian penulis bisa seperti ini. 10. Teman-teman Matematika angkatan 45: Fitriyah, Nurul, Fuka, Maya, Dewi, Ana, Yunda, Tia, Mega, Fina, Haryanto, Hendri, Izzudin, Wulan, Vivi, Isna, Fenny, Hardono, Tiwi, Fikri, Ade, dan lainnya. Terima kasih atas doa, semangat, dan dukungan yang diberikan kepada penulis. 11. Kakak-kakak Matematika angkatan 44. 12. Adik-adik Matematika angkatan 46. 13. Kakak-kakak kost Pondok Sabrina: Mba Ummu, Mba Khusnul, Mba Noja, Mba Chemmy, Mba Devi, Mba Anif. Terima kasih atas doa dan semangat yang diberikan kepada penulis. Penulis menyadari tulisan ini masih memiliki kekurangan dan jauh dari kesempurnaan. Oleh karena itu dibutuhkan kritik dan saran yang membangun dari pembaca. Semoga karya ilmiah ini bermanfaat bagi kita semua, bagi dunia ilmu pengetahuan khususnya Matematika.
Bogor, Desember 2012
Putri Utari
RIWAYAT HIDUP Penulis dilahirkan di Jakarta pada tanggal 9 November 1990 sebagai anak kedua dari dua bersaudara, anak dari Murjono dan Yuliastuti. Pada tahun 1996 penulis menyelesaikan pendidikan di TK Islam Al Fathir Tangerang. Pada tahun 2002 penulis menyelesaikan pendidikan di SD Negeri Poris Gaga 3 Tangerang. Kemudian penulis melanjutkan pendidikan di SMP Negeri 5 Tangerang pada tahun yang sama. Pada tahun 2008 penulis menyelesaikan pendidikan di SMA Negeri 2 Tangerang. Pada tahun 2008 penulis diterima di Institut Pertanian Bogor (IPB) melalui jalur Ujian Seleksi Masuk IPB (USMI) departemen Matematika, Fakultas Matematika dan Ilmu Pengetahuan Alam. Selama mengikuti kegiatan perkuliahan, penulis aktif di berbagai kegiatan mahasiswa yaitu sebagai staf Divisi Keilmuan Gugus Mahasiswa Matematika (GUMATIKA) periode 2009-2010, anggota Divisi Konsumsi pada acara Masa Perkenalan Departemen Matematika (MPD) pada tahun 2010. Selain itu pada tahun yang sama penulis juga menjadi anggota Divisi Kesekretariatan Matematika Ria 2010. Pada tahun 2011, penulis diamanahkan kembali menjadi anggota Divisi Kesekretariatan pada acara Matematika Ria 2011.
DAFTAR ISI Halaman DAFTAR GAMBAR ............................................................................................................. viii DAFTAR TABEL ................................................................................................................. viii DAFTAR LAMPIRAN ......................................................................................................... I.
viii
PENDAHULUAN .......................................................................................................... 1.1 Latar Belakang ........................................................................................................ 1.2 Tujuan ..................................................................................................................... 1.3 Sistematika Penulisan ..............................................................................................
1 1 1 1
II. LANDASAN TEORI ..................................................................................................... 2.1 Pengantar Teori Peluang .......................................................................................... 2.2 Rantai Markov ......................................................................................................... 2.3 Algoritme Expectation Maximization (EM) .............................................................
2 2 3 4
III. MODEL HIDDEN MARKOV DISKRET ...................................................................... 3.1 Algoritme Rabiner ................................................................................................... 3.1.1 Algoritme Forward-Backward ........................................................................ 3.1.2 Algoritme Viterbi ............................................................................................ 3.1.3 Algoritme Baum-Welch .................................................................................. 3.2 Algoritme K-Means Clustering ...............................................................................
5 5 5 6 6 7
IV. APLIKASI DETEKSI PENYALAHGUNAAN KARTU KREDIT ................................. 4.1 Data Input Transaksi Kartu Kredit ........................................................................... 4.2 Pemodelan Transaksi Kartu Kredit .......................................................................... 4.3 Proses Deteksi Kartu Kredit ..................................................................................... 4.4 Hasil Komputasi ...................................................................................................... 4.4.1 Hasil Ketepatan Dugaan Barisan Observasi ..................................................... 4.4.2 Hasil Deteksi ..................................................................................................
8 8 8 9 10 10 10
V. SIMPULAN DAN SARAN ............................................................................................ 5.1 Simpulan ................................................................................................................. 5.2 Saran .......................................................................................................................
12 12 12
DAFTAR PUSTAKA ............................................................................................................
13
DAFTAR GAMBAR Halaman 1
Grafik persentase untuk setiap profil pengeluaran ............................................................
8
2
Proses deteksi penyalahgunaan kartu kredit .....................................................................
9
3
Grafik persentase ketepatan dugaan barisan observasi .....................................................
10
4
Grafik nilai True Positive dan False Positive untuk
11
dan
.............................
DAFTAR TABEL 1
Halaman Tabel persentase ketepatan dugaan barisan observasi ....................................................... 10
2
Tabel akurasi deteksi untuk
dan
.................................................................
11
DAFTAR LAMPIRAN Halaman 1
Data transaksi kartu kredit ..............................................................................................
15
2
Hasil deteksi penyalahgunaan kartu kredit untuk
...............................
18
3
Program untuk mencari parameter dengan menggunakan Mathematica 7.0 ....................
20
4
Output nilai parameter awal dan parameter duga ............................................................
30
5
Program deteksi penyalagunaan kartu kredit menggunakan Mathematica 7.0 .................
33
dan
viii
1
I. PENDAHULUAN 1.1 Latar Belakang Dalam kehidupan sehari-hari terdapat banyak kejadian yang dapat dimodelkan dengan suatu proses stokastik. Setiap kejadian itu terkait erat dengan penyebab kejadiannya. Jika penyebab kejadian tersebut tidak diamati secara langsung dan membentuk rantai Markov, maka pasangan kejadian dan penyebabnya dapat dimodelkan dengan model Hidden Markov. (Hidden Markov Model, HMM) Karakteristik dari model Hidden Markov dicirikan oleh parameter-parameternya yaitu sebaran state awal, matriks transisi dari rantai Markov, dan sebaran peluang dari proses observasi. Parameter-parameter tersebut dapat diduga dengan menggunakan algoritme Rabiner (Rabiner 1989) yang meliputi algoritme forward-backward, algoritme Viterbi, dan algoritme Baum-Welch. Model Hidden Markov juga dapat memodelkan beberapa aplikasi dengan baik. Contohnya pada Kajian Model Hidden Markov Diskret dengan Algoritme Rabiner dan Aplikasinya pada DNA (Wijayanti 2010). Dalam karya ilmiah ini akan dibahas aplikasi model Hidden Markov untuk mendeteksi penyalahgunaan kartu kredit yang diambil dari jurnal berjudul Credit Card Fraud Detection Using Hidden Markov Model (Srivastava et al. 2008). Kemajuan teknologi saat ini telah jauh meningkat. Hal tersebut dapat dilihat dari meningkatnya penggunaan kartu kredit secara dramatis. Saat ini kartu kredit telah banyak digunakan oleh masyarakat dalam proses pembayaran untuk suatu pembelian baik secara online maupun tidak online. Dengan meningkatnya penggunaan kartu kredit di seluruh dunia, maka ada peluang bagi seseorang untuk melakukan penyalahgunaan kartu kredit. Hal tersebut menyebabkan penyalahgunaan terhadap kartu kredit juga meningkat. Jika rincian transaksi kartu kredit merupakan proses observasi dan penyebab besarnya transaksi kartu kredit tersebut diasumsikan membentuk rantai Markov, maka pasangan dari proses observasi dan penyebabnya dapat dimodelkan dengan model Hidden Markov. Oleh karena itu, model Hidden Markov dapat digunakan untuk mendeteksi penyalahgunaan terhadap kartu kredit. Dengan menggunakan data rincian transaksi suatu kartu kredit, maka dapat
diduga parameter modelnya. Algoritme KMeans clustering digunakan untuk menentukan kelompok (cluster) dari data transaksi kartu kredit. Selanjutnya, dengan algoritme forward-backward akan ditentukan nilai peluang munculnya proses observasi. Kemudian dengan algoritme Viterbi akan dipilih hidden state yang memaksimumkan peluang bersama (joint probability) dari hidden state dan proses observasi. Akhirnya dengan algoritme Baum-Welch akan diperoleh nilai dugaan parameter model hidden Markov sedemikian sehingga dapat dengan baik mendeskripsikan rangkaian proses observasi yang terjadi (Wijayanti 2010). Kemudian dilanjutkan dengan proses deteksi penyalahgunaan kartu kredit dengan menghitung peluang observasi untuk setiap transaksi baru. Jika perbandingan selisih peluang observasi transaksi baru dengan peluang observasi transaksi sebelumnya lebih besar atau sama dengan nilai ambang batas, maka transaksi baru tersebut terdeteksi sebagai transaksi penyalahgunaan. Dalam perkembangan lebih lanjut, dibuat suatu program komputasi untuk menyelesaikan masalah model Hidden Markov. Software yang digunakan adalah Mathematica 7.0. Dalam karya ilmiah ini, program tersebut digunakan untuk membantu dalam penyelesaian masalah deteksi penyalahgunaan kartu kredit. 1.2 Tujuan Tujuan karya ilmiah ini adalah: 1. Mempelajari model Hidden Markov diskret dan pendugaan parameternya. 2. Mengimplementasikan model Hidden Markov diskret untuk mendeteksi penyalahgunaan terhadap kartu kredit. 1.3 Sistematika Penulisan Dalam karya ilmiah ini, terlebih dahulu pada Bab I dijelaskan latar belakang dan tujuan penulisan karya ilmiah ini. Pada Bab II dijelaskan beberapa definisi dan teorema yang digunakan pada pembahasan selanjutnya. Kemudian pada Bab III dijelaskan tentang model Hidden Markov diskret dengan algoritme Rabiner dan algoritme K-Means clustering. Setelah itu pada Bab IV dijelaskan aplikasi deteksi penyalahgunaan terhadap kartu kredit, serta Bab V berisi simpulan dan saran.
2
II. LANDASAN TEORI Dalam bab ini dijelaskan beberapa definisi serta teorema yang digunakan dalam pembahasan selanjutnya. 2.1 Pengantar Teori Peluang Definisi 2.1.1 (Percobaan Acak) Dalam suatu percobaan seringkali dilakukan pengulangan yang biasanya dilakukan dalam kondisi yang sama. Semua kemungkinan hasil yang akan muncul dapat diketahui, namun hasil pada percobaan berikutnya tidak dapat diketahui dengan tepat. Percobaan semacam ini, yang dapat diulang dalam kondisi yang sama, disebut percobaan acak. (Ross 1996) Definisi 2.1.2 (Ruang Contoh dan Kejadian) Himpunan dari semua kemungkinan hasil dari suatu percobaan acak disebut ruang contoh, di notasikan dengan Suatu kejadian A adalah himpunan bagian dari (Ghahramani 2005) Definisi 2.1.3 ( Medan- ) adalah suatu himpunan yang anggotanya terdiri atas himpunan bagian dari yang memenuhi kondisi berikut: 1. , 2. Jika . 3. Jika . (Ghahramani 2005) Definisi 2.1.4 (Ukuran Peluang) Misalkan adalah dari ruang contoh . Ukuran peluang adalah suatu fungsi yang memenuhi: 1. 2. Jika adalah himpunan yang untuk saling lepas yaitu setiap pasangan maka
Definisi 2.1.6 (Peluang Bersyarat) Misalkan adalah ruang peluang dan maka peluang A dengan syarat B didefinisikan sebagai (Ghahramani 2005) Teorema 2.1.1 (Teorema Bayes) Misalkan adalah ruang peluang dan Misalkan kejadian terjadi hanya dengan salah satu kejadian , maka peluang bersyarat dari setelah diketahui adalah
(Hogg et al. 2005) Definisi 2.1.7 (Peubah Acak Diskret) Misalkan adalah ruang contoh, adalah dari dan adalah himpunan berhingga. Suatu fungsi disebut peubah acak diskret jika memenuhi sifat berlaku . (Ghahramani 2005) Definisi 2.1.8 (Fungsi Kerapatan Peluang) Misalkan adalah ruang peluang dan adalah himpunan berhingga. Fungsi kerapatan peluang dari peubah acak diskret adalah fungsi didefinisikan oleh . (Grimmet dan Stirzaker 2001) Definisi 2.1.9 (Fungsi Kerapatan Peluang Bersama Dua Peubah Acak Diskret) Misalkan adalah ruang peluang dan adalah himpunan berhingga. Fungsi kerapatan peluang dari peubah acak diskret dan adalah fungsi didefinisikan oleh (Grimmet dan Stirzaker 2001)
Pasangan
disebut ruang peluang. (Ghahramani 2005)
Definisi 2.1.5 (Kejadian Saling Bebas) Kejadian dan dikatakan saling bebas jika Secara umum, himpunan kejadian dikatakan saling bebas jika untuk setiap himpunan bagian berhingga dari . (Grimmet dan Stirzaker 2001)
Definisi 2.1.10 (Fungsi Kerapatan Peluang Bersyarat) Jika dan merupakan peubah acak diskret, maka fungsi kerapatan peluang bersyarat dari jika diberikan terdefinisi untuk setiap sedemikian sehingga 0 . adalah (Ross 1996)
3
Definisi 2.1.11 (Fungsi Kerapatan Marginal) Misalkan adalah fungsi kerapatan peluang bersama dari peubah acak diskret dan Misalkan adalah himpunan nilai yang mungkin dari , dan adalah himpunan nilai yang mungkin dari . Selanjutnya fungsi dan masing-masing disebut fungsi kerapatan marginal dari dan . (Ghahramani 2005)
untuk
semua
kemungkinan
nilai
dari
(Grimmet dan Stirzaker 2001) Definisi 2.2.4 (Matriks Transisi) Misalkan adalah rantai Markov dengan ruang state yang berukuran Matriks transisi berukuran adalah matriks dari peluang transisi untuk (Grimmet dan Stirzaker 2001)
Definisi 2.1.12 (Bebas Stokastik Identik) Misalkan adalah peubah acak yang memiliki fungsi kepekatan yang sama yaitu sehingga
Definisi 2.2.5 (Rantai Markov Homogen) Rantai Markov dengan ruang state disebut homogen jika
(Grimmet dan Stirzaker 2001) dan fungsi kepekatan bersamanya adalah Definisi Peubah acak stokastik identik.
disebut bebas (Hogg et al. 2005)
Definisi 2.1.13 (Nilai Harapan Peubah Acak Diskret) Jika adalah peubah acak diskret dengan fungsi kerapatan peluang maka nilai harapan dari adalah (Ghahramani 2005) 2.2 Rantai Markov Definisi 2.2.1 (Ruang State) Misalkan adalah himpunan nilai dari barisan peubah acak, maka disebut ruang state. (Grimmet dan Stirzaker 2001) Definisi 2.2.2 (Proses Stokastik) Proses Stokastik yang terdefinisi pada ruang peluang adalah suatu himpunan dari peubah acak yang memetakan ruang contoh ke suatu state . (Ross 1996) Definisi 2.2.3 (Rantai Markov dengan Waktu Diskret) Misalkan adalah ruang peluang dan ruang state. Proses stokastik dengan ruang state , disebut rantai Markov dengan waktu diskret jika untuk setiap berlaku:
2.2.6 )
(Peluang
Transisi
Peluang Transisi dari rantai adalah peluang proses Markov berpindah dari ke langkah yang didefinisikan sebagai berikut: (Ross 1996) Definisi 2.2.7 (Terakses) Suatu dari rantai Markov disebut terakses (accessible) dari ditulis jika ada minimal sebuah bilangan (Ross 1996) Definisi 2.2.8 (Berkomunikasi) Dua state dari rantai Markov disebut berkomunikasi (communicate), ditulis , jika state dapat diakses dari state dan state dapat diakses dari . (Ross 1996) Definisi 2.2.9 (Kelas State) Suatu kelas state dari suatu rantai Markov, adalah suatu himpunan tak kosong sehingga semua pasangan state yang merupakan anggota dari berkomunikasi satu dengan yang lainnya, serta tidak ada state yang merupakan anggota yang berkomunikasi dengan suatu state yang bukan anggota dari . (Ross 1996)
4
Definisi 2.2.10 (Tak Tereduksi) Suatu rantai Markov disebut tak tereduksi (irreducible) jika hanya terdapat satu kelas state, yaitu jika semua state-nya berkomunikasi satu dengan yang lainnya. (Ross 1996)
Definisi 2.2.16 (Himpunan Konveks dan Fungsi Konveks) adalah himpunan vektor. Misalkan disebut sebagai himpunan konveks jika untuk semua terdapat sehingga
Definisi 2.2.11 (Berulang) dari suatu rantai Markov disebut State berulang (recurrent) jika (Ross 1996)
Misalkan merupakan fungsi dengan peubah yang terdefinisi pada himpunan konveks . Maka disebut sebagai fungsi konveks jika memenuhi persamaan . (Osborne 1997)
Definisi 2.2.12 (Periode, Periodik, dan Aperiodik) adalah rantai Markov Misalnya yang terdefinisi pada ( ) dengan ruang state . Suatu state disebut periode ditulis jikad adalah persekutuan terbesar (the greatest common divisor) bagi sehingga dinotasikan 0. Suatu state disebut periodik jika ( )>1 dan aperiodik jika (Ross 1996)
2.3 Algoritme EM (Expectation Maximization) Definisi 2.3.1 (Fungsi Likelihood) adalah himpunan Misalkan ukuran peluang yang terdefinisi pada dan kontinu absolut terhadap . Misalkan . Fungsi Likelihood yang digunakan untuk menghitung penduga parameter berdasarkan informasi adalah
Definisi 2.2.13 (Positive Recurrent dan Null Recurrent) Misalnya adalah rantai Markov yang terdefinisi pada ( ) dengan ruang state . Suatu state disebut positive recurrent jika state tersebut adalah recurrent dan berlaku jika proses dimulai dari state maka nilai harapan dari waktu sampai proses tersebut kembali ke state adalah bilangan berhingga (finite). State recurrent yang tidak positive recurrent disebut null recurrent. (Ross 1996) Definisi 2.2.14 (Ergodic) Rantai Markov yang positive recurrent dan aperiodik disebut ergodic. (Ross 1996) Definisi 2.2.15 (Nilai Harapan Rantai Markov Homogen) adalah rantai Markov Misalkan yang ergodic dengan ruang state berukuran dan misalkan merupakan matriks peluang transisi dan maka nilai harapan dari dinotasikan dengan memenuhi dan (Ross 1996)
(Elliott et al. 1995) Definisi 2.3.2 (Penduga Maksimum Likelihood) Penduga Maksimum Likelihood (MLE) didefinisikan oleh (Elliott et al. 1995) Definisi 2.3.3 (Algoritme EM) Umumnya MLE sulit dihitung secara langsung, oleh karena itu diperlukan suatu metode pendugaan berulang, yakni algoritme Expectation Maximization (EM). Langkahlangkah dalam metode tersebut adalah: 1. Set nilai awal parameter dengan . 2. Set dengan . 3. 4.
Cari . Ganti dengan dan ulangi langkah 2 sampai 4 hingga kriteria hentinya tercapai. Menurut Teorema Ketaksamaan Jensen, barisan yang dibangkitkan memberikan nilai fungsi likelihood yang tak turun yaitu: Bentuk likelihood bersyarat.
disebut
pseudo-log-
(Elliot et al. 1995)
5
III. MODEL HIDDEN MARKOV DISKRET Pada bab ini dibahas model Hidden Markov beserta karakteristiknya. Model ini terdiri atas pasangan merupakan penyebab kejadian yang tidak diamati secara langsung dan merupakan rantai Markov. Sedangkan adalah proses observasinya Karakteristik model Hidden Markov dapat dicirikan sebagai berikut: 1) adalah jumlah state rantai Markov yang tersembunyi dalam model. Ruang state yang ditunjukkan yaitu . 2) adalah jumlah simbol observasi yang berbeda. Ruang state yang ditunjukkan yaitu . 3) Matriks peluang state transisi , dimana
4) Matriks peluang dari proses observasi , dimana
5) Vektor peluang state awal dimana
3.1 Algoritme Rabiner Algoritme Rabiner terdiri atas algoritme forward-backward, algoritme Viterbi, dan algoritme Baum-Welch (Rabiner 1989). 3.1.1 Algoritme Forward-Backward Dari ciri-ciri di atas, diperoleh karakteristik model Hidden Markov berupa parameter . Barisan observasi dapat diperoleh dari banyaknya barisan state yang mungkin terjadi. Salah satu barisannya seperti dimana adalah state awal. Peluang observasi yang dihasilkan dari state ini yaitu: (1)
Bukti :
dan dapat diperluas menjadi, (2) Bukti :
. Peluang dari barisan state
yaitu: (3)
Bukti :
Dari persamaan (2) dan (3)
6
(4) Maka dengan demikian, peluang dari barisan observasi model Hidden Markov ditentukan oleh yang dapat ditulis sebagai berikut: (5) Bukti :
Algoritme forward-backward digunakan untuk menentukan peluang munculnya barisan observasi untuk suatu tertentu, yaitu
(9) Pada algoritme ini juga didefinisikan peubah backward,
Pada algoritme ini didefinisikan peubah forward, yaitu peluang bersama observasi dan berada pada state di waktu , jika diberikan model yaitu:
yaitu
(6) Prosedur algoritme forward adalah sebagai berikut: 1. Diberikan nilai awal untuk
(10) barisan observasi parsial dan berada pada state di waktu jika diberikan model Prosedur algoritme backward adalah sebagai berikut: 1. Diberikan nilai awal untuk peluang
(11) 2.
Dengan cara induksi akan diperoleh
(12) (7) (Bukti lihat Wijayanti 2010) 2.
Dengan cara induksi akan diperoleh
(8) 3.
Sehingga diperoleh
3.1.2 Algoritme Viterbi Algoritme Viterbi digunakan untuk menentukan pendugaan barisan state yang memiliki peluang maksimum yang selanjutnya digunakan untuk menduga barisan observasi . Pada masalah ini, akan dipilih sehingga barisan state adalah maksimum.
7
Didefinisikan
sebagai berikut: (13)
Barisan yang menghasilkan oleh
ditunjukkan
(14) Prosedur algoritme Viterbi adalah sebagai berikut: 1. Diberikan nilai awal untuk (15)
2.
. (16) Dengan cara rekursif akan diperoleh
3.
untuk Sehingga diperoleh
terjadi. Prosedur algoritme Baum-Welch adalah sebagai berikut. Didefinisikan yaitu peluang state pada waktu r, dan state j pada waktu , jika diberikan model dan rangkaian pengamatan.
Didefinisikan yaitu peluang berada di state i pada waktu r, jika diberikan rangkaian pengamatan 1,…, dan model .
(17)
(18) Kemudian dilakukan penelusuran balik barisan state optimal. Untuk (19) (Bukti lihat Wijayanti 2010) 3.1.3 Algoritme Baum-Welch Algoritme Baum-Welch digunakan untuk menentukan parameter dugaan model Hidden Markov. Pada masalah ini akan dimaksimumkan nilai peluang observasi untuk memperoleh nilai parameter model Hidden Markov yang dapat dengan baik mendeskripsikan rangkaian observasi yang
Dengan algoritme Expectation Maximization (EM) akan diperoleh yang memaksimumkan nilai peluang observasi Langkah-langkah algoritme EM sebagai berikut: 1. Diketahui sebuah penduga 2.
Hitung nilai fungsi Baum-Welch
3.
Tentukan
4.
Jika kriteria penghentian dipenuhi, maka keluar. Jika tidak terpenuhi, maka , dan ulangi terus langkah 2.
Dengan menggunakan metode Lagrange, dapat diperoleh formula re-estimasi parameter model Hidden Markov, yaitu: (23) (24)
(25) 3.2 Algoritme K-Means Clustering Algoritme K-Means clustering digunakan untuk menentukan kelompok (cluster). Pengelompokan dilakukan dengan
meminimumkan jumlah kuadrat dari jarak antara setiap data dengan pusat atau rata-rata kelompok yang sudah ditentukan.
8
IV. APLIKASI DETEKSI PENYALAHGUNAAN KARTU KREDIT Pada bab ini dibahas aplikasi deteksi penyalahgunaan kartu kredit. Berikut ini terlebih dahulu dijelaskan mengenai data input yang digunakan sebagai data observasi pada model. 4.1 Data Input Transaksi Kartu Kredit Dalam karya ilmiah ini, data yang digunakan merupakan data transaksi kartu kredit seseorang dari tanggal 17 Agustus 2011 sampai dengan 17 Februari 2012. Ada sebanyak 185 data observasi yang digunakan dalam kasus deteksi penyalahgunaan kartu kredit. Data tersebut dibagi menjadi 2 yaitu data untuk pelatihan sebanyak 150 data yang merupakan data transaksi asli dan data untuk pendeteksian sebanyak 35 data yang merupakan data campuran dari data transaksi asli dan data transaksi penyalahgunaan. Penyalahgunaan terhadap kartu kredit ada dua jenis, yaitu penyalahgunaan legal dan penyalahgunaan ilegal. Penyalahgunaan legal adalah penggunaan kartu kredit yang dilakukan oleh orang terdekat atau yang sudah dikenal oleh pemegang kartu kredit atas izin pemegang kartu kredit. Sedangkan penyalahgunaan ilegal adalah penggunaan kartu kredit yang dilakukan oleh orang yang tidak dikenal oleh pemegang kartu kredit tanpa izin dari pemegang kartu kredit. Dalam karya ilmiah ini, data penyalahgunaan yang digunakan merupakan data penyalahgunaan legal karena kartu kredit telah digunakan oleh salah satu keluarga dari pemegang kartu kredit. Data tersebut kemudian dikelompokkan menjadi 3 kelompok profil pengeluaran dengan menggunakan algoritme K-Means clustering, yaitu dan tinggi Dalam komputasi, data input diubah menjadi rendah = 1, sedang = 2, dan tinggi = 3. Data input dapat dilihat pada Lampiran 1. Persentase jumlah data transaksi untuk setiap profil pengeluaran pemegang kartu dapat dilihat pada Gambar 1.
Rendah Sedang Tinggi
88%
10% 2%
Gambar 1 Grafik persentase untuk setiap profil pengeluaran. 4.2 Pemodelan Transaksi Kartu Kredit Saat ini kemajuan teknologi sudah jauh meningkat. Hal tersebut dapat terlihat dari penggunaan kartu kredit yang juga meningkat. Dengan meningkatnya penggunaan kartu kredit, maka peluang bagi seseorang untuk melakukan penyalahgunaan kartu kredit juga semakin tinggi. Besar pengeluaran kartu kredit seseorang terjadi karena beberapa penyebab antara lain adanya discount pada suatu pusat perbelanjaan sehingga pemegang kartu kredit ingin banyak berbelanja sehingga menyebabkan pegeluarannya menjadi tinggi. Selain itu suasana hati pemegang kartu kredit juga dapat menyebabkan keinginan untuk banyak berbelanja sehingga pengeluarannya pun menjadi tinggi. Jenis barang pembelian juga bisa menjadi penyebab dari besarnya pengeluaran kartu kredit, misalnya apabila pemegang kartu kredit ingin membelanjakan barang-barang seperti elektronik atau perhiasan, maka pengeluaran kartu kredit akan tinggi juga. Penyebab kejadian besarnya pengeluaran kartu kredit seseorang diasumsikan bersifat Markov. Artinya meskipun di waktu yang lalu pernah terjadi banyak kejadian yang mempengaruhi besarnya pengeluaran kartu kredit seseorang tetapi penyebab pengeluaran kartu kredit seseorang saat ini cukup dipengaruhi oleh penyebab kejadian saat itu saja. Jadi karena penyebab kejadian pengeluaran kartu kredit seseorang membentuk suatu rantai Markov dan diasumsikan tidak diamati secara langsung, maka untuk mendeteksi kartu kredit dapat dimodelkan menggunakan model Hidden Markov. Proses observasi yang digunakan adalah besarnya pengeluaran transaksi kartu kredit seseorang per hari. Banyaknya data adalah 185, sedangkan penyebab kejadian yang tidak diamati secara langsung pada model adalah
9
penyebab terjadinya transaksi kartu kredit.
besar
pengeluaran
4.3 Proses Deteksi Kartu Kredit Proses pertama merupakan proses pelatihan untuk membangun model Hidden Markov. Kemudian proses selanjutnya merupakan proses pendeteksian transaksi kartu kredit. Setelah mendapatkan parameter model Hidden Markov, kemudian dilanjutkan dengan proses deteksi penyalahgunaan kartu kredit dengan menghitung peluang dari observasi. Misal barisan observasi sampai waktu maka peluang observasi menjadi Misal merupakan observasi baru dari transaksi yang baru pada waktu , maka peluang observasi menjadi
Kemudian didefinisikan
Jika berarti barisan baru yang diterima oleh model Hidden Markov mempunyai peluang yang rendah, maka itu bisa menjadi suatu penyalahgunaan. Transaksi baru yang ditambahkan diperkirakan sebagai penyalahgunaan jika persentase peluangnya di atas nilai ambang batas, yaitu Nilai ambang batas adalah nilai batas deteksi transaksi yang berbahaya (penyalahgunaan) yang dipertimbangkan secara empiris (Srivastava et al. 2008). Jika berbahaya, maka tidak ditambahkan dalam barisan observasi. Jika sebaliknya, maka akan ditambahkan dalam barisan secara permanen. Kemudian barisan baru akan digunakan sebagai barisan dasar untuk menentukan validitas transaksi berikutnya. Proses deteksi penyalahgunaan kartu kredit dapat dilihat pada Gambar 2.
.
Buat cluster (rendah, sedang, tinggi)
Data Transaksi
Identifikasi pengeluaran pemegang kartu kredit
Menyusun barisan data pelatihan
Membangun model Proses Pelatihan Proses Deteksi Transaksi
Masukkan simbol observasi
Tambahkan ke barisan yang tersedia
Tambahkan dalam barisan yang tersedia ke dalam barisan baru
Normal Test
Terdeteksi penyalahgunaan (berbahaya)
Menerima kedua barisan lama dan baru
Tidak Normal
Gambar 2 Proses deteksi penyalahgunaan kartu kredit.
Menghitung
10
4.4 Hasil Komputasi 4.4.1 Hasil Ketepatan Dugaan Barisan Observasi
Berikut adalah hasil persentase ketepatan dugaan yang lebih besar dari 60% dengan menggunakan software Mathematica 7.0.
Tabel 1 Persentase ketepatan dugaan barisan observasi Persentase Ketepatan Dugaan (%)
2
Frekuensi Persentase Ketepatan di atas 60% 9
3
41
95.56
61.54
100.00
92.67
31231
4
65
96.30
61.54
100.00
93.33
7053
5
100
98.52
69.23
100.00
96.00
3157
6
220
97.78
100.00
100.00
98.00
18854
Rendah
Sedang
Tinggi
Barisan Observasi
SeedRandom
87.41
61.54
100.00
85.33
31653
Persentase Ketepatan Dugaan Barisan Observasi (%)
100 90 80
92.67
93.33
96.00
98.00
85.33
70 60 50 40 30 20 10 0 2
3
4 5 6 N Gambar 3 Grafik ketepatan dugaan barisan observasi. Tabel persentase ketepatan di atas menunjukkan bahwa persentase ketepatan dugaan barisan observasi mencapai lebih dari 80%. Dari tabel tersebut dapat dilihat untuk persentase ketepatan barisan observasi mencapai 98.00%. Ini merupakan ketepatan dugaan yang terbaik untuk hasil yang diperoleh. Parameter-parameter model Hidden Markov dapat dilihat pada lampiran. 4.4.2 Hasil Deteksi Setelah mendapatkan tiga parameter model Hidden Markov dengan nilai yang berbeda-beda dan ketepatan dugaan barisan observasi yang mencapai lebih dari 80%, maka selanjutnya akan dilakukan proses deteksi dengan mencari peluang observasi
untuk setiap transaksi baru yang merupakan transaksi campuran antara transaksi asli dan transaksi yang disalahgunakan oleh orang lain. Pertama, nilai yang dipilih untuk pendeteksian adalah dan karena persentase ketepatan dugaan barisan obeservasi yang diperoleh lebih dari 95% dan merupakan hasil ketepatan barisan observasi yang terbaik. Kemudian menentukan nilai ambang batas yang digunakan untuk pendeteksian. Jadi, jika perbandingan selisih pada peluang observasi lebih dari nilai ambang batas, maka transaksi tersebut diduga sebagai transaksi berbahaya, walaupun itu merupakan transaksi asli. Sehingga transaksi baru tersebut tidak dapat dimasukkan dalam transaksi lama. Pendeteksian ini
11
menggunakan standar metrik True Positive (TP), False Positive (FP) dan metrik akurasi untuk mengukur keefektipan sistem. True Positive (TP) merepresentasikan sebagian data transaksi penyalahgunaan yang benar-benar terdeteksi sebagai penyalahgunaan. Sedangkan False Positive
(FP) merepresentasikan sebagian data transaksi asli diduga sebagai transaksi penyalahgunaan. Akurasi merepresentasikan sebagian total banyaknya transaksi (transaksi asli dan transaksi penyalahgunaan) yang terdeteksi dengan benar.
Hasil deteksi dapat dilihat pada grafik berikut:
N=5
100%
N=6 True Positive
80%
False Positive
60% 40%
33% 21%
20%
33% 14%
33%
33% 21%
17%
17% 17%
14%
7%
0% 2%
10% - 90%
91%
2%
10% - 90%
94%
Nilai Ambang Batas Gambar 4 Grafik Nilai True Positive dan False Positive untuk Dari hasil deteksi di atas dapat dilihat dan dengan nilai ambang untuk batas dari 10% sampai dengan 90%, nilai TP (True Positive) mencapai 33% yang artinya peluang data transaksi penyalahgunaan yang benar-benar terdeteksi sebagai penyalahgunaan sebesar 33%. Rendahnya nilai True Positive (33%) yang diperoleh dapat disebabkan oleh data transaksi yang digunakan merupakan data transaksi penyalahgunaan legal dengan pola data transaksi penyalahgunaan menyerupai pola data transaksi asli, sehingga model Hidden Markov sulit membedakan transaksi asli dan transaksi yang disalahgunakan. Nilai akurasi yang diperoleh dari hasil deteksi di atas mencapai 77%.
dan
Data
.
Asli
Penyalahgunaan
Jumlah data
Asli
25 (86%)
4 (67%)
29
Penyalahgunaan
4 (14%)
2 (33%)
6
Jumlah data
29
6
35
Deteksi
Tabel 2 Akurasi deteksi untuk dan dengan nilai ambang batas 10% sampai dengan 90% (banyaknya data yang terdeteksi dan persentasenya).
12
V. SIMPULAN DAN SARAN 5.1 Simpulan Dari pembahasan dapat ditarik simpulan: 1. Transaksi kartu kredit dapat dikaji dengan menggunakan model Hidden Markov diskret. 2. Dari data transaksi kartu kredit yang tersedia dapat dimodelkan dengan model Hidden Markov diskret dengan cukup baik. Hal ini dapat ditunjukkan dengan persentase ketepatan dari dugaan barisan observasi yang lebih dari 80%. 3. Hasil deteksi penyalahgunaan kartu kredit dan yang diperoleh untuk dengan nilai ambang batas dari 10% sampai dengan 90% yaitu nilai TP (True Positive) mencapai 33% dan akurasi mencapai 77%.
5.2 Saran 1. Hasil dalam karya ilmiah ini masih belum sempurna karena data yang digunakan data yang digunakan merupakan data penyalahgunaan legal yang jumlahnya terbatas, sehingga dimungkinkan untuk melakukan penelitian kembali dengan menggunakan data penyalahgunaan ilegal. 2. Program yang digunakan untuk mendeteksi penyalahgunaan masih menggunakan program manual atau menghitung peluang observasi satu per satu untuk setiap transaksi baru, sehingga membutuhkan waktu yang lama. Oleh karena itu, program untuk mendeteksi masih dapat diperbaiki agar hasil yang diperoleh lebih maksimal dan waktu yang dibutuhkan tidak terlalu lama.
13
DAFTAR PUSTAKA Elliot RJ, Aggoun L, Moore JB. 1995. Hidden Markov Models Estimation and Control. Springer Verlag. New York. Ghahramani S. 2005. Fundamental of Probability. Ed. ke-2. Prentice Hall. New Jersey. Grimmet GR, Stirzaker DR. 2001. Probability and Random Processes. Ed. ke-3. Clarendon Press. Oxford. Hogg RV, Craig AT, McKean JW. 2005. Introduction to Mathematical Statistics. Ed. ke-6. Prentice Hall. Englewood Cliffs. New Jersey. Osborne MJ. 1997. Concave and Convex Function of Many Variables.
(http://www.economics.about.com/Convex .htm ). Rabiner LR. 1989. A Tutorial on Hidden Markov Models and Selected Applications in Speech Recognition. Proc. IEEE, Vol 77 No.2, pg.257-286. Ross SM. 1996. Stochastic Processes. Ed. ke2. John Wiley & Sons. New York Srivastava A, Kundu A, Sural S, Majumdar AK. 2008. Credit Card Fraud Detection Using Hidden Markov Model. Proc. IEEE, Vol 5 No.1, pg.37-48. Wijayanti H. 2010. Kajian Model Hidden Markov Diskret dengan Algoritme Rabiner dan Aplikasinya pada DNA. [Tesis] IPB.
14
LAMPIRAN
15
Lampiran 1. Tabel data transaksi kartu kredit
No.
Tanggal Transaksi
Besar Transaksi
cluster
40
25-09-2011
0
1
41
26-09-2011
0
1
1
17-08-2011
1102100
3
42
27-09-2011
0
1
2
18-08-2011
522000
2
43
28-09-2011
0
1
3
19-08-2011
0
1
44
29-09-2011
0
1
4
20-08-2011
0
1
45
30-09-2011
0
1
5
21-08-2011
0
1
46
01-10-2011
0
1
6
22-08-2011
0
1
47
02-10-2011
0
1
03-10-2011
0
1
7
23-08-2011
0
1
48
8
24-08-2011
0
1
49
04-10-2011
0
1
9
25-08-2011
0
1
50
05-10-2011
2007500
3
10
26-08-2011
0
1
51
06-10-2011
184000
1
11
27-08-2011
0
1
52
07-10-2011
28700
1
12
28-08-2011
0
1
53
08-10-2011
0
1
13
29-08-2011
0
1
54
09-10-2011
0
1
14
30-08-2011
0
1
55
10-10-2011
0
1
11-10-2011
0
1
15
31-08-2011
0
1
56
16
01-09-2011
0
1
57
12-10-2011
0
1
17
02-09-2011
0
1
58
13-10-2011
0
1
18
03-09-2011
313120
2
59
14-10-2011
9830
1
19
04-09-2011
0
1
60
15-10-2011
0
1
20
05-09-2011
0
1
61
16-10-2011
101340
1
21
06-09-2011
307500
2
62
17-10-2011
271740
2
22
07-09-2011
0
1
63
18-10-2011
0
1
19-10-2011
0
1
23
08-09-2011
28700
1
64
24
09-09-2011
184000
1
65
20-10-2011
0
1
25
10-09-2011
0
1
66
21-10-2011
0
1
26
11-09-2011
0
1
67
22-10-2011
0
1
27
12-09-2011
0
1
68
23-10-2011
0
1
28
13-09-2011
0
1
69
24-10-2011
0
1
29
14-09-2011
0
1
70
25-10-2011
0
1
30
15-09-2011
0
1
71
26-10-2011
0
1
27-10-2011
0
1
31
16-09-2011
14966
1
72
32
17-09-2011
0
1
73
28-10-2011
0
1
33
18-09-2011
17455
1
74
29-10-2011
0
1
34
19-09-2011
0
1
75
30-10-2011
0
1
35
20-09-2011
0
1
76
31-10-2011
0
1
36
21-09-2011
0
1
77
01-11-2011
0
1
37
22-09-2011
0
1
78
02-11-2011
515000
2
38
23-09-2011
0
1
79
03-11-2011
337500
2
39
24-09-2011
0
1
80
04-11-2011
0
1
16
81
05-11-2011
0
1
124
18-12-2011
74600
1
82
06-11-2011
0
1
125
19-12-2011
0
1
83
07-11-2011
304000
2
126
20-12-2011
0
1
84
08-11-2011
28700
1
127
21-12-2011
0
1
85
09-11-2011
0
1
128
22-12-2011
0
1
86
10-11-2011
0
1
129
23-12-2011
0
1
87
11-11-2011
0
1
130
24-12-2011
0
1
88
12-11-2011
0
1
131
25-12-2011
0
1
89
13-11-2011
233000
2
132
26-12-2011
0
1
90
14-11-2011
0
1
133
27-12-2011
0
1
91
15-11-2011
0
1
134
28-12-2011
0
1
92
16-11-2011
383062
2
135
29-12-2011
0
1
93
17-11-2011
0
1
136
30-12-2011
0
1
94
18-11-2011
0
1
137
31-12-2011
0
1
95
19-11-2011
0
1
138
01-01-2012
0
1
96
20-11-2011
0
1
139
02-01-2012
507500
2
97
21-11-2011
257975
2
140
03-01-2012
0
1
98
22-11-2011
0
1
141
04-01-2012
0
1
99
23-11-2011
0
1
142
05-01-2012
184000
1
100
24-11-2011
0
1
143
06-01-2012
0
1
101
25-11-2011
0
1
144
07-01-2012
0
1
102
26-11-2011
0
1
145
08-01-2012
0
1
103
27-11-2011
0
1
146
09-01-2012
0
1
104
28-11-2011
0
1
147
10-01-2012
0
1
105
29-11-2011
0
1
148
11-01-2012
28700
1
106
30-11-2011
0
1
149
12-01-2012
0
1
107
01-12-2011
0
1
150
13-01-2012
0
1
108
02-12-2011
0
1
151
14-01-2012
0
1
109
03-12-2011
0
1
152
15-01-2012
0
1
110
04-12-2011
0
1
153
16-01-2012
438620
2
111
05-12-2011
741500
2
154
17-01-2012
0
1
112
06-12-2011
0
1
155
18-01-2012
0
1
113
07-12-2011
28700
1
156
19-01-2012
0
1
114
08-12-2011
0
1
157
20-01-2012
0
1
115
09-12-2011
0
1
158
21-01-2012
0
1
116
10-12-2011
0
1
159
22-01-2012
0
1
117
11-12-2011
0
1
160
23-01-2012
655000
2
118
12-12-2011
0
1
161
24-01-2012
0
1
119
13-12-2011
0
1
162
25-01-2012
347000
2
120
14-12-2011
0
1
163
26-01-2012
1002500
3
121
15-12-2011
0
1
164
27-01-2012
0
1
122
16-12-2011
363919
2
165
28-01-2012
0
1
123
17-12-2011
0
1
166
29-01-2012
0
1
17
167
30-01-2012
0
1
168
31-01-2012
195700
1
169
01-02-2012
0
1
170
02-02-2012
1007500
3
171
03-02-2012
0
1
172
04-02-2012
0
1
173
05-02-2012
0
1
174
06-02-2012
184000
1
175
07-02-2012
28700
1
176
08-02-2012
0
1
177
09-02-2012
404000
2
178
10-02-2012
589575
2
179
11-02-2012
0
1
180
12-02-2012
0
1
181
13-02-2012
0
1
182
14-02-2012
0
1
183
15-02-2012
125000
1
184
16-02-2012
425158
2
185
17-02-2012
0
1
Keterangan: Data yang bercetak tebal merupakan data penyalahgunaan.
18
Lampiran 2. Tabel hasil deteksi transaksi kartu kredit untuk
No. 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36
Perbandingan selisih peluang observasi (
Peluang Observasi 1.86899×10-24 1.24589×10-22 1.18287×10-21 1.18287×10-21 1.08105×10-22 1.19905×10-22 1.20033×10-22 1.20007×10-22 1.08201×10-22 1.08200×10-22 1.08201×10-22 9.88870×10-24 1.09681×10-23 9.94152×10-25 1.10290×10-24 1.10256×10-24 1.10912×10-24 1.07940×10-24 1.26906×10-24 1.21296×10-23 1.18021×10-23 2.46774×10-25 2.59863×10-24 2.60289×10-24 2.60210×10-24 2.60224×10-24 2.60222×10-24 2.60222×10-24 2.37823×10-25 2.63782×10-25 2.64064×10-25 2.64006×10-25 2.64017×10-25 2.64015×10-25 2.64015×10-25 2.41290×10-26
-1.22720×10-22 -1.05828×10-21 1.11787×10-32 1.07476×10-21 -1.17999×10-23 -1.27989×10-25 2.61840×10-26 1.18053×10-23 8.43529×10-28 -1.59704×10-28 9.83119×10-23 -1.07944×10-24 9.97399×10-24 -1.08751×10-25 3.42894×10-28 -6.55584×10-27 2.97121×10-26 -1.89657×10-25 -1.08606×10-23 3.27509×10-25 1.15554×10-23 -2.35186×10-24 -4.26386×10-27 7.94996×10-28 -1.42295×10-28 2.39374×10-29 -3.75252×10-30 2.36440×10-24 -2.59589×10-26 -2.81572×10-28 5.76035×10-29 -1.11888×10-29 2.05492×10-30 -3.57271×10-31 2.39886×10-25
9.45048×10-12 0.90861
dan
Terdeteksi berbahaya
Terdeteksi
0.00022 0.09837 7.79592×10-6 0.90861
Terdeteksi
0.90936
Terdeteksi
0.00031 0.02679
0.02700 0.97909
Terdeteksi
0.00031 9.19875×10-6 0.90861
Terdeteksi
0.00022 7.78327×10-6 0.90861
Terdeteksi
19
No. 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36
Peluang Observasi 8.06664×10-25 6.24452×10-23 5.63868×10-22 5.63868×10-22 3.09336×10-23 3.29730×10-23 3.27482×10-23 3.26724×10-23 3.06281×10-23 3.06214×10-23 3.06220×10-23 1.68120×10-24 1.79375×10-24 1.09646×10-25 1.17890×10-25 1.19632×10-25 1.28356×10-25 1.55113×10-25 2.42120×10-25 2.31251×10-24 2.76860×10-24 4.69934×10-26 4.52317×10-25 4.22473×10-25 4.21846×10-25 4.21641×10-25 4.21573×10-25 4.21551×10-25 2.31327×10-26 2.46578×10-26 2.44893×10-26 2.44322×10-26 2.44135×10-26 2.44072×10-26 2.44052×10-26 1.33946×10-27
-6.16386×10-23 -5.01423×10-22 -7.08149×10-29 5.32935×10-22 -2.03940×10-24 2.24761×10-25 7.58214×10-26 2.04426×10-24 6.73031×10-27 -6.28028×10-28 2.89408×10-23 -1.12548×10-25 1.68410×10-24 -8.24304×10-27 -1.74207×10-27 -8.72436×10-27 -2.67570×10-26 -8.70074×10-26 -2.07039×10-24 -4.56093×10-25 2.72161×10-24 -4.05323×10-25 2.98437×10-26 6.26328×10-28 2.04979×10-28 6.80476×10-29 2.25108×10-29 3.98418×10-25 -1.52507×10-27 1.68467×10-28 5.70784×10-29 1.87566×10-29 6.22138×10-30 2.05803×10-30 2.30657×10-26
Perbandingan selisih peluang observasi (
Terdeteksi berbahaya
0.94514
Terdeteksi
0.00682 0.00232 0.06257 0.00022 0.94510
Terdeteksi
0.93887
Terdeteksi
0.98303
Terdeteksi
0.06598 0.00148 0.00049 0.00016 0.00005 0.94513
Terdeteksi
0.00683 0.00233 0.00077 0.00025 0.00008 0.94512
Terdeteksi
Keterangan: Hasil yang bercetak tebal merupakan hasil deteksi untuk data penyalahgunaan.
20
Lampiran 3. Program untuk mencari parameter menggunakan Mathematica 7.0 Untuk
dengan batas Seedrandom
.
Data Input Transaksi Kartu Kredit
data 3, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1;
n 2; n adalah ruang state untuk memodelkan transaksi kartu kredit z0 30 000;SeedRandom untuk mencatat nilai awal yang dibangkitkan secara acak
Program Membangkitkan Nilai Awal Parameter HMM
Menentukan nilai awal A, B, yang dibangkitkan secara acak dan dicatat dengan SeedRandom
transAn_ : Blocka, b, c, a RandomReal0, 1, n , n ; b Plus a; c Transposea b;
transBn_ : Blocka, b, c, a RandomReal0, 1, 3, n ; b Plus a; c Transposea b;
phiAn_ : Blocka, b, c, a RandomReal0, 1, n , 1; b Plus a; c Transposea b;
HasilAkhir ;
Forn 2, n 3, frek 0; PresBar ; zSeedRan ; Forz z0, z 40 000, SeedRandom z; A transAn; B transBn; phi phiAn;
ALGORITME FORWARD
Hitung alpha peluang observasi untuk r1 sampai dengan R tB TransposeB;
bor_ : TransposetBdata1, r ; alpha1 Transposephi bo1;
Alpha NestListTransposeTranspose 1.A bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata 1;
alphaR AlphaLengthTransposedata 1, 1; PeluangObservasi Sum alphaRi, i, n; ALGORITME VITERBI
beta NestListA. 1 bo 2, 2 1 &,
ConstantArray1, n, LengthTransposedata, LengthTransposedata 1;
21
Delta NestListTableMax 1 TransposeAi, i, n bo 2 1,
2 1 &, alpha1, 1, LengthTransposedata 1;
deltaR DeltaLengthTransposedata, 1; PeluangStateMaks MaxdeltaR; BST TablePositionDeltai, 1, MaxDeltai, 1, i, LengthTransposedata;
BarisanStateTerbaik TableBSTi, 11, i, LengthTransposedata; BarisanObservasiDuga
TableRandomChoiceBBarisanStateTerbaiki 1, 2, 3, 1, i, LengthBarisanStateTerbaik Flatten;
banding Flattendata, BarisanObservasiDuga Transpose;
JmldtL_ : Countbanding, L, _;
JmldtdgR_, P_ : Countbanding, R , P ;
PresTepatdgK_ : JmldtdgK , K JmldtK 100 N; PresTepatdg1 PresTepatdg1; PresTepatdg2 PresTepatdg2; PresTepatdg3 PresTepatdg3;
PresTepDgBar Sum Jmldtdgi, i, i, 3 LengthTransposedata 100 N; IfPresTepatdg1 60 && PresTepatdg2 60 && PresTepatdg3 60, frek frek 1; PresBar AppendPresBar, PresTepDgBar; zSeedRan AppendzSeedRan, z; Printz, PresTepatdg1, PresTepatdg2, PresTepatdg3, PresTepDgBar, 0; z ; Printn; Printfrek; PrintMaxPresBar;
PrintzSeedRanPositionPresBar, MaxPresBar Flatten;
Hasil
n, frek, MaxPresBar,
zSeedRanPositionPresBar, MaxPresBar Flatten Flatten;
HasilAkhir AppendHasilAkhir, Hasil; n ; PrintHasilAkhir;
Pendugaan Parameter Dengan Algoritme Baum Welch THA TransposeHasilAkhir; X THA1; Y THA2; Z THA3; M THA4; Fori 1, i 3, n Xi; z Mi;
22
Program Pembangkit Matriks Peluang
transAn_ : Blocka, b, c, a RandomReal0, 1, n , n ; b Plus a; c Transposea b;
transBn_ : Blocka, b, c, a RandomReal0, 1, 3, n ; b Plus a; c Transposea b;
phiAn_ : Blocka, b, c, a RandomReal0, 1, n , 1; b Plus a; c Transposea b;
SeedRandom z; A transAn; B transBn;
phi phiAn; tB TransposeB;
bor_ : TransposetBdata1, r ;
alpha1 Transposephi bo1;
Alpha NestListTransposeTranspose 1.A bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata 1;
alphaR AlphaLengthTransposedata 1, 1;
beta NestListA. 1 bo 2, 2 1 &,
ConstantArray1, n, LengthTransposedata, LengthTransposedata 1;
albetb TableAlphai, 1 betaLengthTransposedata i, 1 boi 1, i, LengthTransposedata 1;
Epsi TableTableFlattenalbetbj Ai, i, n, j, LengthTransposedata 1;
SumEpsi TableSum Epsili, j, i, n, j, n, l, LengthTransposedata 1;
Epsilon TableEpsii SumEpsii, i, LengthTransposedata 1;
Gama TableTotalEpsiloni, j, i, LengthTransposedata 1, j, n;
SumEpsilon TableSum Epsiloni, jl, i, LengthTransposedata 1, j, n, l, n;
SumGama TableSum Gamai, j, i, LengthTransposedata 1, j, n; Adugaan TableSumEpsiloni, j SumGamai, i, n, j, n N; phidugaan ListGama1 N;
databedug DeleteFlattendata, LengthFlattendata; Bedug Transpose
TableSum GamaPositiondatabedug, 1i, 2, j, i, LengthPositiondatabedug, 1, j, n,
TableSum GamaPositiondatabedug, 2i, 2, j, i, LengthPositiondatabedug, 2, j, n,
TableSum GamaPositiondatabedug, 3i, 2, j, i, LengthPositiondatabedug, 3, j, n;
Bdugaan TableBedugi, j SumGamai, i, n, j, 3 N;
23
Program Output Nilai Dugaan Parameter Print Grid
" N", n,
"Seed Random yang mempunyai Ketepatan Dugaan Barisan Observasi Maksimal", z, " Ketepatan Dugaan Barisan Observasi",
Zi Grid" Nilai Awal",
"A", "", MatrixForm A,
" B", "", MatrixForm B,
" phi", "", MatrixForm phi Grid" Nilai Dugaan",
"A", "", MatrixForm Adugaan,
" B", "", MatrixForm Bdugaan,
" phi", "", MatrixForm phidugaan; i ;
Untuk
dengan batas Seedrandom
.
Data Input Transaksi Kartu Kredit
data 3, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1;
n 4; n adalah ruang state untuk memodelkan transaksi kartu kredit z0 3000; SeedRandom untuk mencatat nilai awal yang dibangkitkan secara acak
Program Membangkitkan Nilai Awal Parameter HMM
Menentukan nilai awal A, B, yang dibangkitkan secara acak dan dicatat dengan SeedRandom transAn_ : Blocka, b, c, a RandomReal0, 1, n , n ; b Plus a; c Transposea b;
transBn_ : Blocka, b, c, a RandomReal0, 1, 3, n ; b Plus a; c Transposea b;
phiAn_ : Blocka, b, c, a RandomReal0, 1, n , 1; b Plus a; c Transposea b;
HasilAkhir ;
24
Forn 4, n 5, frek 0; PresBar ; zSeedRan ; Forz z0, z 10 000, SeedRandom z; A transAn; B transBn; phi phiAn;
ALGORITME FORWARD
Hitung alpha peluang observasi untuk r1 sampai dengan R tB TransposeB;
bor_ : TransposetBdata1, r ; alpha1 Transposephi bo1;
Alpha NestListTransposeTranspose 1.A bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata 1;
alphaR AlphaLengthTransposedata 1, 1; PeluangObservasi Sum alphaRi, i, n; ALGORITME VITERBI
beta NestListA. 1 bo 2, 2 1 &,
ConstantArray1, n, LengthTransposedata, LengthTransposedata 1;
Delta
NestListTableMax 1 TransposeAi, i, n bo 2 1,
2 1 &, alpha1, 1, LengthTransposedata 1;
deltaR DeltaLengthTransposedata, 1; PeluangStateMaks MaxdeltaR; BST TablePositionDeltai, 1, MaxDeltai, 1, i, LengthTransposedata;
BarisanStateTerbaik TableBSTi, 11, i, LengthTransposedata; BarisanObservasiDuga
TableRandomChoiceBBarisanStateTerbaiki 1, 2, 3, 1, i, LengthBarisanStateTerbaik Flatten;
banding Flattendata, BarisanObservasiDuga Transpose;
JmldtL_ : Countbanding, L, _;
JmldtdgR_, P_ : Countbanding, R , P ;
PresTepatdgK_ : JmldtdgK , K JmldtK 100 N; PresTepatdg1 PresTepatdg1; PresTepatdg2 PresTepatdg2; PresTepatdg3 PresTepatdg3;
PresTepDgBar Sum Jmldtdgi, i, i, 3 LengthTransposedata 100 N; IfPresTepatdg1 60 && PresTepatdg2 60 && PresTepatdg3 60, frek frek 1; PresBar AppendPresBar, PresTepDgBar; zSeedRan AppendzSeedRan, z; Printz, PresTepatdg1, PresTepatdg2, PresTepatdg3, PresTepDgBar, 0; z ; Printn; Printfrek; PrintMaxPresBar;
PrintzSeedRanPositionPresBar, MaxPresBar Flatten;
25
Hasil
n, frek, MaxPresBar,
zSeedRanPositionPresBar, MaxPresBar Flatten Flatten;
HasilAkhir AppendHasilAkhir, Hasil; n ; PrintHasilAkhir; Pendugaan Parameter Dengan Algoritme Baum Welch THA TransposeHasilAkhir; X THA1; Y THA2; Z THA3; M THA4; Fori 1, i 3, n Xi; z Mi;
Program Pembangkit Matriks Peluang
transAn_ : Blocka, b, c, a RandomReal0, 1, n , n ; b Plus a; c Transposea b;
transBn_ : Blocka, b, c, a RandomReal0, 1, 3, n ; b Plus a; c Transposea b;
phiAn_ : Blocka, b, c, a RandomReal0, 1, n , 1; b Plus a; c Transposea b;
SeedRandom z; A transAn; B transBn;
phi phiAn; tB TransposeB;
bor_ : TransposetBdata1, r ;
alpha1 Transposephi bo1;
Alpha NestListTransposeTranspose 1.A bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata 1;
alphaR AlphaLengthTransposedata 1, 1;
beta NestListA. 1 bo 2, 2 1 &,
ConstantArray1, n, LengthTransposedata, LengthTransposedata 1;
albetb TableAlphai, 1 betaLengthTransposedata i, 1 boi 1, i, LengthTransposedata 1;
Epsi TableTableFlattenalbetbj Ai, i, n, j, LengthTransposedata 1;
SumEpsi TableSum Epsili, j, i, n, j, n, l, LengthTransposedata 1;
Epsilon TableEpsii SumEpsii, i, LengthTransposedata 1;
Gama TableTotalEpsiloni, j, i, LengthTransposedata 1, j, n;
SumEpsilon TableSum Epsiloni, jl, i, LengthTransposedata 1, j, n, l, n;
26
SumGama TableSum Gamai, j, i, LengthTransposedata 1, j, n; Adugaan TableSumEpsiloni, j SumGamai, i, n, j, n N; phidugaan ListGama1 N;
databedug DeleteFlattendata, LengthFlattendata; Bedug Transpose
TableSum GamaPositiondatabedug, 1i, 2, j, i, LengthPositiondatabedug, 1, j, n,
TableSum GamaPositiondatabedug, 2i, 2, j, i, LengthPositiondatabedug, 2, j, n,
TableSum GamaPositiondatabedug, 3i, 2, j, i, LengthPositiondatabedug, 3, j, n;
Bdugaan TableBedugi, j SumGamai, i, n, j, 3 N;
Program Output Nilai Dugaan Parameter Print Grid
" N", n,
"Seed Random yang mempunyai Ketepatan Dugaan Barisan Observasi Maksimal", z, " Ketepatan Dugaan Barisan Observasi",
Zi Grid" Nilai Awal",
"A", "", MatrixForm A,
" B", "", MatrixForm B,
" phi", "", MatrixForm phi Grid" Nilai Dugaan",
"A", "", MatrixForm Adugaan,
" B", "", MatrixForm Bdugaan,
" phi", "", MatrixForm phidugaan; i ;
Untuk
dengan batas Seedrandom
.
Data Input Transaksi Kartu Kredit
data 3, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1;
n 6; n adalah ruang state untuk memodelkan transaksi kartu kredit z0 10 000; SeedRandom untuk mencatat nilai awal yang dibangkitkan secara acak
27
Program Membangkitkan Nilai Awal Parameter HMM
Menentukan nilai awal A, B, yang dibangkitkan secara acak dan dicatat dengan SeedRandom transAn_ : Blocka, b, c, a RandomReal0, 1, n , n ; b Plus a; c Transposea b;
transBn_ : Blocka, b, c, a RandomReal0, 1, 3, n ; b Plus a; c Transposea b;
phiAn_ : Blocka, b, c, a RandomReal0, 1, n , 1; b Plus a; c Transposea b;
HasilAkhir ;
Forn 6, n 6, frek 0; PresBar ; zSeedRan ; Forz z0, z 20 000, SeedRandom z; A transAn; B transBn; phi phiAn;
ALGORITME FORWARD
Hitung alpha peluang observasi untuk r1 sampai dengan R tB TransposeB;
bor_ : TransposetBdata1, r ; alpha1 Transposephi bo1;
Alpha NestListTransposeTranspose 1.A bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata 1;
alphaR AlphaLengthTransposedata 1, 1; PeluangObservasi Sum alphaRi, i, n; ALGORITME VITERBI
beta NestListA. 1 bo 2, 2 1 &,
ConstantArray1, n, LengthTransposedata, LengthTransposedata 1;
Delta
NestListTableMax 1 TransposeAi, i, n bo 2 1,
2 1 &, alpha1, 1, LengthTransposedata 1;
deltaR DeltaLengthTransposedata, 1; PeluangStateMaks MaxdeltaR; BST TablePositionDeltai, 1, MaxDeltai, 1,
i, LengthTransposedata;
BarisanStateTerbaik TableBSTi, 11, i, LengthTransposedata; BarisanObservasiDuga
TableRandomChoiceBBarisanStateTerbaiki 1, 2, 3, 1, i, LengthBarisanStateTerbaik Flatten;
banding Flattendata, BarisanObservasiDuga Transpose;
JmldtL_ : Countbanding, L, _;
JmldtdgR_, P_ : Countbanding, R , P ;
28
PresTepatdgK_ : JmldtdgK , K JmldtK 100 N; PresTepatdg1 PresTepatdg1; PresTepatdg2 PresTepatdg2; PresTepatdg3 PresTepatdg3;
PresTepDgBar Sum Jmldtdgi, i, i, 3 LengthTransposedata 100 N; IfPresTepatdg1 60 && PresTepatdg2 60 && PresTepatdg3 60, frek frek 1; PresBar AppendPresBar, PresTepDgBar; zSeedRan AppendzSeedRan, z; Printz, PresTepatdg1, PresTepatdg2, PresTepatdg3, PresTepDgBar, 0; z ; Printn; Printfrek; PrintMaxPresBar;
PrintzSeedRanPositionPresBar, MaxPresBar Flatten;
Hasil
n, frek, MaxPresBar,
zSeedRanPositionPresBar, MaxPresBar Flatten Flatten;
HasilAkhir AppendHasilAkhir, Hasil; n ; PrintHasilAkhir; Pendugaan Parameter Dengan Algoritme Baum Welch THA TransposeHasilAkhir; X THA1; Y THA2; Z THA3; M THA4; Fori 1, i 1, n Xi; z Mi;
Program Pembangkit Matriks Peluang
transAn_ : Blocka, b, c, a RandomReal0, 1, n , n ; b Plus a; c Transposea b;
transBn_ : Blocka, b, c, a RandomReal0, 1, 3, n ; b Plus a; c Transposea b;
phiAn_ : Blocka, b, c, a RandomReal0, 1, n , 1; b Plus a; c Transposea b;
SeedRandom z; A transAn; B transBn;
phi phiAn; tB TransposeB;
bor_ : TransposetBdata1, r ;
alpha1 Transposephi bo1;
Alpha NestListTransposeTranspose 1.A bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata 1;
alphaR AlphaLengthTransposedata 1, 1;
29
beta NestListA. 1 bo 2, 2 1 &,
ConstantArray1, n, LengthTransposedata, LengthTransposedata 1;
albetb TableAlphai, 1 betaLengthTransposedata i, 1 boi 1, i, LengthTransposedata 1;
Epsi TableTableFlattenalbetbj Ai, i, n, j, LengthTransposedata 1;
SumEpsi TableSum Epsili, j, i, n, j, n, l, LengthTransposedata 1;
Epsilon TableEpsii SumEpsii, i, LengthTransposedata 1;
Gama TableTotalEpsiloni, j, i, LengthTransposedata 1, j, n;
SumEpsilon TableSum Epsiloni, jl, i, LengthTransposedata 1, j, n, l, n;
SumGama TableSum Gamai, j, i, LengthTransposedata 1, j, n; Adugaan TableSumEpsiloni, j SumGamai, i, n, j, n N; phidugaan ListGama1 N;
databedug DeleteFlattendata, LengthFlattendata; Bedug Transpose
TableSum GamaPositiondatabedug, 1i, 2, j, i, LengthPositiondatabedug, 1, j, n,
TableSum GamaPositiondatabedug, 2i, 2, j, i, LengthPositiondatabedug, 2, j, n,
TableSum GamaPositiondatabedug, 3i, 2, j, i, LengthPositiondatabedug, 3, j, n;
Bdugaan TableBedugi, j SumGamai, i, n, j, 3 N; Program Output Nilai Dugaan Parameter Print Grid
" N", n,
"Seed Random yang mempunyai Ketepatan Dugaan Barisan Observasi Maksimal", z, " Ketepatan Dugaan Barisan Observasi",
Zi Grid" Nilai Awal",
"A", "", MatrixForm A,
" B", "", MatrixForm B,
" phi", "", MatrixForm phi Grid" Nilai Dugaan",
"A", "", MatrixForm Adugaan,
" B", "", MatrixForm Bdugaan,
" phi", "", MatrixForm phidugaan; i ;
30
Lampiran 4. Output nilai parameter awal dan parameter duga N 2 Seed Random yang mempunyai Ketepatan Dugaan Barisan Observasi Maksimal 31 653 Ketepatan Dugaan Barisan Observasi 85.3333 Nilai Awal A B phi Nilai Dugaan A B phi
0.39331 0.60669 0.403884 0.596116 0.564087 0.391925 0.0439879 0.847429 0.0782065 0.0743648 0.0192293 0.980771
0.231007 0.768993 0.238936 0.761064 0.900029 0.0865064 0.0134647 0.898624 0.0879956 0.0133807 0.503985 0.496015
N 3 Seed Random yang mempunyai Ketepatan Dugaan Barisan Observasi Maksimal 31 231 Ketepatan Dugaan Barisan Observasi 92.6667 Nilai Awal A
B phi
0.138727 0.56445 0.511885 0.329679 0.340517 0.376118 0.952675 0.0165709 0.129403 0.675667 0.358725 0.174008 0.12018 0.310284
0.296823 0.158436 0.283366 0.0307538 0.19493 0.467267 0.569535
Nilai Dugaan A
B phi
0.414649 0.353072 0.82249 0.110859 0.690118 0.159525 0.815094 0.161045 0.938671 0.0535091 0.906822 0.0797661 0.428916 0.243442
0.232279 0.0666509 0.150357 0.0238612 0.00782024 0.0134117 0.327641
31
N 4 Seed Random yang mempunyai Ketepatan Dugaan Barisan Observasi Maksimal 7053 Ketepatan Dugaan Barisan Observasi 93.3333 Nilai Awal A
B phi
0.285066 0.33441 0.216447 0.164077 0.315302 0.337722 0.345903 0.00107366 0.280062 0.279989 0.0762525 0.363696 0.292378 0.29212 0.310122 0.105381 0.401679 0.0392994 0.559021 0.953992 0.0232288 0.0227796 0.37982 0.145581 0.474599 0.406329 0.485464 0.108207 0.319199 0.232895 0.0698628 0.378043
Nilai Dugaan A
B phi
0.125772 0.684398 0.109794 0.0800364 0.138245 0.686868 0.174367 0.000520464 0.135387 0.627848 0.0423803 0.194384 0.137883 0.639026 0.168146 0.0549452 0.90347 0.0838957 0.0126346 0.922063 0.0663569 0.0115806 0.868575 0.115359 0.0160663 0.900452 0.0858796 0.013668 0.248186 0.156038 0.373259 0.222516
N 5 Seed Random yang mempunyai Ketepatan Dugaan Barisan Observasi Maksimal 3157 Ketepatan Dugaan Barisan Observasi 96. Nilai Awal 0.079903 0.428402 0.188425 0.388075 0.272801
A
B
phi
0.388816 0.368145 0.0737536 0.0893823 0.199567 0.232335 0.0530097 0.0866862 0.28943 0.156885 0.251607 0.113653 0.156711 0.387102 0.0198785 0.0482338 0.223915 0.358247 0.0827858 0.0622508 0.310376 0.188719 0.500905 0.984584 0.0135389 0.00187681 0.39929 0.544795 0.0559157 0.537615 0.164098 0.298287 0.342201 0.41238 0.245418 0.262975 0.305094 0.175311 0.246236 0.010384
Nilai Dugaan
A
B
phi
0.0313804 0.239256 0.0951883 0.218934 0.13888
0.698871 0.237732 0.0192338 0.012783 0.510102 0.213353 0.0196586 0.0176297 0.669186 0.130316 0.084402 0.0209078 0.404626 0.359084 0.00744678 0.00990909 0.521715 0.299879 0.0279856 0.0115404 0.923218 0.0717466 0.00503574 0.888041 0.0912436 0.0207158 0.929881 0.0585486 0.0115707 0.857408 0.122958 0.0196339 0.885693 0.100398 0.0139085 0.0993446 0.268758 0.172116 0.256009 0.203773
32
N 6 Seed Random yang mempunyai Ketepatan Dugaan Barisan Observasi Maksimal 18 854 Ketepatan Dugaan Barisan Observasi 98. Nilai Awal
A
B
phi
0.14544 0.106327 0.283051 0.239758 0.130804 0.0946192 0.157999 0.0793258 0.319338 0.136935 0.147378 0.159025 0.154031 0.248953 0.101834 0.0674307 0.240047 0.187705 0.23784 0.0846836 0.0274142 0.298484 0.273154 0.078425 0.071521 0.195624 0.231019 0.156598 0.208749 0.136489 0.156938 0.28523 0.11371 0.0549448 0.380211 0.0089663 0.00134162 0.117761 0.880898 0.426909 0.256665 0.316426 0.375343 0.604857 0.0198 0.963948 0.0253693 0.0106826 0.0286946 0.753613 0.217692 0.380068 0.274274 0.345658 0.259436 0.055281 0.173323 0.230043 0.0632001 0.218717
Nilai Dugaan
A
B
phi
0.00261408 0.0602434 0.125473 0.771455 0.0140175 0.0261975 0.00411701 0.0651592 0.205224 0.63877 0.0228969 0.063832 0.00572443 0.291659 0.0933397 0.448626 0.0531906 0.10746 0.00397361 0.0445998 0.011296 0.892737 0.0272096 0.0201837 0.00165109 0.142362 0.131534 0.647183 0.0287328 0.048538 0.00648661 0.371637 0.115915 0.406555 0.0936978 0.00570886 0.921725 0.0702043 0.00807064 0.882474 0.105034 0.0124917 0.845064 0.129125 0.0258102 0.950961 0.0399153 0.00912395 0.892153 0.0956748 0.0121723 0.820026 0.153669 0.0263052 0.139107 0.162146 0.189087 0.203784 0.12666 0.179217
33
Lampiran 5. Program deteksi transaksi kartu kredit menggunakan Mathematica 7.0 Data Input Transaksi Kartu Kredit Data Pelatihan
data 3, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1; tB TransposeBdugaan; bor_ : TransposetBdata1, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata 1;
alphaR AlphaLengthTransposedata 1, 1; PeluangObservasi Sum alphaRi, i, n; PrintPeluangObservasi;
Pendeteksian Transaksi Kartu Kredit
Menghitung Peluang Observasi untuk Setiap Transaksi Baru
Data transaksi yang dideteksi = {1,1,2,1,1,1,1,1,1,2,1,2,3,1,1,1,1,1,1,3,1,1,1,1,1,1,2,2,1,1,1,1,1 ,2,1} data2 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1; tB TransposeBdugaan;
bor_ : TransposetBdata21, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata2 1;
alphaR AlphaLengthTransposedata2 1, 1; PeluangObservasi2 Sum alphaRi, i, n; PrintPeluangObservasi2;
deltaAlpha1 PeluangObservasi PeluangObservasi2 deltaAlpha1 PeluangObservasi data3 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1;
34
tB TransposeBdugaan;
bor_ : TransposetBdata31, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata3 1;
alphaR AlphaLengthTransposedata3 1, 1; PeluangObservasi3 Sum alphaRi, i, n; PrintPeluangObservasi3;
deltaAlpha2 PeluangObservasi2 PeluangObservasi3 deltaAlpha2 PeluangObservasi2 data4 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2; tB TransposeBdugaan;
bor_ : TransposetBdata41, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata4 1;
alphaR AlphaLengthTransposedata4 1, 1; PeluangObservasi4 Sum alphaRi, i, n; PrintPeluangObservasi4;
deltaAlpha3 PeluangObservasi3 PeluangObservasi4 deltaAlpha3 PeluangObservasi3 data5 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1; tB TransposeBdugaan;
bor_ : TransposetBdata51, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata5 1;
alphaR AlphaLengthTransposedata5 1, 1; PeluangObservasi5 Sum alphaRi, i, n; PrintPeluangObservasi5;
deltaAlpha4 PeluangObservasi4 PeluangObservasi5 deltaAlpha4 PeluangObservasi4
35
data6 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1; tB TransposeBdugaan;
bor_ : TransposetBdata61, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata6 1;
alphaR AlphaLengthTransposedata6 1, 1; PeluangObservasi6 Sum alphaRi, i, n; PrintPeluangObservasi6;
deltaAlpha5 PeluangObservasi5 PeluangObservasi6 deltaAlpha5 PeluangObservasi5 data7 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1; tB TransposeBdugaan;
bor_ : TransposetBdata71, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata7 1;
alphaR AlphaLengthTransposedata7 1, 1; PeluangObservasi7 Sum alphaRi, i, n; PrintPeluangObservasi7;
deltaAlpha6 PeluangObservasi6 PeluangObservasi7 deltaAlpha6 PeluangObservasi6 data8 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1;
36
tB TransposeBdugaan;
bor_ : TransposetBdata81, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata8 1;
alphaR AlphaLengthTransposedata8 1, 1; PeluangObservasi8 Sum alphaRi, i, n; PrintPeluangObservasi8;
deltaAlpha7 PeluangObservasi7 PeluangObservasi8 deltaAlpha7 PeluangObservasi7 data9 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1; tB TransposeBdugaan;
bor_ : TransposetBdata91, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata9 1;
alphaR AlphaLengthTransposedata9 1, 1; PeluangObservasi9 Sum alphaRi, i, n; PrintPeluangObservasi9;
deltaAlpha8 PeluangObservasi8 PeluangObservasi9 deltaAlpha8 PeluangObservasi8 data10 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1; tB TransposeBdugaan;
bor_ : TransposetBdata101, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata10 1;
alphaR AlphaLengthTransposedata10 1, 1; PeluangObservasi10 Sum alphaRi, i, n; PrintPeluangObservasi10;
deltaAlpha9 PeluangObservasi9 PeluangObservasi10 deltaAlpha9 PeluangObservasi9
37
data11 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2; tB TransposeBdugaan;
bor_ : TransposetBdata111, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata11 1;
alphaR AlphaLengthTransposedata11 1, 1; PeluangObservasi11 Sum alphaRi, i, n; PrintPeluangObservasi11;
deltaAlpha10 PeluangObservasi10 PeluangObservasi11 deltaAlpha10 PeluangObservasi10 data12 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1; tB TransposeBdugaan;
bor_ : TransposetBdata121, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata12 1;
alphaR AlphaLengthTransposedata12 1, 1; PeluangObservasi12 Sum alphaRi, i, n; PrintPeluangObservasi12;
deltaAlpha11 PeluangObservasi11 PeluangObservasi12 deltaAlpha11 PeluangObservasi11 data13 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 2;
38
tB TransposeBdugaan;
bor_ : TransposetBdata131, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata13 1;
alphaR AlphaLengthTransposedata13 1, 1; PeluangObservasi13 Sum alphaRi, i, n; PrintPeluangObservasi13;
deltaAlpha12 PeluangObservasi12 PeluangObservasi13 deltaAlpha12 PeluangObservasi12 data14 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 2, 3; tB TransposeBdugaan;
bor_ : TransposetBdata141, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata14 1;
alphaR AlphaLengthTransposedata14 1, 1; PeluangObservasi14 Sum alphaRi, i, n; PrintPeluangObservasi14;
deltaAlpha13 PeluangObservasi13 PeluangObservasi14 deltaAlpha13 PeluangObservasi13 data15 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 2, 1; tB TransposeBdugaan;
bor_ : TransposetBdata151, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata15 1;
alphaR AlphaLengthTransposedata15 1, 1; PeluangObservasi15 Sum alphaRi, i, n; PrintPeluangObservasi15;
deltaAlpha14 PeluangObservasi14 PeluangObservasi15 deltaAlpha14 PeluangObservasi14
39
data16 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 2, 1, 1; tB TransposeBdugaan;
bor_ : TransposetBdata161, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata16 1;
alphaR AlphaLengthTransposedata16 1, 1; PeluangObservasi16 Sum alphaRi, i, n; PrintPeluangObservasi16;
deltaAlpha15 PeluangObservasi15 PeluangObservasi16 deltaAlpha15 PeluangObservasi15 data17 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1; tB TransposeBdugaan;
bor_ : TransposetBdata171, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata17 1;
alphaR AlphaLengthTransposedata17 1, 1; PeluangObservasi17 Sum alphaRi, i, n; PrintPeluangObservasi17;
deltaAlpha16 PeluangObservasi16 PeluangObservasi17 deltaAlpha16 PeluangObservasi16 data18 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1;
40
tB TransposeBdugaan;
bor_ : TransposetBdata181, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata18 1;
alphaR AlphaLengthTransposedata18 1, 1; PeluangObservasi18 Sum alphaRi, i, n; PrintPeluangObservasi18;
deltaAlpha17 PeluangObservasi17 PeluangObservasi18 deltaAlpha17 PeluangObservasi17 data19 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1; tB TransposeBdugaan;
bor_ : TransposetBdata191, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata19 1;
alphaR AlphaLengthTransposedata19 1, 1; PeluangObservasi19 Sum alphaRi, i, n; PrintPeluangObservasi19;
deltaAlpha18 PeluangObservasi18 PeluangObservasi19 deltaAlpha18 PeluangObservasi18 data20 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1; tB TransposeBdugaan;
bor_ : TransposetBdata201, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata20 1;
alphaR AlphaLengthTransposedata20 1, 1; PeluangObservasi20 Sum alphaRi, i, n; PrintPeluangObservasi20;
deltaAlpha19 PeluangObservasi19 PeluangObservasi20 deltaAlpha19 PeluangObservasi19
41
data21 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 3; tB TransposeBdugaan;
bor_ : TransposetBdata211, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata21 1;
alphaR AlphaLengthTransposedata21 1, 1; PeluangObservasi21 Sum alphaRi, i, n; PrintPeluangObservasi21;
deltaAlpha20 PeluangObservasi20 PeluangObservasi21 deltaAlpha20 PeluangObservasi20 data22 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 3, 1; tB TransposeBdugaan;
bor_ : TransposetBdata221, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata22 1;
alphaR AlphaLengthTransposedata22 1, 1; PeluangObservasi22 Sum alphaRi, i, n; PrintPeluangObservasi22;
deltaAlpha21 PeluangObservasi21 PeluangObservasi22 deltaAlpha21 PeluangObservasi21 data23 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 3, 1;
42
tB TransposeBdugaan;
bor_ : TransposetBdata231, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata23 1;
alphaR AlphaLengthTransposedata23 1, 1; PeluangObservasi23 Sum alphaRi, i, n; PrintPeluangObservasi23;
deltaAlpha22 PeluangObservasi22 PeluangObservasi23 deltaAlpha22 PeluangObservasi22 data24 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 3, 1, 1; tB TransposeBdugaan;
bor_ : TransposetBdata241, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata24 1;
alphaR AlphaLengthTransposedata24 1, 1; PeluangObservasi24 Sum alphaRi, i, n; PrintPeluangObservasi24;
deltaAlpha23 PeluangObservasi23 PeluangObservasi24 deltaAlpha23 PeluangObservasi23 data25 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1; tB TransposeBdugaan;
bor_ : TransposetBdata251, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata25 1;
alphaR AlphaLengthTransposedata25 1, 1; PeluangObservasi25 Sum alphaRi, i, n; PrintPeluangObservasi25;
deltaAlpha24 PeluangObservasi24 PeluangObservasi25 deltaAlpha24 PeluangObservasi24
43
data26 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1; tB TransposeBdugaan;
bor_ : TransposetBdata261, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata26 1;
alphaR AlphaLengthTransposedata26 1, 1; PeluangObservasi26 Sum alphaRi, i, n; PrintPeluangObservasi26;
deltaAlpha25 PeluangObservasi25 PeluangObservasi26 deltaAlpha25 PeluangObservasi25 data27 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1; tB TransposeBdugaan;
bor_ : TransposetBdata271, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata27 1;
alphaR AlphaLengthTransposedata27 1, 1; PeluangObservasi27 Sum alphaRi, i, n; PrintPeluangObservasi27;
deltaAlpha26 PeluangObservasi26 PeluangObservasi27 deltaAlpha26 PeluangObservasi26 data28 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 2;
44
tB TransposeBdugaan;
bor_ : TransposetBdata281, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata28 1;
alphaR AlphaLengthTransposedata28 1, 1; PeluangObservasi28 Sum alphaRi, i, n; PrintPeluangObservasi28;
deltaAlpha27 PeluangObservasi27 PeluangObservasi28 deltaAlpha27 PeluangObservasi27 data29 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 2, 2; tB TransposeBdugaan;
bor_ : TransposetBdata291, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata29 1;
alphaR AlphaLengthTransposedata29 1, 1; PeluangObservasi29 Sum alphaRi, i, n; PrintPeluangObservasi29;
deltaAlpha28 PeluangObservasi28 PeluangObservasi29 deltaAlpha28 PeluangObservasi28 data30 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 2, 1; tB TransposeBdugaan;
bor_ : TransposetBdata301, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata30 1;
alphaR AlphaLengthTransposedata30 1, 1; PeluangObservasi30 Sum alphaRi, i, n; PrintPeluangObservasi30;
deltaAlpha29 PeluangObservasi29 PeluangObservasi30 deltaAlpha29 PeluangObservasi29
45
data31 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 2, 1, 1; tB TransposeBdugaan;
bor_ : TransposetBdata311, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata31 1;
alphaR AlphaLengthTransposedata31 1, 1; PeluangObservasi31 Sum alphaRi, i, n; PrintPeluangObservasi31;
deltaAlpha30 PeluangObservasi30 PeluangObservasi31 deltaAlpha30 PeluangObservasi30 data32 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 2, 1, 1, 1; tB TransposeBdugaan;
bor_ : TransposetBdata321, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata32 1;
alphaR AlphaLengthTransposedata32 1, 1; PeluangObservasi32 Sum alphaRi, i, n; PrintPeluangObservasi32;
deltaAlpha31 PeluangObservasi31 PeluangObservasi32 deltaAlpha31 PeluangObservasi31 data33 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1;
46
tB TransposeBdugaan;
bor_ : TransposetBdata331, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata33 1;
alphaR AlphaLengthTransposedata33 1, 1; PeluangObservasi33 Sum alphaRi, i, n; PrintPeluangObservasi33;
deltaAlpha32 PeluangObservasi32 PeluangObservasi33 deltaAlpha32 PeluangObservasi32 data34 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1; tB TransposeBdugaan;
bor_ : TransposetBdata341, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata34 1;
alphaR AlphaLengthTransposedata34 1, 1; PeluangObservasi34 Sum alphaRi, i, n; PrintPeluangObservasi34;
deltaAlpha33 PeluangObservasi33 PeluangObservasi34 deltaAlpha33 PeluangObservasi33 data35 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2; tB TransposeBdugaan;
bor_ : TransposetBdata351, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata35 1;
alphaR AlphaLengthTransposedata35 1, 1; PeluangObservasi35 Sum alphaRi, i, n; PrintPeluangObservasi35;
deltaAlpha34 PeluangObservasi34 PeluangObservasi35 deltaAlpha34 PeluangObservasi34
47
data36 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1; tB TransposeBdugaan;
bor_ : TransposetBdata361, r ;
alpha1 Transposephidugaan bo1; Alpha
NestListTransposeTranspose 1.Adugaan bo 2 1, 2 1 &, alpha1, 1, LengthTransposedata36 1;
alphaR AlphaLengthTransposedata36 1, 1; PeluangObservasi36 Sum alphaRi, i, n; PrintPeluangObservasi36;
deltaAlpha35 PeluangObservasi35 PeluangObservasi36 deltaAlpha35 PeluangObservasi35
